Changeset 915 in lmdz_wrf for trunk/tools

Timestamp:

Jun 20, 2016, 8:08:32 PM (8 years ago)

Author:

lfita

Message:

Working until `pinterp'

• A new function `cleaning_varsfile': Function to keep only a list of variables from a file
• If a variable can be computed directly from the model of via a 'diagnostic' the later one is preferred

File:

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

trunk/tools/model_graphics.bash

@@ -144 +144 @@
       Spyout=`echo ${pyout} | tr '\n' '#' | tr ' ' '!'`
       ferrmsgF ${pyn} ${fname} "python!'varaddattr'!'units'!failed#"${Spyout} ${wrff}
+    fi
+  done
+}
+
+function cleaning_varsfile() {
+# Function to keep only a list of variables from a file
+#   filen= file to clean
+#   keepvars= ':' separated list of variables to keep
+  fname='cleaning_varsfile'
+
+  filen=$1
+  keepvars=$2
+
+  fvars=`${pyHOME}/nc_var.py -o ivars -f ${filen} | grep allvars | tr '=' ' ' |      \
+    awk '{print $3}'`
+  fvs=`echo ${fvars} | tr ':' ' '`
+
+  for fv in ${fvs}; do
+    if ! $(isInlist ${keepvars} ${fv}); then
+      pyout=`python ${pyHOME}/nc_var.py -o varrm -v ${fv} -f ${filen}`
+      pyn=$?
+      Spyout=`echo ${pyout} | tr '\n' '#' | tr ' ' '!'`
+      ferrmsg ${pyn} ${fname} "python!'varrm'!'${fv}'!failed#"${Spyout}
     fi
   done
@@ -220 +243 @@
 
     if test ! -f ${filen} && test ! -f ${odir}/${CFvarn}_${fhead}.nc; then
-      if test ! ${modvar} = 'None'; then
+# Since model direct values are retrieved from `variables_valules.dat' which was initially coincived
+#   as a way to only give variable attributes, range and color bars, if a variable has a diagnostic
+#   way to be computed, the later one will be preferred
+
+      if test ! ${modvar} = 'None' && test ${diagvar} = 'None'; then
         # model variable
         values=${modvar}',0,-1,-1'
@@ -251 +278 @@
           idiagv=`expr ${idiagv} + 1`
         done
-
         pyout=`python ${pyHOME}/diagnostics.py -f ${cf} -d ${dims} -v ${diagn}'|'${diagc}`
         pyn=$?
@@ -276 +302 @@
       fi
     fi
-
 
     ifile=`expr ${ifile} + 1`
@@ -463 +488 @@
     esac
   fi
+
+  cleaning_varsfile ${filen} ${CFvarn}':lon:lat:pres:time'
 #  exit
 # LLUIS