source: trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_S.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: 3.7 KB
Line 
1pro std_ts_S, masknp, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
2
3  compile_opt idl2, strictarrsubs
4 
5@common
6@std_common
7
8; get exp1 info
9  vsal = getenv('VAR1_S')      &   prefix = getenv('V1S_PREF')         &   suffix = getenv('V1S_SUFF')
10  vssh = getenv('VAR1_SSH')    &   sshprefix = getenv('V1SSH_PREF')    &   sshsuffix = getenv('V1SSH_SUFF')
11; get exp2 info
12  vsal2 = getenv('VAR2_S')     &   prefix2 = getenv('V2S_PREF')        &   suffix2 = getenv('V2S_SUFF')
13  vssh2 = getenv('VAR2_SSH')   &   sshprefix2 = getenv('V2SSH_PREF')   &   sshsuffix2 = getenv('V2SSH_SUFF')
14;
15  cdti3 = string(cnt, format = '(i3.3)')
16  print, cdti3 + ') ' + blabla
17  filename = cdti3 + '_ts_S_'+prefix
18  if prefix NE prefix2 then filename = filename + '_'+prefix2
19  if KEYWORD_SET(postscript) then openps, filename+'_1.ps', portrait = 1
20
21  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')'
22
23; read levitus data
24  std_file_Levitus_S =  isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato)
25  Lev = read_ncdf(getenv('VAR_SAL_3D'), filename = std_file_Levitus_S, /nostruct )
26  Lev_xyz = moyenne(Lev, 'xyz', mask2d = masknp)
27  levz = moyenne(temporary(Lev), 'xy', mask2d = masknp, /KEEPBOTTOM)
28;
29  iodir = std_iodir_data
30;
31
32; read exp1 data
33  std_ts_read, vsal, date1, date2, prefix, suffix, ts_Sal, ts_z $
34               , WITHSSH = vssh, SSHPREFIX = sshprefix, SSHSUFFIX = sshsuffix, LEVZ = levz
35
36  title = prefix+'!C'+blabla+' '+d1_d2
37  pltt, ts_Sal, 't', 34.54, 34.76, date1, date2, /REMPLI $
38        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex
39
40  title = prefix+' - Levitus!C'+blabla+' '+d1_d2
41  pltt, ts_Sal - Lev_xyz, 't', -.05, .05, date1, date2, /REMPLI $
42        , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /NOERASE, _extra = ex
43
44  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_1.png  />  ' ]
45  if KEYWORD_SET(postscript) then closeps
46  if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1
47
48  title = prefix+' - Levitus!C ZT-plot (0-1500m) '+blabla+' '+d1_d2
49  pltt, ts_z, 'zt', -.5, .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so' $
50        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex
51 
52  title = prefix+' - Levitus!C ZT-plot (0-6000m) '+blabla+' '+d1_d2
53  pltt, ts_z, 'zt', -.5, .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so' $
54        , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex
55 
56  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ]
57  if KEYWORD_SET(postscript) then closeps
58
59  IF prefix NE prefix2 THEN BEGIN
60
61    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')'
62; read exp2 data
63    tsave = time
64    std_ts_read, vsal2, date1_2, date2_2, prefix2, suffix2, ts_Sal2, ts_z2 $
65                 , WITHSSH = vssh2, SSHPREFIX = sshprefix2, SSHSUFFIX = sshsuffix2, LEVZ = levz
66    time = tsave   &   IF n_elements(time) NE jpt THEN stop
67
68    if KEYWORD_SET(postscript) then openps, filename+'_3.ps', portrait = 1
69
70    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla
71    pltt, ts_Sal - ts_Sal2, 't', -0.1, 0.1, date1, date2, /REMPLI $
72          , COLOR = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, window = 2, _extra = ex
73
74    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C ZT-plot (0-6000m) '+blabla
75    pltt, ts_z - ts_z2, 'zt', -0.1, 0.1, INTER = .01, date1, date2, /REMPLI, style = 'so0so' $
76          , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex
77
78    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_3.png  />  ' ]
79    if KEYWORD_SET(postscript) then closeps
80
81  ENDIF
82
83  return
84end
85
Note: See TracBrowser for help on using the repository browser.