source: branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_ArcSal.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.2 KB
RevLine 
[2751]1pro std_plot_ArcSal, S1, S2in, SLevin, Z100 = z100, POSTSCRIPT = postscript, _extra = ex
[2709]2
[2751]3  compile_opt idl2, strictarrsubs
[2709]4
[2751]5@common 
6@std_common
[2709]7
[2751]8  CASE n_params() OF
9    2:BEGIN
10      Slev = S2in
11    END
12    3:BEGIN
13      IF S2in.arr[0] EQ -1 THEN return
14      S2 = S2in
15      Slev = Slevin
16    END
17  ENDCASE
18;
19  cdti3 = string(cnt, format = '(i3.3)')
20  print, cdti3 + ') ' + blabla
21;
22  IF keyword_set(z100) THEN tmp = min(abs(100 - gdept), ind) ELSE ind = 0
23  sdepref = strtrim(round(gdept[ind]), 1)+'m'
24  filename = cdti3 + '_Arctic_Sal'+sdepref+'_'+std_file1_T
25  IF keyword_set(S2) THEN filename = filename + '_' + std_file2_T
26  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
27;
[4320]28  domdef, 20, 380, 50, 90
[2751]29;
30  varunit = S1.unit
31  titleorg = 'Salinity ('+sdepref+')!C'
32;
33  IF keyword_set(S2) THEN BEGIN
34    title = titleorg+std_file1_T+ ' - '+std_file2_T
35    plt, S1.arr - S2.arr, MIN = -4., MAX = 4., INTER = 0.2, CELL_FILL = 2,  STYLE = 'so0so', format = '(f4.1)'  $     
36         , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $     
37         , /ORTHO, MAP = [90, 0, 0], LATDEL = 5, boxzoom = [ind, ind], /zindex, /portrait, _extra = ex
38  ENDIF ELSE BEGIN
39    title = titleorg+std_file1_T
40    plt, S1, MIN = 30.0, MAX = 36., INTER = 0.2, CELL_FILL = 2, format = '(f4.1)' $
41         , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $
42         , /ORTHO, MAP = [90, 0, 0], LATDEL = 5, boxzoom = [ind, ind], /zindex, /PORTRAIT, _extra = ex
43  ENDELSE
44;
45  IF keyword_set(S2) THEN BEGIN
46    title = titleorg+std_file2_T+ ' - Levitus'
47    tmp = S2.arr - SLev.arr
48  ENDIF ELSE BEGIN
49    title = titleorg+std_file1_T+ ' - Levitus'
50    tmp = S1.arr - SLev.arr
51  ENDELSE
52  plt, temporary(tmp), MIN = -4., MAX = 4., INTER = 0.2, CELL_FILL = 2,  STYLE = 'so0so', format = '(f4.1)'  $     
53       , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $     
54       , /ORTHO, MAP = [90, 0, 0], LATDEL = 5, boxzoom = [ind, ind], /zindex, /NOERASE, _extra = ex
55;
56  domdef
57;
58  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
59  if KEYWORD_SET(postscript) then closeps
60 
61  return
[2709]62end
Note: See TracBrowser for help on using the repository browser.