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_IceThick.pro in branches/UKMO/r5518_INGV1_WAVE-coupling/NEMOGCM/CONFIG/ORCA2_SAS_LIM/IDL_scripts – NEMO

source: branches/UKMO/r5518_INGV1_WAVE-coupling/NEMOGCM/CONFIG/ORCA2_SAS_LIM/IDL_scripts/std_plot_IceThick.pro @ 7152

Last change on this file since 7152 was 7152, checked in by jcastill, 7 years ago

Initial implementation of wave coupling branch - INGV wave branch + UKMO wave coupling branch

File size: 1.9 KB
Line 
1pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, 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 = 'IceThi'
12  IF keyword_set(arc) THEN var = var+'_Arc_'
13  IF keyword_set(ant) THEN var = var+'_Ant_'
14  IF keyword_set(march) THEN var = var+'March'
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 = 10.
25    vint = 0.5
26    fmt = '(i2)'
27    vmind = -3.
28    vmaxd = 3.
29    vintd = 0.5
30    fmtd = '(i2)'
31    map = [90, 0, 0]
32  ENDIF
33  IF keyword_set(ant) THEN BEGIN
34    domdef, 20, 380, -90, -50
35    vmin = 0.
36    vmax = 3.
37    vint = 0.2
38    fmt = '(f4.1)'
39    vmind = -1.
40    vmaxd = 1.
41    vintd = 0.2
42    fmtd = '(f4.1)'
43    map = [-90, 0, 0]
44  ENDIF
45;
46  varunit = Ithi1.unit
47;
48  title = var+'!C'+std_file1_I
49  plt, (Ithi1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $
50       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $
51       , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex
52;                                ;
53  if std_file1_I NE std_file2_I then begin           
54    title = title + std_file2_I
55    plt, Ithi1.arr - Ithi2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $
56         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $
57         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex
58  endif                       
59 
60  domdef
61
62  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
63  if KEYWORD_SET(postscript) then closeps
64
65  return
66end
Note: See TracBrowser for help on using the repository browser.