source: trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_S100m.pro @ 2760

Last change on this file since 2760 was 2760, checked in by smasson, 10 years ago

bugfix in IDL_scripts

File size: 1.5 KB
Line 
1pro std_plot_S100m, S1, S2, SLev, 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  tmp = min(abs(100 - gdept), ind)
11  sdepref = strtrim(round(gdept[ind]), 1)+'m'
12  filename = cdti3 + '_S'+sdepref+'_'+std_file1_T
13  if std_file1_T NE std_file2_T then filename = filename + '_'+std_file2_T
14  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
15;
16  title = 'Salinity ('+sdepref+')!C'+std_file1_T
17  plt, S1, MIN = 33., MAX = 38.2, INTER = .2, format = '(f4.1)' $
18       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, boxzoom = [floor(gdept[ind]), ceil(gdept[ind])], /PORTRAIT, _extra = ex
19;
20  if std_file1_T NE std_file2_T then begin
21    title = title+ ' - '+std_file2_T
22    plt, S1.arr[*, *, 0] - S2.arr[*, *, 0], MIN = -1., MAX = 1., INTER = .1, STYLE = 'so0so', format = '(f4.1)' $
23         , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, boxzoom = [floor(gdept[ind]), ceil(gdept[ind])], /noerase, _extra = ex   
24  endif else begin
25    title = title+ ' - Levitus'
26    plt, S1.arr-SLev.arr, MIN = -1., MAX = 1., INTER = 0.1, STYLE = 'so0so', format = '(f4.1)' $
27         , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, boxzoom = [floor(gdept[ind]), ceil(gdept[ind])], /noerase, _extra = ex
28  endelse
29
30  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
31  if KEYWORD_SET(postscript) then closeps
32
33  return
34end
Note: See TracBrowser for help on using the repository browser.