Changeset 340 in lmdz_wrf for trunk

Timestamp:

Mar 1, 2015, 11:12:17 PM (10 years ago)

Author:

lfita

Message:

Adding check when simulated variable has dimensions not given as reference [X,Y,Z,T]

File:

• trunk/tools/validation_sim.py (modified) (6 diffs)

trunk/tools/validation_sim.py

@@ -3 +3 @@
 ## e.g. sfcEneAvigon # validation_sim.py -d X@west_east@None,Y@south_north@None,T@Time@time -D X@XLONG@longitude,Y@XLAT@latitude,T@time@time -k single-station -l 4.878773,43.915876,12. -o /home/lluis/DATA/obs/HyMeX/IOP15/sfcEnergyBalance_Avignon/OBSnetcdf.nc -s /home/lluis/PY/wrfout_d01_2012-10-18_00:00:00.tests -v HFX@H,LH@LE,GRDFLX@G
 ## e.g. AIREP # validation_sim.py -d X@west_east@lon2D,Y@south_north@lat2D,Z@bottom_top@z2D,T@Time@time -D X@XLONG@longitude,Y@XLAT@latitude,Z@WRFz@alti,T@time@time -k trajectory -o /home/lluis/DATA/obs/HyMeX/IOP15/AIREP/2012/10/AIREP_121018.nc -s /home/lluis/PY/wrfout_d01_2012-10-18_00:00:00.tests -v WRFt@t,WRFtd@td,WRFws@u,WRFwd@dd
-## e.g. ATRCore # validation_sim.py -d X@west_east@lon2D,Y@south_north@lat2D,Z@bottom_top@z2D,T@Time@CFtime -D X@XLONG@longitude,Y@XLAT@latitude,Z@WRFz@altitude,T@time@time -k trajectory -o /home/lluis/DATA/obs/HyMeX/IOP15/ATRCore/V3/ATR_1Hz-HYMEXBDD-SOP1-v3_20121018_as120051.nc -s /home/lluis/PY/wrfout_d01_2012-10-18_00:00:00.tests -v WRFT@air_temperature@subc@273.15
-## e.g. BAMED # validation_sim.py -d X@west_east@lon2D,Y@south_north@lat2D,Z@bottom_top@z2D,T@Time@CFtime -D X@XLONG@longitude,Y@XLAT@latitude,Z@WRFz@altitude,T@time@time -k trajectory -o /home/lluis/DATA/obs/HyMeX/IOP15/BAMED/BAMED_SOP1_B12_TOT5.nc -s /home/lluis/PY/wrfout_d01_2012-10-18_00:00:00.tests -v WRFT@tas_north
+## e.g. ATRCore # validation_sim.py -d X@west_east@lon2D,Y@south_north@lat2D,Z@bottom_top@z2D,T@Time@CFtime -D X@XLONG@longitude,Y@XLAT@latitude,Z@WRFz@altitude,T@time@time -k trajectory -o /home/lluis/DATA/obs/HyMeX/IOP15/ATRCore/V3/ATR_1Hz-HYMEXBDD-SOP1-v3_20121018_as120051.nc -s /home/lluis/PY/wrfout_d01_2012-10-18_00:00:00.tests -v WRFt@air_temperature@subc@273.15,WRFp@air_pressure,WRFrh@relative_humidity,WRFrh@relative_humidity_Rosemount,WRFwd@wind_from_direction,WRFws@wind_speed
+## e.g. BAMED # validation_sim.py -d X@west_east@lon2D,Y@south_north@lat2D,Z@bottom_top@z2D,T@Time@CFtime -D X@XLONG@longitude,Y@XLAT@latitude,Z@WRFz@altitude,T@time@time -k trajectory -o /home/lluis/DATA/obs/HyMeX/IOP15/BAMED/BAMED_SOP1_B12_TOT5.nc -s /home/lluis/PY/wrfout_d01_2012-10-18_00:00:00.tests -v WRFt@tas_north,WRFp@pressure,WRFrh@hus,U@uas,V@vas
 
 import numpy as np
@@ -796 +796 @@
     quit(-1)
 else:
+    simdims = {}
+    obsdims = {}
     print main +': couple of dimensions _______'
     dims = {}
@@ -807 +809 @@
             quit(-1)
         dims[dsecs[0]] = [dsecs[1], dsecs[2]]
+        simdims[dsecs[0]] = dsecs[1]
+        obsdims[dsecs[0]] = dsecs[2]
+
         print '  ',dsecs[0],':',dsecs[1],',',dsecs[2]
         
@@ -1104 +1109 @@
 for ivar in range(Nvars):
     simobsvalues = []
-# Values spatially around the point (+/- [Ngrid] points)
-    simobsSvalues = []
 
     varsimobs = valvars[ivar][0] + '_' + valvars[ivar][1]
+    print '  ' + varsimobs + '... .. .'
 
     if not oobs.variables.has_key(valvars[ivar][1]):
@@ -1125 +1129 @@
     else:
         ovsim = osim.variables[valvars[ivar][0]]
+
+    for idn in ovsim.dimensions:
+        if not searchInlist(simdims.values(),idn):
+            print errormsg
+            print '  ' + main + ": dimension '" + idn + "' of variable '" +          \
+              valvars[ivar][0] + "' not provided as reference coordinate [X,Y,Z,T] !!"
+            quit(-1)
 
     ovobs = oobs.variables[valvars[ivar][1]]
@@ -1190 +1201 @@
                     sliceS.append(slice(0,maxy-miny))
                     sliceS.append(slice(0,maxx-minx))
-                    
+
                     simobsSvalues[tuple(sliceS)] = slicevar
                     if ivar == 0: