Ignore:
Timestamp:
2014-02-25T16:20:43+01:00 (7 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_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.