Changeset 912 in lmdz_wrf for trunk/tools

Timestamp:

Jun 20, 2016, 1:41:24 PM (8 years ago)

Author:

lfita

Message:

Simly working `pinterp'

File:

• trunk/tools/model_graphics.bash (modified) (6 diffs)

trunk/tools/model_graphics.bash

@@ -68 +68 @@
 #  msg: message to write ('!' for spaces, '#' for end of line)
 #  FileName: name of the file to remove in case of error
-  comdexit=$1
+  comdexit=`expr $1 + 0`
   ref=$2
   msg=$3
@@ -75 +75 @@
   if test ${comdexit} -ne 0; then
     echo "ERROR -- error -- ERROR -- error"
-    echo "  "${ref}": "$(echo ${msg} | tr '!' ' ' | tr '#' '\n')" !!"
+    echo "  "${ref}": "$(echo ${msg} | tr '!' ' ' | tr '\#' '\n')" !!"
     rm ${FileName} >& /dev/null
     exit
@@ -205 +205 @@
     ifS=`printf "%0${Ntotfiles}d" ${ifile}`
 
-    # Computing variable
+# Computing variable
+    # Changing file head when it is a pressure-interpolated variable
     if test ${vark} == 'pinterp'; then
       fhead=${headerf}p
@@ -351 +352 @@
   cd ${odir}
 
-  filen=${CFvarn}_${headerf}_${vark}.nc
+  if test ${vark} == 'pinterp'; then
+    fhead=${headerf}p
+  else
+    fhead=${headerf}
+  fi
+  filen=${CFvarn}_${fhead}_${vark}.nc
   if ${scratch}; then rm ${filen} >& /dev/null; fi
 
@@ -394 +400 @@
       'pinterp')
         vals=${plevels}',1,1'
-        python $pyHOME/nc_var.py -o pinterp -f ${cfiles} -S ${vals} -v ${CFvarn}
+        echo "python $pyHOME/nc_var.py -o pinterp -f ${cfiles} -S ${vals} -v ${CFvarn}"
+        pyout=`python $pyHOME/nc_var.py -o pinterp -f ${cfiles} -S ${vals}           \
+          -v ${CFvarn}`
+        pyn=$?
+        Spyout=`echo ${pyout} | tr '\n' '#' | tr ' ' '!'`
+        ferrmsg ${pyn} ${fname} "python!'pinterp'!failed#"${Spyout}
+        cp pinterp.nc ${filen}
+
+        # adding CF lon,lat,time in WRF files
+        if test ${headerf:0:3} = 'wrf'; then
+          WRF_toCF ${filen} ${vdnx} ${vdny}
+        fi
       ;;
       # temporal mean values
@@ -731 +748 @@
 
         if test ! ${vark} = 'diff'; then
-          compute_statistics ${iwdir} ${CFv}_${fileh}.nc ${owdir} ${cvar}            \
+          if test ${vark} = 'pinterp'; then
+            fhead=${fileh}'p'
+          else
+            fhead=${fileh}
+          fi
+          compute_statistics ${iwdir} ${CFv}_${fhead}.nc ${owdir} ${cvar}            \
             ${moddims} ${modvdims} ${scratch}
           isc=`expr ${isc} + 1`