| 1 | #/usr/bin/ksh
|
|---|
| 2 | #
|
|---|
| 3 | # Script to rerun rose suite on it-pp-06, to compare against IDL version.
|
|---|
| 4 |
|
|---|
| 5 | ulimit -S -s unlimited
|
|---|
| 6 |
|
|---|
| 7 | #IFILE0=data/CASE_07_ref1d.nc
|
|---|
| 8 | IFILE0=data/CASE_12_ref1d.nc
|
|---|
| 9 | #IFILE0=data/CASE_05_ref1d.nc
|
|---|
| 10 | #IFILE0=data/CASE_51_ref1d.nc
|
|---|
| 11 |
|
|---|
| 12 | MFILE=data/MSIS_coeff.nc
|
|---|
| 13 |
|
|---|
| 14 | CFILE=data/default_wopt.nml
|
|---|
| 15 |
|
|---|
| 16 | THINFILE=data/ropp_thin_default.dat
|
|---|
| 17 |
|
|---|
| 18 | # These two might need commenting out, as they were originally missing from the rose suite version.
|
|---|
| 19 | GEOPOT_COEF=data/egm96.dat
|
|---|
| 20 | GEOPOT_CORR=data/corrcoef.dat
|
|---|
| 21 |
|
|---|
| 22 |
|
|---|
| 23 | #------------------------------------------------------
|
|---|
| 24 | #0. Cntl: Rose
|
|---|
| 25 | #------------------------------------------------------
|
|---|
| 26 |
|
|---|
| 27 | #if false ; then # bypass the control
|
|---|
| 28 |
|
|---|
| 29 | EXEC=/data/local/idculv/ropp_test_folder_02042019/gfortran/ropp_pp/tools/ropp_pp_wopt_tool
|
|---|
| 30 |
|
|---|
| 31 | IFILE=$IFILE0
|
|---|
| 32 |
|
|---|
| 33 | OFILE1=`echo $IFILE0 |sed -es/'.nc'/'_1.1.nc'/`
|
|---|
| 34 |
|
|---|
| 35 | TFILE=`echo $IFILE0 |sed -es/'.nc'/'_rose.out'/`
|
|---|
| 36 |
|
|---|
| 37 | echo "$EXEC $IFILE -o $OFILE1 -c $CFILE -f > $TFILE 2>&1"
|
|---|
| 38 | $EXEC $IFILE -o $OFILE1 -c $CFILE -f > $TFILE 2>&1
|
|---|
| 39 |
|
|---|
| 40 | # Now some messy preparation for invert tool
|
|---|
| 41 | ncrename -vBangle_FSI,bangle -ddim_009,dim_lev1b $OFILE1
|
|---|
| 42 | cmd="impact=Impact_height_FSI+6378137.0; lat_tp=impact*0.0; lon_tp=impact*0.0; azimuth_tp=impact*0.0; \
|
|---|
| 43 | impact_L1=impact; impact_L2=impact; bangle_L1=bangle; bangle_L2=bangle; \
|
|---|
| 44 | bangle_L1_sigma=impact*0.0; bangle_L2_sigma=impact*0.0; bangle_sigma=impact*0.0; \
|
|---|
| 45 | bangle_L1_qual=impact*0.0; bangle_L2_qual=impact*0.0; bangle_qual=impact*0.0"
|
|---|
| 46 | ncap2 -s "$cmd" $OFILE1 -O temp.nc
|
|---|
| 47 | ncatted -O -a valid_range,impact,o,d,'-10000000., 10000000.' temp.nc
|
|---|
| 48 |
|
|---|
| 49 | EXEC=/data/local/idculv/ropp_test_folder_02042019/gfortran/ropp_io/tools/ropp2ropp
|
|---|
| 50 |
|
|---|
| 51 | OFILE2=`echo $IFILE0 |sed -es/'.nc'/'_2.1.nc'/`
|
|---|
| 52 |
|
|---|
| 53 | echo "$EXEC temp.nc -o $OFILE2 -d -p $THINFILE >> $TFILE 2>&1"
|
|---|
| 54 | $EXEC temp.nc -o $OFILE2 -d -p $THINFILE >> $TFILE 2>&1
|
|---|
| 55 |
|
|---|
| 56 |
|
|---|
| 57 | # Now pass through invert to generate refracs and bangles
|
|---|
| 58 | EXEC=/data/local/idculv/ropp_test_folder_02042019/gfortran/ropp_pp/tools/ropp_pp_invert_tool
|
|---|
| 59 |
|
|---|
| 60 | OFILE=`echo $IFILE0 |sed -es/'.nc'/'.1.nc'/`
|
|---|
| 61 |
|
|---|
| 62 | echo "$EXEC $OFILE2 -m MSIS -o $OFILE >> $TFILE 2>&1"
|
|---|
| 63 | $EXEC $OFILE2 -m MSIS -o $OFILE >> $TFILE 2>&1
|
|---|
| 64 |
|
|---|
| 65 |
|
|---|
| 66 |
|
|---|
| 67 | #fi # bypass the control
|
|---|
