-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathgenit
executable file
·111 lines (96 loc) · 2.58 KB
/
genit
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
#!/bin/bash
case_dir=/work2/noaa/gsienkf/weihuang/jedi/case_study
run_dir=run_80.40t1n_36p
output=1
#caselist=(amsua iasi)
caselist=(amsua iasi)
for i in ${!caselist[@]}
do
gridfile=regrid/${caselist[$i]}
if [ ! -f ${gridfile} ]
then
cd regrid
rm -f fv3latlon.nc
#interp.sh ${case_dir}/sfc-letkf/analysis/increment/
interp.sh ${case_dir}/${caselist[$i]}/${run_dir}/analysis/increment/
cd ..
mv regrid/fv3latlon.nc ${gridfile}
fi
if [ "${caselist[$i]}" == "iasi" ]
then
obsfile=${case_dir}/${caselist[$i]}/manmade-iasi-obs/iasi_metop-a_obs_2020121500_m.nc4
else
obsfile=${case_dir}/${caselist[$i]}/manmade-amsua-obs/amsua_n19_obs_2020121500_m.nc4
fi
if [ -f ${obsfile} ]
then
python new-obs-jedi.py \
--output=${output} \
--casename=${caselist[$i]} \
--varname=T \
--gridfile=${gridfile} \
--obsfile=${obsfile}
fi
done
exit 0
caselist=(aircraft sfcship sondes)
#var_list=(tsen tv uv q)
var_list=(tsen)
for i in ${!caselist[@]}
do
gridfile=regrid/${caselist[$i]}
if [ ! -f ${gridfile} ]
then
cd regrid
rm -f fv3latlon.nc
#interp.sh ${case_dir}/sfc-letkf/analysis/increment/
interp.sh ${case_dir}/${caselist[$i]}/${run_dir}/analysis/increment/
cd ..
mv regrid/fv3latlon.nc ${gridfile}
fi
for j in ${!var_list[@]}
do
echo "i = $i, j = $j"
obsfile=${case_dir}/${caselist[$i]}/ioda_v2_data/${caselist[$i]}_${var_list[$j]}_obs_2020011006.nc4
if [ -f ${obsfile} ]
then
python new-obs-jedi.py \
--output=${output} \
--casename=${caselist[$i]}_${var_list[$j]} \
--varname=T \
--gridfile=${gridfile} \
--obsfile=${obsfile}
fi
done
done
caselist=(satwind scatwind vadwind windprof)
var_list=(uv)
run_dir=run_80.40t1n_36p
output=1
for i in ${!caselist[@]}
do
gridfile=regrid/${caselist[$i]}
if [ ! -f ${gridfile} ]
then
cd regrid
rm -f fv3latlon.nc
#interp.sh ${case_dir}/sfc-letkf/analysis/increment/
interp.sh ${case_dir}/${caselist[$i]}/${run_dir}/analysis/increment/
cd ..
mv regrid/fv3latlon.nc ${gridfile}
fi
for j in ${!var_list[@]}
do
echo "i = $i, j = $j"
obsfile=${case_dir}/${caselist[$i]}/ioda_v2_data/${caselist[$i]}_obs_2020011006.nc4
if [ -f ${obsfile} ]
then
python new-obs-jedi.py \
--output=${output} \
--casename=${caselist[$i]}_${var_list[$j]} \
--varname=T \
--gridfile=${gridfile} \
--obsfile=${obsfile}
fi
done
done