source: branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_mld.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: 1.8 KB
Line 
1pro std_plot_mld, MLD1, MLD2in, MLDin, POSTSCRIPT = postscript, _extra = ex
2
3  compile_opt idl2, strictarrsubs
4
5@common 
6@std_common
7
8  CASE n_params() OF
9    2:BEGIN
10      MLD = MLD2in
11    END
12    3:BEGIN
13      IF MLD2in.arr[0] EQ -1 THEN return
14      MLD2 = MLD2in
15      MLD = MLDin
16    END
17  ENDCASE
18;
19  cdti3 = string(cnt, format = '(i3.3)')
20  print, cdti3 + ') ' + blabla
21  filename = cdti3 + '_MLD_'+std_file1_T
22  IF keyword_set(MLD2) THEN filename = filename + '_' + std_file2_T
23  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
24;
25  varunit = MLD1.unit
26  titleorg = 'MLD!C'
27;
28  IF keyword_set(MLD2) THEN BEGIN
29    title = titleorg+std_file1_T+ ' - '+std_file2_T
30    plt, MLD1.arr - MLD2.arr, MIN = -80., MAX = 80., INTER = 10., FORMAT = '(I3)'  $
31         , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, /NOCONTOUR, /PORTRAIT, _extra = ex
32  ENDIF ELSE BEGIN
33    title = titleorg+std_file1_T
34    plt, MLD1,label=4, cb_label=[0, 15, 25, 50, 75, 100, 125, 150, 250, 400, 600, 850],/NOCONTOUR,FORMAT = '(I3)' $
35            , small = [1, 2, 1],COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex
36;ORI SF    plt, MLD1, MIN = 0., MAX = 500., INTER = 25., /NOCONTOUR, FORMAT = '(I3)' $
37;ORI SF         , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex
38  ENDELSE
39;
40  IF keyword_set(MLD2) THEN BEGIN
41    title = titleorg+std_file2_T+ ' - DeBoyer'
42    tmp = MLD2.arr - MLD.arr
43  ENDIF ELSE BEGIN
44    title = titleorg+std_file1_T+ ' - DeBoyer'
45    tmp = MLD1.arr - MLD.arr
46  ENDELSE
47  plt, temporary(tmp), MIN = -80., MAX = 80., INTER = 10., FORMAT = '(I3)'  $
48       , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, /NOCONTOUR, /NOERASE, _extra = ex
49
50  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
51  if KEYWORD_SET(postscript) then closeps
52;
53  return
54end
55
Note: See TracBrowser for help on using the repository browser.