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.
Changeset 4495 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro – NEMO

Ignore:
Timestamp:
2014-02-10T16:51:42+01:00 (10 years ago)
Author:
flavoni
Message:

commit FRAM strait export for idl plots, see ticket #724

File:
1 edited

Legend:

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

    r4395 r4495  
    7171  ICE_area_export = reform(ICE_area_export, 12, nyr) 
    7272  ICE_area_export = total(ICE_area_export,2)/nyr 
    73   ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'} 
     73  ; ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'}  ; annual mean 
     74  ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'}    ; monthly mean 
     75   
    7476  ; 
    7577  ICE_vol_export = reform(ICE_vol_export, 12, nyr) 
    7678  ICE_vol_export = total(ICE_vol_export,2)/nyr 
    77   ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'} 
     79  ; ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'}      ; annual mean 
     80  ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 30 , unit : '10^3 Km^3/month'}      ; monthly mean 
    7881 
    7982  ; 
     
    8285  jpt=12 
    8386  time=julday(1,15,1900)+30*lindgen(12) 
    84   pltt, ICE_area_export, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 0., MAX = 2. $ 
    85         , small = [1, 2, 1],YTITLE = '10^6 Km^2/year',  TITLE = title, _extra = ex 
     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     
    8691  ; 
    8792  ;title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 
    8893  title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 
    89   pltt, ICE_vol_export, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 6. $ 
    90         , small = [1, 2, 2], YTITLE = '10^3 Km^3/year', TITLE = title, _extra = ex 
     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 
    9198  ; 
    9299  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
     
    119126    ICE_area_export_2 = reform(ICE_area_export_2, 12, nyr) 
    120127    ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 
    121     ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'} 
     128    ; ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'}   ; annual mean 
     129    ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'}     ; monthly mean 
    122130    ; 
    123131    ICE_vol_export_2 = reform(ICE_vol_export_2, 12, nyr) 
    124132    ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 
    125     ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'} 
     133    ; ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'}     ; annual mean 
     134    ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 30 , unit : '10^3 Km^3/month'}       ; monthly mean 
    126135    ; 
    127136    ; 
     
    131140    jpt=12 
    132141    time=julday(1,15,1900)+30*lindgen(12) 
    133     pltt, ICE_area_export.arr, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, MIN = 0., MAX = 2., window = 2 $ 
    134         , small = [1, 2, 1],YTITLE = '10^6 Km^2/year',  TITLE = title, _extra = ex 
    135     pltt, ICE_area_export_2.arr ,'t', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ 
    136         , /ov1d, COLOR = 250, small = [1, 2, 1],YTITLE = '10^6 Km^2/year',  TITLE = title, _extra = ex 
    137     ; 
     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 
    138150    ; 
    139151    title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
    140     pltt, ICE_vol_export.arr, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 6. $ 
    141         , small = [1, 2, 2], YTITLE = '10^3 Km^3/year', TITLE = title, _extra = ex 
    142     pltt, ICE_vol_export_2.arr, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE  $ 
    143         , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^3 Km^3/year', TITLE = title, _extra = ex 
     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 
    144160    ; 
    145161 
Note: See TracChangeset for help on using the changeset viewer.