Changeset 4497


Ignore:
Timestamp:
2014-02-17T15:42:08+01:00 (7 years ago)
Author:
flavoni
Message:

update idl scripts, see ticket #724

Location:
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
Files:
10 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_AMOC.pro

    r4495 r4497  
    9595 
    9696  title = prefix+' '+d1_d2+'!C'+'Max Atlantic MOC between 40N and 50N' 
    97   pltt, t45, 't', 0., 30., date1, date2, /REMPLI, /PORTRAIT $ 
     97  pltt, t45, 't', MIN = 0., MAX = 30., date1, date2, /REMPLI, /PORTRAIT, XGRIDSTYLE = 1 $ 
    9898        ,  small = [1, 2, 1],  TITLE = title, YTITLE = varunit, _extra = ex 
    9999  IF prefix NE prefix2 then begin 
    100100    title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Max Atlantic MOC between 40N and 50N' 
    101     pltt, t45, 't', 0., 30., date1, date2, /REMPLI, /NOERASE $ 
     101    pltt, t45, 't', MIN = 0., MAX = 30., date1, date2, /REMPLI , /NOERASE, XGRIDSTYLE = 1 $ 
    102102          ,  small = [1, 2, 2],  TITLE = title, YTITLE = varunit, _extra = ex ; BLACK 
    103     pltt, t45_2, 't', 0., 30., date1, date2, /REMPLI, /NOERASE $ 
     103    pltt, t45_2, 't', date1, date2, /REMPLI, /NOERASE $ 
    104104          ,  /ov1d, COLOR = 250, small = [1, 2, 2],  TITLE = title, YTITLE = varunit, _extra = ex ; RED 
    105105  endif 
     
    110110 
    111111  title = prefix+' '+d1_d2+'!C'+'Max AntArctic Bottom Water between 80S and 65S '+d1_d2 
    112   pltt, -t70, 't', 0., 20., date1, date2, /REMPLI, /PORTRAIT, window = 2 $ 
     112  pltt, -t70, 't', MIN = 0., MAX = 20., date1, date2, /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2 $ 
    113113        , small = [1, 2, 1],  TITLE = title, YTITLE = varunit, _extra = ex 
    114114  if prefix NE prefix2 then begin  
    115115    title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Max AntArctic Bottom Water between 80S and 65S' 
    116     pltt, -t70, 't', 0., 20., date1, date2, /REMPLI, /NOERASE $ 
     116    pltt, -t70, 't', MIN = 0.,MAX =  20., date1, date2, /REMPLI, /NOERASE, XGRIDSTYLE = 1 $ 
    117117          ,  small = [1, 2, 2],  TITLE = title, YTITLE = varunit, _extra = ex ; BLACK 
    118     pltt, -t70_2, 't', 0., 20., date1, date2, /REMPLI, /NOERASE $ 
     118    pltt, -t70_2, 't', date1, date2, /REMPLI, /NOERASE $ 
    119119          , /ov1d, COLOR = 250, small = [1, 2, 2],  TITLE = title, YTITLE = varunit, _extra = ex ; RED 
    120120  endif 
     
    125125 
    126126  title = prefix+' '+d1_d2+'!C '+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2 
    127   pltt, -t50, 't', 5., 30., date1, date2, /REMPLI, /PORTRAIT, window = 3 $ 
     127  pltt, -t50, 't', MIN = 5., MAX = 30., date1, date2, /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 3 $ 
    128128        ,  small = [1, 2, 1],  TITLE = title, YTITLE = varunit, _extra = ex 
    129129  if prefix NE prefix2 then begin 
    130130    title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Max AntArctic Bottom Cell between 65S and 30N' 
    131     pltt, -t50 , 't', 5., 30., date1, date2, /REMPLI, /NOERASE $ 
     131    pltt, -t50 , 't', MIN = 5., MAX = 30., date1, date2, /REMPLI, /NOERASE, XGRIDSTYLE = 1 $ 
    132132          ,  small = [1, 2, 2],  TITLE = title, YTITLE = varunit, _extra = ex ; BLACK 
    133     pltt, -t50_2, 't', 5., 30., date1, date2, /REMPLI, /NOERASE  $ 
     133    pltt, -t50_2, 't', date1, date2, /REMPLI, /NOERASE  $ 
    134134          , /ov1d, COLOR = 250, small = [1, 2, 2],  TITLE = title, YTITLE = varunit, _extra = ex ; RED 
    135135  endif 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_Drake.pro

    r4495 r4497  
    4141; 
    4242  title = prefix+' '+d1_d2+'!C'+blabla 
    43   pltt, drk1, 't', 0., 200., date1, date2, /REMPLI $ 
     43  pltt, drk1, 't', 0., 200., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    4444        , small = [1, 2, 1],  TITLE = title, YTITLE = varunit, /PORTRAIT, _extra = ex 
    4545   
     
    6060 
    6161    title = prefix+' (BLACK) - '+prefix2+' (RED)  '+d1_d2_2+'!C'+blabla 
    62     pltt, drk1 , 't', 0., 200., date1, date2, /REMPLI, /NOERASE $ 
     62    pltt, drk1 , 't', 0., 200., date1, date2, /REMPLI, /NOERASE, XGRIDSTYLE = 1 $ 
    6363          , COLOR = 000, small = [1, 2, 2],  TITLE = title, YTITLE = varunit, _extra = ex ; BLACK 
    6464    pltt, drk2, 't', 0., 200., date1, date2, /REMPLI, /NOERASE $ 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_EMP.pro

    r4320 r4497  
    2727  ts_EMP.arr = ts_EMP.arr * ( 1.E-09 * surf_oce )   &   ts_EMP.unit = 'Sv' 
    2828  title = prefix+' '+d1_d2+'!C'+blabla 
    29   pltt, ts_EMP, 't', -2., 2., date1, date2, /REMPLI $ 
     29  pltt, ts_EMP, 't', MIN = -2.,MAX = 2., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    3030        , small = [1, 2, 1], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex 
    3131   
     
    3838    time = tsave   &   IF n_elements(time) NE jpt THEN stop 
    3939 
    40     title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 
    41     pltt, ts_EMP.arr - ts_EMP2.arr, 't', -2., 2., date1, date2, /REMPLI $ 
     40    title = prefix+' (BLACK) - ' +prefix2+' (RED) '+d1_d2_2+'!C'+blabla 
     41    pltt, ts_EMP.arr - ts_EMP2.arr, 't', MIN = -2., MAX = 2., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    4242          , color = 250, small = [1, 2, 2], YTITLE = 'Sv', TITLE = title, /noerase, _extra = ex 
    4343 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro

    r4495 r4497  
    9191  jpt=12 
    9292  time=julday(1,15,1900)+30*lindgen(12) 
    93   pltt, ICE_N, 't', 0., 15., /REMPLI, /PORTRAIT, MIN = 4., MAX = 16. $ 
     93  pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    9494       , COLOR = 000 , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    95   pltt, ICE_N_15, 't', 0., 15., /REMPLI, /PORTRAIT, LINESTYLE=2 $       ;;; dashed lines is LINESTYLE=2  $ 
     95  pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $       ;;; dashed lines is LINESTYLE=2  $ 
    9696        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    97   pltt, vICE_area_NH, 't', 0., 15.,  /REMPLI, /PORTRAIT $  
     97  pltt, vICE_area_NH, 't',  /REMPLI, /PORTRAIT $  
    9898        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    99   pltt, vICE_ext_NH, 't', 0., 15.,  /REMPLI, /PORTRAIT, LINESTYLE=2 $   ;;; dashed lines is LINESTYLE=2  $ 
     99  pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $   ;;; dashed lines is LINESTYLE=2  $ 
    100100        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    101101; 
    102102  title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
    103   pltt, ICE_S, 't', 0., 15., /REMPLI, /NOERASE , MIN = 0., MAX = 20. $ 
     103  pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
    104104        ,COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    105   pltt, ICE_S_15, 't', 0., 15., /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     105  pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    106106        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    107   pltt,  vICE_area_SH, 't', 0., 15., /REMPLI, /PORTRAIT $  
     107  pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT $  
    108108         , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    109   pltt,  vICE_ext_SH, 't', 0., 15., /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     109  pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    110110        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    111111; 
     
    162162  jpt=12 
    163163  time=julday(1,15,1900)+30*lindgen(12) 
    164   pltt, ICE_N, 't', -.2, .2,  /REMPLI, /PORTRAIT, MIN = 4., MAX = 16., window = 2 $ 
     164  pltt, ICE_N, 't', MIN = 4, MAX = 16,  /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 
    165165        , COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; BLACK 
    166   pltt, ICE_N2, 't', -.2, .2,  /REMPLI, /PORTRAIT $ 
     166  pltt, ICE_N2, 't', /REMPLI, /PORTRAIT $ 
    167167        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; RED 
    168   pltt, ICE_N_15, 't', -.2, .2,  /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     168  pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
    169169        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    170   pltt, ICE_N2_15, 't', -.2, .2, /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     170  pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
    171171        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    172   pltt, vICE_area_NH, 't', -2., 2.,  /REMPLI, /PORTRAIT $  
     172  pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT $  
    173173        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; light blue 
    174   pltt, vICE_ext_NH, 't', -2., 2.,  /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     174  pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    175175        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; blu scuro 
    176176; 
    177177  title ='Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+'Global Annual Mean Ice Area (CONTINUS)'+'!C'+ 'and Extend minus 15% (DASHED)' 
    178178;  title ='Southern Hemisphere'+'!C' 
    179   pltt, ICE_S.arr, 't', -2., 2.,  /REMPLI, /NOERASE, MIN = 0., MAX = 20. $ 
     179  pltt, ICE_S, 't', MIN = 0., MAX = 20.,  /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    180180         , COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    181   pltt, ICE_S2.arr, 't', -2., 2., /REMPLI, /NOERASE $ 
     181  pltt, ICE_S2, 't', /REMPLI, /NOERASE $ 
    182182        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    183   pltt, ICE_S_15.arr , 't', -2., 2., /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     183  pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
    184184        , /ov1d, COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    185   pltt, ICE_S2_15.arr , 't', -2., 2., /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     185  pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
    186186        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    187   pltt,  vICE_area_SH, 't', -2., 2.,  /REMPLI, /PORTRAIT $  
     187  pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT $  
    188188          , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    189   pltt,  vICE_ext_SH, 't', -2., 2.,  /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     189  pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    190190        , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    191191; 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro

    r4495 r4497  
    4949  indy2= indy1  
    5050 ; 
     51 ;OBSERVATIONS : mean seasonal cycle/month 
     52  vol_obs = [0.261625, 0.230750, 0.325375, 0.252000, 0.172500, 0.0805000, 0.0805000, 0.0805000, 0.0805000, 0.176500, 0.148500, 0.235000] 
     53  area_obs = [0.103292, 0.0997500, 0.107625, 0.0944167, 0.0612083, 0.0262500, 0.0262500, 0.0262500, 0.0262500, 0.0843750, 0.0914583, 0.104083] 
     54 ; 
    5155  domdef, indx1, indx2, indy1, indy2, /xindex, /yindex,/memeindices 
    5256  ICE = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, /nostruct) 
     
    8589  jpt=12 
    8690  time=julday(1,15,1900)+30*lindgen(12) 
    87   ;; pltt, ICE_area_export, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 0., MAX = .5 $ 
    88   pltt, ICE_area_export, 't', 0., 15., /REMPLI, /PORTRAIT,MIN = 0., MAX = .5 $ 
    89   ;;      , small = [1, 2, 1], YTITLE = varunit, TITLE = title, _extra = ex 
    90         , small = [1, 2, 1],YTITLE = '10^6 Km^2/month',  TITLE = title, _extra = ex     
    91   ; 
     91  pltt, ICE_area_export, 't', /REMPLI, /PORTRAIT, MIN = 0., MAX = .5 , XGRIDSTYLE = 1 $ 
     92        , small = [1, 2, 1],YTITLE = '10^6 Km^2/month',  TITLE = title, DATE_FORMAT = '%M', _extra = ex     
     93  pltt, area_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4  $    ; light blue 
     94         , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^6 Km^2/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
     95  ; 
     96  tot_area_expo = total(ICE_area_export.arr) 
     97  ; 
     98  xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     99  xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     100  xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    92101  ;title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 
    93102  title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 
    94   ;; pltt, ICE_vol_export, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = .5 $ 
    95   pltt, ICE_vol_export, 't', 0., 15., /REMPLI, /NOERASE , MIN = 0., MAX = .5 $ 
    96   ;;       , small = [1, 2, 2], YTITLE = varunit, TITLE = title, _extra = ex 
    97        , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, _extra = ex 
    98   ; 
     103  pltt, ICE_vol_export, 't',  /REMPLI, MIN = 0., MAX = .5 , /NOERASE, XGRIDSTYLE = 1 $ 
     104       , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
     105  pltt, vol_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4  $    ; light blue 
     106         , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^6 Km^2/month',DATE_FORMAT = '%M', TITLE = title, _extra = ex 
     107  ; 
     108  tot_vol_expo = total(ICE_vol_export.arr) 
     109  tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 
     110  xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     111  xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_vol_expo, 1)+' Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     112  xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     113  xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
     114 
    99115  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
    100116  if KEYWORD_SET(postscript) then closeps 
     
    140156    jpt=12 
    141157    time=julday(1,15,1900)+30*lindgen(12) 
    142     ;; pltt, ICE_area_export.arr, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, MIN = 0., MAX = .5, window = 2 $ 
    143     pltt, ICE_area_export, 't', 0., 15., /REMPLI, /PORTRAIT, MIN = 0., MAX = .5, window = 2 $ 
    144     ;;    , small = [1, 2, 1], YTITLE = varunit, TITLE = title, _extra = ex 
    145         , small = [1, 2, 1], YTITLE = '10^6 Km^2/month',  TITLE = title, _extra = ex     
    146     ;; pltt, ICE_area_export_2.arr ,'t', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ 
    147     pltt, ICE_area_export_2 ,'t', 0., 15., /REMPLI, /PORTRAIT $ 
    148     ;;    , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = varunit,  TITLE = title, _extra = ex 
    149         , /ov1d, COLOR = 250, small = [1, 2, 1],YTITLE = '10^6 Km^2/month',  TITLE = title, _extra = ex 
    150     ; 
     158    pltt, ICE_area_export, 't', /REMPLI, /PORTRAIT, MIN = 0., MAX = .5, XGRIDSTYLE = 1, window = 2 $ 
     159        , small = [1, 2, 1], YTITLE = '10^6 Km^2/month',  TITLE = title, DATE_FORMAT = '%M', _extra = ex     
     160    pltt, ICE_area_export_2 ,'t', /REMPLI, /PORTRAIT, /NOERASE $ 
     161        , /ov1d, COLOR = 250, small = [1, 2, 1],YTITLE = '10^6 Km^2/month',  TITLE = title, DATE_FORMAT = '%M', _extra = ex 
     162    pltt, area_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4  $    ; light blue 
     163         , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^6 Km^2/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
     164    ; 
     165    tot_area_expo = total(ICE_area_export.arr) 
     166    tot_area_expo_2 = total(ICE_area_export_2.arr) 
     167    ; 
     168    xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     169    xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model 1= '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     170    xyouts, julday(6,15,1900), 0.44, 'Tot. Annual Export Model 2= '+strtrim(tot_area_expo_2, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     171    xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    151172    title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
    152     ;; pltt, ICE_vol_export.arr, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = .5 $ 
    153     pltt, ICE_vol_export, 't', 0., 15., /REMPLI, /NOERASE , MIN = 0., MAX = .5 $ 
    154     ;;    , small = [1, 2, 2], YTITLE = varunit, TITLE = title, _extra = ex 
    155         , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, _extra = ex 
    156     ;; pltt, ICE_vol_export_2.arr, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE  $ 
    157     pltt, ICE_vol_export_2, 't', 0., 15.,  /REMPLI, /NOERASE  $ 
    158     ;;    , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = varunit, TITLE = title, _extra = ex 
    159         , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, _extra = ex 
    160     ; 
     173    pltt, ICE_vol_export, 't', /REMPLI,  MIN = 0., MAX = .5, /NOERASE, XGRIDSTYLE = 1 $ 
     174        , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
     175    pltt, ICE_vol_export_2, 't', /REMPLI, /NOERASE  $ 
     176        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M',  _extra = ex 
     177    pltt, vol_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4  $    ; light blue 
     178         , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^6 Km^2/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
     179    ; 
     180    tot_vol_expo = total(ICE_vol_export.arr) 
     181    tot_vol_expo_2 = total(ICE_vol_export_2.arr) 
     182    tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 
     183    tot_vol_expo_2_Sv = tot_vol_expo_2 * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 
     184    xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     185    xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model 1 = '+strtrim(tot_vol_expo, 1)+' Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     186    xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     187    xyouts, julday(6,15,1900), 0.42, 'Tot. Annual Export Model 2 = '+strtrim(tot_vol_expo_2, 1)+' Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     188    xyouts, julday(9,15,1900), 0.40, 'in Sv = '+strtrim(tot_vol_expo_2_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     189    xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    161190 
    162191    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ] 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro

    r4495 r4497  
    6565  jpt=12 
    6666  time=julday(1,15,1900)+30*lindgen(12) 
    67   pltt, ICE_vol_N, 't', 0., 15., /REMPLI, /PORTRAIT,MIN = 0., MAX = 30000. $ 
     67  pltt, ICE_vol_N, 't', MIN = 0., MAX = 30000., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    6868        , small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
    6969; 
    7070  title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (BLACK)' 
    71   pltt, ICE_vol_S, 't', 0., 15., /REMPLI, /NOERASE , MIN = 0., MAX = 11000. $ 
     71  pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
    7272        , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
    7373; 
     
    112112    jpt=12 
    113113    time=julday(1,15,1900)+30*lindgen(12) 
    114     pltt, ICE_vol_N, 't', -125., -60., /REMPLI, /PORTRAIT, MIN = 0., MAX = 30000., window = 2 $ 
     114    pltt, ICE_vol_N, 't', MIN = 0., MAX = 30000., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 
    115115          , small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
    116     pltt, ICE_vol_N2 , 't', -.2, .2, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
     116    pltt, ICE_vol_N2 , 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    117117        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
    118118    ; 
    119119    title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Global Annual Mean Ice Volume' 
    120     pltt, ICE_vol_S, 't', -620., -10., /REMPLI, /NOERASE, MIN = 0., MAX = 11000. $ 
     120    pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    121121          , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex    
    122     pltt, ICE_vol_S2, 't', -620., -10., /REMPLI, /NOERASE $ 
     122    pltt, ICE_vol_S2, 't', /REMPLI, /NOERASE $ 
    123123          , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex    
    124124 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_Q.pro

    r4495 r4497  
    2323  ts_Q = rseries_ncdf(vq, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 
    2424  title = prefix+' '+d1_d2+'!C'+blabla 
    25   pltt, ts_Q, 't', -4, 4, date1, date2, /REMPLI $ 
     25  pltt, ts_Q, 't', MIN = -4, MAX = 4, date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    2626        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex 
    2727   
     
    3434 
    3535    title = prefix+' (BLACK)  - '+prefix2+' (RED) '+d1_d2_2+'!C'+blabla 
    36     pltt, ts_Q.arr , 't', -4., 4., date1, date2, /REMPLI $ 
     36    pltt, ts_Q.arr , 't', MIN = -4., MAX = 4., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    3737          , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 
    38     pltt, ts_Q2.arr, 't', -4., 4., date1, date2, /REMPLI $ 
     38    pltt, ts_Q2.arr, 't', date1, date2, /REMPLI $ 
    3939          , /ov1d, color = 250, small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 
    4040 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_S.pro

    r4495 r4497  
    3535 
    3636  title = prefix+' '+d1_d2+'!C'+blabla 
    37   pltt, ts_Sal, 't', 34.54, 34.76, date1, date2, /REMPLI $ 
     37  pltt, ts_Sal, 't', MIN = 34.54, MAX = 34.76, date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    3838        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex 
    3939 
    4040  title = prefix+' '+d1_d2+' - Levitus!C'+blabla 
    41   pltt, ts_Sal - Lev_xyz, 't', -.05, .05, date1, date2, /REMPLI $ 
     41  pltt, ts_Sal - Lev_xyz, 't', MIN = -.05, MAX = .05, date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    4242        , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /NOERASE, _extra = ex 
    4343 
     
    4747 
    4848  title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-1500m) '+blabla 
    49   pltt, ts_z, 'zt', -.5, .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so' $  
     49  pltt, ts_z, 'zt', MIN = -.5, MAX = .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $  
    5050        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex 
    5151   
    5252  title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-6000m) '+blabla 
    53   pltt, ts_z, 'zt', -.5, .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so' $  
     53  pltt, ts_z, 'zt', MIN = -.5, MAX = .5, INTER = .05, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $  
    5454        , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 
    5555   
     
    6969 
    7070    title = prefix+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 
    71     pltt, ts_Sal - ts_Sal2, 't', -0.1, 0.1, date1, date2, /REMPLI $ 
     71    pltt, ts_Sal - ts_Sal2, 't', MIN = -0.1, MAX = 0.1, date1, date2, /REMPLI , XGRIDSTYLE = 1 $ 
    7272          , COLOR = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, window = 2, _extra = ex 
    7373 
    7474    title = prefix+' - '+prefix2+' '+d1_d2_2+'!C ZT-plot (0-6000m) '+blabla 
    75     pltt, ts_z - ts_z2, 'zt', -0.1, 0.1, INTER = .01, date1, date2, /REMPLI, style = 'so0so' $  
     75    pltt, ts_z - ts_z2, 'zt', MIN = -0.1, MAX = 0.1, INTER = .01, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $  
    7676          , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 
    7777 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_SSH.pro

    r4495 r4497  
    2222; 
    2323  ts_SSH = rseries_ncdf(vssh, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 
    24   ;SF ts_SSH.arr = ts_SSH.arr * 100   &   ts_SSH.unit = 'cm' 
    2524  ts_SSH.arr = ts_SSH.arr   &   ts_SSH.unit = 'm' 
    2625  title = prefix+' '+d1_d2+'!C'+blabla 
    27   pltt, ts_SSH, 't', -12., 12., date1, date2, /REMPLI $ 
    28   ;SF       , small = [1, 2, 1], YTITLE = 'cm', TITLE = title, /PORTRAIT, _extra = ex 
     26  pltt, ts_SSH, 't', MIN = -12., MAX = 12., date1, date2, /REMPLI , XGRIDSTYLE = 1 $ 
    2927        , small = [1, 2, 1], YTITLE = 'm', TITLE = title, /PORTRAIT, _extra = ex 
    3028   
     
    3432    tsave = time 
    3533    ts_SSH2 = rseries_ncdf(vssh2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', mask2d = masknp) 
    36     ;SF ts_SSH2.arr = ts_SSH2.arr * 100   &   ts_SSH2.unit = 'cm' 
    3734    ts_SSH2.arr = ts_SSH2.arr   &   ts_SSH2.unit = 'm' 
    3835    time = tsave   &   IF n_elements(time) NE jpt THEN stop 
    3936 
    4037    title = prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+blabla 
    41     pltt, ts_SSH.arr - ts_SSH2.arr, 't', -10., 10., date1, date2, /REMPLI $ 
    42           ;SF , color = 250, small = [1, 2, 2], YTITLE = 'cm', TITLE = title, /noerase, _extra = ex 
     38    pltt, ts_SSH.arr - ts_SSH2.arr, 't', MIN = -10., MAX = 10., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    4339          , color = 250, small = [1, 2, 2], YTITLE = 'm', TITLE = title, /noerase, _extra = ex 
    4440 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_T.pro

    r4495 r4497  
    3535 
    3636  title = prefix+' '+d1_d2+'!C'+blabla 
    37   pltt, ts_Temp, 't', 1., 4.5, date1, date2, /REMPLI $ 
     37  pltt, ts_Temp, 't', MIN = 1., MAX = 4.5, date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    3838        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex 
    3939 
    4040  title = prefix+' '+d1_d2+' - Levitus!C'+blabla 
    41   pltt, ts_Temp - Lev_xyz, 't', -1., 1., date1, date2, /REMPLI $ 
     41  pltt, ts_Temp - Lev_xyz, 't', MIN = -1., MAX = 1., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    4242        , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /NOERASE, _extra = ex 
    4343 
     
    4747 
    4848  title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-1500m) '+blabla 
    49   pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so' $  
     49  pltt, ts_z, 'zt', MIN = -2.,MAX = 2., inter = .1, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $  
    5050        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex 
    5151   
    5252  title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-6000m) '+blabla 
    53   pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so' $  
     53  pltt, ts_z, 'zt', MIN = -2., MAX = 2., inter = .1, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $  
    5454        , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 
    5555   
     
    6969 
    7070    title = prefix+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 
    71     pltt, ts_Temp - ts_Temp2, 't', -1., 1., date1, date2, /REMPLI $ 
     71    pltt, ts_Temp - ts_Temp2, 't', MIN = -1., MAX = 1., date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
    7272          , COLOR = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, window = 2, _extra = ex 
    7373 
    7474    title = prefix+' - '+prefix2+' '+d1_d2_2+'!C ZT-plot (0-6000m) '+blabla 
    75     pltt, ts_z - ts_z2, 'zt', -1., 1., inter = .1, date1, date2, /REMPLI, style = 'so0so' $  
     75    pltt, ts_z - ts_z2, 'zt', MIN = -1., MAX = 1., inter = .1, date1, date2, /REMPLI, style = 'so0so', XGRIDSTYLE = 1 $  
    7676          , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 
    7777 
Note: See TracChangeset for help on using the changeset viewer.