Changeset 908 in lmdz_wrf for trunk/tools

Timestamp:

Jun 19, 2016, 10:05:46 PM (8 years ago)

Author:

lfita

Message:

Adding vertical interpolation

Location:

trunk/tools

Files:

• model_graphics.bash (modified) (10 diffs)
• model_graphics_template.dat (modified) (4 diffs)

trunk/tools/model_graphics.bash

@@ -162 +162 @@
 #        lsec: longitudinal section (longitudinal value must be given, [var]@[lat])
 #        lmean: longitudinal mean values
+#        pinterp: pressure interpolation (to the given $plevels)
 #        tmean: temporal mean values
 #        zsum: vertical aggregated values
@@ -205 +206 @@
 
     # Computing variable
-    filen=${odir}/${CFvarn}_${headerf}_${ifile}-${Ntotfiles}.nc
+    if test ${vark} == 'pinterp'; then
+      fhead=${headerf}p
+    else
+      fhead=${headerf}
+    fi
+
+    filen=${odir}/${CFvarn}_${fhead}_${ifile}-${Ntotfiles}.nc
     if ${scratch}; then 
       rm ${filen} >& /dev/null
-      rm ${odir}/${CFvarn}_${headerf}.nc >& /dev/null
+      rm ${odir}/${CFvarn}_${fhead}.nc >& /dev/null
     fi
 
-    if test ! -f ${filen} && test ! -f ${odir}/${CFvarn}_${headerf}.nc; then
+    if test ! -f ${filen} && test ! -f ${odir}/${CFvarn}_${fhead}.nc; then
       if test ! ${modvar} = 'None'; then
         # model variable
@@ -255 +262 @@
         WRF_toCF ${filen} ${vdnx} ${vdny}
       fi
+      # Attaching necessary variables for the pressure interpolation
+      if test ${vark} == 'pinterp'; then
+        requiredinterpvars='P:PB:PSFC:PH:PHB:HGT:T:QVAPOR:'
+        rqvs=`echo ${requiredinterpvars} | tr ':' ' '`
+        echo "  "${fname}": adding variables: "${rqvs}" to allow pressure interpolation"
+        for rqv in ${rqvs}; do
+          pyout=`${pyHOME}/nc_var.py -o fvaradd -S ${cf},${rqv} -f ${filen}`
+          pyn=$?
+          Spyout=`echo ${pyout} | tr '\n' '#' | tr ' ' '!'`
+          ferrmsgF ${pyn} ${fname} "python!'fvaradd'!failed#"${Spyout} ${filen}
+        done
+      fi
     fi
+
 
     ifile=`expr ${ifile} + 1`
@@ -262 +282 @@
 
   # Joining variable files
-  filen=${odir}/${CFvarn}_${headerf}.nc
+  filen=${odir}/${CFvarn}_${fhead}.nc
   if ${scratch}; then rm ${filen} >& /dev/null; fi
 
   if test ! -f ${filen}; then
-    pyout=`python ${pyHOME}/nc_var.py -f ${CFvarn}_${headerf}_                       \
-      -o netcdf_fold_concatenation -S ./,time -v all`
+    pyout=`python ${pyHOME}/nc_var.py -f ${CFvarn}'_'${fhead}'_,-,.nc'               \
+      -o netcdf_fold_concatenation_HMT -S ./,time -v all`
     pyn=$?
     Spyout=`echo ${pyout} | tr '\n' '#' | tr ' ' '!'`
-    ferrmsg ${pyn} ${fname} "python!'netcdf_fold_concatenation'!failed#"${Spyout}
-
-    mv netcdf_fold_concatenated.nc ${filen}
+    ferrmsg ${pyn} ${fname} "python!'netcdf_fold_concatenation_HMT'!failed#"${Spyout}
+    mv netcdf_fold_concatenated_HMT.nc ${filen}
     if test -f ${filen}; then
-      rm ${CFvarn}_${headerf}_*
+      rm ${CFvarn}_${fhead}_*-*.nc
     fi
   fi
@@ -294 +313 @@
 #        lmean: longitudinal mean values
 #        lsec: longitudinal section (longitudinal value must be given, [var]@[lat])
+#        pinterp: pressure interpolation (to the given $plevels)
 #        tmean: temporal mean values
 #        xmean: x-axis mean values
@@ -371 +391 @@
         exit
       ;;
+      # pinterp: pressure interpolation (to the given $plevels)
+      'pinterp')
+        vals=${plevels}',1,1'
+        python $pyHOME/nc_var.py -o pinterp -f ${cfiles} -S ${vals} -v ${CFvarn}
+      ;;
       # temporal mean values
       'tmean')
@@ -441 +466 @@
 else
   scratch=false
+fi
+
+if test ${debug} = 'true'; then
+  dbg=true
+else
+  dbg=false
 fi
 
@@ -574 +605 @@
           'lsec')
             varvks=${varlsec}
+          ;;
+          'pinterp')
+            varvks=${varpinterp}
           ;;
           'tmean')
@@ -689 +723 @@
       modv=`echo ${cvar} | tr '|' ' ' | awk '{print $4}'`
       diagv=`echo ${cvar} | tr '|' ' ' | awk '{print $5}'`
-      echo "        "${CFv}"; "${modv}" "${diagv}
+      if ${dbg}; then echo "        "${CFv}"; "${modv}" "${diagv}; fi
 
       if test ! ${modv} = 'None' || test ! ${diagv} = 'None'; then
@@ -705 +739 @@
         fi
       else
-        echo "        not computing '"${CFv}"' for model '"${mod}"' !!"
+        if ${dbg}; then echo "        not '"${CFv}"' for model '"${mod}"' !!"; fi
       fi
 

trunk/tools/model_graphics_template.dat

@@ -20 +20 @@
 # Srcatch
 scratch=false
+
+# Debug
+debug=false
 
 # Folder with the files
@@ -54 +57 @@
 #   lmean: longitudinal mean values
 #   lsec: longitudinal section (longitudinal value must be given, [var]@[lat])
+#   pinterp: pressure interpolation (to the given $plevels, see below)
 #   tmean: temporal mean values
 #   xmean: x-axis mean values
 #   ymean: y-axis mean values
 #   zsum: vertical aggregated values
-varkinds = direct:tmean:xmean:diff
+varkinds = direct:tmean:xmean:diff:pinterp
 
 # variables (providing CF-name, script must provide equivalent variable for each model)
@@ -67 +71 @@
 varxmean = ta:ua:va:hur
 vardiff = ualmean:valmean:talmean:huslmean:uas:vas:ps:pr
+varpinterp = ta:ua:va:hur
 
+# Pressure levels
+#   Common vertical pressure levels to compute vertical differences
+plevels = 100000.:97500.:95000.:92500.:90000.:85000.:80000.:75000.:70000.:65000.:60000.:55000.:50000.:45000.:40000.:35000.:30000.:25000.:20000.:15000.:10000.:5000.:2500.:1000.:500.:250.
 # Kind of plots
 #   diffmap2Dsfc: 2D map of surface differences values of 1 variable
@@ -79 +87 @@
 diffmap2Dz = ualmean:valmean:talmean:huslmean
 
-# Pressure levels
-#   Common vertical pressure levels to compute vertical differences
-plevels = 100000.,97500.,95000.,92500.,90000.,85000.,80000.,75000.,70000.,65000.,60000.,55000.,50000.,45000.,40000.,35000.,30000.,25000.,20000.,15000.,10000.,5000.,2500.,1000.,500.,250.
-
 # Figures output 
 kindfig = pdf