Ignore:
Timestamp:
2014-02-25T12:41:51+01:00 (7 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_last_year_5D.pro

    r4320 r4514  
    203203  ENDELSE 
    204204; 
     205  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
     206  caldat, time, mm 
     207  march = where(mm EQ 3, cnt) 
     208  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     209  sept = where(mm EQ 9, cnt) 
     210  Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit}  
     211  undefine, Iage_1 
     212; 
     213  IF std_file2_I NE std_file1_I THEN BEGIN 
     214     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
     215     caldat, time, mm 
     216     march = where(mm EQ 3, cnt) 
     217     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     218     sept = where(mm EQ 9, cnt) 
     219     Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     220     undefine, Iage_2 
     221  ENDIF ELSE BEGIN  
     222     Iage_march_2 = {arr:-1} 
     223     Iage_sept_2 = {arr:-1} 
     224  ENDELSE 
     225; 
    205226  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    206227  ;; help, Ifra_1.arr, jpi,jpj 
     
    225246  ENDELSE 
    226247; 
    227 ;; SF commentato perché non funziona snow 
    228 ;; SF   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
    229 ;; SF   help, Isnow_1.arr, jpi,jpj 
    230 ;; SF   caldat, time, mm 
    231 ;; SF   april = where(mm EQ 4, cnt) 
    232 ;; SF   Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
    233 ;; SF   sept = where(mm EQ 9, cnt) 
    234 ;; SF   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
    235 ;; SF   ;undefine, Isnow_1 
    236 ;; SF ; 
    237 ;; SF   IF std_file2_I NE std_file1_I THEN BEGIN 
    238 ;; SF     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
    239 ;; SF     caldat, time, mm 
    240 ;; SF     april = where(mm EQ 4, cnt) 
    241 ;; SF     Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    242 ;; SF     sept = where(mm EQ 9, cnt) 
    243 ;; SF     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    244 ;; SF   ;  undefine, Isnow_2 
    245 ;; SF   ENDIF ELSE BEGIN  
    246 ;; SF     Isnow_april_2 = {arr:-1} 
    247 ;; SF     Isnow_sept_2 = {arr:-1} 
    248 ;; SF   ENDELSE 
     248   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
     249   caldat, time, mm 
     250   march = where(mm EQ 3, cnt) 
     251   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
     252   sept = where(mm EQ 9, cnt) 
     253   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
     254   ;undefine, Isnow_1 
     255 ; 
     256   IF std_file2_I NE std_file1_I THEN BEGIN 
     257     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
     258     caldat, time, mm 
     259     march = where(mm EQ 3, cnt) 
     260     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     261     sept = where(mm EQ 9, cnt) 
     262     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     263   ;  undefine, Isnow_2 
     264   ENDIF ELSE BEGIN  
     265     Isnow_march_2 = {arr:-1} 
     266     Isnow_sept_2 = {arr:-1} 
     267   ENDELSE 
     268; 
     269   Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I )  
     270   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
     271   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations  
     272   help, Ifra_1 
     273   help, mask 
     274   caldat, time, mm 
     275   march = where(mm EQ 3, cnt) 
     276   febr = where(mm EQ 2, cnt) 
     277   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     278   Isal_1.arr = Isal_1.arr * msk 
     279   help, Isal_1 
     280   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     281   help, Isal_march_1 
     282   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     283   sept = where(mm EQ 9, cnt) 
     284   Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
     285   ;undefine, Isal_1 
     286 ; 
     287   Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     288   Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 
     289   msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations  
     290   IF std_file2_I NE std_file1_I THEN BEGIN 
     291     Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     292     Isal_2.arr = Isal_2.arr * msk 
     293     caldat, time, mm 
     294     march = where(mm EQ 3, cnt) 
     295     febr = where(mm EQ 2, cnt) 
     296     Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     297     Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     298     sept = where(mm EQ 9, cnt) 
     299     Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     300   ;  undefine, Isal_2 
     301   ENDIF ELSE BEGIN  
     302     Isal_febr_2 = {arr:-1} 
     303     Isal_march_2 = {arr:-1} 
     304     Isal_sept_2 = {arr:-1} 
     305   ENDELSE 
    249306; 
    250307  jpt = 1 
     
    282339  cnt = cnt+1   &   blabla = 'Erp salinity damping term' 
    283340  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 
    284  
     341; 
    285342  cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term' 
    286343  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex 
    287  
     344; 
    288345  cnt = cnt+1   &   blabla = 'Net heat flux' 
    289346  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex 
    290  
     347; 
    291348  cnt = cnt+1   &   blabla = 'Meridionnal Heat Transport' 
    292349  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex 
    293  
     350; 
    294351  cnt = cnt+1   &   blabla = 'Global Barotropic stream Function' 
    295352  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex 
    296  
     353; 
    297354  cnt = cnt+1   &   blabla = 'mean Temperature diff with New Reynolds' 
    298355  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex 
    299  
     356; 
    300357  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus' 
    301358  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex 
    302  
     359; 
    303360  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus' 
    304361  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex        
    305  
     362; 
    306363  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2' 
    307364  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex 
    308  
     365; 
    309366  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters' 
    310367  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex        
    311   
     368;  
    312369  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters' 
    313370  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex 
    314  
     371; 
    315372  cnt = cnt+1   &   blabla = 'mean Temperature diff with Levitus at z=100 meters' 
    316373  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex 
    317    
     374;   
    318375  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus at z=100 meters' 
    319376  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex 
    320   
    321 ;;SF need monthly climatology to plot difference of MAX of MLD   
    322   cnt = cnt+1   &   blabla = 'Mixed layer depth & difference with DeBoyer' 
    323   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD, _extra = ex 
    324     
     377 
     378  cnt = cnt+1   &   blabla = 'Mixed layer depth' 
     379  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD, _extra = ex 
     380; 
    325381  cnt = cnt+1   &   blabla = 'Mixed layer depth differences' 
    326382  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD_2, MLD, _extra = ex 
     
    339395  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Global' 
    340396  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex 
    341  
     397; 
    342398  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Atlantic' 
    343399  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex 
    344    
     400;   
    345401  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Indian' 
    346402  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex 
    347    
     403;   
    348404  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Pacific' 
    349405  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex 
    350  
     406; 
    351407  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Global' 
    352408  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex 
    353  
     409; 
    354410  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Atlantic' 
    355411  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex 
    356  
     412; 
    357413  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Indian' 
    358414  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex  
    359  
     415; 
    360416  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Pacific' 
    361417  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex  
    362    
     418;   
    363419;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
    364420;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
    365  
     421; 
    366422  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
    367423  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 
    368  
     424; 
    369425  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT' 
    370426  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 
    371    
     427;   
    372428;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
    373429;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 
    374  
     430;  
    375431  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
    376432  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 
    377  
     433; 
    378434  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT' 
    379435  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 
    380  
     436; 
     437  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
     438  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 
     439 
     440  cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT' 
     441  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 
     442 
     443  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH' 
     444  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 
     445; 
     446  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
     447  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 
     448; 
    381449  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH' 
    382450  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
    383    
     451;   
    384452  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: SEPT' 
    385453  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 
    386    
     454;   
    387455  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH' 
    388456  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex 
    389  
     457; 
    390458  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT' 
    391459  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
    392  
    393 ;; SF non funziona, dice che deve essere una struttura e non e!!!! 
    394 ;;SF  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL' 
    395 ;;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 
    396 ;; SF  
    397 ;; SF   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
    398 ;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
    399 ;; SF   
    400 ;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL' 
    401 ;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 
    402 ;; SF   
    403 ;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
    404 ;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
    405 ;  
    406 ;  cnt = cnt+1   &   blabla = 'Arctic Ice Extend: ' 
    407 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex 
    408    
    409 ;  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: ' 
    410 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex 
    411  
     460; 
     461  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: MARCH' 
     462  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 
     463;  
     464  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
     465  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
     466 
     467  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH' 
     468  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 
     469;  
     470  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
     471  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
     472;  
     473  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: MARCH' 
     474  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 
     475;  
     476  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT' 
     477  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 
     478 
     479  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: FEBRUARY' 
     480  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 
     481;  
     482  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT' 
     483  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex 
     484;  
    412485  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)' 
    413486  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex 
    414  
     487 
    415488  cnt = cnt+1   &   blabla = 'Meridional stream Function: Atlantic' 
    416489  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex 
    417  
     490 ; 
    418491  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indian' 
    419492  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex 
    420  
     493 ; 
    421494  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indo-Pacific' 
    422495  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex 
    423  
     496; 
    424497  cnt = cnt+1   &   blabla = 'Equatorial Temperature' 
    425498  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex 
    426  
     499; 
    427500  cnt = cnt+1   &   blabla = 'Equatorial Salinity' 
    428501  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex 
    429  
     502; 
    430503  cnt = cnt+1   &   blabla = 'Equatorial zonal velocity' 
    431504  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex 
    432  
     505; 
    433506  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=700' 
    434507  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex 
    435  
     508; 
    436509  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=1000' 
    437510  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex 
    438  
     511; 
    439512  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40°N' 
    440513  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex 
    441  
     514; 
    442515  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38°N' 
    443516  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 
    444  
     517; 
    445518  cnt = cnt+1   &   blabla = 'Vertical Global mean T & S' 
    446519  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex 
    447  
     520; 
    448521  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*] 
    449522   
Note: See TracChangeset for help on using the changeset viewer.