source: trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro @ 2751

Last change on this file since 2751 was 2751, checked in by smasson, 10 years ago

large update of IDL_scripts for Is-ENES

File size: 17.4 KB
Line 
1pro std_plot_all, 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  PRINT, ''
10  PRINT, '  ############################################'
11  PRINT, ''
12  PRINT, '                    LAUNCH of std_plots'
13  PRINT, ''
14  PRINT, '  ############################################'
15  PRINT, ''
16;
17  std_iodir_data    = isadirectory(getenv('DIR_DATA'), title = 'path of data in NetCdf format')
18  std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'), title = 'path of climatological data')
19  std_iodir_mask    = isadirectory(getenv('DIR_MASK'), title = 'path of mask files (ex: subbasins)')
20; meshmask
21  std_file_mesh = isafile(getenv('FILE_MESH_MASK'), title = 'mesh_mask', iodir = std_iodir_mask)
22  std_file_msksub = isafile(getenv('FILE_MASK_SUBDOMAIN'), title = 'sub-bassin masks', iodir = std_iodir_mask)
23
24; Levitus 98
25  std_file_Levitus_T =  isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato)
26  std_file_Levitus_S =  isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato)
27  std_file_reynolds =  isafile(getenv('FILE_SST'), title = 'Reynolds', iodir = std_iodir_climato)
28  std_file_oaflux =  isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato)
29  std_file_mld =  isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato)
30
31  IF strlowcase(getenv('FILE_GEOHEAT')) EQ 'no' THEN std_file_geoheat = 'no' $
32  ELSE std_file_geoheat =  isafile(getenv('FILE_GEOHEAT'), title = 'Geothermal heating', iodir = std_iodir_climato)
33;
34  allrec =  1 - keyword_set(long(getenv('READ_ONLY_FIRST_RECORD')))
35; Output run experience1
36  std_file1_T     = isafile(getenv('FILE1_T'), title = 'exp1 grid T input file', iodir = std_iodir_data)
37  std_file1_U     = isafile(getenv('FILE1_U'), title = 'exp1 grid U input file', iodir = std_iodir_data)
38  std_file1_V     = isafile(getenv('FILE1_V'), title = 'exp1 grid V input file', iodir = std_iodir_data)
39  std_file1_I     = isafile(getenv('FILE1_I'), title = 'exp1 ice    input file', iodir = std_iodir_data)
40 
41; Output run experience2
42  std_file2_T     = isafile(getenv('FILE2_T'), title = 'exp2 grid T input file', iodir = std_iodir_data)
43  std_file2_U     = isafile(getenv('FILE2_U'), title = 'exp2 grid U input file', iodir = std_iodir_data)
44  std_file2_V     = isafile(getenv('FILE2_V'), title = 'exp2 grid V input file', iodir = std_iodir_data)
45  std_file2_I     = isafile(getenv('FILE2_I'), title = 'exp2 ice    input file', iodir = std_iodir_data)
46
47  PRINT, ''
48  PRINT, '  std_iodir_data : ' + std_iodir_data
49  PRINT, '  std_file1T : ' + std_file1_T
50  PRINT, '  std_file1U : ' + std_file1_U
51  PRINT, '  std_file1V : ' + std_file1_V
52;  PRINT, ' std_file1W : ' + std_file1_W
53  PRINT, '  std_file2I : ' + std_file1_I
54  PRINT, '  std_file2T : ' + std_file2_T
55  PRINT, '  std_file2U : ' + std_file2_U
56  PRINT, '  std_file2V : ' + std_file2_V
57;  PRINT, ' std_file2W : ' + std_file2_W
58  PRINT, '  std_file2I : ' + std_file2_I
59  PRINT, ''
60
61;#########################################################################
62;##########################  Load Grids   ################################
63;#########################################################################
64; load the grid
65  load_orca, std_file_mesh 
66; reading variables
67  masknp = read_ncdf('tmaskutil', file = std_file_mesh, /nostruct, /cont_nofill)
68;#########################################################################
69;############################  Read Data  ################################
70;#########################################################################
71;
72  allrec =  1; - keyword_set(long(getenv('READ_ONLY_FIRST_RECORD')))
73;
74;;; 3D ;;;
75; temperature
76  T1 = read_ncdf(getenv('VAR1_T'), allrecords = allrec, direc = 't', filename = std_file1_T )
77  IF std_file2_T NE std_file1_T THEN BEGIN
78    T2 = read_ncdf(getenv('VAR2_T'), allrecords = allrec, direc = 't', filename = std_file2_T )
79  ENDIF ELSE T2 = {arr:-1}
80  TLev = read_ncdf(getenv('VAR_TEMP_3D'), filename = std_file_Levitus_T )
81  TRey = read_ncdf(getenv('VAR_SST'), filename = std_file_reynolds )
82
83; salinity
84  S1 = read_ncdf(getenv('VAR1_S'), allrecords = allrec, direc = 't', filename = std_file1_T )
85  IF std_file2_T NE std_file1_T THEN BEGIN
86    S2 = read_ncdf(getenv('VAR2_S'), allrecords = allrec, direc = 't', filename = std_file2_T )
87  ENDIF ELSE S2 = {arr:-1}
88  SLev = read_ncdf(getenv('VAR_SAL_3D'), filename = std_file_Levitus_S )
89
90;;; 2D ;;;
91; Net Downward heat flux
92  Q1 = read_ncdf(getenv('VAR1_QNET'), allrecords = allrec, direc = 't', filename = std_file1_T )
93  IF std_file2_T NE std_file1_T THEN BEGIN
94    Q2 = read_ncdf(getenv('VAR2_QNET'), allrecords = allrec, direc = 't', filename = std_file2_T )
95  ENDIF ELSE Q2 = {arr:-1}
96; Geothermal heating
97  IF std_file_geoheat EQ 'no' THEN geo = {arr:float(getenv('VAR_GEOHEAT'))} $
98  ELSE geo = read_ncdf(getenv('VAR_GEOHEAT'), filename =  std_file_geoheat )
99  geo = geo.arr*1.e-3          ; convert into W/m2
100;climatology
101  QNET = read_ncdf(getenv('VAR_FLUX'), filename = std_file_oaflux )
102
103; erp (evaporation damping)
104  ERP1 = read_ncdf(getenv('VAR1_ERP'), allrecords = allrec, direc = 't', filename = std_file1_T )
105  ERP1 = {arr:ERP1.arr * 86400., unit:'mm/day', grid:'T'}
106  IF std_file2_T NE std_file1_T THEN BEGIN
107    ERP2 = read_ncdf(getenv('VAR2_ERP'), allrecords = allrec, direc = 't', filename = std_file2_T )
108    ERP2 = {arr:ERP2.arr * 86400., unit:'mm/day', grid:'T'}
109  ENDIF ELSE ERP2 = {arr:-1}
110
111; emp (evaporation minus precipitation)
112  EMP1 = read_ncdf(getenv('VAR1_EMP'), allrecords = allrec, direc = 't', filename = std_file1_T )
113  EMP1 = {arr:EMP1.arr * 86400., unit:'mm/day', grid:'T'}
114  IF std_file2_T NE std_file1_T THEN BEGIN
115    EMP2 = read_ncdf(getenv('VAR2_EMP'), allrecords = allrec, direc = 't', filename = std_file2_T )
116    EMP2 = {arr:EMP2.arr * 86400., unit:'mm/day', grid:'T'}
117  ENDIF ELSE EMP2 = {arr:-1}
118 
119;mixed layer depth
120  MLD1 = read_ncdf(getenv('VAR1_MLD'), allrecords = allrec, direc = 't', filename = std_file1_T ) ; 10 m
121  IF std_file2_T NE std_file1_T THEN BEGIN
122    MLD2 = read_ncdf(getenv('VAR2_MLD'), allrecords = allrec, direc = 't', filename = std_file2_T ) ; 10 m
123  ENDIF ELSE MLD2 = {arr:-1}
124;climatology
125  MLD = read_ncdf(getenv('VAR_MLD'), filename = std_file_mld )
126
127; velocities
128  U1 = read_ncdf(getenv('VAR1_U'), allrecords = allrec, direc = 't', filename = std_file1_U )
129  IF strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' THEN BEGIN
130    U1.arr = U1.arr / e3u_3d(/e2) * umask()
131    U1.unit = 'm/s'
132  ENDIF
133  IF std_file2_U NE std_file1_U THEN BEGIN
134    U2 = read_ncdf(getenv('VAR2_U'), allrecords = allrec, direc = 't', filename = std_file2_U )
135    IF strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' THEN BEGIN
136      U2.arr = U2.arr / e3u_3d(/e2) * umask()
137      U2.unit = 'm/s'
138    ENDIF
139  ENDIF ELSE U2 = {arr:-1}
140;
141  V1 = read_ncdf(getenv('VAR1_V'), allrecords = allrec, direc = 't', filename = std_file1_V )
142  IF strlowcase(getenv('VAR1_V')) EQ 'vocetr_eff' THEN BEGIN
143    V1.arr = V1.arr / e3v_3d(/e1) * vmask()
144    V1.unit = 'm/s'
145  ENDIF
146  IF std_file2_V NE std_file1_V THEN BEGIN
147    V2 = read_ncdf(getenv('VAR2_V'), allrecords = allrec, direc = 't', filename = std_file2_V )
148    IF strlowcase(getenv('VAR2_V')) EQ 'vocetr_eff' THEN BEGIN
149      V2.arr = V2.arr / e3v_3d(/e1) * vmask()
150      V2.unit = 'm/s'
151    ENDIF
152  ENDIF ELSE V2 = {arr:-1}
153
154; ice
155  Ithi_1 = read_ncdf(getenv('VAR1_Ithick'), allrecords = allrec, filename = std_file1_I )
156  caldat, time, mm
157  march = where(mm EQ 3, cnt)
158  Ithi_march_1 = {arr:1./float(cnt) * total(Ithi_1.arr[*, *, temporary(march)], 3), unit:Ithi_1.unit}
159  sept = where(mm EQ 3, cnt)
160  Ithi_sept_1 = {arr:1./float(cnt) * total(Ithi_1.arr[*, *, temporary(sept)], 3), unit:Ithi_1.unit}
161  undefine, Ithi_1
162;
163  IF std_file2_I NE std_file1_I THEN BEGIN
164    Ithi_2 = read_ncdf(getenv('VAR2_Ithick'), allrecords = allrec, filename = std_file2_I )
165    caldat, time, mm
166    march = where(mm EQ 3, cnt)
167    Ithi_march_2 = {arr:1./float(cnt) * total(Ithi_2.arr[*, *, temporary(march)], 3), unit:Ithi_2.unit}
168    sept = where(mm EQ 3, cnt)
169    Ithi_sept_2 = {arr:1./float(cnt) * total(Ithi_2.arr[*, *, temporary(sept)], 3), unit:Ithi_2.unit}
170    undefine, Ithi_2
171  ENDIF ELSE BEGIN
172    Ithi_march_2 = {arr:-1}
173    Ithi_sept_2 = {arr:-1}
174  ENDELSE
175;
176  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )
177  caldat, time, mm
178  march = where(mm EQ 3, cnt)
179  Ifra_march_1 = {arr:1./float(cnt) * total(Ifra_1.arr[*, *, temporary(march)], 3), unit:Ifra_1.unit}
180  sept = where(mm EQ 9, cnt)
181  Ifra_sept_1 = {arr:1./float(cnt) * total(Ifra_1.arr[*, *, temporary(sept)], 3), unit:Ifra_1.unit}
182  undefine, Ifra_1
183;
184  IF std_file2_I NE std_file1_I THEN BEGIN
185    Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )
186    caldat, time, mm
187    march = where(mm EQ 3, cnt)
188    Ifra_march_2 = {arr:1./float(cnt) * total(Ifra_2.arr[*, *, temporary(march)], 3), unit:Ifra_2.unit}
189    sept = where(mm EQ 9, cnt)
190    Ifra_sept_2 = {arr:1./float(cnt) * total(Ifra_2.arr[*, *, temporary(sept)], 3), unit:Ifra_2.unit}
191    undefine, Ifra_2
192  ENDIF ELSE BEGIN
193    Ifra_march_2 = {arr:-1}
194    Ifra_sept_2 = {arr:-1}
195  ENDELSE
196;
197  jpt = 1
198;
199; shorter file names for legends...
200;
201  std_file1_T = file_basename(std_file1_T,'.nc')
202  std_file1_T = (strsplit(std_file1_T,'_grid_T',/extract,/regex))[0]
203  std_file2_T = file_basename(std_file2_T,'.nc')
204  std_file2_T = (strsplit(std_file2_T,'_grid_T',/extract,/regex))[0]
205  std_file1_U = file_basename(std_file1_U,'.nc')
206  std_file1_U = (strsplit(std_file1_U,'_grid_U',/extract,/regex))[0]
207  std_file2_U = file_basename(std_file2_U,'.nc')
208  std_file2_U = (strsplit(std_file2_U,'_grid_U',/extract,/regex))[0]
209  std_file1_V = file_basename(std_file1_V,'.nc')
210  std_file1_V = (strsplit(std_file1_V,'_grid_V',/extract,/regex))[0]
211  std_file2_V = file_basename(std_file2_V,'.nc')
212  std_file2_V = (strsplit(std_file2_V,'_grid_V',/extract,/regex))[0]
213  std_file1_I = file_basename(std_file1_I,'.nc')
214  std_file1_I = (strsplit(std_file1_I,'_icemod',/extract,/regex))[0]
215  std_file2_I = file_basename(std_file2_I,'.nc')
216  std_file2_I = (strsplit(std_file2_I,'_icemod',/extract,/regex))[0]
217
218;#########################################################################
219;######################  STANDARD PLOTS   ################################
220;#########################################################################
221
222  IF keyword_set(doplot) EQ 0 THEN doplot = 0
223
224; fixed color tabled
225  lct, 64
226  cnt = 0
227  htmltxt = ''
228;
229  cnt = cnt+1   &   blabla = 'Erp salinity damping term'
230  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex
231
232  cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term'
233  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex
234
235  cnt = cnt+1   &   blabla = 'Net heat flux'
236  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex
237
238  cnt = cnt+1   &   blabla = 'Meridionnal Heat Transport'
239  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex
240
241  cnt = cnt+1   &   blabla = 'Global Barotropic stream Function'
242  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex
243
244  cnt = cnt+1   &   blabla = 'mean Temperature diff with New Reynolds'
245  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex
246
247  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus'
248  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex
249
250  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus'
251  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex       
252
253  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2'
254  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex
255
256  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters'
257  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex       
258 
259  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters'
260  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex
261
262  cnt = cnt+1   &   blabla = 'mean Temperature diff with Levitus at z=100 meters'
263  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex
264 
265  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus at z=100 meters'
266  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex
267 
268  cnt = cnt+1   &   blabla = 'Mixed layer depth'
269  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD, _extra = ex
270
271  cnt = cnt+1   &   blabla = 'Mixed layer depth differences'
272  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD2, MLD, _extra = ex
273
274  cnt = cnt+1   &   blabla = 'Zonal mean Mixed layer depth'
275  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ZonMld, MLD1, MLD2, MLD, _extra = ex
276 
277  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Global'
278  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex
279
280  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Atlantic'
281  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex
282 
283  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Indian'
284  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex
285 
286  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Pacific'
287  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex
288
289  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Global'
290  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex
291
292  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Atlantic'
293  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex
294
295  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Indian'
296  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex
297
298  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Pacific'
299  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex
300 
301  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH'
302  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex
303 
304  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT'
305  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex
306 
307  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH'
308  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex
309
310  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT'
311  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex
312
313  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH'
314  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex
315 
316  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: SEPT'
317  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex
318 
319  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH'
320  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex
321
322  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT'
323  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex
324
325  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)'
326  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex
327
328  cnt = cnt+1   &   blabla = 'Meridional stream Function: Atlantic'
329  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex
330
331  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indian'
332  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex
333
334  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indo-Pacific'
335  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex
336
337  cnt = cnt+1   &   blabla = 'Equatorial Temperature'
338  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex
339
340  cnt = cnt+1   &   blabla = 'Equatorial Salinity'
341  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex
342
343  cnt = cnt+1   &   blabla = 'Equatorial zonal velocity'
344  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex
345
346  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=700'
347  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex
348
349  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=1000'
350  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex
351
352  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40°N'
353  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex
354
355  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38°N'
356  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex
357
358  cnt = cnt+1   &   blabla = 'Vertical Global mean T & S'
359  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex
360
361  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*]
362 
363  return
364END
Note: See TracBrowser for help on using the repository browser.