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_IceSal.pro in branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts – NEMO

source: branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceSal.pro @ 4514

Last change on this file since 4514 was 4514, checked in by flavoni, 9 years ago

add ice pots & time series, see ticket #724

File size: 1.7 KB
Line 
1pro std_plot_IceSal, Isal1, Isal2, 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 = 'IceSal'
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     map = [90, 0, 0]
25  ENDIF
26  IF keyword_set(ant) THEN BEGIN
27     domdef, 20, 380, -90, -50
28     map = [-90, 0, 0]
29  ENDIF
30  ;
31  varunit = Isal1.unit
32  ;
33   title = var+'!C'+std_file1_I
34   plt, Isal1.arr, MIN = 0., MAX = 12., INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $
35        , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $
36        , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex
37 ;                                ;
38   if std_file1_I NE std_file2_I then begin
39     title = title + std_file2_I
40     plt, Isal1.arr - Isal2.arr, MIN = -4, MAX = 4, INTER = 1., STYLE = 'so0so', format = '(f4.1)' $
41          , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $
42          , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex
43 endif
44 
45   domdef
46 
47   htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
48   if KEYWORD_SET(postscript) then closeps
49 
50   return
51 end
52
Note: See TracBrowser for help on using the repository browser.