source: branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro @ 4394

Last change on this file since 4394 was 4394, checked in by flavoni, 7 years ago

add time series of FRAM STRAIT, and update IDL plots, see ticket: #724

File size: 4.0 KB
Line 
1pro std_ts_all, doplot = doplot, _extra=ex
2
3  compile_opt idl2, strictarrsubs
4
5@common
6@std_common
7
8  PRINT, ''
9  PRINT, '  ############################################'
10  PRINT, ''
11  PRINT, '                   plot of all TIME SERIES'
12  PRINT, ''
13  PRINT, '  ############################################'
14  PRINT, ''
15;
16  std_iodir_data    = isadirectory(getenv('DIR_DATA'), title = 'path of data in NetCdf format')
17  std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'), title = 'path of climatological data')
18  std_iodir_mask    = isadirectory(getenv('DIR_MASK'), title = 'path of mask files (ex: subbasins)')
19; meshmask
20  std_file_mesh = isafile(getenv('FILE_MESH_MASK'), title = 'mesh_mask', iodir = std_iodir_mask)
21  std_file_msksub = isafile(getenv('FILE_MASK_SUBDOMAIN'), title = 'sub-bassin masks', iodir = std_iodir_mask)
22; load the grid
23  load_orca, std_file_mesh 
24; reading variables
25  masknp = read_ncdf('tmaskutil', file = std_file_mesh, /nostruct, /cont_nofill)
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
37  date1   = long(getenv('DATE1'))     &   date2   = long(getenv('DATE2'))
38  date1_2 = long(getenv('DATE1_2'))   &   date2_2 = long(getenv('DATE2_2'))
39
40  allrec =  1 - keyword_set(long(getenv('READ_ONLY_FIRST_RECORD')))
41
42;#########################################################################
43;######################  STANDARD PLOTS   ################################
44;#########################################################################
45
46  IF keyword_set(doplot) EQ 0 THEN doplot = 0
47
48; fixed color tabled
49  lct, 64
50  cnt = 0
51  htmltxt = ''
52;
53  cnt = cnt+1   &   blabla = 'Global Mean Temperature'
54  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_T, masknp, POSTSCRIPT = postscript, _extra = ex
55;
56  cnt = cnt+1   &   blabla = 'Global Mean Salinity'
57  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_S, masknp, POSTSCRIPT = postscript, _extra = ex
58;
59  cnt = cnt+1   &   blabla = 'Global Mean SSH'
60  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SSH, masknp, POSTSCRIPT = postscript, _extra = ex
61;
62  cnt = cnt+1   &   blabla = 'Global Mean Q net'
63  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Q, masknp, POSTSCRIPT = postscript, _extra = ex
64;
65  cnt = cnt+1   &   blabla = 'Global Mean EMP'
66  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_EMP, masknp, POSTSCRIPT = postscript, _extra = ex
67;
68  cnt = cnt+1   &   blabla = 'Drake Transport'
69  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Drake, masknp, POSTSCRIPT = postscript, _extra = ex
70;
71  cnt = cnt+1   &   blabla = 'Max AMOC'
72  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_AMOC, masknp, POSTSCRIPT = postscript, _extra = ex
73;
74  cnt = cnt+1   &   blabla = 'Sea-Ice cover'
75  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE, masknp, POSTSCRIPT = postscript, _extra = ex
76;
77  cnt = cnt+1   &   blabla = 'Sea-Ice Volume'
78  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex
79
80;SF  cnt = cnt+1   &   blabla = 'SNOW Volume'
81;SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SNOW_Vol, masknp, POSTSCRIPT = postscript, _extra = ex
82;
83  cnt = cnt+1   &   blabla = 'Sea-Ice Volume Export at Fram Strait'
84  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_FRAM, masknp, POSTSCRIPT = postscript, _extra = ex
85;
86  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_ts_html_body.txt', htmltxt[1:*]
87
88  return
89END
Note: See TracBrowser for help on using the repository browser.