New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
std_plot_all_last_year_5D.pro in branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts – NEMO

source: branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro @ 4184

Last change on this file since 4184 was 4184, checked in by flavoni, 10 years ago

update IDL_scripts for ORCA2_LIM and ORCA2_LIM3, add plot of max_mld output 5days, see ticket: #724

File size: 22.0 KB
Line 
1pro std_plot_all_last_year_5D, doplot = doplot, _extra = ex
2
3  compile_opt idl2, strictarrsubs
4
5@common
6@std_common
7                                ; scripts for nemo v3_2 and v3_3
8
9;; script for max MLD, output 5 days
10
11  PRINT, ''
12  PRINT, '  ############################################'
13  PRINT, ''
14  PRINT, '                    LAUNCH of std_plots_last_year_5D'
15  PRINT, ''
16  PRINT, '  ############################################'
17  PRINT, ''
18;
19  std_iodir_data    = isadirectory(getenv('DIR_DATA'),     title = 'path of data in NetCdf format')
20  std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'),  title = 'path of climatological data')
21  std_iodir_mask    = isadirectory(getenv('DIR_MASK'),     title = 'path of mask files (ex: subbasins)')
22; meshmask
23  std_file_mesh = isafile(getenv('FILE_MESH_MASK'),        title = 'mesh_mask', iodir = std_iodir_mask)
24  std_file_msksub = isafile(getenv('FILE_MASK_SUBDOMAIN'), title = 'sub-bassin masks', iodir = std_iodir_mask)
25
26; climatologies
27  std_file_Levitus_T =  isafile(getenv('FILE_TEMP_3D'),    title = 'Levitus_T', iodir = std_iodir_climato)
28  std_file_Levitus_S =  isafile(getenv('FILE_SAL_3D'),     title = 'Levitus_S', iodir = std_iodir_climato)
29  std_file_reynolds  =  isafile(getenv('FILE_SST'),        title = 'Reynolds', iodir = std_iodir_climato)
30  std_file_oaflux    =  isafile(getenv('FILE_FLUX'),       title = 'oaflux', iodir = std_iodir_climato)
31  std_file_mld       =  isafile(getenv('FILE_MLD'),        title = 'Mixed layer depth', iodir = std_iodir_climato)
32  std_file_ice       =  isafile(getenv('FILE_ICE'),        title = 'ICE', iodir = std_iodir_climato)
33  std_file_snow_arc  =  isafile(getenv('FILE_SNOW_ARC'),   title = 'SNOW_ARC', iodir = std_iodir_climato)
34  std_file_snow_ant  =  isafile(getenv('FILE_SNOW_ANT'),   title = 'SNOW_ANT', iodir = std_iodir_climato)
35
36  IF strlowcase(getenv('FILE_GEOHEAT')) EQ 'no' THEN std_file_geoheat = 'no' $
37  ELSE std_file_geoheat =  isafile(getenv('FILE_GEOHEAT'), title = 'Geothermal heating', iodir = std_iodir_climato)
38;
39  allrec =  1 - keyword_set(long(getenv('READ_ONLY_FIRST_RECORD')))
40; Output run experience1
41  std_file1_T     = isafile(getenv('FILE1_T'), title = 'exp1 grid T input file', iodir = std_iodir_data)
42  std_file1_U     = isafile(getenv('FILE1_U'), title = 'exp1 grid U input file', iodir = std_iodir_data)
43  std_file1_V     = isafile(getenv('FILE1_V'), title = 'exp1 grid V input file', iodir = std_iodir_data)
44  std_file1_I     = isafile(getenv('FILE1_I'), title = 'exp1 ice    input file', iodir = std_iodir_data)
45 
46; Output run experience2
47  std_file2_T     = isafile(getenv('FILE2_T'), title = 'exp2 grid T input file', iodir = std_iodir_data)
48  std_file2_U     = isafile(getenv('FILE2_U'), title = 'exp2 grid U input file', iodir = std_iodir_data)
49  std_file2_V     = isafile(getenv('FILE2_V'), title = 'exp2 grid V input file', iodir = std_iodir_data)
50  std_file2_I     = isafile(getenv('FILE2_I'), title = 'exp2 ice    input file', iodir = std_iodir_data)
51
52  PRINT, ''
53  PRINT, '  std_iodir_data : ' + std_iodir_data
54  PRINT, '  std_file1T : ' + std_file1_T
55  PRINT, '  std_file1U : ' + std_file1_U
56  PRINT, '  std_file1V : ' + std_file1_V
57;  PRINT, ' std_file1W : ' + std_file1_W
58  PRINT, '  std_file2I : ' + std_file1_I
59  PRINT, '  std_file2T : ' + std_file2_T
60  PRINT, '  std_file2U : ' + std_file2_U
61  PRINT, '  std_file2V : ' + std_file2_V
62;  PRINT, ' std_file2W : ' + std_file2_W
63  PRINT, '  std_file2I : ' + std_file2_I
64  PRINT, ''
65
66;#########################################################################
67;##########################  Load Grids   ################################
68;#########################################################################
69; load the grid
70  load_orca, std_file_mesh 
71; reading variables
72  masknp = read_ncdf('tmaskutil', file = std_file_mesh, /nostruct, /cont_nofill)
73;#########################################################################
74;############################  Read Data  ################################
75;#########################################################################
76;
77  allrec =  1; - keyword_set(long(getenv('READ_ONLY_FIRST_RECORD')))
78;
79;;; 3D ;;;
80; temperature
81  T1 = read_ncdf(getenv('VAR1_T'), allrecords = allrec, direc = 't', filename = std_file1_T )
82  IF std_file2_T NE std_file1_T THEN BEGIN
83    T2 = read_ncdf(getenv('VAR2_T'), allrecords = allrec, direc = 't', filename = std_file2_T )
84  ENDIF ELSE T2 = {arr:-1}
85  TLev = read_ncdf(getenv('VAR_TEMP_3D'), filename = std_file_Levitus_T )
86  TRey = read_ncdf(getenv('VAR_SST'), filename = std_file_reynolds )
87
88; salinity
89  S1 = read_ncdf(getenv('VAR1_S'), allrecords = allrec, direc = 't', filename = std_file1_T )
90  IF std_file2_T NE std_file1_T THEN BEGIN
91    S2 = read_ncdf(getenv('VAR2_S'), allrecords = allrec, direc = 't', filename = std_file2_T )
92  ENDIF ELSE S2 = {arr:-1}
93  SLev = read_ncdf(getenv('VAR_SAL_3D'), filename = std_file_Levitus_S )
94
95;;; 2D ;;;
96; Net Downward heat flux
97  Q1 = read_ncdf(getenv('VAR1_QNET'), allrecords = allrec, direc = 't', filename = std_file1_T )
98  IF std_file2_T NE std_file1_T THEN BEGIN
99    Q2 = read_ncdf(getenv('VAR2_QNET'), allrecords = allrec, direc = 't', filename = std_file2_T )
100  ENDIF ELSE Q2 = {arr:-1}
101; Geothermal heating
102  IF std_file_geoheat EQ 'no' THEN geo = {arr:float(getenv('VAR_GEOHEAT'))} $
103  ELSE geo = read_ncdf(getenv('VAR_GEOHEAT'), filename =  std_file_geoheat )
104  geo = geo.arr*1.e-3          ; convert into W/m2
105;climatology
106  QNET = read_ncdf(getenv('VAR_FLUX'), filename = std_file_oaflux )
107
108; erp (evaporation damping)
109  ERP1 = read_ncdf(getenv('VAR1_ERP'), allrecords = allrec, direc = 't', filename = std_file1_T )
110  ERP1 = {arr:ERP1.arr * 86400., unit:'mm/day', grid:'T'}
111  IF std_file2_T NE std_file1_T THEN BEGIN
112    ERP2 = read_ncdf(getenv('VAR2_ERP'), allrecords = allrec, direc = 't', filename = std_file2_T )
113    ERP2 = {arr:ERP2.arr * 86400., unit:'mm/day', grid:'T'}
114  ENDIF ELSE ERP2 = {arr:-1}
115
116; emp (evaporation minus precipitation)
117  EMP1 = read_ncdf(getenv('VAR1_EMP'), allrecords = allrec, direc = 't', filename = std_file1_T )
118  EMP1 = {arr:EMP1.arr * 86400., unit:'mm/day', grid:'T'}
119  IF std_file2_T NE std_file1_T THEN BEGIN
120    EMP2 = read_ncdf(getenv('VAR2_EMP'), allrecords = allrec, direc = 't', filename = std_file2_T )
121    EMP2 = {arr:EMP2.arr * 86400., unit:'mm/day', grid:'T'}
122  ENDIF ELSE EMP2 = {arr:-1}
123 
124  ;mixed layer depth
125  ;; nyear=jpt/73
126  MLD_1 = read_ncdf(getenv('VAR1_MLD'), allrecords = allrec, filename = std_file1_T )
127  caldat, time, mm
128  march = where(mm EQ 3, cnt)
129  MLD_march_1 = max(MLD_1.arr[*,*,march], dim=3)
130  MLD_march_1 = {arr: MLD_march_1, unit: MLD_1.unit}
131  sept = where(mm EQ 9, cnt)
132  MLD_sept_1 = max(MLD_1.arr[*,*,sept], dim=3)
133  MLD_sept_1 = {arr: MLD_sept_1, unit: MLD_1.unit}
134  MLD_1 = {arr: max(MLD_1.arr[*,*,*], dim=3), unit: MLD_1.unit}
135  IF std_file2_T NE std_file1_T THEN BEGIN
136    MLD_2 = read_ncdf(getenv('VAR2_MLD'), allrecords = allrec, filename = std_file2_T ) ; 10 m
137    caldat, time, mm
138    march = where(mm EQ 3, cnt)
139    MLD_march_2 = max(MLD_2.arr[*,*,march], dim=3)
140    MLD_march_2 = {arr: MLD_march_2, unit: MLD_2.unit}
141    sept = where(mm EQ 9, cnt)
142    MLD_sept_2 = max(MLD_2.arr[*,*,sept], dim=3)
143    MLD_sept_2 = {arr:MLD_sept_2, unit: MLD_2.unit}
144    MLD_2 = {arr: max(MLD_2.arr[*,*,*], dim=3), unit: MLD_2.unit}
145  ENDIF ELSE BEGIN
146    MLD_2 = {arr:-1}
147    MLD_march_2 = {arr:-1}
148    MLD_sept_2 = {arr:-1}
149  ENDELSE
150 
151  ;climatology
152  MLD = read_ncdf(getenv('VAR_MLD'), filename = std_file_mld )
153; if monthly climatology
154;  MLD = {arr: max(MLD.arr[*,*,*], dim=3), unit: MLD.unit}
155
156; velocities
157  U1 = read_ncdf(getenv('VAR1_U'), allrecords = allrec, direc = 't', filename = std_file1_U )
158  IF strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' THEN BEGIN
159    U1.arr = U1.arr / e3u_3d(/e2) * umask()
160    U1.unit = 'm/s'
161  ENDIF
162  IF std_file2_U NE std_file1_U THEN BEGIN
163    U2 = read_ncdf(getenv('VAR2_U'), allrecords = allrec, direc = 't', filename = std_file2_U )
164    IF strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' THEN BEGIN
165      U2.arr = U2.arr / e3u_3d(/e2) * umask()
166      U2.unit = 'm/s'
167    ENDIF
168  ENDIF ELSE U2 = {arr:-1}
169;
170  V1 = read_ncdf(getenv('VAR1_V'), allrecords = allrec, direc = 't', filename = std_file1_V )
171  IF strlowcase(getenv('VAR1_V')) EQ 'vocetr_eff' THEN BEGIN
172    V1.arr = V1.arr / e3v_3d(/e1) * vmask()
173    V1.unit = 'm/s'
174  ENDIF
175  IF std_file2_V NE std_file1_V THEN BEGIN
176    V2 = read_ncdf(getenv('VAR2_V'), allrecords = allrec, direc = 't', filename = std_file2_V )
177    IF strlowcase(getenv('VAR2_V')) EQ 'vocetr_eff' THEN BEGIN
178      V2.arr = V2.arr / e3v_3d(/e1) * vmask()
179      V2.unit = 'm/s'
180    ENDIF
181  ENDIF ELSE V2 = {arr:-1}
182
183; ice
184  Ithi_1 = read_ncdf(getenv('VAR1_Ithick'), allrecords = allrec, filename = std_file1_I )
185  caldat, time, mm
186  april = where(mm EQ 4, cnt)
187  Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}
188  sept = where(mm EQ 9, cnt)
189  Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}
190  undefine, Ithi_1
191;
192  IF std_file2_I NE std_file1_I THEN BEGIN
193     Ithi_2 = read_ncdf(getenv('VAR2_Ithick'), allrecords = allrec, filename = std_file2_I )
194     caldat, time, mm
195     april = where(mm EQ 4, cnt)
196     Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}
197     sept = where(mm EQ 9, cnt)
198     Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}
199     undefine, Ithi_2
200  ENDIF ELSE BEGIN
201     Ithi_april_2 = {arr:-1}
202     Ithi_sept_2 = {arr:-1}
203  ENDELSE
204;
205  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )
206  ;; help, Ifra_1.arr, jpi,jpj
207  caldat, time, mm
208  march = where(mm EQ 3, cnt)
209  Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}
210  sept = where(mm EQ 9, cnt)
211  Ifra_sept_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}
212  undefine, Ifra_1
213;
214  IF std_file2_I NE std_file1_I THEN BEGIN
215    Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )
216    caldat, time, mm
217    march = where(mm EQ 3, cnt)
218    Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit}
219    sept = where(mm EQ 9, cnt)
220    Ifra_sept_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_2.unit}
221    undefine, Ifra_2
222  ENDIF ELSE BEGIN
223    Ifra_march_2 = {arr:-1}
224    Ifra_sept_2 = {arr:-1}
225  ENDELSE
226;
227;; SF commentato perché non funziona snow
228;; SF   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )
229;; SF   help, Isnow_1.arr, jpi,jpj
230;; SF   caldat, time, mm
231;; SF   april = where(mm EQ 4, cnt)
232;; SF   Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}
233;; SF   sept = where(mm EQ 9, cnt)
234;; SF   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}
235;; SF   ;undefine, Isnow_1
236;; SF ;
237;; SF   IF std_file2_I NE std_file1_I THEN BEGIN
238;; SF     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )
239;; SF     caldat, time, mm
240;; SF     april = where(mm EQ 4, cnt)
241;; SF     Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
242;; SF     sept = where(mm EQ 9, cnt)
243;; SF     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
244;; SF   ;  undefine, Isnow_2
245;; SF   ENDIF ELSE BEGIN
246;; SF     Isnow_april_2 = {arr:-1}
247;; SF     Isnow_sept_2 = {arr:-1}
248;; SF   ENDELSE
249;
250  jpt = 1
251;
252; shorter file names for legends...
253;
254  std_file1_T = file_basename(std_file1_T,'.nc')
255  std_file1_T = (strsplit(std_file1_T,'_grid_T',/extract,/regex))[0]
256  std_file2_T = file_basename(std_file2_T,'.nc')
257  std_file2_T = (strsplit(std_file2_T,'_grid_T',/extract,/regex))[0]
258  std_file1_U = file_basename(std_file1_U,'.nc')
259  std_file1_U = (strsplit(std_file1_U,'_grid_U',/extract,/regex))[0]
260  std_file2_U = file_basename(std_file2_U,'.nc')
261  std_file2_U = (strsplit(std_file2_U,'_grid_U',/extract,/regex))[0]
262  std_file1_V = file_basename(std_file1_V,'.nc')
263  std_file1_V = (strsplit(std_file1_V,'_grid_V',/extract,/regex))[0]
264  std_file2_V = file_basename(std_file2_V,'.nc')
265  std_file2_V = (strsplit(std_file2_V,'_grid_V',/extract,/regex))[0]
266  std_file1_I = file_basename(std_file1_I,'.nc')
267  std_file1_I = (strsplit(std_file1_I,'_icemod',/extract,/regex))[0]
268  std_file2_I = file_basename(std_file2_I,'.nc')
269  std_file2_I = (strsplit(std_file2_I,'_icemod',/extract,/regex))[0]
270
271;#########################################################################
272;######################  STANDARD PLOTS   ################################
273;#########################################################################
274
275  IF keyword_set(doplot) EQ 0 THEN doplot = 0
276
277; fixed color tabled
278  lct, 64
279  cnt = 0
280  htmltxt = ''
281;
282  cnt = cnt+1   &   blabla = 'Erp salinity damping term'
283  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex
284
285  cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term'
286  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex
287
288  cnt = cnt+1   &   blabla = 'Net heat flux'
289  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex
290
291  cnt = cnt+1   &   blabla = 'Meridionnal Heat Transport'
292  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex
293
294  cnt = cnt+1   &   blabla = 'Global Barotropic stream Function'
295  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex
296
297  cnt = cnt+1   &   blabla = 'mean Temperature diff with New Reynolds'
298  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex
299
300  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus'
301  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex
302
303  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus'
304  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex       
305
306  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2'
307  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex
308
309  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters'
310  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex       
311 
312  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters'
313  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex
314
315  cnt = cnt+1   &   blabla = 'mean Temperature diff with Levitus at z=100 meters'
316  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex
317 
318  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus at z=100 meters'
319  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex
320 
321;;SF need monthly climatology to plot difference of MAX of MLD 
322  cnt = cnt+1   &   blabla = 'Mixed layer depth & difference with DeBoyer'
323  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD, _extra = ex
324   
325  cnt = cnt+1   &   blabla = 'Mixed layer depth differences'
326  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD_2, MLD, _extra = ex
327
328  ;SF: add 2 cases for max of Mixed Layer depth on march & septembre
329  cnt = cnt+1   &   blabla = 'Mixed layer depth on March'
330  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_max_mld, MLD_march_1, MLD_march_2, /MARCH, _extra = ex
331
332  cnt = cnt+1   &   blabla = 'Mixed layer depth on September'
333  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_max_mld, MLD_sept_1, MLD_sept_2, /SEPT, _extra = ex
334
335;; SF: commented because not ok for max MLD, output 5 days
336;; cnt = cnt+1   &   blabla = 'Zonal mean Mixed layer depth'
337;; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ZonMld, MLD_1, MLD_2, MLD, _extra = ex
338 
339  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Global'
340  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex
341
342  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Atlantic'
343  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex
344 
345  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Indian'
346  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex
347 
348  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Pacific'
349  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex
350
351  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Global'
352  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex
353
354  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Atlantic'
355  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex
356
357  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Indian'
358  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex
359
360  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Pacific'
361  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex
362 
363;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH'
364;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex
365
366  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL'
367  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex
368
369  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT'
370  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex
371 
372;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH'
373;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex
374
375  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL'
376  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex
377
378  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT'
379  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex
380
381  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH'
382  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex
383 
384  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: SEPT'
385  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex
386 
387  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH'
388  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex
389
390  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT'
391  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex
392
393;; SF non funziona, dice che deve essere una struttura e non e!!!!
394;;SF  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL'
395;;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex
396;; SF
397;; SF   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT'
398;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex
399;; SF 
400;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL'
401;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex
402;; SF 
403;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT'
404;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex
405;
406;  cnt = cnt+1   &   blabla = 'Arctic Ice Extend: '
407;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex
408 
409;  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: '
410;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex
411
412  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)'
413  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex
414
415  cnt = cnt+1   &   blabla = 'Meridional stream Function: Atlantic'
416  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex
417
418  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indian'
419  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex
420
421  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indo-Pacific'
422  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex
423
424  cnt = cnt+1   &   blabla = 'Equatorial Temperature'
425  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex
426
427  cnt = cnt+1   &   blabla = 'Equatorial Salinity'
428  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex
429
430  cnt = cnt+1   &   blabla = 'Equatorial zonal velocity'
431  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex
432
433  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=700'
434  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex
435
436  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=1000'
437  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex
438
439  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40°N'
440  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex
441
442  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38°N'
443  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex
444
445  cnt = cnt+1   &   blabla = 'Vertical Global mean T & S'
446  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex
447
448  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*]
449 
450  return
451END
Note: See TracBrowser for help on using the repository browser.