Changeset 2028 in lmdz_wrf

Timestamp:

Aug 6, 2018, 4:23:09 PM (7 years ago)

Author:

lfita

Message:

Adding `inf_operSlist'
Improving `inf_operSlist' checking of sign preservation of interval

Location:

trunk/tools

Files:

• generic_tools.py (modified) (1 diff)
• obs-sim_Comparison.bash (modified) (2 diffs)

trunk/tools/generic_tools.py

@@ -13672 +13672 @@
         Nval = len(listv)
         dlist = np.array(listv[1:Nval]) - np.array(listv[0:Nval-1])
-        if not np.all(dlist > 0.) or np.all(dlist < 0.): 
+        if dlist[0] > 0. and not np.all(dlist > 0.): 
             print errormsg
-            print '  ' + fname + ": operation '" + oper + "' requires sorted list !!"
+            print '  ' + fname + ": operation '" + oper + "' requires sort list !!"
             print "    provided one '" + listS + "' is not!!"
+            if dlist[0] > 0.:
+                for iv in range(1,Nval-2): 
+                    if dlist[iv] < 0.: 
+                        print '    change of sign of interval at value:', iv,        \
+                      listv[iv], listv[iv+1]
+        if dlist[0] < 0. and not np.all(dlist < 0.): 
+            print errormsg
+            print '  ' + fname + ": operation '" + oper + "' requires sort list !!"
+            print "    provided one '" + listS + "' is not!!"
+            if dlist[0] < 0.:
+                for iv in range(1,Nval-2): 
+                    if dlist[iv] > 0.: 
+                        print '    change of sign of interval at value:', iv,        \
+                          listv[iv], listv[iv+1]
             quit(-1)
         sign = listv[1] - listv[0]

trunk/tools/obs-sim_Comparison.bash

@@ -1470 +1470 @@
         obsuavaf=${odir}/sounding_uava_${stid}.nc
         obspres=`python $pyHOME/nc_var.py -o varout -f ${obsuavaf} -S time:${it}     \
-          -v pres | awk '{print $2}' | tr '\n' ' '`
-        kobs=0
-        for obsp0 in ${obspres}; do
-          obsp=`echo ${obsp0} | awk '{printf ("%d",$1*1000)}'`
-          if test ${obsp} -gt ${nsimpres}; then kobs=`expr ${kobs} + 1`
-          else kobs=`expr ${kobs} + 1`; break; fi
-        done
+          -v pres | awk '{print $2}' | tr '\n' ':'`
+        vals=${obspres}',threshold,:,'${nsimpres}
+        kobs=`python ${pyHOME}/generic.py -o inf_operSlist -S ${vals}`
+        if test $? -ne 0; then
+            echo ${errormsg}
+            echo "  python failed !!"
+            echo python ${pyHOME}/generic.py -o inf_operSlist -S ${vals}
+            exit
+        fi
+        echo "      obs level of sim npres "${nsimpres}": "${kobs}
 
         mWRcomplabs='obs,'${expn}
@@ -1538 +1541 @@
 
     multi_WindRose_plot ${fivals}
+    exit
 
     it=`expr ${it} + 1`