source: branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro @ 4320

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

update IDL_scripts for ORCA2_LIM and ORCA2_LIM3, add plot of max_mld output 5days, see ticket: #724

File size: 2.0 KB
Line 
1pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex
2
3  compile_opt idl2, strictarrsubs
4
5@common 
6@std_common
7
8  cdti3 = string(cnt, format = '(i3.3)')
9  print, cdti3 + ') ' + blabla
10
11  var = 'SnowThi'
12  IF keyword_set(arc) THEN var = var+'_Arc_'
13  IF keyword_set(ant) THEN var = var+'_Ant_'
14  IF keyword_set(april) THEN var = var+'April'
15  IF keyword_set(sept) THEN var = var+'Sept'
16
17  filename = cdti3 + '_'+var+'_'+std_file1_I
18  if std_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I
19  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
20 
21  IF keyword_set(arc) THEN BEGIN
22    domdef, 20, 380, 50, 90
23    vmin = 0.
24    vmax = 2.
25    vint = 0.1
26    ;fmt = '(i2)'
27    fmt = '(f4.1)'
28    vmind = -1.
29    vmaxd = 1.
30    vintd = 0.1
31    ;fmtd = '(i2)'
32    fmtd = '(f4.1)'
33    map = [90, 0, 0]
34  ENDIF
35  IF keyword_set(ant) THEN BEGIN
36    domdef, 20, 380, -90, -50
37    vmin = 0.
38    vmax = 2.
39    vint = 0.1
40    fmt = '(f4.1)'
41    vmind = -1.
42    vmaxd = 1.
43    vintd = 0.1
44    fmtd = '(f4.1)'
45    map = [-90, 0, 0]
46  ENDIF
47;
48  varunit = Isnow1.unit
49;
50  title = var+'!C'+std_file1_I
51  plt, (Isnow1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $
52       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $
53       , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex
54;                                ;
55  if std_file1_I NE std_file2_I then begin           
56    title = title + std_file2_I
57    plt, Isnow1.arr - Isnow2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $
58         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $
59         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex
60  endif                       
61 
62  domdef
63
64  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
65  if KEYWORD_SET(postscript) then closeps
66
67  return
68end
Note: See TracBrowser for help on using the repository browser.