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

Ignore:
Timestamp:
2014-02-25T12:41:51+01:00 (10 years ago)
Author:
flavoni
Message:

add ice pots & time series, 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

    r4507 r4514  
    187187  ENDELSE 
    188188; 
     189  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
     190  caldat, time, mm 
     191  march = where(mm EQ 3, cnt) 
     192  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     193  sept = where(mm EQ 9, cnt) 
     194  Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit}  
     195  undefine, Iage_1 
     196; 
     197  IF std_file2_I NE std_file1_I THEN BEGIN 
     198     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
     199     caldat, time, mm 
     200     march = where(mm EQ 3, cnt) 
     201     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     202     sept = where(mm EQ 9, cnt) 
     203     Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     204     undefine, Iage_2 
     205  ENDIF ELSE BEGIN  
     206     Iage_march_2 = {arr:-1} 
     207     Iage_sept_2 = {arr:-1} 
     208  ENDELSE 
     209; 
    189210  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    190   ;help, Ifra_1.arr, jpi,jpj 
    191211  caldat, time, mm 
    192212  march = where(mm EQ 3, cnt) 
    … …  
    210230; 
    211231   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
    212    ;help, Isnow_1.arr, jpi,jpj 
    213232   caldat, time, mm 
    214    april = where(mm EQ 4, cnt) 
    215    Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
     233   march = where(mm EQ 3, cnt) 
     234   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
    216235   sept = where(mm EQ 9, cnt) 
    217236   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
    … …  
    221240     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
    222241     caldat, time, mm 
    223      april = where(mm EQ 4, cnt) 
    224      Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     242     march = where(mm EQ 3, cnt) 
     243     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    225244     sept = where(mm EQ 9, cnt) 
    226245     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    227246   ;  undefine, Isnow_2 
    228247   ENDIF ELSE BEGIN  
    229      Isnow_april_2 = {arr:-1} 
     248     Isnow_march_2 = {arr:-1} 
    230249     Isnow_sept_2 = {arr:-1} 
    231250   ENDELSE 
    232251; 
    233252   Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I )  
    234    ;help, Isal_1.arr, jpi,jpj 
     253   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
     254   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations  
     255   help, Ifra_1 
     256   help, mask 
    235257   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}  
     258   march = where(mm EQ 3, cnt) 
     259   febr = where(mm EQ 2, cnt) 
     260   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     261   Isal_1.arr = Isal_1.arr * msk 
     262   help, Isal_1 
     263   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     264   help, Isal_march_1 
     265   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    238266   sept = where(mm EQ 9, cnt) 
    239267   Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
    240268   ;undefine, Isal_1 
    241269 ; 
     270   Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     271   Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 
     272   msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations  
    242273   IF std_file2_I NE std_file1_I THEN BEGIN 
    243274     Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     275     Isal_2.arr = Isal_2.arr * msk 
    244276     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} 
     277     march = where(mm EQ 3, cnt) 
     278     febr = where(mm EQ 2, cnt) 
     279     Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     280     Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
    247281     sept = where(mm EQ 9, cnt) 
    248282     Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
    249283   ;  undefine, Isal_2 
    250284   ENDIF ELSE BEGIN  
    251      Isal_april_2 = {arr:-1} 
     285     Isal_febr_2 = {arr:-1} 
     286     Isal_march_2 = {arr:-1} 
    252287     Isal_sept_2 = {arr:-1} 
    253288   ENDELSE 
    … …  
    375410  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 
    376411; 
     412  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
     413  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 
     414;   
     415  cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT' 
     416  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 
     417;   
     418  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH' 
     419  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 
     420; 
     421  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
     422  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 
     423; 
    377424  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH' 
    378425  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
    … …  
    387434  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
    388435; 
    389   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL' 
    390   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 
     436  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: MARCH' 
     437  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 
    391438;  
    392439  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
    393440  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
    394441;   
    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 
     442  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH' 
     443  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 
    397444;  
    398445  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
    399446  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
    400447;  
    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 
     448  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: MARCH' 
     449  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 
    403450;  
    404451  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT' 
    405452  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 
    406453;   
    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 
     454  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: FEBRUARY' 
     455  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 
    409456;  
    410457  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT' 
Note: See TracChangeset for help on using the changeset viewer.