Changeset 128

Timestamp:

05/29/09 11:55:47 (15 years ago)

Author:

ericg

Message:

Added mean and rms info in plt xy map requires new moyenne.pro

File:

• trunk/procs/plt_map.pro (modified) (3 diffs)

trunk/procs/plt_map.pro

@@ -795 +795 @@
          ;      
                print, '     '+cmd.var+' data to plot min and max :            ', $
-                min(fld(where (fld NE valmask)), /nan), max(fld(where (fld NE valmask)), /nan)
+                min(fld(where (fld LE valmask/10)), /nan), max(fld(where (fld LE valmask/10)), /nan)
                print, ''
 
@@ -801 +801 @@
 
                niveau = 1
+               
+               ; compute mean and rms of field
+               zmean = moyenne(fld, 'xy', boite = box_plt, /NAN)
+               zrms = moyenne(fld, 'xy', boite = box_plt, /RMS, /NAN)
+               filleg = ' [region mean/rms='+string(strtrim(zmean, 2))+'/'+string(strtrim(zrms, 2))+']'
 
                IF vecplot EQ 0 THEN BEGIN 
@@ -808 +813 @@
                      fldzm = replicate(1, (size(fld))(1))#fldzm
                      fld = fld-fldzm
-                     filleg = ' Eddy '
+                     edd_txt = ' Eddy '
                   ENDIF ELSE BEGIN 
-                     filleg = ''
+                     edd_txt = '' 
                   ENDELSE 
+                  
 
                   IF hotyp EQ 'xyt' THEN date_txt = 'monthly'
 
-                  pltcmd = 'plt,{a:fld, d:'''+date_txt+''', n:'''+filleg+varlegend+' '+legbox+legsurf+''', u:'''+field.units+''', g:'''+vargrid+'''}'+mincmd+maxcmd+intcmd+com_strplt+strcont+map_cmd
+                  pltcmd = 'plt,{a:fld, d:'''+date_txt+''', n:'''+edd_txt+varlegend+' '+legbox+legsurf+''', u:'''+field.units+filleg+''', g:'''+vargrid+'''}'+mincmd+maxcmd+intcmd+com_strplt+strcont+map_cmd
                   ;; For 2D plots with overlays (correlation and proba for instance)
                   ;; Scratch the titles and define contours style for the second plot