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_IceSal.pro in branches/2014/dev_r4650_UKMO14.11_SETTE_OBSASM/NEMOGCM/CONFIG/ORCA2_LIM_OBS/IDL_scripts – NEMO

source: branches/2014/dev_r4650_UKMO14.11_SETTE_OBSASM/NEMOGCM/CONFIG/ORCA2_LIM_OBS/IDL_scripts/std_plot_IceSal.pro @ 4751

Last change on this file since 4751 was 4751, checked in by djlea, 10 years ago

Changes to include an OBS test in SETTE. At the moment this uses an example profile observation.

File size: 1.9 KB
Line 
1pro std_plot_IceSal, Isal1, Isal2, Ifra1, Ifra2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, 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;
11  var = 'IceSal'
12  IF keyword_set(arc) THEN var = var+'_Arc_'
13  IF keyword_set(ant) THEN var = var+'_Ant_'
14  IF keyword_set(febr) THEN var = var+'Febr'
15  IF keyword_set(march) THEN var = var+'March'
16  IF keyword_set(sept) THEN var = var+'Sept'
17
18  filename = cdti3 + '_'+var+'_'+std_file1_I
19  if std_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I
20  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
21 
22  IF keyword_set(arc) THEN BEGIN
23     domdef, 20, 380, 50, 90
24     map = [90, 0, 0]
25  ENDIF
26  IF keyword_set(ant) THEN BEGIN
27     domdef, 20, 380, -90, -50
28     map = [-90, 0, 0]
29  ENDIF
30  ;
31  varunit = Isal1.unit
32  ;
33   title = var+'!C'+std_file1_I
34
35   Isal1.arr = Isal1.arr * ( Ifra1.arr gt 0.15 )
36   plt, Isal1.arr - 1.E-04, MIN = 0., MAX = 10., INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $
37        , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $
38        , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex
39 ;                                ;
40   if std_file1_I NE std_file2_I then begin
41     title = title + std_file2_I
42     
43     Isal2.arr = Isal2.arr * ( Ifra2.arr gt 0.15 )
44     
45     plt, Isal1.arr - Isal2.arr, MIN = -2., MAX = 2., INTER = 1., STYLE = 'so0so', format = '(f4.1)' $
46          , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $
47          , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex
48 endif
49 
50   domdef
51 
52   htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
53   if KEYWORD_SET(postscript) then closeps
54 
55   return
56 end
57
Note: See TracBrowser for help on using the repository browser.