Changeset 1573 in lmdz_wrf for trunk/tools

Timestamp:

May 12, 2017, 9:25:55 PM (8 years ago)

Author:

lfita

Message:

Adding use of 'auto'
Changing example

Location:

trunk/tools

Files:

• documentation/plotting/gallery.html (modified) (1 diff)
• documentation/plotting/time_lag.html (modified) (2 diffs)
• drawing_tools.py (modified) (3 diffs)

trunk/tools/documentation/plotting/gallery.html

@@ -38 +38 @@
       <A HREF="subbasin.html"><IMG SRC="figs/Caceres_subbasin.png" ALT="plot sub-basin from ORCHIDEE's routing" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>     
       <A HREF="Taylor.html"><IMG SRC="figs/Taylor.png" ALT="Taylor diagram" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
-      <A HREF="time_lag.html"><IMG SRC="figs/mon1_timelag_pr-evap_SACZ.png" ALT="Time lag" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
+      <A HREF="time_lag.html"><IMG SRC="figs/seas1_timelag_pr-evap_SESA_anom.png" ALT="Time lag" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="time_lag.html"><IMG SRC="figs/mon6_timelag_pr-evap_SACZ_shadline.png" ALT="Time lag" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>
       <A HREF="time_lag.html"><IMG SRC="figs/mon1_timelag_pr-evap_SACZ_repeatPeriod.png" ALT="Time lag" STYLE="float:left;width:200px;height:200px" TARGET="value"></IMG></A>

trunk/tools/documentation/plotting/time_lag.html

@@ -13 +13 @@
       Plot a time-lag figure with multiple sources (x, previous values; y, future values)
       <TABLE>
-        <TD><TR><IMG WIDTH="100%" SRC="figs/mon1_timelag_pr-evap_SACZ.png"></IMG></TR></TD>
+        <TD><TR><IMG WIDTH="100%" SRC="figs/seas1_timelag_pr-evap_SESA_anom.png"></IMG></TR></TD>
         <TD><TR><IMG WIDTH="100%" SRC="figs/mon6_timelag_pr-evap_SACZ_shadline.png"></IMG></TR></TD>
         <TD><TR><IMG WIDTH="100%" SRC="figs/mon1_timelag_pr-evap_SACZ_repeatPeriod.png"></IMG></TR></TD>
@@ -117 +117 @@
     </DIV>
     <DIV CLASS="valins">
-      * Plotting with cloud of points <BR>
-      &nbsp;&nbsp;$ python ${pyHOME}/drawing.py -o draw_time_lag -f 'Forcing_rainstats_OKstomate_SESA_norm_meanstd_stat2.nc;time_centered,-1;time_counter;rainnorm%DiPolo_evapstats_SESA_norm_meanstd_OKstomate.nc;time_counter,-1|stats,2;time_counter;evapnorm#Forcing_rainstats_OKstomate_SACZ_norm_meanstd_stat2.nc;time_centered,-1;time_counter;rainnorm%DiPolo_evapstats_SACZ_norm_meanstd_OKstomate.nc;time_counter,-1|stats,2;time_counter;evapnorm' -S 'auto;centered;pr;1;evap;1;-5,5;-2,2;year|seconds!since!1958-01-01!00:00:00|exct,5,y|%Y;SESA,SACZ;simplepts,1,#AA0000@#0000AA,x,2.;1-month!time-lag!at!SACZ!between!pr!&!evap!normalized!by!mean/stddev;0|12;png;mon1_timelag_pr-evap_SACZ;yes' <BR>
+      * Plotting with cloud of points seasonal time-lag<BR>
+      &nbsp;&nbsp;$ python ${pyHOME}/drawing.py -o draw_time_lag -f 'Forcing_rainstats_OKstomate_CRUNCEP_spinup_SESA_seasmean_anom.nc;time_centered,0@-9@4;time_centered;rainanom%DiPolo_evapstats_OKstomate_CRUNCEP_SESA_seasmean_anom.nc;time_centered,0@-9@4;time_centered;evapanom#Forcing_rainstats_OKstomate_CRUNCEP_spinup_SESA_seasmean_anom.nc;time_centered,1@-9@4;time_centered;rainanom%DiPolo_evapstats_OKstomate_CRUNCEP_SESA_seasmean_anom.nc;time_centered,1@-9@4;time_centered;evapanom#Forcing_rainstats_OKstomate_CRUNCEP_spinup_SESA_seasmean_anom.nc;time_centered,2@-9@4;time_centered;rainanom%DiPolo_evapstats_OKstomate_CRUNCEP_SESA_seasmean_anom.nc;time_centered,2@-9@4;time_centered;evapanom#Forcing_rainstats_OKstomate_CRUNCEP_spinup_SESA_seasmean_anom.nc;time_centered,3@-9@4;time_centered;rainanom%DiPolo_evapstats_OKstomate_CRUNCEP_SESA_seasmean_anom.nc;time_centered,3@-9@4;time_centered;evapanom' -S 'auto;centered;pr;mm/d;evap;mm/d;extremescentered;extremescentered;year|seconds!since!1958-01-01!00:00:00|exct,5,y|%Y;DJF,MAM,JJA,SON;simplepts,1,#AA0000@#0000AA@#00AA00@#AAAA00,x,2.;1-seas!time-lag!at!SESA!between!pr!&!evap!;0|12;png;seas1_timelag_pr-evap_SESA_anom;yes' <BR>
       * Plotting with connected lines <BR>
       &nbsp;&nbsp;$ python ${pyHOME}/drawing.py -o draw_time_lag -f 'Forcing_rainstats_OKstomate_SESA_norm_meanstd_stat2.nc;time_centered,-1;time_counter;rainnorm%DiPolo_evapstats_SESA_norm_meanstd_OKstomate.nc;time_counter,-1|stats,2;time_counter;evapnorm#Forcing_rainstats_OKstomate_SACZ_norm_meanstd_stat2.nc;time_centered,-1;time_counter;rainnorm%DiPolo_evapstats_SACZ_norm_meanstd_OKstomate.nc;time_counter,-1|stats,2;time_counter;evapnorm' -S 'auto;centered;pr;1;evap;1;-5,5;-2,2;year|seconds!since!1958-01-01!00:00:00|exct,5,y|%Y;SESA,SACZ;shadTline,1,spectral,%s,horizontal,-@-.,2.,x,2.;6-month!time-lag!at!SACZ!between!pr!&!evap!normalized!by!mean/stddev;0|12;png;mon6_timelag_pr-evap_SACZ_shadline;yes' <BR>

trunk/tools/drawing_tools.py

@@ -9897 +9897 @@
       ifclose= whether figure should be closed
     """
+    from scipy import stats as sts
     fname = 'plot_time_lag'
 
     Nvalues = len(prevals)
 
-    dxn = gen.fillValue
-    dxx = -gen.fillValue
-    dyn = gen.fillValue
-    dyx = -gen.fillValue
-
-    for iv in range(Nvalues):
-        minv = np.min(prevals[iv])
-        maxv = np.max(prevals[iv])
-        if minv < dxn: dxn = minv
-        if maxv > dxx: dxx = maxv
-        minv = np.min(futvals[iv])
-        maxv = np.max(futvals[iv])
-        if minv < dyn: dyn = minv
-        if maxv > dyx: dyx = maxv
+    dxn = prerange[0]
+    dxx = prerange[1]
+    dyn = futrange[0]
+    dyx = futrange[1]
 
     axstyles = ['pretty', 'Nfix', 'Stringlist', 'Vfix']
@@ -10073 +10064 @@
 
             Nlag = int(figv[1])
-            col = gen.str_list_rep(figv[2],'@',Nvalues)
+            if figv[2] == 'auto':
+                col = colorsauto
+            else:
+                col = gen.str_list_rep(figv[2],'@',Nvalues)
             mark = gen.str_list_rep(figv[3],'@',Nvalues)
             size = gen.str_list_rep(figv[4],'@',Nvalues)
@@ -10083 +10077 @@
                 plt.plot(xvs[0:dimt-Nlag], yvs[Nlag:dimt], ',', color=col[iv],       \
                   marker=mark[iv], markersize=size[iv])
+            # Computing correlation
+            corr, p_value = sts.pearsonr(xvs[0:dimt-Nlag], yvs[Nlag:dimt])
+            txt='corr= ' + '{0:.3g}'.format(corr) + ' p: ' + '{0:.2g}'.format(p_value)
+            plt.annotate(txt,xy=(0.95,0.05*(Nvalues-iv)), xycoords='figure fraction',\
+              color=col[iv], fontsize=lsiz, horizontalalignment='right')
         else:
             print errormsg