Changeset 2033 in lmdz_wrf for trunk

Timestamp:

Aug 7, 2018, 1:34:32 AM (6 years ago)

Author:

lfita

Message:

Adding `uavaFROMobswswd', computing woind components from observations (wind origin)

Location:

trunk/tools

Files:

• diagnostics.inf (modified) (1 diff)
• diagnostics.py (modified) (1 diff)
• obs-sim_Comparison.bash (modified) (3 diffs)

trunk/tools/diagnostics.inf

@@ -42 +42 @@
 va, WRFva, U@V@SINALPHA@COSALPHA
 uavaFwswd, uavaFROMwswd, ws@wd
+uavaFOwswd, uavaFROMobswswd, ws@wd
 uavaz, WRFzwind, U@V@WRFz@U10@V10@SINALPHA@COSALPHA@z=100.
 uavaz, WRFzwind_log, U@V@WRFz@U10@V10@SINALPHA@COSALPHA@z=100.

trunk/tools/diagnostics.py

@@ -911 +911 @@
         ncvar.insert_variable(ncobj, 'va', va, dnamesvar, dvnamesvar, newnc)
 
+# ua va from obs ws wd (deg)
+    elif diagn == 'uavaFROMobswswd':
+            
+        var0 = ncobj.variables[depvars[0]][:]
+        var1 = ncobj.variables[depvars[1]][:]
+
+        ua = var0*np.cos((var1+180.)*np.pi/180.)
+        va = var0*np.sin((var1+180.)*np.pi/180.)
+
+        dnamesvar = ncobj.variables[depvars[0]].dimensions
+        dvnamesvar = ncvar.var_dim_dimv(dnamesvar,dnames,dvnames)
+
+        ncvar.insert_variable(ncobj, 'ua', ua, dnamesvar, dvnamesvar, newnc)
+        ncvar.insert_variable(ncobj, 'va', va, dnamesvar, dvnamesvar, newnc)
 
 # WRFbils fom WRF as HFX + LH

trunk/tools/obs-sim_Comparison.bash

@@ -913 +913 @@
     figheader=`echo ${valuesfig} | tr '%' ' ' | awk '{print $9}'`
 
-    Nfiles=`echo ${labs} | tr ':' ' ' | wc -w | awk '{print $1-1}'`
+    Nfiles=`echo ${labs} | tr ',' ' ' | wc -w | awk '{print $1-1}'`
     frqm='1'
     iff=1
-    while test ${iff} -le ${Nfiles}; do frqm=${frqm}'@10'; done
-
-    values='linepoint;multicol;pres;-;'${lmark}':'${frqm}';spectral;auto:'${labs}':'
+    while test ${iff} -le ${Nfiles}; do frqm=${frqm}'@10'; iff=`expr ${iff} + 1`; done
+
+    values='linepoint;multicol;pres;-;'${lmark}';'${frqm}';spectral;auto:'${labs}':'
     values=${values}'WindRose!obs!,!sim!comparison!at!'${stid}'!on!'${timeS}'!UTC:'
     values=${values}${kfig}':auto:''True'
@@ -1268 +1268 @@
   if ${fscratch}; then rm ${ouavasndf}; fi
   if test ! -f ${ouavasndf}; then
-    python ${pyHOME}/diagnostics.py -f ${sndf} -d 'time@time,pres@pres' -v 'uavaFROMwswd|ws@wd'
+    python ${pyHOME}/diagnostics.py -f ${sndf} -d 'time@time,pres@pres'              \
+      -v 'uavaFROMobswswd|ws@wd'
     if test $? -ne 0; then
       echo ${errmsg}
       echo "  python failed!!"
-      echo python ${pyHOME}/diagnostics.py -f ${sndf} -d 'time@time,pres@pres' -v 'uavaFROMwswd|ws@wd'
+      echo python ${pyHOME}/diagnostics.py -f ${sndf} -d 'time@time,pres@pres'       \
+        -v 'uavaFROMobswswd|ws@wd'
       exit
     fi
     echo " * Getting ua,va winds from station "${stid}": " >> ${ofileins}
     echo python ${pyHOME}/diagnostics.py -f ${sndf} -d "'"time@time,pres@pres"'"     \
-      -v "'"uavaFROMwswd|ws@wd"'" >> ${ofileins}
+      -v "'"uavaFROMobswswd|ws@wd"'" >> ${ofileins}
     echo " " >> ${ofileins}    
     mv diagnostics.nc ${ouavasndf}
@@ -1549 +1551 @@
 
     multi_WindRose_plot ${fivals}
-    exit
 
     it=`expr ${it} + 1`