source: trunk/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_S.pro @ 2150

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

commit idl tools for post-processing, for ORCA2_LIM configuration, see ticket #724

File size: 3.9 KB
Line 
1pro ts_mean_S, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
2
3  compile_opt idl2, strictarrsubs
4 
5  @common
6 
7  @initorca2
8 
9  std_iodir_climato = '/Users/sflod/idl_PLOTS/CLIMATOLOGIES/'
10  std_iodir_interan = '/Users/sflod/idl_PLOTS/INTERAN/'
11
12
13  if KEYWORD_SET(POSTSCRIPT) then begin
14     openps, sEXP1+'_Salinity_mean.ps', portrait = 1
15  endif
16 
17  filename = std_iodir_climato+'data_1y_salinity_nomask.nc'
18  Lev = read_ncdf('vosaline', filename = filename, /nostruct )
19  Lev_xyz = moyenne(Lev, 'xyz')
20  Lev_z = moyenne(Lev, 'xy',/KEEPBOTTOM)
21
22  ts_Sal = 0.
23  ts_Sal_2 = 0.
24  ts_Sal_Lev = 0.
25  ts_Sal_Lev_2 = 0.
26  ts_Time = 0.
27  ts_z = fltarr(jpk)
28  ts_z_2 = fltarr(jpk)
29
30  for year=year1,year2 do begin
31     
32     date1 = year * 10000L + 101
33     date2 = year * 10000L + 1231
34     filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1)
35     filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1)
36
37     if filename eq '' then stop
38     Sal   = read_ncdf('so',date1 ,date2 , filename = filename, /nostruct)
39     Sal_2 = read_ncdf('so',date1 ,date2 , filename = filename2, /nostruct)
40     SSH   =  read_ncdf('zos',date1 ,date2 , filename = filename, /nostruct)
41     SSH_2 =  read_ncdf('zos',date1 ,date2 , filename = filename2, /nostruct)
42
43     ts_Sal = [ ts_Sal, moyenne(Sal, 'xyz', ssh = SSH) ]
44     ts_Sal_2 = [ ts_Sal_2, moyenne(Sal_2, 'xyz', ssh = SSH_2) ]
45     ts_Sal_Lev = [ ts_Sal_Lev, moyenne(Sal, 'xyz', ssh = SSH) - Lev_xyz ]
46     ts_Sal_Lev_2 = [ ts_Sal_Lev_2, moyenne(Sal_2, 'xyz', ssh = SSH_2) - Lev_xyz ]
47     ts_Time = [ ts_Time, Time] ; Time = date
48
49     ts_z = [ [ts_z], [moyenne(Sal, 'xy',/KEEPBOTTOM, ssh = SSH)- Lev_z] ]
50     ts_z_2 = [ [ts_z_2], [moyenne(Sal_2, 'xy',/KEEPBOTTOM, ssh = SSH_2)- Lev_z] ]
51
52  endfor
53 
54  ts_Sal = ts_Sal[1:*]             ; remove first record of 0
55  ts_Sal_2 = ts_Sal_2[1:*]         ; remove first record of 0
56  ts_Sal_Lev = ts_Sal_Lev[1:*]     ; remove first record of 0
57  ts_Sal_Lev_2 = ts_Sal_Lev_2[1:*] ; remove first record of 0
58 
59  ts_z = ts_z[*,1:*]            ; remove first record of 0
60  ts_z_2 = ts_z_2[*,1:*]            ; remove first record of 0
61 
62  Time = ts_Time[1:*]           ; remove first record of 0
63  jpt = n_elements(Time)
64
65
66  title = sEXP1+' Global Annual Mean Salinity --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
67  pltt, ts_Sal, 't', 34.54, 34.76, INTER = .2, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
68       , small = [1,2,1], YTITLE='psu', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
69
70  title = sEXP1+' Global Annual Mean Salinity - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
71  pltt, ts_Sal_Lev, 't', -.05, .05, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
72       , small = [1,2,2], YTITLE='psu',  TITLE = title, FORMAT = '(I3)', /NOERASE, _extra = ex
73
74;;;;;;;;;;;;;;;;; new page ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
75ERASE
76 
77  title = sEXP1+' ZT-plot Salinity Mean (0-1500m) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
78  pltt,ts_z ,'zt', -.5, .5, INTER = .025, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
79       , small = [1,2,1], YTITLE='psu', STYLE = 'so0so',  TITLE = title, boxzoom = 1500., /PORTRAIT, _extra = ex
80
81  title = sEXP1+' ZT-plot Salinity Mean (0-6000) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
82  pltt,ts_z,'zt', -.5, .5, INTER = .025, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
83       , small = [1,2,2], YTITLE='psu', STYLE = 'so0so',  TITLE = title, boxzoom = 6000., zoom = 1500.,  /NOERASE, _extra = ex
84 
85
86if sEXP2 ne sEXP1 then begin
87
88ERASE
89
90   title = sEXP1+' - '+sEXP2+' Global Annual Mean Salinity --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
91   pltt, ts_Sal - ts_Sal_2, 't', -1., 1., INTER = .2, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
92         , COLOR = 250, small = [1,2,1], YTITLE='psu', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
93   
94endif
95
96 
97  if KEYWORD_SET(POSTSCRIPT) then begin
98     closeps
99  endif
100 
101 
102  return
103end
Note: See TracBrowser for help on using the repository browser.