source: trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_T.pro @ 2709

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

update idl scripts, see ticket 724

File size: 4.2 KB
Line 
1pro ts_mean_T, 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  std_iodir_data = '/Users/sflod/idl_PLOTS/INTERAN/mercure/'
12
13
14  if KEYWORD_SET(POSTSCRIPT) then begin
15     openps, sEXP1+'_Temperature_mean.ps', portrait = 1
16  endif
17 
18  masknp = read_ncdf('tmaskutil', file = 'mesh_mask_partial.nc', /nostruct)
19
20  filename = std_iodir_climato+'data_1y_potential_temperature_nomask.nc'
21  Lev = read_ncdf('votemper', filename = filename, /nostruct )
22  Lev_xyz = moyenne(Lev, 'xyz', mask2d = masknp)
23  Lev_z = moyenne(Lev, 'xy', mask2d = masknp,/KEEPBOTTOM)
24
25  ts_Temp = 0.
26  ts_Temp_2 = 0.
27  ts_Temp_Lev = 0.
28  ts_Temp_Lev_2 = 0.
29  ts_Time = 0.
30  ts_z = fltarr(jpk)
31  ts_z_2 = fltarr(jpk)
32 
33  for year=year1,year2 do begin
34     
35     date1 = year * 10000L + 101L
36     date2 = year * 10000L + 1231L
37 
38print,s_iodir_data
39
40     filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1)
41     filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1)
42
43     if filename eq '' then stop
44     Temp   = read_ncdf('thetao',date1 ,date2 , filename = filename, /nostruct)
45     Temp_2 = read_ncdf('thetao',date1 ,date2 , filename = filename2, /nostruct)
46     SSH    =  read_ncdf('zos',date1 ,date2 , filename = filename, /nostruct)
47     SSH_2  =  read_ncdf('zos',date1 ,date2 , filename = filename2, /nostruct)
48
49     ts_Temp = [ ts_Temp, moyenne(Temp, 'xyz', mask2d = masknp,ssh = SSH) ]
50     ts_Temp_2 = [ ts_Temp_2, moyenne(Temp_2, 'xyz', mask2d = masknp,ssh = SSH_2) ]
51     ts_Temp_Lev = [ ts_Temp_Lev, moyenne(Temp, 'xyz', mask2d = masknp,ssh = SSH) - Lev_xyz ]
52     ts_Temp_Lev_2 = [ ts_Temp_Lev_2, moyenne(Temp_2, 'xyz', mask2d = masknp,ssh = SSH_2) - Lev_xyz ]
53     ts_Time = [ ts_Time, Time] ; Time = date
54
55     ts_z = [ [ts_z], [moyenne(Temp, 'xy', mask2d = masknp,/KEEPBOTTOM, ssh = SSH)- Lev_z] ]
56     ts_z_2 = [ [ts_z_2], [moyenne(Temp_2, 'xy', mask2d = masknp,/KEEPBOTTOM, ssh = SSH_2)- Lev_z] ]
57
58  endfor
59
60  ts_Temp = ts_Temp[1:*]        ; remove first record of 0
61  ts_Temp_2 = ts_Temp_2[1:*]        ; remove first record of 0
62  ts_Temp_Lev = ts_Temp_Lev[1:*] ; remove first record of 0
63  ts_Temp_Lev_2 = ts_Temp_Lev_2[1:*] ; remove first record of 0
64  ts_z = ts_z[*,1:*]             ; remove first record of 0
65  ts_z_2 = ts_z_2[*,1:*]             ; remove first record of 0
66 
67  Time = ts_Time[1:*]           ; remove first record of 0
68  jpt = n_elements(Time)
69 
70  title = sEXP1+' Global Annual Mean Temperature --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
71  pltt, ts_Temp, 't', 1., 4.5, INTER = .5, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
72       , small = [1,2,1], YTITLE='Deg C', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
73
74  title = sEXP1+' Global Annual Mean Temperature - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
75  pltt,ts_Temp_Lev, 't', -1., 1., INTER = .25,  year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
76       , small = [1,2,2], YTITLE='Deg C', TITLE = title, FORMAT = '(I3)', /NOERASE, _extra = ex
77
78  ;;;;;;;;;;;;;;;;; new page  ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
79
80ERASE
81
82  title = sEXP1+' ZT-plot Temperature Mean (0-1500m) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
83  pltt,ts_z ,'zt', -2., 2., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
84       , small = [1,2,1], YTITLE='Deg C', TITLE = title, boxzoom = 1500., /NOERASE, _extra = ex
85 
86  title = sEXP1+' ZT-plot Temperature Mean (0-6000m) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
87  pltt,ts_z ,'zt', -2., 2., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
88       , small = [1,2,2], YTITLE='Deg C', TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex
89 
90if sEXP2 ne sEXP1 then begin
91 
92ERASE
93 
94   title = sEXP1+' - '+sEXP2+' Global Annual Mean Temperature --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
95   pltt, ts_Temp - ts_Temp_2, 't', -1., 1., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
96         , COLOR = 250, small = [1,2,1], YTITLE='Deg C', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
97   
98endif
99
100 
101  if KEYWORD_SET(POSTSCRIPT) then begin
102     closeps
103  endif
104 
105 
106  return
107end
Note: See TracBrowser for help on using the repository browser.