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