source: trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.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.2 KB
Line 
1pro std_ts_ICE, 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  vICE1 = getenv('VAR1_ICE')   &   prefix = getenv('V1ICE_PREF')    &   suffix = getenv('V1ICE_SUFF')
10; get exp2 info
11  vICE2 = getenv('VAR2_ICE')   &   prefix2 = getenv('V2ICE_PREF')   &   suffix2 = getenv('V2ICE_SUFF')
12;
13  cdti3 = string(cnt, format = '(i3.3)')
14  print, cdti3 + ') ' + blabla
15  filename = cdti3 + '_ts_AMOC_'+prefix
16  if prefix NE prefix2 then filename = filename + '_'+prefix2
17  if KEYWORD_SET(postscript) then openps, filename+'_1.ps', portrait = 1
18;
19  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')'
20;
21  iodir = std_iodir_data
22
23  domdef, 0, jpi-1, 30, 90, /xindex
24  ICE_N = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp)
25  ICE_N.arr = ICE_N.arr * 1.e-12   &   ICE_N.unit = '10^12 m^2'
26  domdef, 0, jpi-1, -90, -30, /xindex
27  ICE_S = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp)
28  ICE_S.arr = ICE_S.arr * 1.e-12   &   ICE_S.unit = '10^12 m^2'
29
30  title = prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.)'
31  pltt, ICE_N, 't', 0., 15., date1, date2, /REMPLI, /PORTRAIT $
32        , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex
33 
34  title = prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)'
35  pltt, ICE_S, 't', 0., 15., date1, date2, /REMPLI, /NOERASE $
36        , small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex
37 
38  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_1.png  />  ' ]
39  if KEYWORD_SET(postscript) then closeps
40
41  if prefix NE prefix2 then BEGIN
42
43    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')'
44
45    tsave = time
46    domdef, 0, jpi-1, 30, 90, /xindex
47    ICE_N2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp)
48    ICE_N2.arr = ICE_N2.arr * 1.e-12   &   ICE_N2.unit = '10^12 m^2'
49    domdef, 0, jpi-1, -90, -30, /xindex
50    ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp)
51    ICE_S2.arr = ICE_S2.arr * 1.e-12   &   ICE_S2.unit = '10^12 m^2'
52    time = tsave   &   IF n_elements(time) NE jpt THEN stop
53
54    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1
55
56    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.)'
57    pltt, ICE_N.arr - ICE_N2.arr, 't', -2., 2., date1, date2, /REMPLI, /PORTRAIT, window = 2 $
58          , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex
59   
60    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)'
61    pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., date1, date2, /REMPLI, /NOERASE $
62          , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex
63
64    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ]
65    if KEYWORD_SET(postscript) then closeps
66
67  endif
68
69  domdef
70 
71
72  return
73end
Note: See TracBrowser for help on using the repository browser.