Changeset 2049 in lmdz_wrf

Timestamp:

Aug 7, 2018, 7:05:26 PM (7 years ago)

Author:

lfita

Message:

Adding case when there are no diagnostics in sfc map

File:

• trunk/tools/get_WRFsounding-surface_data.bash (modified) (2 diffs)

trunk/tools/get_WRFsounding-surface_data.bash

@@ -314 +314 @@
 
   # Computing sfc diagnostics
-  sfcmapdiagd='Time@WRFtime:south_north@XLAT:west_east@XLONG'
-  diagns=`echo ${sfcdiagnames} | tr ':' ' '`
-
-  # getting that variables only required for the diagnostics
-  vardiags=''
-  idiag=1
-  for diagn in ${diagns}; do
-    diagv=`cat $pyHOME/diagnostics.inf | grep WRF | grep ${diagn}',' | tr ',' ' ' |  \
-      awk '{print $2"|"$3}'`
-    if test ${idiag} -eq 1; then
-      sfcdiagvals=${diagv}
-      vardiags=`echo ${diagv} | tr '|' ' ' | awk '{print $2}' | tr '@' ':'`
-    else
-      sfcdiagvals=${sfcdiagvals}','${diagv}
-      vardiag=`echo ${diagv} | tr '|' ' ' | awk '{print $2}' | tr '@' ' '`
-      for vd in ${vardiag}; do
-        if ! $(isin_list ${vardiags} ${vd}); then vardiags=${vardiags}':'${vd}; fi
+  if test ! ${sfcdiagnames} = 'None'; then
+    sfcmapdiagd='Time@WRFtime:south_north@XLAT:west_east@XLONG'
+    diagns=`echo ${sfcdiagnames} | tr ':' ' '`
+
+    # getting that variables only required for the diagnostics
+    vardiags=''
+    idiag=1
+    for diagn in ${diagns}; do
+      diagv=`cat $pyHOME/diagnostics.inf | grep WRF | grep ${diagn}',' | tr ',' ' ' |  \
+        awk '{print $2"|"$3}'`
+      if test ${idiag} -eq 1; then
+        sfcdiagvals=${diagv}
+        vardiags=`echo ${diagv} | tr '|' ' ' | awk '{print $2}' | tr '@' ':'`
+      else
+        sfcdiagvals=${sfcdiagvals}','${diagv}
+        vardiag=`echo ${diagv} | tr '|' ' ' | awk '{print $2}' | tr '@' ' '`
+        for vd in ${vardiag}; do
+          if ! $(isin_list ${vardiags} ${vd}); then vardiags=${vardiags}':'${vd}; fi
+        done
+      fi
+      idiag=`expr ${idiag} + 1`
+    done
+
+    # Inlcuding non-diagnostics
+    if test ! ${nondiagvs} = 'None'; then
+      nvns=`echo ${nondiagvs} | tr ':' ' '`
+      for vn in ${nvns}; do
+        if ! $(isin_list ${vardiags} ${vn}); then vardiags=${vardiags}':'${vn}; fi
       done
     fi
-    idiag=`expr ${idiag} + 1`
-  done
-
-  # Inlcuding non-diagnostics
-  if test ! ${nondiagvs} = 'None'; then
-    nvns=`echo ${nondiagvs} | tr ':' ' '`
-    for vn in ${nvns}; do
-      if ! $(isin_list ${vardiags} ${vn}); then vardiags=${vardiags}':'${vn}; fi
-    done
-  fi
+  else
+    vardiags=${nondiagvs}
+  fi
+
   vardiag=`echo ${vardiags} | tr ':' ','`
   
@@ -369 +374 @@
 
   if ${fscratch}; then rm ${simsfcmapdiagsf}; fi
-  if test ! -f ${simsfcmapdiagsf}; then
-    python ${pyHOME}/diagnostics.py -f ${simsfcmapf} -d ${sfcmapdiagd}               \
-      -v ${sfcdiagvals}
-    if test $? -ne 0; then
-      echo ${errmsg}
-      echo "  python failed!!"
+    if test ! ${sfcdiagnames} = 'None'; then
+      if test ! -f ${simsfcmapdiagsf}; then
+      python ${pyHOME}/diagnostics.py -f ${simsfcmapf} -d ${sfcmapdiagd}             \
+        -v ${sfcdiagvals}
+      if test $? -ne 0; then
+        echo ${errmsg}
+        echo "  python failed!!"
+        echo python ${pyHOME}/diagnostics.py -f ${simsfcmapf} -d "'"${sfcmapdiagd}"'"\
+          -v "'"${sfcdiagvals}"'"
+        exit
+      fi
+      echo "* diagnostics for sfc 2D map: " >> ${ofileins}
       echo python ${pyHOME}/diagnostics.py -f ${simsfcmapf} -d "'"${sfcmapdiagd}"'"  \
-        -v "'"${sfcdiagvals}"'"
-      exit
-    fi
-    echo "* diagnostics for sfc 2D map: " >> ${ofileins}
-    echo python ${pyHOME}/diagnostics.py -f ${simsfcmapf} -d "'"${sfcmapdiagd}"'"    \
-      -v "'"${sfcdiagvals}"'" >> ${ofileins}
-    echo " " >> ${ofileins}
-    mv diagnostics.nc ${simsfcmapdiagsf}
-    if $(isin_list ${sfcdiagnames} wds); then
-      #To rad
-      python ${pyHOME}/nc_var.py -o valmod -S mulc,0.0174532925199                   \
-        -f ${simsfcmapdiagsf} -v wds
-      python ${pyHOME}/nc_var.py -o varaddattr -S 'units|rad' -f ${simsfcmapdiagsf}  \
-        -v wds
-    fi
-
+        -v "'"${sfcdiagvals}"'" >> ${ofileins}
+      echo " " >> ${ofileins}
+      mv diagnostics.nc ${simsfcmapdiagsf}
+      if $(isin_list ${sfcdiagnames} wds); then
+        #To rad
+        python ${pyHOME}/nc_var.py -o valmod -S mulc,0.0174532925199                 \
+          -f ${simsfcmapdiagsf} -v wds
+        python ${pyHOME}/nc_var.py -o varaddattr -S 'units|rad' -f ${simsfcmapdiagsf}\
+          -v wds
+      fi
+    fi
+    # If there was not diagnostics
+    if test ! -f ${simsfcmapdiagsf}; then
+      dddv='Time@time'
+      selvars='XLONG,XLAT,time'
+      python ${pyHOME}/nc_var.py -o selvar -S ${dddv} -f ${simsfcmapf} -v ${selvars}
+      if test $? -ne 0; then
+        echo ${errmsg}
+        echo "  python failed!!"
+        echo python ${pyHOME}/nc_var.py -o selvar -S "'"${ddv}"'" -f ${simsfmapf}    \
+          -v ${vars}
+        exit
+        rm ${simsfcmapf}
+      fi
+      echo "* creation of diagnostics file for 2D map:" >>  ${ofileins}
+      echo python ${pyHOME}/nc_var.py -o selvar -S "'"${dddv}"'" -f ${simsfcmapf}    \
+          -v ${selvars} >> ${ofileins}
+      echo " " >> ${ofileins}
+      mv selvar_new.nc ${simsfcmapdiagsf}
+    fi 
     #Adding non-diagnostic variables
     if test ! ${nondiagvs} = 'None'; then