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 4395 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM – NEMO

Ignore:
Timestamp:
2014-02-04T10:01:26+01:00 (10 years ago)
Author:
flavoni
Message:

update IDL scripts, see ticket: #724

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

Legend:

Unmodified
Added
Removed

  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r4394 r4395  
    338338;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
    339339;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
    340  
     340; 
    341341  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
    342342  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 
    … …  
    347347;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
    348348;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 
    349  
     349;  
    350350  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
    351351  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 
    … …  
    385385;SF  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: ' 
    386386;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex 
    387  
     387; 
    388388  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)' 
    389389  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex 
    390   
     390;  
    391391  cnt = cnt+1   &   blabla = 'Meridional stream Function: Atlantic' 
    392392  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex 
    393   
     393 ; 
    394394  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indian' 
    395395  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex 
    396   
     396 ; 
    397397  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indo-Pacific' 
    398398  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex 
    399  
     399; 
    400400  cnt = cnt+1   &   blabla = 'Equatorial Temperature' 
    401401  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex 
    402  
     402; 
    403403  cnt = cnt+1   &   blabla = 'Equatorial Salinity' 
    404404  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex 
    405  
     405; 
    406406  cnt = cnt+1   &   blabla = 'Equatorial zonal velocity' 
    407407  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex 
    408  
     408; 
    409409  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=700' 
    410410  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex 
    411  
     411; 
    412412  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=1000' 
    413413  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex 
    414  
     414; 
    415415  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40°N' 
    416416  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex 
    417  
     417; 
    418418  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38°N' 
    419419  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 
    420  
     420; 
    421421  cnt = cnt+1   &   blabla = 'Vertical Global mean T & S' 
    422422  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex 
    423  
     423; 
    424424  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*] 
    425425   

  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro

    r4394 r4395  
    8888  ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 
    8989  ; 
    90   title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+' Global Annual Mean Ice Area (Black SOLID simulation)'+'!C'+ 'and Extend minus 15% (Blue SOLID simulation)'+'!C'+'Observation (dashed)' 
     90  title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+' Global Annual Mean Ice Area (BLACK)'+'!C'+ 'and Extend minus 15% (BLUE)'+'!C'+'Observation (DASHED)' 
    9191  jpt=12 
    9292  time=julday(1,15,1900)+30*lindgen(12) 
    93   pltt, ICE_N, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 4., MAX = 16. $ 
    94         , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    95   pltt, ICE_N_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; dashed lines is LINESTYLE=2  $ 
    96         , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    97   pltt, vICE_area_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    98          , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    99   pltt, vICE_ext_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     93  pltt, ICE_N, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, MIN = 4., MAX = 16. $ 
     94       , COLOR = 000 , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     95  pltt, ICE_N_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $       ;;; dashed lines is LINESTYLE=2  $ 
     96        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
     97  pltt, vICE_area_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $  
    10098        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    101 ; 
    102   title ='Southern Hemisphere !C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Area (Black SOLID simulation)'+'!C'+ 'and Extend minus 15% (Blue SOLID simulation)'+'!C'+'Observation (dashed)' 
     99  pltt, vICE_ext_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $   ;;; dashed lines is LINESTYLE=2  $ 
     100        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     101; 
     102  title ='Southern Hemisphere !C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Area (BLACK)'+'!C'+ 'and Extend minus 15% (BLUE)'+'!C'+'Observation (DASHED)' 
    103103  pltt, ICE_S, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 20. $ 
    104         , small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    105   pltt, ICE_S_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $  
    106         , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    107   pltt,  vICE_area_SH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    108          , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     104        ,COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     105  pltt, ICE_S_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     106        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
     107  pltt,  vICE_area_SH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $  
     108         , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    109109  pltt,  vICE_ext_SH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    110110        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    … …  
    159159 
    160160 
    161   title = 'Northern Hemisphere'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2+'!C'+' Global Annual Mean Ice Area (BLUE)'+'!C'+ 'and Extend minus 15% (RED)' 
     161  title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
    162162  jpt=12 
    163163  time=julday(1,15,1900)+30*lindgen(12) 
    164   pltt, ICE_N.arr - ICE_N2.arr, 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT, window = 2 $ 
    165         , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    166   pltt, ICE_N_15.arr -  ICE_N2_15.arr , 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    167         , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    168 ;SF  pltt, vICE_area_NH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    169 ;SF        , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    170 ;SF  pltt, vICE_ext_NH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    171 ;SF        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    172 ; 
    173 ;  title ='Southern Hemisphere'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2+' - '+'!C'+'Global Annual Mean Ice Area (BLUE)'+'!C'+ 'and Extend minus 15% (RED)' 
    174   title ='Southern Hemisphere'+'!C' 
    175   pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., 19000101, 19001231, /REMPLI, /NOERASE $ 
    176          , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    177   pltt, ICE_S_15.arr -  ICE_S2_15.arr , 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    178         , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    179 ;SF  pltt,  vICE_area_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    180 ;SF         , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    181 ;SF  pltt,  vICE_ext_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    182 ;SF        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     164  pltt, ICE_N.arr, 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 4., MAX = 16., window = 2 $ 
     165        , COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; BLACK 
     166  pltt, ICE_N2.arr, 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT, window = 2 $ 
     167        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; RED 
     168  pltt, ICE_N_15.arr, 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     169        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
     170  pltt, ICE_N2_15.arr , 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     171        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     172  pltt, vICE_area_NH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT $  
     173        , /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., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     175        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; blu scuro 
     176; 
     177  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)' 
     178;  title ='Southern Hemisphere'+'!C' 
     179  pltt, ICE_S.arr, 't', -2., 2., 19000101, 19001231, /REMPLI, /NOERASE, MIN = 0., MAX = 20. $ 
     180         , COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     181  pltt, ICE_S2.arr, 't', -2., 2., 19000101, 19001231, /REMPLI, /NOERASE $ 
     182        , /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., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     184        , /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., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     186        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     187  pltt,  vICE_area_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT $  
     188          , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     189  pltt,  vICE_ext_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     190        , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    183191; 
    184192 

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

    r4394 r4395  
    128128    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
    129129 
    130     title = 'Fram Strait Areal Export'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2 
     130    title = 'Fram Strait Areal Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
    131131    jpt=12 
    132132    time=julday(1,15,1900)+30*lindgen(12) 
    133     pltt, ICE_area_export - ICE_area_export_2, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 0., MAX = 2. $ 
     133    pltt, ICE_area_export.arr, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, MIN = 0., MAX = 2., window = 2 $ 
    134134        , 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 
    135137    ; 
    136     title = 'Fram Strait Volume Export'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2 
    137     pltt, ICE_vol_export - ICE_vol_export_2, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 6. $ 
     138    ; 
     139    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. $ 
    138141        , 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 
    139144    ; 
    140145 

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

    r4320 r4395  
    6262  ICE_vol_S = {arr:ICE_vol_S * 1.e-9, unit : '10^9 Km^3'} 
    6363  ; 
    64   title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 
     64  title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Volume (BLACK) ' 
    6565  jpt=12 
    6666  time=julday(1,15,1900)+30*lindgen(12) 
    … …  
    6868        , small = [1, 2, 1], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 
    6969; 
    70   title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 
     70  title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (BLACK)' 
    7171  pltt, ICE_vol_S, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 11000. $ 
    7272        , small = [1, 2, 2], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 
    … …  
    109109    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
    110110 
    111     title = 'Northern Hemisphere'+'!C'+'Global Annual Mean Ice Volume'+prefix+' - '+prefix2+' '+d1_d2_2 
     111    title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 +'!C'+'Global Annual Mean Ice Volume (BLACK) ' 
    112112    jpt=12 
    113113    time=julday(1,15,1900)+30*lindgen(12) 
    114     pltt, ICE_vol_N.arr - ICE_vol_N2.arr, 't', -125., -60., 19000101, 19001231, /REMPLI, /PORTRAIT, window = 2 $ 
     114    pltt, ICE_vol_N.arr, 't', -125., -60., 19000101, 19001231, /REMPLI, /PORTRAIT, MIN = 0., MAX = 30000., window = 2 $ 
    115115          , small = [1, 2, 1], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 
     116    pltt, ICE_vol_N2.arr , 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
     117        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 
    116118    ; 
    117     title = 'Southern Hemisphere'+'!C'+'Global Annual Mean Ice Volume'+prefix+' - '+prefix2+' '+d1_d2_2 
    118     pltt, ICE_vol_S.arr - ICE_vol_S2.arr, 't', -620., -10., 19000101, 19001231, /REMPLI, /NOERASE $ 
     119    title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Global Annual Mean Ice Volume' 
     120    pltt, ICE_vol_S.arr, 't', -620., -10., 19000101, 19001231, /REMPLI, /NOERASE, MIN = 0., MAX = 11000. $ 
    119121          , small = [1, 2, 2], YTITLE = 'Km^3 ', TITLE = title, _extra = ex    
     122    pltt, ICE_vol_S2.arr, 't', -620., -10., 19000101, 19001231, /REMPLI, /NOERASE $ 
     123          , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = 'Km^3 ', TITLE = title, _extra = ex    
    120124 
    121125    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_vardef.sh

    r4394 r4395  
    1818idl_command=/Applications/itt/idl64/bin/idl 
    1919# 
    20 PS_DIR=$( pwd )/ts_500yfwb2_ps 
    21 PDF_DIR=$( pwd )/ts_500yfwb2_pdf 
     20PS_DIR=$( pwd )/ts_sbcmodMV_l334plus_ps 
     21PDF_DIR=$( pwd )/ts_sbcmodMV_l334plus_pdf 
    2222HTML_DIR=$( pwd )/html 
    2323SAXO_DIR=/Users/sflod/SAXO_DIR 
    … …  
    5252READ_ONLY_FIRST_RECORD=0   # if 0 then read all records in files else read only the first reacord in each file 
    5353# 
    54 DATE1=20010101   ;   DATE2=25001231 
     54#DATE1=20010101   ;   DATE2=25001231 
     55DATE1=20410101   ;   DATE2=20501231 
    5556#DATE1=00010101      ;   DATE2=00101231 
    5657# 
    57 VAR1_T=votemper         ;   V1T_PREF=500yfwb2      ;   V1T_SUFF=_1Y_grid_T.nc 
    58 VAR1_S=vosaline         ;   V1S_PREF=500yfwb2      ;   V1S_SUFF=_1Y_grid_T.nc 
    59 VAR1_SSH=sossheig       ;   V1SSH_PREF=500yfwb2    ;   V1SSH_SUFF=_1Y_grid_T.nc 
    60 VAR1_Q=sohflxdo         ;   V1Q_PREF=500yfwb2      ;   V1Q_SUFF=_1Y_grid_T.nc 
    61 VAR1_EMP=sowaflup       ;   V1EMP_PREF=500yfwb2    ;   V1EMP_SUFF=_1Y_grid_T.nc 
    62 VAR1_U=vozoeftr         ;   V1U_PREF=500yfwb2      ;   V1U_SUFF=_1Y_grid_U.nc 
    63 VAR1_V=vomeeftr         ;   V1V_PREF=500yfwb2      ;   V1V_SUFF=_1Y_grid_V.nc 
    64 VAR1_ICE=iiceconc       ;   V1ICE_PREF=500yfwb2    ;   V1ICE_SUFF=_1M_icemod.nc 
    65 VAR1_Ithick=iicethic    ;   V1It_PREF=500yfwb2     ;   V1It_SUFF=_1M_icemod.nc 
    66 VAR1_SNOW=isnowthi      ;   V1SNOW_PREF=500yfwb2   ;   V1SNOW_SUFF=_1M_icemod.nc 
    67 VAR1_IvelV=iicevelv     ;   V1IvV_PREF=500yfwb2    ;   V1IvV_SUFF=_1M_icemod.nc 
     58VAR1_T=thetao        ;   V1T_PREF=sbcmodMV      ;   V1T_SUFF=_1Y_grid_T.nc 
     59VAR1_S=so            ;   V1S_PREF=sbcmodMV      ;   V1S_SUFF=_1Y_grid_T.nc 
     60VAR1_SSH=zos         ;   V1SSH_PREF=sbcmodMV    ;   V1SSH_SUFF=_1Y_grid_T.nc 
     61VAR1_Q=tohfls        ;   V1Q_PREF=sbcmodMV      ;   V1Q_SUFF=_1Y_grid_T.nc 
     62VAR1_EMP=wfo         ;   V1EMP_PREF=sbcmodMV    ;   V1EMP_SUFF=_1Y_grid_T.nc 
     63VAR1_U=uocetr_eff    ;   V1U_PREF=sbcmodMV      ;   V1U_SUFF=_1Y_grid_U.nc 
     64VAR1_V=vocetr_eff    ;   V1V_PREF=sbcmodMV      ;   V1V_SUFF=_1Y_grid_V.nc 
     65VAR1_ICE=siconc      ;   V1ICE_PREF=sbcmodMV    ;   V1ICE_SUFF=_1M_icemod.nc 
     66VAR1_Ithick=sithic   ;   V1It_PREF=sbcmodMV     ;   V1It_SUFF=_1M_icemod.nc 
     67VAR1_SNOW=sndept     ;   V1SNOW_PREF=sbcmodMV   ;   V1SNOW_SUFF=_1M_icemod.nc 
     68VAR1_IvelV=sivelv    ;   V1IvV_PREF=sbcmodMV    ;   V1IvV_SUFF=_1M_icemod.nc 
    6869# 
    6970#===================== EXP2 ===================== 
    7071# 
    7172#DATE1_2=20010101   ;   DATE2_2=20101231 
    72 DATE1_2=20010101   ;   DATE2_2=25001231 
     73DATE1_2=20410101   ;   DATE2_2=20501231 
     74#DATE1_2=20010101   ;   DATE2_2=25001231 
    7375# 
    74 VAR2_T=votemper         ;   V2T_PREF=500yfwb2      ;   V2T_SUFF=_1Y_grid_T.nc 
    75 VAR2_S=vosaline         ;   V2S_PREF=500yfwb2      ;   V2S_SUFF=_1Y_grid_T.nc 
    76 VAR2_SSH=sossheig       ;   V2SSH_PREF=500yfwb2    ;   V2SSH_SUFF=_1Y_grid_T.nc 
    77 VAR2_Q=sohflxdo         ;   V2Q_PREF=500yfwb2      ;   V2Q_SUFF=_1Y_grid_T.nc 
    78 VAR2_EMP=sowaflup       ;   V2EMP_PREF=500yfwb2    ;   V2EMP_SUFF=_1Y_grid_T.nc 
    79 VAR2_U=vozoeftr         ;   V2U_PREF=500yfwb2      ;   V2U_SUFF=_1Y_grid_U.nc 
    80 VAR2_V=vomeeftr         ;   V2V_PREF=500yfwb2      ;   V2V_SUFF=_1Y_grid_V.nc 
    81 VAR2_ICE=iiceconc       ;   V2ICE_PREF=500yfwb2    ;   V2ICE_SUFF=_1M_icemod.nc 
    82 VAR2_Ithick=iicethic    ;   V2It_PREF=500yfwb2     ;   V2It_SUFF=_1M_icemod.nc 
    83 VAR2_SNOW=isnowthi      ;   V2SNOW_PREF=500yfwb2   ;   V2SNOW_SUFF=_1M_icemod.nc 
    84 VAR2_IvelV=iicevelv     ;   V2IvV_PREF=500yfwb2    ;   V2IvV_SUFF=_1M_icemod.nc 
     76VAR2_T=thetao      ;   V2T_PREF=l334plus      ;   V2T_SUFF=_1Y_grid_T.nc 
     77VAR2_S=so          ;   V2S_PREF=l334plus      ;   V2S_SUFF=_1Y_grid_T.nc 
     78VAR2_SSH=zos       ;   V2SSH_PREF=l334plus    ;   V2SSH_SUFF=_1Y_grid_T.nc 
     79VAR2_Q=tohfls      ;   V2Q_PREF=l334plus      ;   V2Q_SUFF=_1Y_grid_T.nc 
     80VAR2_EMP=wfo       ;   V2EMP_PREF=l334plus    ;   V2EMP_SUFF=_1Y_grid_T.nc 
     81VAR2_U=uocetr_eff  ;   V2U_PREF=l334plus      ;   V2U_SUFF=_1Y_grid_U.nc 
     82VAR2_V=vocetr_eff  ;   V2V_PREF=l334plus      ;   V2V_SUFF=_1Y_grid_V.nc 
     83VAR2_ICE=siconc    ;   V2ICE_PREF=l334plus    ;   V2ICE_SUFF=_1M_icemod.nc 
     84VAR2_Ithick=sithic ;   V2It_PREF=l334plus     ;   V2It_SUFF=_1M_icemod.nc 
     85VAR2_SNOW=sndept   ;   V2SNOW_PREF=l334plus   ;   V2SNOW_SUFF=_1M_icemod.nc 
     86VAR2_IvelV=sivelv  ;   V2IvV_PREF=sbcmodMV    ;   V2IvV_SUFF=_1M_icemod.nc 
    8587# 
    8688######################### Export Variables ############################### 
    … …  
    115117export VAR1_SNOW    V1SNOW_PREF  V1SNOW_SUFF 
    116118export VAR1_IvelV   V1IvV_PREF   V1IvV_PREF 
     119export VAR2_IvelV   V2IvV_PREF  V2IvV_PREF 
     120 
    117121#===================== EXP2 ===================== 
    118122export DATE1_2      DATE2_2 
    119 export VAR2_T       V2T_PREF     V2T_SUFF 
    120 export VAR2_S       V2S_PREF     V2S_SUFF 
    121 export VAR2_SSH     V2SSH_PREF   V2SSH_SUFF 
    122 export VAR2_Q       V2Q_PREF     V2Q_SUFF 
    123 export VAR2_EMP     V2EMP_PREF   V2EMP_SUFF 
    124 export VAR2_U       V2U_PREF     V2U_SUFF 
    125 export VAR2_V       V2V_PREF     V2V_SUFF 
    126 export VAR2_ICE     V2ICE_PREF   V2ICE_SUFF 
    127 export VAR2_Ithick  V2It_PREF    V2It_SUFF 
    128 export VAR2_SNOW    V2SNOW_PREF  V2SNOW_SUFF 
    129 export VAR2_IvelV   V2IvV_PREF   V2IvV_PREF 
     123export VAR2_T       V2T_PREF    V2T_SUFF 
     124export VAR2_S       V2S_PREF    V2S_SUFF 
     125export VAR2_SSH     V2SSH_PREF  V2SSH_SUFF 
     126export VAR2_Q       V2Q_PREF    V2Q_SUFF 
     127export VAR2_EMP     V2EMP_PREF  V2EMP_SUFF 
     128export VAR2_U       V2U_PREF    V2U_SUFF 
     129export VAR2_V       V2V_PREF    V2V_SUFF 
     130export VAR2_ICE     V2ICE_PREF  V2ICE_SUFF 
     131export VAR2_Ithick  V2It_PREF   V2It_SUFF 
     132export VAR2_SNOW    V2SNOW_PREF V2SNOW_SUFF 
     133export VAR2_IvelV   V2IvV_PREF  V2IvV_PREF 
     134 
    130135# 
Note: See TracChangeset for help on using the changeset viewer.