Changeset 4755


Ignore:
Timestamp:
2014-09-02T15:55:01+02:00 (6 years ago)
Author:
flavoni
Message:

add new idl plots; mld new map for example, see ticket #724

Location:
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
Files:
2 added
10 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/README

    r4394 r4755  
    2424 
    2525 
    26  
    27  
     26#---------------------------------------------------------------------- 
     27# HOW TO USE  
     28#---------------------------------------------------------------------- 
     29./std_main.sh -plot -pdf 
     30or 
     31./std_main.sh -ts -pdf 
    2832 
    2933 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro

    r4559 r4755  
    2121   
    2222  IF keyword_set(arc) THEN BEGIN  
    23     domdef, 20, 380, 50, 90 
     23    domdef, 20, 380, 40, 90 
    2424    map = [90, 0, 0] 
    2525  ENDIF 
    2626  IF keyword_set(ant) THEN BEGIN  
    27     domdef, 20, 380, -90, -50 
     27    domdef, 20, 380, -90, -40 
    2828    map = [-90, 0, 0] 
    2929  ENDIF 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r4559 r4755  
    172172  Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 
    173173  jan = where(mm EQ 1, cnt) 
     174  Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    174175  sept = where(mm EQ 9, cnt) 
    175   Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    176176  Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    177177  undefine, Ithi_1 
     
    288288     Isal_sept_2 = {arr:-1} 
    289289   ENDELSE 
    290 ; 
     290 
    291291  IvelU_1 = read_ncdf(getenv('VAR1_IvelU'), allrecords = allrec, filename = std_file1_I )  
    292292  IvelV_1 = read_ncdf(getenv('VAR1_IvelV'), allrecords = allrec, filename = std_file1_I )  
     
    335335    Ivelo_sept_2 = {arr:-1} 
    336336  ENDELSE 
    337 ; 
     337 
    338338  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
    339339  caldat, time, mm 
     
    396396  cnt = cnt+1   &   blabla = 'Erp salinity damping term' 
    397397  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 
    398 ; 
     398;; 
    399399  cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term' 
    400400  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex 
     
    439439  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD2, MLD, _extra = ex 
    440440; 
     441  cnt = cnt+1   &   blabla = 'Mixed layer depth ortho plan ARCTIC' 
     442  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_MLD_ortho, MLD1, MLD2, /ARC, _extra = ex 
     443;sf 
     444  cnt = cnt+1   &   blabla = 'Mixed layer depth ortho plan ANTARTIC' 
     445  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_MLD_ortho, MLD1, MLD2, /ANT, _extra = ex 
     446;sf 
    441447  cnt = cnt+1   &   blabla = 'Zonal mean Mixed layer depth' 
    442448  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ZonMld, MLD1, MLD2, MLD, _extra = ex 
     
    564570  cnt = cnt+1   &   blabla = 'Antartic Ice Velocity: SEPT' 
    565571  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_sept_1, IvelU_sept_2, IvelV_sept_1, IvelV_sept_2, Ivelo_sept_1, Ivelo_sept_2, /ANT,  /SEPT, _extra = ex 
    566 ; 
     572 
    567573  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
    568574  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh

    r4562 r4755  
    2121#idl_command=idl71 
    2222# 
    23 PS_DIR=$( pwd )/OR2LPv36_ps_plot 
    24 PDF_DIR=$( pwd )/OR2LPv36_pdf_plot 
     23PS_DIR=$( pwd )/p4H25a50-testht_ps_plot 
     24PDF_DIR=$( pwd )/p4H25a50-testht_pdf_plot 
    2525HTML_DIR=$( pwd )/html_plot 
    2626SAXO_DIR=/Users/sflod/SAXO_DIR 
    2727# 
    28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/500y_pisces     # path of data in NetCDF format 
     28DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
    2929DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    3030DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
     
    5454# 
    5555#===================== EXP1 ===================== 
    56 # 23610101_23701231 
    57 FILE1_T=OR2LPv36_23610101_23701231_1Y_grid_T.nc   # exp1 grid T input file 
    58 FILE1_U=OR2LPv36_23610101_23701231_1Y_grid_U.nc   # exp1 grid U input file 
    59 FILE1_V=OR2LPv36_23610101_23701231_1Y_grid_V.nc   # exp1 grid V input file 
    60 FILE1_I=OR2LPv36_23610101_23701231_1M_icemod.nc   # exp1 ice    input file 
     56# 
     57FILE1_T=p4H25a50_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file 
     58FILE1_U=p4H25a50_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file 
     59FILE1_V=p4H25a50_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file 
     60FILE1_I=p4H25a50_20410101_20501231_1M_icemod.nc   # exp1 ice    input file 
    6161VAR1_T=thetao 
    6262VAR1_S=so         
    63 VAR1_QNET=tohfls 
    64 VAR1_ERP=sowafldp       
     63VAR1_QNET=qt 
     64VAR1_ERP=wfcorr 
    6565VAR1_EMP=wfo 
    6666VAR1_MLD=mldr10_1 
    6767VAR1_U=uocetr_eff 
    6868VAR1_V=vocetr_eff 
    69 VAR1_Ithick=iicethic 
    70 VAR1_Ifrac=soicecov 
    71 VAR1_Isnow=isnowthi 
    72 #VAR1_Isal=sisali 
    73 #VAR1_Iage=siages 
    74 VAR1_IvelU=iicevelu 
    75 VAR1_IvelV=iicevelv 
    76 #VAR1_Ivelo=sivelo 
     69VAR1_Ithick=sithic 
     70VAR1_Ifrac=siconc 
     71VAR1_Isnow=snvolu 
     72VAR1_Isal=sisali 
     73VAR1_Iage=siages 
     74VAR1_IvelU=sivelu 
     75VAR1_IvelV=sivelv 
     76VAR1_Ivelo=sivelo 
    7777# 
    7878#===================== EXP2 ===================== 
    7979# 
    80 FILE2_T=OR2LPv36_23610101_23701231_1Y_grid_T.nc   # exp1 grid T input file 
    81 FILE2_U=OR2LPv36_23610101_23701231_1Y_grid_U.nc   # exp1 grid U input file 
    82 FILE2_V=OR2LPv36_23610101_23701231_1Y_grid_V.nc   # exp1 grid V input file 
    83 FILE2_I=OR2LPv36_23610101_23701231_1M_icemod.nc   # exp1 ice    input file 
     80FILE2_T=testht_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file 
     81FILE2_U=testht_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file 
     82FILE2_V=testht_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file 
     83FILE2_I=testht_20410101_20501231_1M_icemod.nc   # exp1 ice    input file 
    8484VAR2_T=thetao 
    8585VAR2_S=so 
    86 VAR2_QNET=tohfls 
    87 VAR2_ERP=sowafldp 
     86VAR2_QNET=qt 
     87VAR2_ERP=wfcorr 
    8888VAR2_EMP=wfo 
    8989VAR2_MLD=mldr10_1 
    9090VAR2_U=uocetr_eff 
    9191VAR2_V=vocetr_eff 
    92 VAR2_Ithick=iicethic 
    93 VAR2_Ifrac=soicecov 
    94 VAR2_Isnow=isnowthi 
    95 #VAR2_Isal=sisali 
    96 #VAR2_Iage=siages 
    97 VAR2_IvelU=iicevelu 
    98 VAR2_IvelV=iicevelv 
    99 #VAR2_Ivelo=sivelo 
     92VAR2_Ithick=sithic 
     93VAR2_Ifrac=siconc 
     94VAR2_Isnow=snvolu 
     95VAR2_Isal=sisali 
     96VAR2_Iage=siages 
     97VAR2_IvelU=sivelu 
     98VAR2_IvelV=sivelv 
     99VAR2_Ivelo=sivelo 
    100100# 
    101101######################### Export Variables ############################### 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names

    r4559 r4755  
    7474VAR1_IvelU=iicevelu 
    7575VAR1_IvelV=iicevelv 
    76 VAR1_Ivelo=iicevelo 
     76#VAR1_Ivelo=iicevelo 
    7777# 
    7878#===================== EXP2 ===================== 
     
    9999VAR2_IvelU=iicevelu 
    100100VAR2_IvelV=iicevelv 
    101 VAR2_Ivelo=iicevelo 
     101#VAR2_Ivelo=iicevelo 
    102102# 
    103103######################### Export Variables ############################### 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro

    r4497 r4755  
    8888  ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 
    8989  ; 
    90   title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     90  ;;title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     91  title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
    9192  jpt=12 
    9293  time=julday(1,15,1900)+30*lindgen(12) 
    93   pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
     94  pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, LINESTYLE=2, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    9495       , COLOR = 000 , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    95   pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $       ;;; dashed lines is LINESTYLE=2  $ 
     96  pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT $       ;;; dashed lines is LINESTYLE=2  $ 
    9697        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    97   pltt, vICE_area_NH, 't',  /REMPLI, /PORTRAIT  $  
     98  pltt, vICE_area_NH, 't',  /REMPLI, /PORTRAIT, LINESTYLE=2  $  
    9899        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    99   pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $   ;;; dashed lines is LINESTYLE=2  $ 
     100  pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT $   ;;; dashed lines is LINESTYLE=2  $ 
    100101        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    101102; 
    102   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', MIN = 0., MAX = 20., /REMPLI, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
     103  title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
     104  pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, LINESTYLE=2, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
    104105        ,COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    105   pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     106  pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT $  
    106107        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    107   pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT  $  
     108  pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2  $  
    108109         , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    109   pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     110  pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT $  
    110111        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    111112; 
     
    159160 
    160161 
    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)' 
     162  ;;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)' 
     163  title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
    162164  jpt=12 
    163165  time=julday(1,15,1900)+30*lindgen(12) 
    164   pltt, ICE_N, 't', MIN = 4, MAX = 16,  /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 
     166  pltt, ICE_N, 't', MIN = 4, MAX = 16,  /REMPLI, /PORTRAIT, LINESTYLE=2, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 
    165167        , COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; BLACK 
    166   pltt, ICE_N2, 't', /REMPLI, /PORTRAIT $ 
     168  pltt, ICE_N2, 't', /REMPLI, /PORTRAIT , LINESTYLE=2 $ 
    167169        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; RED 
    168   pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     170  pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    169171        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    170   pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     172  pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    171173        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    172   pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT  $  
     174  pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2  $  
    173175        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; light blue 
    174   pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     176  pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT $  
    175177        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; blu scuro 
    176178; 
    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)' 
     179  title ='Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+'Global Annual Mean Ice Area (DASHED)'+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
    178180;  title ='Southern Hemisphere'+'!C' 
    179   pltt, ICE_S, 't', MIN = 0., MAX = 20.,  /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
     181  pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, LINESTYLE=2, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    180182         , COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    181   pltt, ICE_S2, 't', /REMPLI, /NOERASE $ 
     183  pltt, ICE_S2, 't', /REMPLI, /NOERASE, LINESTYLE=2 $ 
    182184        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    183   pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     185  pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    184186        , /ov1d, COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    185   pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     187  pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    186188        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    187   pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT $  
     189  pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    188190          , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    189   pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     191  pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT $  
    190192        , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    191193; 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro

    r4514 r4755  
    55arr2d = arr2d[t1:t2] 
    66arr2d = reform(arr2d,12,nyear) ; put in 2D array 
    7 arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 
     7;arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 
     8arr2d = arr2d[*, nyear-1] ;  select last year 
    89 
    910return, arr2d 
     
    3435; 
    3536  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 
     37  d2 = '('+strtrim(date2, 1)+')' 
    3638; 
    3739  iodir = std_iodir_data 
     
    7476  ;; AREA 
    7577  ICE_area_export = reform(ICE_area_export, 12, nyr) 
    76   ICE_area_export = total(ICE_area_export,2)/nyr 
     78  ;ICE_area_export = total(ICE_area_export,2)/nyr  ; old version monthly mean average  over all years 
     79  ICE_area_export = ICE_area_export[*,nyr-1]  ; dim= 12, index 0-11; last year choosen 
    7780  ; ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'}  ; annual mean 
    7881  ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'}    ; monthly mean 
     
    8083  ; 
    8184  ICE_vol_export = reform(ICE_vol_export, 12, nyr) 
    82   ICE_vol_export = total(ICE_vol_export,2)/nyr 
     85  ;ICE_vol_export = total(ICE_vol_export,2)/nyr 
     86  ICE_vol_export = ICE_vol_export[*,nyr-1] 
    8387  ; ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 365 , unit : '10^3 Km^3/year'}      ; annual mean 
    8488  ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}      ; monthly mean 
    8589 
    8690  ; 
    87   title = 'Fram Strait Areal Export'+'!C'+prefix+' '+d1_d2 
     91  ;title = 'Fram Strait Areal Export: LAST YEAR'+'!C'+prefix+' '+d1_d2 
     92  title = 'Fram Strait Areal Export: LAST YEAR'+'!C'+prefix+' '+d2 
    8893  jpt=12 
    8994  time=julday(1,15,1900)+30*lindgen(12) 
     
    99104  xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    100105  ; 
    101   title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 
     106  ;title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' '+d1_d2 
     107  title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' '+d2 
    102108  pltt, ICE_vol_export, 't',  /REMPLI, MIN = 60., MAX = 500. , /NOERASE, XGRIDSTYLE = 1 $ 
    103109       , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
     
    118124 
    119125    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 
     126    d2_2 = '('+strtrim(date2_2, 1)+')' 
    120127    tsave = time 
    121128    domdef, indx1, indx2, indy1, indy2, /xindex, /yindex,/memeindices 
     
    140147    ;; AREA 
    141148    ICE_area_export_2 = reform(ICE_area_export_2, 12, nyr) 
    142     ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 
     149    ;ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 
     150    ICE_area_export_2 = ICE_area_export_2[*,nyr-1] 
    143151    ; ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'}   ; annual mean 
    144152    ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'}     ; monthly mean 
    145153    ; 
    146154    ICE_vol_export_2 = reform(ICE_vol_export_2, 12, nyr) 
    147     ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 
     155    ;ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 
     156    ICE_vol_export_2 = ICE_vol_export_2[*,nyr-1] 
    148157    ; ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'}     ; annual mean 
    149158    ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}       ; monthly mean 
     
    152161    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
    153162 
    154     title = 'Fram Strait Areal Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     163    ;title = 'Fram Strait Areal Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     164    title = 'Fram Strait Areal Export LAST YEAR'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d2_2 
    155165    jpt=12 
    156166    time=julday(1,15,1900)+30*lindgen(12) 
     
    170180    xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    171181    ; 
    172     title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     182    ;title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     183    title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d2_2 
    173184    pltt, ICE_vol_export, 't', /REMPLI,  MIN = 60., MAX = 500., /NOERASE, XGRIDSTYLE = 1 $ 
    174185        , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro

    r4504 r4755  
    7070; 
    7171  title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (BLACK)' 
    72   pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
     72  pltt, ICE_vol_S, 't', MIN = 0., MAX = 12000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
    7373        , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
    7474; 
     
    119119    ; 
    120120    title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Global Annual Mean Ice Volume' 
    121     pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
     121    pltt, ICE_vol_S, 't', MIN = 0., MAX = 12000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    122122          , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex    
    123123    pltt, ICE_vol_S2, 't',  /REMPLI, /NOERASE $ 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro

    r4559 r4755  
    6060  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SSH, masknp, POSTSCRIPT = postscript, _extra = ex 
    6161; 
    62   cnt = cnt+1   &   blabla = 'Global Mean Q net' 
    63   IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Q, masknp, POSTSCRIPT = postscript, _extra = ex 
     62;;;;  cnt = cnt+1   &   blabla = 'Global Mean Q net' 
     63;;;;  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Q, masknp, POSTSCRIPT = postscript, _extra = ex 
    6464; 
    6565  cnt = cnt+1   &   blabla = 'Global Mean EMP' 
     
    7777  cnt = cnt+1   &   blabla = 'Sea-Ice Volume' 
    7878  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 
    79 ; 
    80   cnt = cnt+1   &   blabla = 'ICE Velocity ' 
    81   IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 
    82 ; 
     79;SF; and probably useless 
     80;SF  cnt = cnt+1   &   blabla = 'ICE Velocity ' 
     81;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 
     82;SF; 
    8383  cnt = cnt+1   &   blabla = 'Sea-Ice Volume Export at Fram Strait' 
    8484  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_FRAM, masknp, POSTSCRIPT = postscript, _extra = ex 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh

    r4559 r4755  
    1818idl_command=/Applications/itt/idl64/bin/idl 
    1919# 
    20 PS_DIR=$( pwd )/ts_heatbug_ps 
    21 PDF_DIR=$( pwd )/ts_heatbug_pdf 
     20PS_DIR=$( pwd )/ts_p4H25a50-testht_ps 
     21PDF_DIR=$( pwd )/ts_p4H25a50-testht_pdf 
    2222HTML_DIR=$( pwd )/html 
    2323SAXO_DIR=/Users/sflod/SAXO_DIR 
    2424# 
    25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada      # path of data in NetCDF format 
     25DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
    2626DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    2727DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
     
    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 
    55 DATE1=20310101   ;   DATE2=20401231 
    56 #DATE1=00010101      ;   DATE2=00101231 
     54DATE1=20010101         ;   DATE2=20501231 
    5755# 
    58 VAR1_T=thetao        ;   V1T_PREF=heatbug       ;   V1T_SUFF=_1Y_grid_T.nc 
    59 VAR1_S=so            ;   V1S_PREF=heatbug      ;   V1S_SUFF=_1Y_grid_T.nc 
    60 VAR1_SSH=zos         ;   V1SSH_PREF=heatbug    ;   V1SSH_SUFF=_1Y_grid_T.nc 
    61 VAR1_Q=tohfls        ;   V1Q_PREF=heatbug      ;   V1Q_SUFF=_1Y_grid_T.nc 
    62 VAR1_EMP=wfo         ;   V1EMP_PREF=heatbug    ;   V1EMP_SUFF=_1Y_grid_T.nc 
    63 VAR1_U=uocetr_eff    ;   V1U_PREF=heatbug      ;   V1U_SUFF=_1Y_grid_U.nc 
    64 VAR1_V=vocetr_eff    ;   V1V_PREF=heatbug      ;   V1V_SUFF=_1Y_grid_V.nc 
    65 VAR1_ICE=siconc      ;   V1ICE_PREF=heatbug    ;   V1ICE_SUFF=_1M_icemod.nc 
    66 VAR1_Ithick=sithic   ;   V1It_PREF=heatbug     ;   V1It_SUFF=_1M_icemod.nc 
    67 VAR1_SNOW=sndept     ;   V1SNOW_PREF=heatbug   ;   V1SNOW_SUFF=_1M_icemod.nc 
    68 VAR1_IvelV=sivelv    ;   V1IvV_PREF=heatbug    ;   V1IvV_SUFF=_1M_icemod.nc 
    69 VAR1_Ivel=sivelo     ;   V1Iv_PREF=heatbug     ;   V1Iv_SUFF=_1M_icemod.nc 
     56VAR1_T=thetao          ;   V1T_PREF=p4H25a50      ;   V1T_SUFF=_1Y_grid_T.nc 
     57VAR1_S=so              ;   V1S_PREF=p4H25a50      ;   V1S_SUFF=_1Y_grid_T.nc 
     58VAR1_SSH=zos         ;   V1SSH_PREF=p4H25a50    ;   V1SSH_SUFF=_1Y_grid_T.nc 
     59##VAR1_Q=qt              ;   V1Q_PREF=p4H25a50      ;   V1Q_SUFF=_1Y_grid_T.nc 
     60VAR1_EMP=wfo         ;   V1EMP_PREF=p4H25a50    ;   V1EMP_SUFF=_1Y_grid_T.nc 
     61VAR1_U=uocetr_eff      ;   V1U_PREF=p4H25a50      ;   V1U_SUFF=_1Y_grid_U.nc 
     62VAR1_V=vocetr_eff      ;   V1V_PREF=p4H25a50      ;   V1V_SUFF=_1Y_grid_V.nc 
     63VAR1_ICE=siconc      ;   V1ICE_PREF=p4H25a50    ;   V1ICE_SUFF=_1M_icemod.nc 
     64VAR1_Ithick=sithic    ;   V1It_PREF=p4H25a50     ;   V1It_SUFF=_1M_icemod.nc 
     65VAR1_SNOW=snthic    ;   V1SNOW_PREF=p4H25a50   ;   V1SNOW_SUFF=_1M_icemod.nc 
     66VAR1_IvelV=sivelv    ;   V1IvV_PREF=p4H25a50    ;   V1IvV_SUFF=_1M_icemod.nc  
     67VAR1_Ivel=sivelo      ;   V1Iv_PREF=p4H25a50     ;   V1Iv_SUFF=_1M_icemod.nc 
    7068# 
    7169#===================== EXP2 ===================== 
    7270# 
    73 #DATE1_2=20010101   ;   DATE2_2=20101231 
    74 DATE1_2=20310101   ;   DATE2_2=20401231 
    75 #DATE1_2=20010101   ;   DATE2_2=25001231 
     71DATE1_2=20010101       ;   DATE2_2=20501231 
    7672# 
    77 VAR2_T=thetao      ;   V2T_PREF=heatbug      ;   V2T_SUFF=_1Y_grid_T.nc 
    78 VAR2_S=so          ;   V2S_PREF=heatbug      ;   V2S_SUFF=_1Y_grid_T.nc 
    79 VAR2_SSH=zos       ;   V2SSH_PREF=heatbug    ;   V2SSH_SUFF=_1Y_grid_T.nc 
    80 VAR2_Q=tohfls      ;   V2Q_PREF=heatbug      ;   V2Q_SUFF=_1Y_grid_T.nc 
    81 VAR2_EMP=wfo       ;   V2EMP_PREF=heatbug    ;   V2EMP_SUFF=_1Y_grid_T.nc 
    82 VAR2_U=uocetr_eff  ;   V2U_PREF=heatbug      ;   V2U_SUFF=_1Y_grid_U.nc 
    83 VAR2_V=vocetr_eff  ;   V2V_PREF=heatbug      ;   V2V_SUFF=_1Y_grid_V.nc 
    84 VAR2_ICE=siconc    ;   V2ICE_PREF=heatbug    ;   V2ICE_SUFF=_1M_icemod.nc 
    85 VAR2_Ithick=sithic ;   V2It_PREF=heatbug     ;   V2It_SUFF=_1M_icemod.nc 
    86 VAR2_SNOW=sndept   ;   V2SNOW_PREF=heatbug   ;   V2SNOW_SUFF=_1M_icemod.nc 
    87 VAR2_IvelV=sivelv  ;   V2IvV_PREF=heatbug    ;   V2IvV_SUFF=_1M_icemod.nc 
    88 VAR2_Ivel=sivelo   ;   V2Iv_PREF=heatbug     ;   V2Iv_SUFF=_1M_icemod.nc 
     73VAR2_T=thetao          ;   V2T_PREF=testht     ;   V2T_SUFF=_1Y_grid_T.nc 
     74VAR2_S=so           ;   V2S_PREF=testht     ;   V2S_SUFF=_1Y_grid_T.nc 
     75VAR2_SSH=zos         ;   V2SSH_PREF=testht   ;   V2SSH_SUFF=_1Y_grid_T.nc 
     76##VAR2_Q=qt            ;   V2Q_PREF=testht     ;   V2Q_SUFF=_1Y_grid_T.nc 
     77VAR2_EMP=wfo            ;   V2EMP_PREF=testht   ;   V2EMP_SUFF=_1Y_grid_T.nc 
     78VAR2_U=uocetr_eff      ;   V2U_PREF=testht     ;   V2U_SUFF=_1Y_grid_U.nc 
     79VAR2_V=vocetr_eff      ;   V2V_PREF=testht     ;   V2V_SUFF=_1Y_grid_V.nc 
     80VAR2_ICE=siconc         ;   V2ICE_PREF=testht   ;   V2ICE_SUFF=_1M_icemod.nc 
     81VAR2_Ithick=sithic       ;   V2It_PREF=testht    ;   V2It_SUFF=_1M_icemod.nc 
     82VAR2_SNOW=snthic       ;   V2SNOW_PREF=testht  ;   V2SNOW_SUFF=_1M_icemod.nc 
     83VAR2_IvelV=sivelv       ;   V2IvV_PREF=testht   ;   V2IvV_SUFF=_1M_icemod.nc  
     84VAR2_Ivel=sivelo        ;   V2IvV_PREF=testht   ;   V2IvV_SUFF=_1M_icemod.nc  
    8985# 
    9086######################### Export Variables ############################### 
Note: See TracChangeset for help on using the changeset viewer.