New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
std_plot_SnowThick.pro in NEMO/trunk/cfgs/ORCA2_ICE_PISCES/IDL_scripts – NEMO

source: NEMO/trunk/cfgs/ORCA2_ICE_PISCES/IDL_scripts/std_plot_SnowThick.pro @ 11536

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

commit idl scripts deleted

File size: 2.2 KB
Line 
1pro std_plot_SnowThick, Isnow1, Isnow2, Ifra1, Ifra2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, 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(febr) THEN var = var+'Febr'
15  IF keyword_set(march) THEN var = var+'March'
16  IF keyword_set(sept) THEN var = var+'Sept'
17
18  filename = cdti3 + '_'+var+'_'+std_file1_I
19  if std_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I
20  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
21 
22  IF keyword_set(arc) THEN BEGIN
23    domdef, 20, 380, 50, 90
24    vmin = 0.
25    vmax = .6
26    vint = 0.1
27    ;fmt = '(i2)'
28    fmt = '(f4.1)'
29    vmind = -1.
30    vmaxd = 1.
31    vintd = 0.1
32    ;fmtd = '(i2)'
33    fmtd = '(f4.1)'
34    map = [90, 0, 0]
35  ENDIF
36  IF keyword_set(ant) THEN BEGIN
37    domdef, 20, 380, -90, -50
38    vmin = 0.
39    vmax = 1.
40    vint = 0.1
41    fmt = '(f4.1)'
42    vmind = -1.
43    vmaxd = 1.
44    vintd = 0.1
45    fmtd = '(f4.1)'
46    map = [-90, 0, 0]
47  ENDIF
48;
49  varunit = Isnow1.unit
50;
51  title = var+'!C'+std_file1_I
52 
53  Isnow1.arr = Isnow1.arr * ( Ifra1.arr gt 0.15 )
54 
55  plt, (Isnow1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $
56       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $
57       , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex
58;                                ;
59  if std_file1_I NE std_file2_I then begin           
60    title = title + std_file2_I
61   
62    Isnow2.arr = Isnow2.arr * ( Ifra2.arr gt 0.15 )
63     
64    plt, Isnow1.arr - Isnow2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $
65         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $
66         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex
67  endif                       
68 
69  domdef
70
71  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
72  if KEYWORD_SET(postscript) then closeps
73
74  return
75end
Note: See TracBrowser for help on using the repository browser.