Changeset 2012 in lmdz_wrf

Timestamp:

Aug 2, 2018, 3:42:26 PM (7 years ago)

Author:

lfita

Message:

Working version with the generic definition of output file (preventing generalization...)

File:

• trunk/tools/obs-sim_Comparison.bash (modified) (6 diffs)

trunk/tools/obs-sim_Comparison.bash

@@ -789 +789 @@
 
 # sounding data
-isnd=1
-for sndf in ${sndfiles}; do
-  stid=`python $pyHOME/nc_var.py -o grattr -f ${sndf} -S Station_number`
-  sndlon=`ncdump -h ${sndf} | grep Station_longitude | awk '{print $3}'`
-  sndlat=`ncdump -h ${sndf} | grep Station_latitude | awk '{print $3}'`
-  sndheight=`ncdump -h ${sndf} | grep Station_elevation | awk '{print $3}'`
-
-  if test ${isnd} -eq 1; then 
-    sndstids=${stid}'@'${sndf}
-    echo "snd,"${stid}","${sndlon}","${sndlat}","${sndheight} > ${wdir}/stations.inf
+iost=1
+obstype='sounding'
+prn='snd'
+listfiles=${sndfiles}
+wrfdiags=${wrfsnddiags}
+simdiagd=${simsnddiagd}
+for ostf in ${listfiles}; do
+  stid=`python $pyHOME/nc_var.py -o grattr -f ${ostf} -S Station_number`
+  ostlon=`ncdump -h ${ostf} | grep Station_longitude | awk '{print $3}'`
+  ostlat=`ncdump -h ${ostf} | grep Station_latitude | awk '{print $3}'`
+  ostheight=`ncdump -h ${ostf} | grep Station_elevation | awk '{print $3}'`
+
+  if test ${iost} -eq 1; then 
+    oststids=${stid}'@'${ostf}
+    echo ${prn}","${stid}","${ostlon}","${ostlat}","${ostheight} > ${wdir}/stations.inf
   else 
-    sndstids=${sndstids}':'${stid}'@'${sndf}
-    echo "snd,"${stid}","${sndlon}","${sndlat}","${sndheight} >> ${wdir}/stations.inf
+    oststids=${oststids}':'${stid}'@'${ostf}
+    echo ${prn}","${stid}","${ostlon}","${ostlat}","${ostheight} >> ${wdir}/stations.inf
   fi
 
@@ -818 +823 @@
         awk '{print $3"_"$4}'`
 
-      simsndptf=${odir}'/simout_vars_sndpt_'${stid}'_'${expl}'_'${fdate}'.nc'
-      simsnddiagsf=${odir}'/simout_snddiags_'${stid}'_'${expl}'_'${fdate}'.nc'
-
-      fvals=${sndlon}'#'${sndlat}'#'${simf}'#XLONG#XLAT#'${wrfsnddiags}'#'
-      fvals=${fvals}${simsnddiagd}
+      simostptf=${odir}'/simout_vars_'${prn}'pt_'${stid}'_'${expl}'_'${fdate}'.nc'
+      simostdiagsf=${odir}'/simout_'${prn}'diags_'${stid}'_'${expl}'_'${fdate}'.nc'
+
+      fvals=${ostlon}'#'${ostlat}'#'${simf}'#XLONG#XLAT#'${wrfdiags}'#'
+      fvals=${fvals}${simdiagd}
       get_WRFsounding_data ${fvals}
 
@@ -829 +834 @@
 
     # Concatenating files
-    ofile=${odir}/simout_sndvars_${stid}_${expn}.nc
+    ofile=${odir}/simout_${prn}vars_${stid}_${expl}.nc
     if ${fscratch}; then rm ${ofile} >& /dev/null; fi
     if test ! -f ${ofile}; then
-      Nconfc=`ls -1 ${odir}/simout_snddiags_*${stid}*nc | wc -l | awk '{print $1}'`
+      Nconfc=`ls -1 ${odir}/simout_${prn}diags_*${stid}*${expl}*.nc | wc -l | awk '{print $1}'`
       if test ${Nconfc} -gt 1; then 
         values=${odir}',Time,time'
-        HMT='simout_snddiags_,'${stid}'_'${expn},'nc'
+        HMT='simout_'${prn}'diags_,'${stid}'_'${expn},'nc'
         python ${pyHOME}/nc_var.py -o netcdf_fold_concatenation_HMT -S ${values}      \
          -f ${HMT} -v all
@@ -845 +850 @@
           exit
         fi
-        echo " * Concatenating sounding data at point: "${sndlon}", "${sndlat} >> ${ofileins}
+        echo " * Concatenating "${obstype}" data at point: "${ostlon}", "${ostlat} >> ${ofileins}
         echo python  ${pyHOME}/nc_var.py -o netcdf_fold_concatenation_HMT -S \
           "'"${values}"'" -f ${HMT} -v all >> ${ofileins}
@@ -851 +856 @@
         mv netcdf_fold_concatenated_HMT.nc ${ofile}
       else
-        ffilen=`ls -1 ${odir}/simout_snddiags_*${stid}*nc`
+        ffilen=`ls -1 ${odir}/simout_${prn}diags_*${stid}*${expl}*.nc`
         cp ${ffilen} ${ofile}
       fi
@@ -859 +864 @@
   done
 
-  isnd=`expr ${isnd} + 1`
+  iost=`expr ${iost} + 1`
 done
 exit