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_mld.pro in branches/2012/dev_r3322_NOCS09_SAS/NEMOGCM/CONFIG/ORCA2_SAS_LIM/IDL_scripts – NEMO

source: branches/2012/dev_r3322_NOCS09_SAS/NEMOGCM/CONFIG/ORCA2_SAS_LIM/IDL_scripts/std_plot_mld.pro @ 3331

Last change on this file since 3331 was 3331, checked in by sga, 12 years ago

NEMO 2012 development branch dev_r3322_NOCS09_SAS
Code changes made for compilation and running of StandAlone? Surface scheme (tinkering still required)

File size: 1.6 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, MIN = 0., MAX = 500., INTER = 25., /NOCONTOUR, FORMAT = '(I3)' $
35         , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex
36  ENDELSE
37;
38  IF keyword_set(MLD2) THEN BEGIN
39    title = titleorg+std_file2_T+ ' - DeBoyer'
40    tmp = MLD2.arr - MLD.arr
41  ENDIF ELSE BEGIN
42    title = titleorg+std_file1_T+ ' - DeBoyer'
43    tmp = MLD1.arr - MLD.arr
44  ENDELSE
45  plt, temporary(tmp), MIN = -80., MAX = 80., INTER = 10., FORMAT = '(I3)'  $
46       , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, /NOCONTOUR, /NOERASE, _extra = ex
47
48  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
49  if KEYWORD_SET(postscript) then closeps
50;
51  return
52end
53
Note: See TracBrowser for help on using the repository browser.