Changeset 4507


Ignore:
Timestamp:
2014-02-21T16:05:20+01:00 (7 years ago)
Author:
flavoni
Message:

commit ICE Salinity plots for ORCA2_LIM, see ticket #724

Location:
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
Files:
2 added
11 edited

Legend:

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

    r4501 r4507  
    188188; 
    189189  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    190   help, Ifra_1.arr, jpi,jpj 
     190  ;help, Ifra_1.arr, jpi,jpj 
    191191  caldat, time, mm 
    192192  march = where(mm EQ 3, cnt) 
     
    210210; 
    211211   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
    212    help, Isnow_1.arr, jpi,jpj 
     212   ;help, Isnow_1.arr, jpi,jpj 
    213213   caldat, time, mm 
    214214   april = where(mm EQ 4, cnt) 
     
    231231   ENDELSE 
    232232; 
     233   Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I )  
     234   ;help, Isal_1.arr, jpi,jpj 
     235   caldat, time, mm 
     236   april = where(mm EQ 4, cnt) 
     237   Isal_april_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     238   sept = where(mm EQ 9, cnt) 
     239   Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
     240   ;undefine, Isal_1 
     241 ; 
     242   IF std_file2_I NE std_file1_I THEN BEGIN 
     243     Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     244     caldat, time, mm 
     245     april = where(mm EQ 4, cnt) 
     246     Isal_april_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     247     sept = where(mm EQ 9, cnt) 
     248     Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     249   ;  undefine, Isal_2 
     250   ENDIF ELSE BEGIN  
     251     Isal_april_2 = {arr:-1} 
     252     Isal_sept_2 = {arr:-1} 
     253   ENDELSE 
     254; 
    233255  jpt = 1 
    234256; 
     
    368390  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 
    369391;  
    370    cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
    371    IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
    372  
    373    cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL' 
    374    IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 
    375 ;  
    376    cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
    377    IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
    378 ;  
    379 ; NOT NEEDED, only for TIME SERIES NOW ICE EXTEND 
    380 ;SF  cnt = cnt+1   &   blabla = 'Arctic Ice Extend: ' 
    381 ;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex 
    382 ;SF   
    383 ;SF  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: ' 
    384 ;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex 
    385 ; 
     392  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
     393  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
     394 
     395  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL' 
     396  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 
     397;  
     398  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
     399  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
     400;  
     401  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: APRIL' 
     402  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_april_1, Isal_april_2, /ARC, /APRIL, _extra = ex 
     403;  
     404  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT' 
     405  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 
     406 
     407  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: APRIL' 
     408  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_april_1, Isal_april_2, /ANT, /APRIL, _extra = ex 
     409;  
     410  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT' 
     411  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex 
     412;  
    386413  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)' 
    387414  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh

    r4504 r4507  
    2121#idl_command=idl71 
    2222# 
    23 PS_DIR=$( pwd )/sbcmodMV_ps_plot 
    24 PDF_DIR=$( pwd )/sbcmodMV_pdf_plot 
     23PS_DIR=$( pwd )/sbcmodMV-l334plus_ps_plot 
     24PDF_DIR=$( pwd )/sbcmodMV-l334plus_pdf_plot 
    2525HTML_DIR=$( pwd )/html_plot 
    2626SAXO_DIR=/Users/sflod/SAXO_DIR 
     
    7070VAR1_Ifrac=siconc 
    7171VAR1_Isnow=sndept 
     72VAR1_Isal=sisali 
    7273# 
    7374#===================== EXP2 ===================== 
    7475# 
    75 FILE2_T=sbcmodMV_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file 
    76 FILE2_U=sbcmodMV_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file 
    77 FILE2_V=sbcmodMV_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file 
    78 FILE2_I=sbcmodMV_20410101_20501231_1M_icemod.nc   # exp1 ice    input file 
     76FILE2_T=l334plus_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file 
     77FILE2_U=l334plus_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file 
     78FILE2_V=l334plus_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file 
     79FILE2_I=l334plus_20410101_20501231_1M_icemod.nc   # exp1 ice    input file 
    7980VAR2_T=thetao 
    8081VAR2_S=so 
     
    8889VAR2_Ifrac=siconc 
    8990VAR2_Isnow=sndept 
     91VAR2_Isal=sisali 
    9092# 
    9193######################### Export Variables ############################### 
     
    111113export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    112114export VAR1_U      VAR1_V 
    113 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     115export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    114116#===================== EXP2 ===================== 
    115117export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    116118export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    117119export VAR2_U      VAR2_V 
    118 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     120export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    119121# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.500yfwb0_100y

    r4504 r4507  
    7878VAR1_Ifrac=iiceconc 
    7979VAR1_Isnow=isnowthi 
     80VAR1_Isal=iicesali 
    8081# 
    8182#===================== EXP2 ===================== 
     
    102103VAR2_Ifrac=iiceconc 
    103104VAR2_Isnow=isnowthi 
     105VAR2_Isal=iicesali 
    104106# 
    105107######################### Export Variables ############################### 
     
    125127export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    126128export VAR1_U      VAR1_V 
    127 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     129export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    128130#===================== EXP2 ===================== 
    129131export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    130132export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    131133export VAR2_U      VAR2_V 
    132 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     134export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    133135# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.500yfwb0_500y

    r4504 r4507  
    7878VAR1_Ifrac=iiceconc 
    7979VAR1_Isnow=isnowthi 
     80VAR1_Isal=iicesali 
    8081# 
    8182#===================== EXP2 ===================== 
     
    102103VAR2_Ifrac=iiceconc 
    103104VAR2_Isnow=isnowthi 
     105VAR2_Isal=iicesali 
    104106# 
    105107######################### Export Variables ############################### 
     
    125127export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    126128export VAR1_U      VAR1_V 
    127 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     129export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    128130#===================== EXP2 ===================== 
    129131export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    130132export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    131133export VAR2_U      VAR2_V 
    132 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     134export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    133135# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.l378sal6

    r4504 r4507  
    7878VAR1_Ifrac=iiceconc 
    7979VAR1_Isnow=isnowthi 
     80VAR1_Isal=iicesali 
    8081# 
    8182#===================== EXP2 ===================== 
     
    102103VAR2_Ifrac=iiceconc 
    103104VAR2_Isnow=isnowthi 
     105VAR2_Isal=iicesali 
    104106# 
    105107######################### Export Variables ############################### 
     
    125127export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    126128export VAR1_U      VAR1_V 
    127 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     129export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    128130#===================== EXP2 ===================== 
    129131export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    130132export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    131133export VAR2_U      VAR2_V 
    132 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     134export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    133135# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.l3fwb178-500yfwb0

    r4504 r4507  
    7878VAR1_Ifrac=iiceconc 
    7979VAR1_Isnow=isnowthi 
     80VAR1_Isal=iicesali 
    8081# 
    8182#===================== EXP2 ===================== 
     
    9899VAR2_Ifrac=iiceconc 
    99100VAR2_Isnow=isnowthi 
     101VAR2_Isal=iicesali 
    100102# 
    101103######################### Export Variables ############################### 
     
    121123export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    122124export VAR1_U      VAR1_V 
    123 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     125export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    124126#===================== EXP2 ===================== 
    125127export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    126128export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    127129export VAR2_U      VAR2_V 
    128 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     130export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    129131# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.last_year_5D

    r4504 r4507  
    8181VAR1_Ifrac=iiceconc 
    8282VAR1_Isnow=isnowthi 
     83VAR1_Isal=iicesali 
    8384# 
    8485#===================== EXP2 ===================== 
     
    112113VAR2_Ifrac=iiceconc 
    113114VAR2_Isnow=isnowthi 
     115VAR2_Isal=iicesali 
    114116# 
    115117######################### Export Variables ############################### 
     
    135137export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    136138export VAR1_U      VAR1_V 
    137 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     139export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    138140#===================== EXP2 ===================== 
    139141export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    140142export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    141143export VAR2_U      VAR2_V 
    142 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     144export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    143145# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.new_names

    r4504 r4507  
    7070VAR1_Ifrac=siconc 
    7171VAR1_Isnow=sndept 
     72VAR1_Isal=sisali 
    7273# 
    7374#===================== EXP2 ===================== 
     
    8889VAR2_Ifrac=siconc 
    8990VAR2_Isnow=sndept 
     91VAR2_Isal=sisali 
    9092# 
    9193######################### Export Variables ############################### 
     
    111113export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    112114export VAR1_U      VAR1_V 
    113 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     115export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    114116#===================== EXP2 ===================== 
    115117export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    116118export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    117119export VAR2_U      VAR2_V 
    118 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     120export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    119121# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.qsriceok

    r4504 r4507  
    7878VAR1_Ifrac=iiceconc 
    7979VAR1_Isnow=isnowthi 
     80VAR1_Isal=iicesali 
    8081# 
    8182#===================== EXP2 ===================== 
     
    102103VAR2_Ifrac=iiceconc 
    103104VAR2_Isnow=isnowthi 
     105VAR2_Isal=iicesali 
    104106# 
    105107######################### Export Variables ############################### 
     
    125127export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    126128export VAR1_U      VAR1_V 
    127 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     129export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    128130#===================== EXP2 ===================== 
    129131export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    130132export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    131133export VAR2_U      VAR2_V 
    132 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     134export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    133135# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1

    r4504 r4507  
    7070VAR1_Ifrac=sic 
    7171VAR1_Isnow=isnowthi 
     72VAR1_Isal=iicesali 
    7273# 
    7374#===================== EXP2 ===================== 
     
    8889VAR2_Ifrac=sic 
    8990VAR2_Isnow=isnowthi 
     91VAR2_Isal=iicesali 
    9092# 
    9193######################### Export Variables ############################### 
     
    111113export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    112114export VAR1_U      VAR1_V 
    113 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     115export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    114116#===================== EXP2 ===================== 
    115117export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    116118export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    117119export VAR2_U      VAR2_V 
    118 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     120export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    119121# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2

    r4504 r4507  
    7070VAR1_Ifrac=sic 
    7171VAR1_Isnow=isnowthi 
     72VAR1_Isal=iicesali 
    7273# 
    7374#===================== EXP2 ===================== 
     
    8889VAR2_Ifrac=sic 
    8990VAR2_Isnow=isnowthi 
     91VAR2_Isal=iicesali 
    9092# 
    9193######################### Export Variables ############################### 
     
    111113export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    112114export VAR1_U      VAR1_V 
    113 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow 
     115export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
    114116#===================== EXP2 ===================== 
    115117export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    116118export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    117119export VAR2_U      VAR2_V 
    118 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow 
     120export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
    119121# 
Note: See TracChangeset for help on using the changeset viewer.