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_IceFrac.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_IceFrac.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: 2.2 KB
Line 
1pro std_plot_IceFrac, 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 = 'IceFra'
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 = Ifra1.unit
32;
33  title = var+'!C'+std_file1_I
34  Ifra1.arr = Ifra1.arr * ( Ifra1.arr gt 0.15 )
35  plt, Ifra1.arr - 1.E-04, MIN = 0., MAX = 1., INTER = .1, /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $
36       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $
37       , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex
38;                                ;
39  if std_file1_I NE std_file2_I then begin           
40    title = title + std_file2_I
41    Ifra2.arr = Ifra2.arr * ( Ifra2.arr gt 0.15 )
42    plt, Ifra1.arr - Ifra2.arr, MIN = -1, MAX = 1, INTER = .1, STYLE = 'so0so', format = '(f4.1)' $
43         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $
44         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex
45;;SF  finire, aggiungere lettura obs
46;;SF  endif else begin
47;;SF    title = title+ ' - Observations'
48;;SF    plt, Ifra1.arr - Ifra_obs.arr, MIN = -1., MAX = 1., INTER = 0.1, STYLE = 'so0so', format = '(f4.1)' $
49;;SF         , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, CELL_FILL = 2, TITLE = title $
50;;SF         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex
51;;SF  endelse 
52endif
53
54  domdef
55
56  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
57  if KEYWORD_SET(postscript) then closeps
58
59  return
60end
Note: See TracBrowser for help on using the repository browser.