Changeset 4498 for branches

Timestamp:

2014-02-18T15:24:03+01:00 (10 years ago)

Author:

flavoni

Message:

fix small improvements in idl plots, see ticket #724

Location:

branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts

Files:

• std_plot_mht.pro (modified) (1 diff)
• std_ts_ICE_FRAM.pro (modified) (4 diffs)
• std_ts_ICE_Vol.pro (modified) (2 diffs)

branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_mht.pro

@@ -67 +67 @@
     plt1d, htr2_atl, 'y', COLOR = 50, /ov1d,  _extra = ex
 
-    varunit = 'Tw'
-    title = 'MHT (TeraW) differences!C'+std_file1_T+' - '+std_file2_T
-    plt1d, (htr1-htr2)*1.e3, 'y', min = -250., max = 250., TITLE = title $
+    title = 'MHT (Black) & Atlantic MHT (Blue) '+std_file1_T+' !C (Red) '+std_file2_T
+    plt1d, htr1, 'y', min = -2., max = 2., TITLE = title $
            , SMALL = [2, 2, 3], TICKLEN = 1, XGRIDSTYLE = 2 $
-           , YTITLE = 'TW', YGRIDSTYLE = 2, /NOERASE, _extra = ex
-    plt1d, (htr1_atl - htr2_atl)*1.e3, 'y', COLOR = 50, /ov1d,  _extra = ex
+           , YTITLE = 'PW', YGRIDSTYLE = 2, /NOERASE, _extra = ex
+    plt1d, htr2, 'y', COLOR = 250, /ov1d,  _extra = ex
+    plt1d, htr1_atl, 'y', COLOR = 50, /ov1d,  _extra = ex
+    plt1d, htr2_atl, 'y', COLOR = 250, /ov1d,  _extra = ex
+    ;plt1d, (htr1-htr2)*1.e3, 'y', min = -250., max = 250., TITLE = title $
+    ;       , SMALL = [2, 2, 3], TICKLEN = 1, XGRIDSTYLE = 2 $
+    ;       , YTITLE = 'TW', YGRIDSTYLE = 2, /NOERASE, _extra = ex
+    ;plt1d, (htr1_atl - htr2_atl)*1.e3, 'y', COLOR = 50, /ov1d,  _extra = ex
     
   endif 

branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro

@@ -10 +10 @@
 end 
 
-;
-;
+;; here start procedure that use function read_arr2d
 pro std_ts_ICE_FRAM, masknp, s_iodir_data,  POSTSCRIPT = postscript, _extra = ex
 
@@ -26 +25 @@
   vICE2     = getenv('VAR2_ICE')     &   prefix2 = getenv('V2ICE_PREF')   &   suffix2 = getenv('V2ICE_SUFF')
   v2_Ithick = getenv('VAR2_Ithick')  &   prefix2 = getenv('V2It_PREF')    &   suffix2 = getenv('V2It_SUFF')
-  v2_IV     = getenv('VAR2_IvelV')   &   prefix = getenv('V2IvV_PREF')    &   suffix = getenv('V2IvV_SUFF')
+  v2_IV     = getenv('VAR2_IvelV')   &   prefix2 = getenv('V2IvV_PREF')   &   suffix2 = getenv('V2IvV_SUFF')
 
   cdti3 = string(cnt, format = '(i3.3)')
@@ -101 +100 @@
   ;title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)'
   title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2
-  pltt, ICE_vol_export, 't',  /REMPLI, MIN = 0., MAX = .5 , /NOERASE, XGRIDSTYLE = 1 $
+  pltt, ICE_vol_export, 't',  /REMPLI, MIN = 0., MAX = .6 , /NOERASE, XGRIDSTYLE = 1 $
        , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex
   pltt, vol_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4  $    ; light blue
@@ -171 +170 @@
     xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2
     title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2
-    pltt, ICE_vol_export, 't', /REMPLI,  MIN = 0., MAX = .5, /NOERASE, XGRIDSTYLE = 1 $
+    pltt, ICE_vol_export, 't', /REMPLI,  MIN = 0., MAX = .6, /NOERASE, XGRIDSTYLE = 1 $
         , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex
     pltt, ICE_vol_export_2, 't', /REMPLI, /NOERASE  $

branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro

@@ -65 +65 @@
   jpt=12
   time=julday(1,15,1900)+30*lindgen(12)
-  pltt, ICE_vol_N, 't', MIN = 0., MAX = 30000., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $
+  pltt, ICE_vol_N, 't', MIN = 0., MAX = 40000., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $
         , small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex
 ;
@@ -112 +112 @@
     jpt=12
     time=julday(1,15,1900)+30*lindgen(12)
-    pltt, ICE_vol_N, 't', MIN = 0., MAX = 30000., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $
+    pltt, ICE_vol_N, 't', MIN = 0., MAX = 40000., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $
           , small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex
     pltt, ICE_vol_N2 , 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $