Changeset 4518


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

change months in ice plot, see ticket #724

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

Legend:

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

    r4516 r4518  
    1 pro std_plot_IceAge, Iage1, Iage2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_IceAge, Iage1, Iage2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    1212  IF keyword_set(arc) THEN var = var+'_Arc_' 
    1313  IF keyword_set(ant) THEN var = var+'_Ant_' 
     14  IF keyword_set(febr) THEN var = var+'Febr' 
    1415  IF keyword_set(march) THEN var = var+'March' 
    1516  IF keyword_set(sept) THEN var = var+'Sept' 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro

    r4320 r4518  
    1 pro std_plot_IceFrac, Ifra1, Ifra2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_IceFrac, Ifra1, Ifra2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    1212  IF keyword_set(arc) THEN var = var+'_Arc_' 
    1313  IF keyword_set(ant) THEN var = var+'_Ant_' 
     14  IF keyword_set(febr) THEN var = var+'Febr' 
    1415  IF keyword_set(march) THEN var = var+'March' 
    1516  IF keyword_set(sept) THEN var = var+'Sept' 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceThick.pro

    r4320 r4518  
    1 pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, APRIL = april, JAN = jan, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    1313  IF keyword_set(ant) THEN var = var+'_Ant_' 
    1414  IF keyword_set(april) THEN var = var+'April' 
     15  IF keyword_set(jan) THEN var = var+'Jan' 
    1516  IF keyword_set(sept) THEN var = var+'Sept' 
    1617 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro

    r4514 r4518  
    1 pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    1212  IF keyword_set(arc) THEN var = var+'_Arc_' 
    1313  IF keyword_set(ant) THEN var = var+'_Ant_' 
     14  IF keyword_set(febr) THEN var = var+'Febr' 
    1415  IF keyword_set(march) THEN var = var+'March' 
    1516  IF keyword_set(sept) THEN var = var+'Sept' 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r4517 r4518  
    103103;climatology  
    104104  QNET = read_ncdf(getenv('VAR_FLUX'), filename = std_file_oaflux ) 
    105   
     105 
    106106; erp (evaporation damping) 
    107107  ERP1 = read_ncdf(getenv('VAR1_ERP'), allrecords = allrec, direc = 't', filename = std_file1_T ) 
     
    170170  april = where(mm EQ 4, cnt) 
    171171  Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 
     172  jan = where(mm EQ 1, cnt) 
    172173  sept = where(mm EQ 9, cnt) 
     174  Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    173175  Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    174176  undefine, Ithi_1 
     
    179181     april = where(mm EQ 4, cnt) 
    180182     Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
     183     jan = where(mm EQ 1, cnt) 
    181184     sept = where(mm EQ 9, cnt) 
    182185     Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
     186     Ithi_jan_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
    183187     undefine, Ithi_2 
    184188  ENDIF ELSE BEGIN  
    185189     Ithi_april_2 = {arr:-1} 
    186190     Ithi_sept_2 = {arr:-1} 
     191     Ithi_jan_2 = {arr:-1} 
    187192  ENDELSE 
    188193; 
    189194  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
    190195  caldat, time, mm 
     196  febr = where(mm EQ 2, cnt) 
     197  Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
    191198  march = where(mm EQ 3, cnt) 
    192199  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     
    198205     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
    199206     caldat, time, mm 
     207     febr = where(mm EQ 2, cnt) 
     208     Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
    200209     march = where(mm EQ 3, cnt) 
    201210     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     
    204213     undefine, Iage_2 
    205214  ENDIF ELSE BEGIN  
     215     Iage_febr_2 = {arr:-1} 
    206216     Iage_march_2 = {arr:-1} 
    207217     Iage_sept_2 = {arr:-1} 
     
    210220  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    211221  caldat, time, mm 
     222  febr = where(mm EQ 2, cnt) 
     223  Ifra_febr_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
    212224  march = where(mm EQ 3, cnt) 
    213225  Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
     
    219231    Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )  
    220232    caldat, time, mm 
     233    febr = where(mm EQ 2, cnt) 
     234    Ifra_febr_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
    221235    march = where(mm EQ 3, cnt) 
    222236    Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
     
    225239    undefine, Ifra_2 
    226240  ENDIF ELSE BEGIN  
     241    Ifra_febr_2 = {arr:-1} 
    227242    Ifra_march_2 = {arr:-1} 
    228243    Ifra_sept_2 = {arr:-1} 
     
    231246   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
    232247   caldat, time, mm 
     248   febr = where(mm EQ 2, cnt) 
     249   Isnow_febr_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
    233250   march = where(mm EQ 3, cnt) 
    234251   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
     
    240257     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
    241258     caldat, time, mm 
     259     febr = where(mm EQ 2, cnt) 
     260     Isnow_febr_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    242261     march = where(mm EQ 3, cnt) 
    243262     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     
    246265   ;  undefine, Isnow_2 
    247266   ENDIF ELSE BEGIN  
     267     Isnow_febr_2 = {arr:-1} 
    248268     Isnow_march_2 = {arr:-1} 
    249269     Isnow_sept_2 = {arr:-1} 
     
    428448  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 
    429449; 
    430 ;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
    431 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
     450  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: JAN' 
     451  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, /ARC, /JAN, _extra = ex 
    432452; 
    433453  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
    434454  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 
    435 ; 
    436   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT' 
    437   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 
    438  
    439 ;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
    440 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 
    441 ;  
     455 
    442456  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
    443457  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 
     
    453467 
    454468  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH' 
    455   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 
     469  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, /ANT, /FEBR, _extra = ex 
    456470; 
    457471  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
     
    465479 
    466480  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH' 
    467   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex 
     481  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 
    468482; 
    469483  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT' 
     
    477491 
    478492  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH' 
    479   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 
     493  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, /ANT, /FEBR, _extra = ex 
    480494;  
    481495  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_1Y.pro

    r4517 r4518  
    1 pro std_plot_all, doplot = doplot, _extra = ex 
     1pro std_plot_all_1Y, doplot = doplot, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    103103;climatology  
    104104  QNET = read_ncdf(getenv('VAR_FLUX'), filename = std_file_oaflux ) 
    105     
     105 
    106106; erp (evaporation damping) 
    107107  ERP1 = read_ncdf(getenv('VAR1_ERP'), allrecords = allrec, direc = 't', filename = std_file1_T ) 
     
    170170  april = where(mm EQ 4, cnt) 
    171171  Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 
     172  jan = where(mm EQ 1, cnt) 
    172173  sept = where(mm EQ 9, cnt) 
     174  Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    173175  Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    174176  undefine, Ithi_1 
     
    179181     april = where(mm EQ 4, cnt) 
    180182     Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
     183     jan = where(mm EQ 1, cnt) 
    181184     sept = where(mm EQ 9, cnt) 
    182185     Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
     186     Ithi_jan_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
    183187     undefine, Ithi_2 
    184188  ENDIF ELSE BEGIN  
    185189     Ithi_april_2 = {arr:-1} 
    186190     Ithi_sept_2 = {arr:-1} 
     191     Ithi_jan_2 = {arr:-1} 
    187192  ENDELSE 
    188193; 
    189194  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
    190195  caldat, time, mm 
     196  febr = where(mm EQ 2, cnt) 
     197  Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
    191198  march = where(mm EQ 3, cnt) 
    192199  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     
    198205     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
    199206     caldat, time, mm 
     207     febr = where(mm EQ 2, cnt) 
     208     Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
    200209     march = where(mm EQ 3, cnt) 
    201210     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     
    204213     undefine, Iage_2 
    205214  ENDIF ELSE BEGIN  
     215     Iage_febr_2 = {arr:-1} 
    206216     Iage_march_2 = {arr:-1} 
    207217     Iage_sept_2 = {arr:-1} 
     
    210220  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    211221  caldat, time, mm 
     222  febr = where(mm EQ 2, cnt) 
     223  Ifra_febr_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
    212224  march = where(mm EQ 3, cnt) 
    213225  Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
     
    219231    Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )  
    220232    caldat, time, mm 
     233    febr = where(mm EQ 2, cnt) 
     234    Ifra_febr_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
    221235    march = where(mm EQ 3, cnt) 
    222236    Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
     
    225239    undefine, Ifra_2 
    226240  ENDIF ELSE BEGIN  
     241    Ifra_febr_2 = {arr:-1} 
    227242    Ifra_march_2 = {arr:-1} 
    228243    Ifra_sept_2 = {arr:-1} 
     
    231246   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
    232247   caldat, time, mm 
     248   febr = where(mm EQ 2, cnt) 
     249   Isnow_febr_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
    233250   march = where(mm EQ 3, cnt) 
    234251   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
     
    240257     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
    241258     caldat, time, mm 
     259     febr = where(mm EQ 2, cnt) 
     260     Isnow_febr_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    242261     march = where(mm EQ 3, cnt) 
    243262     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     
    246265   ;  undefine, Isnow_2 
    247266   ENDIF ELSE BEGIN  
     267     Isnow_febr_2 = {arr:-1} 
    248268     Isnow_march_2 = {arr:-1} 
    249269     Isnow_sept_2 = {arr:-1} 
     
    428448  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 
    429449; 
    430 ;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
    431 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
     450  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: JAN' 
     451  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, /ARC, /JAN, _extra = ex 
    432452; 
    433453  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
    434454  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 
    435 ; 
    436   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT' 
    437   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 
    438  
    439 ;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
    440 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 
    441 ;  
     455 
    442456  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
    443457  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 
     
    453467 
    454468  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH' 
    455   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 
     469  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, /ANT, /FEBR, _extra = ex 
    456470; 
    457471  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
     
    465479 
    466480  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH' 
    467   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex 
     481  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 
    468482; 
    469483  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT' 
     
    477491 
    478492  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH' 
    479   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 
     493  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, /ANT, /FEBR, _extra = ex 
    480494;  
    481495  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro

    r4517 r4518  
    105105;climatology  
    106106  QNET = read_ncdf(getenv('VAR_FLUX'), filename = std_file_oaflux ) 
    107   
     107 
    108108; erp (evaporation damping) 
    109109  ERP1 = read_ncdf(getenv('VAR1_ERP'), allrecords = allrec, direc = 't', filename = std_file1_T ) 
     
    186186  april = where(mm EQ 4, cnt) 
    187187  Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 
     188  jan = where(mm EQ 1, cnt) 
    188189  sept = where(mm EQ 9, cnt) 
     190  Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    189191  Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    190192  undefine, Ithi_1 
     
    195197     april = where(mm EQ 4, cnt) 
    196198     Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
     199     jan = where(mm EQ 1, cnt) 
    197200     sept = where(mm EQ 9, cnt) 
    198201     Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
     202     Ithi_jan_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_2.unit}  
    199203     undefine, Ithi_2 
    200204  ENDIF ELSE BEGIN  
    201205     Ithi_april_2 = {arr:-1} 
    202206     Ithi_sept_2 = {arr:-1} 
     207     Ithi_jan_2 = {arr:-1} 
    203208  ENDELSE 
    204209; 
    205210  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
    206211  caldat, time, mm 
     212  febr = where(mm EQ 2, cnt) 
     213  Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
    207214  march = where(mm EQ 3, cnt) 
    208215  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     
    214221     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
    215222     caldat, time, mm 
     223     febr = where(mm EQ 2, cnt) 
     224     Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
    216225     march = where(mm EQ 3, cnt) 
    217226     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     
    220229     undefine, Iage_2 
    221230  ENDIF ELSE BEGIN  
     231     Iage_febr_2 = {arr:-1} 
    222232     Iage_march_2 = {arr:-1} 
    223233     Iage_sept_2 = {arr:-1} 
     
    225235; 
    226236  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    227   ;; help, Ifra_1.arr, jpi,jpj 
    228237  caldat, time, mm 
     238  febr = where(mm EQ 2, cnt) 
     239  Ifra_febr_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
    229240  march = where(mm EQ 3, cnt) 
    230241  Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
     
    236247    Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )  
    237248    caldat, time, mm 
     249    febr = where(mm EQ 2, cnt) 
     250    Ifra_febr_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
    238251    march = where(mm EQ 3, cnt) 
    239252    Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
     
    242255    undefine, Ifra_2 
    243256  ENDIF ELSE BEGIN  
     257    Ifra_febr_2 = {arr:-1} 
    244258    Ifra_march_2 = {arr:-1} 
    245259    Ifra_sept_2 = {arr:-1} 
     
    248262   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
    249263   caldat, time, mm 
     264   febr = where(mm EQ 2, cnt) 
     265   Isnow_febr_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
    250266   march = where(mm EQ 3, cnt) 
    251267   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
     
    257273     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
    258274     caldat, time, mm 
     275     febr = where(mm EQ 2, cnt) 
     276     Isnow_febr_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    259277     march = where(mm EQ 3, cnt) 
    260278     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     
    263281   ;  undefine, Isnow_2 
    264282   ENDIF ELSE BEGIN  
     283     Isnow_febr_2 = {arr:-1} 
    265284     Isnow_march_2 = {arr:-1} 
    266285     Isnow_sept_2 = {arr:-1} 
     
    450469  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 
    451470; 
    452 ;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
    453 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
     471  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: JAN' 
     472  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, /ARC, /JAN, _extra = ex 
    454473; 
    455474  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
    456475  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 
    457 ; 
    458   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT' 
    459   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 
    460  
    461 ;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
    462 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 
    463 ;  
     476 
    464477  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
    465478  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 
     
    475488 
    476489  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH' 
    477   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 
     490  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, /ANT, /FEBR, _extra = ex 
    478491; 
    479492  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
     
    487500 
    488501  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH' 
    489   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex 
     502  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 
    490503; 
    491504  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT' 
     
    499512 
    500513  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH' 
    501   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 
     514  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, /ANT, /FEBR, _extra = ex 
    502515;  
    503516  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
Note: See TracChangeset for help on using the changeset viewer.