Changeset 1767 in lmdz_wrf

Timestamp:

Jan 26, 2018, 5:40:21 PM (7 years ago)

Author:

lfita

Message:

Working version of `draw_lines_time' with slicing
Introducing mm/day' and mm/month' to the `units_lunits'

Location:

trunk/tools

Files:

• drawing.py (modified) (5 diffs)
• drawing_tools.py (modified) (1 diff)

trunk/tools/drawing.py

@@ -32 +32 @@
 ## e.g. # $ drawing.py -f WRF/current/hurs_wrfout_tturb_xmean_last.nc,WRF/micro1/hurs_wrfout_tturb_xmean_last.nc,WRF/micro2/hurs_wrfout_tturb_xmean_last.nc,WRF_LMDZ/AR40/hurs_wrfout_tturb_xmean_last.nc,WRF_LMDZ/NPv31/hurs_wrfout_tturb_xmean_last.nc,LMDZ/AR40/hurs_reproj-histins_tturb_xmean_last.nc,LMDZ/NPv31/hurs_reproj-histins_tturb_xmean_last.nc -o draw_lines -S 'lat:x:lat ($degrees\ North$):auto:auto:wcurr,wmp1,wmp2,wlmdza,wlmdzb,lmdza,lmdzb:hurs:all model-experiments meridional hurs$_{[tturb\ xmean\ last]}$:0|auto:None:-:,:2.:2.:all:Nlines_hurs_wrfout_tturb-xmean-last:png:True' -v hursturbmean
 ## e.g. # drawing.py -o draw_lines_time -f wrfout_d01_2001-11-11_00:00:00_west_east_B20-E20-I1_south_north_B20-E20-I1.nc,wrfout_d01_2001-11-11_00:00:00_west_east_B25-E25-I1_south_north_B25-E25-I1.nc,wrfout_d01_2001-11-11_00:00:00_west_east_B35-E35-I1_south_north_B35-E35-I1.nc -S 'WRFtime;y;time ([DD]${[HH]}$);auto;we=20$\times$sn=20,we=25$\times$sn=25,we=35$\times$sn=35;tas;tas|evolution|at|3|different|grid|points;None;time|hours!since!1949-12-01_00:00:00|exct,3,h|%d$^{%H}$;0|12;pdf;-;r,g,b;.;2.;2.;all;-1;True' -v T2
+## e.g. # drawing.py -o draw_lines_time -f 'TRMM_1999-2010_prec_mm-year_hydyearsum_sellonlatbox2005_fldmean.nc%time|-1;lon|0;lat|0,CRUNCEP_1999-2010_prec_mm-year_hydyearsum_sellonlatbox2005_fldmean.nc%time|-1;lon|0;lat|0,WFDEIgpcc_1999-2010_prec_mm-year_hydyearsum_sellonlatbox2005_fldmean.nc%time|-1;lon|0;lat|0' -S 'time;y;time;auto;TRMM,CRUNCEP,WFDEIgpcc;prec;prec|hydrological|year;auto,auto;time|exct,1,y|%Y;0|auto;png;-;k,b,g;.;2;2;all;-1;True' -v prec
 ## e.g. # drawing.py -o draw_Neighbourghood_evol -S 'vas:Time|-1|WRFtime,south_north|44|XLAT,west_east|88|XLONG:south_north,west_east:5:auto:time|($[DD]^{[HH]}$),time|($[DD]^{[HH]}$):exct,1,h|exct,3,h:$%d^{%H}$,$%d^{%H}$:5|pts|neighbourghood|temporal|evolution|on|2001|Nov.|at|x=88,|y=44:0.0,20.:rainbow,auto,auto:auto:png:vas_Neigh_evol:True' -f ~/PY/wrfout_d01_2001-11-11_00:00:00 -v V10
 ## e.g. # drawing.py -o draw_points -S 'SuperStorm/tslist.dat,#,3,2,1:SuperStorm|sfc|stations:auto:cyl,i:labelled,10,r:auto:None:0:png:stations_loc:True' -f 'geo_em.d02.nc,XLONG_M,XLAT_M,HGT_M,Time|0@west_east|30;180;1@south_north|175;255;1,height,0.,1500.,terrain,auto,auto,m'
@@ -2372 +2373 @@
           [vmin],[vmax]: minimum and maximum values where [vmNN] can also be:
             'auto': the computed minimumm or maximum of the values  
-        [timevals]: [timen]|[units]|[kind]|[tfmt] time labels characteristics
+        [timevals]: [timen]|[kind]|[tfmt] time labels characteristics
            [timen]; name of the time variable
-           [units]; units string according to CF conventions ([tunits] since 
-             [YYYY]-[MM]-[DD] [[HH]:[MI]:[SS]], '!' for spaces) this is used to 
-             re-compute the time values to be the same reference for all the lines
            [kind]; kind of output
              'Nval': according to a given number of values as 'Nval',[Nval]
@@ -2409 +2407 @@
         [close]: Whether figure should be finished or not
       varnames= ',' list of variable names to plot (assuming only 1 variable per file
-        consecutively)
+        consecutively). With a single value will be used for all files
       values= 'time;y;time ([DD]${[HH]}$);32x32;$wss^{*}$;wss Taylor's turbulence term;time|hours!since!1949-12-01_00:00:00;exct,12,h|%d$^{%H}$;2;pdf'
     """
@@ -2518 +2516 @@
 
     timename = timevals.split('|')[0]
-    timeunit = timevals.split('|')[1].replace('!',' ')
-    timekind = timevals.split('|')[2]
-    timefmt = timevals.split('|')[3]
+    timekind = timevals.split('|')[1]
+    timefmt = timevals.split('|')[2]
 
     if rangevals == 'None':
@@ -2645 +2642 @@
             reftvals = gen.coincident_CFtimes(timevals0, timeunit, tunits0)
         else:
+            timeunit = tunits0
             reftvals = timevals0
 

trunk/tools/drawing_tools.py

@@ -1953 +1953 @@
         elif u == 'm/s': lu='$ms^{-1}$'
         elif u == 'mmh-3': lu='$mmh^{-3}$'
-        elif u == 'mm/d': lu='$mmd^{-1}$'
-        elif u == 'mmd-1': lu='$mmd^{-1}$'
+        elif u == 'mm/d': lu='$mm\ d^{-1}$'
+        elif u == 'mmd-1': lu='$mm\ d^{-1}$'
+        elif u == 'mm/day': lu='$mm\ day^{-1}$'
+        elif u == 'mmday-1': lu='$mm\ day^{-1}$'
+        elif u == 'mmmon-1': lu='$mm\ mon^{-1}$'
+        elif u == 'mm mon-1': lu='$mm\ mon^{-1}$'
+        elif u == 'mmmonth-1': lu='$mm\ month^{-1}$'
+        elif u == 'mm month-1': lu='$mm\ month^{-1}$'
+        elif u == 'mm yr-1': lu='$mm\ yr^{-1}$'
+        elif u == 'mmyear-1': lu='$mm\ year^{-1}$'
+        elif u == 'mm year-1': lu='$mm\ year^{-1}$'
         elif u == 'ms-1': lu='$ms^{-1}$'
         elif u == 'm s-1': lu='$ms^{-1}$'