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

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

change months in ice plot, see ticket #724

File:
1 edited

Legend:

Unmodified
Added
Removed

  • 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' 
Note: See TracChangeset for help on using the changeset viewer.