Changeset 43


Ignore:
Timestamp:
04/11/11 18:05:51 (13 years ago)
Author:
pinsard
Message:

add minimal header in paper01 IDL files

Location:
trunk/src/paper01
Files:
41 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/paper01/fig1/fig1_no_obs_map.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _fig1_no_obs_map.pro 
     3; 
     4; =================== 
     5; fig1_no_obs_map.pro 
     6; =================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> fig1_no_obs_map 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    234pro fig1_no_obs_map 
    335@common 
    4 ;------------------------------------------------------------ 
    536 
    637reinitplt, /z,/invert 
    … …  
    839 
    940openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1141close,/all 
    1242 
    … …  
    136166plots, 202.46, 22.45, psym=2  ;; WHOTS 
    137167 
    138 ;----------------------------------------------------------- 
    139168closeps 
    140169 
    … …  
    144173return 
    145174end 
    146 ;-------------------------------------------------------------------------- 
     175 
    147176function x_site_location, site 
    148177    n1=strpos(site, 's') 
    … …  
    159188return, float(x) 
    160189end 
    161 ;-------------------------------------------------------------------------- 
     190 
    162191function y_site_location, site 
    163192    n1=strpos(site, 'e') 
    … …  
    185214return,float(y) 
    186215end 
    187  
    188 ;-------------------------------------------------------------------------- 
    189  

  • trunk/src/paper01/fig10/fig10_swr_correction.pro

    r41 r43  
     1;+ 
     2; .. _fig10_swr_correction.pro 
     3; 
     4; ======================== 
     5; fig10_swr_correction.pro 
     6; ======================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> fig10_swr_correction 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
     34 
    135pro fig10_swr_correction 
    236@common 
    3 ;------------------------------------------------------------ 
    437reinitplt, /z,/invert 
    538key_portrait = 1 
    639 
    740openps, FILENAME = 'idl.ps' 
    8 ;------------------------------------------------------------ 
    941close,/all 
    1042 
    … …  
    4880significance_test,mean,std 
    4981 
    50 ;;------------------------------------------------------------------------ 
    5182closeps 
    5283 

  • trunk/src/paper01/fig11/swr_isccp_tropflux_correlation_fig10.pro

    r41 r43  
     1;+ 
     2; .. _swr_isccp_tropflux_correlation_fig10.pro: 
     3; 
     4; ======================================== 
     5; swr_isccp_tropflux_correlation_fig10.pro 
     6; ======================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> swr_isccp_tropflux_correlation_fig10 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro swr_isccp_tropflux_correlation_fig10 
    235@common 
    3 ;-------------------------------------------------------------------------- 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    638 
    739openps, FILENAME = 'idl.ps' 
    8 ;-------------------------------------------------------------------------- 
    940st=19890101 & en=20071231 
    1041 
    … …  
    2455     contour=sst-273.15, contmin=24, contmax=29, contint=1.5 
    2556 
    26 ;-------------------------------------------------------------------------- 
    2757closeps 
    2858fig='correlation_rmsd_olra_swra_smooth.ps' 
    … …  
    3262stop 
    3363end 
    34 ;-------------------------------------------------------------------------- 

  • trunk/src/paper01/fig11/swr_isccp_tropflux_new_v1.pro

    r41 r43  
     1;+ 
     2; 
     3; .. _swr_isccp_tropflux_new_v1.pro: 
     4; 
     5; ============================= 
     6; swr_isccp_tropflux_new_v1.pro 
     7; ============================= 
     8; 
     9; DESCRIPTION 
     10; =========== 
     11; 
     12; SEE ALSO 
     13; ======== 
     14; 
     15; :ref:`tropflux_profile.sh` 
     16; 
     17; EXAMPLES 
     18; ======== 
     19; 
     20; :: 
     21; 
     22;  IDL> @tropflux_init 
     23;  IDL> swr_isccp_tropflux_new_v1 
     24; 
     25; EVOLUTIONS 
     26; ========== 
     27; 
     28; $Id$ 
     29; 
     30; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     31; 
     32;   * minimal header 
     33; 
     34;- 
    135pro swr_isccp_tropflux_new_v1 
    236@common 
    3 ;-------------------------------------------------------------------------- 
    437reinitplt, /z,/invert 
    538key_portrait = 1 
    639 
    740openps, FILENAME = 'idl.ps' 
    8 ;-------------------------------------------------------------------------- 
    941st=19890101 & en=20071231 
    1042 
    … …  
    5082help, swi, olri 
    5183 
    52 ;-------------------------------------------------------------------------- 
    5384caldat, time,mon,day,yea 
    5485swm=olr*0. & olrm=olr*0. 
    … …  
    6192swa=sw_isccp-swm        ; swr anomaly 
    6293olra=olr-olrm     ; olr anomaly 
    63 ;-------------------------------------------------------------------------- 
    6494help, swa, olra 
    6595m=reform(sw_trop(*,*,0)) 
    … …  
    81111  endfor 
    82112endfor 
    83   
     113 
    84114cor=fltarr(jpi,jpj)+!values.f_nan & rms=fltarr(jpi,jpj)+!values.f_nan 
    85115for jj=1,jpj-2 do begin 
    … …  
    92122 
    93123save, cor, cor_olr, filename='correlation_isccp_olr.idl' 
    94   
     124 
    95125lct=64 & cs=0.9 
    96126marge1=[0,0,-5,0] 
    … …  
    118148     subt='',title='b) RMS Difference (OLRA, SWRA)',xminor=1,yminor=1 
    119149 
    120 ;-------------------------------------------------------------------------- 
    121150closeps 
    122151fig='correlation_rmsd_olra_swra_smooth.ps' 
    … …  
    127156 
    128157end 
    129 ;-------------------------------------------------------------------------- 

  • trunk/src/paper01/fig12/swr_statistics_map_2000_2009_v50.pro

    r41 r43  
     1;+ 
     2; .. _swr_statistics_map_2000_2009_v50.pro 
     3; 
     4; ==================================== 
     5; swr_statistics_map_2000_2009_v50.pro 
     6; ==================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> swr_statistics_map_2000_2009_v50 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro swr_statistics_map_2000_2009_v50 
    235@common 
    3 ;------------------------------------------------------------ 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242bias_mi=-20 & bias_ma=20 & bias_int=2 
    1343std_mi=0.7   & std_ma=1.31  & std_int=0.05 
    … …  
    90120for n=0,NN-1 do begin 
    91121  x=lon(n) 
    92   y=lat(n)  
    93   c=cor_trop(n)  
    94    
    95   plots, x,y,psym=8,symsize=1.5,color=0 
    96 endfor 
    97 ;------------------------------------------------------------ 
     122  y=lat(n) 
     123  c=cor_trop(n) 
     124 
     125  plots, x,y,psym=8,symsize=1.5,color=0 
     126endfor 
    98127plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $ 
    99128   title='2) SWR Correlation -  ERAI', subtitle='', small=[1,4,2],/noer,/rempl, marge=marge 
    … …  
    116145for n=0,NN-1 do begin 
    117146  x=lon(n) 
    118   y=lat(n)  
    119   c=cor_era(n)  
    120    
    121   plots, x,y,psym=8,symsize=1.5,color=0 
    122 endfor 
    123 ;------------------------------------------------------------ 
     147  y=lat(n) 
     148  c=cor_era(n) 
     149 
     150  plots, x,y,psym=8,symsize=1.5,color=0 
     151endfor 
    124152plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $ 
    125153   title='3) SWR Correlation - OAFlux', subtitle='', small=[1,4,3],/noer,/rempl, marge=marge 
    … …  
    144172  y=lat(n) 
    145173  c=cor_oaflx(n) 
    146   
    147   plots, x,y,psym=8,symsize=1.5,color=0 
    148 endfor 
    149 ;------------------------------------------------------------ 
     174 
     175  plots, x,y,psym=8,symsize=1.5,color=0 
     176endfor 
    150177plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int, $ 
    151178   title='4) SWR Correlation - SWR_olr', subtitle='', small=[1,4,4],/noer,/rempl, marge=marge1 
    … …  
    170197  y=lat(n) 
    171198  c=cor_olr(n) 
    172   
     199 
    173200  plots, x,y,psym=8,symsize=1.5,color=0 
    174201endfor 
    175202 
    176203erase 
    177 ;------------------------------------------------------------ 
    178204 
    179205plt, -5+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int, marge=marge,/nocolorb, $ 
    … …  
    197223for n=0,NN-1 do begin 
    198224  x=lon(n) 
    199   y=lat(n)  
    200   c=bias_trop(n)  
    201    
    202   plots, x,y,psym=8,symsize=1.5,color=0 
    203 endfor 
    204 ;------------------------------------------------------------ 
     225  y=lat(n) 
     226  c=bias_trop(n) 
     227 
     228  plots, x,y,psym=8,symsize=1.5,color=0 
     229endfor 
    205230plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int,/noer, marge=marge, $ 
    206231   title='2) SWR Mean bias - ERAI', subtitle='', small=[1,4,2],/rempl,format='(i3)',/nocolorb 
    … …  
    225250  y=lat(n) 
    226251  c=bias_era(n) 
    227    
    228   plots, x,y,psym=8,symsize=1.5,color=0 
    229 endfor 
    230  
    231 ;------------------------------------------------------------ 
     252 
     253  plots, x,y,psym=8,symsize=1.5,color=0 
     254endfor 
     255 
    232256 
    233257plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int,/noer, marge=marge, $ 
    … …  
    253277  y=lat(n) 
    254278  c=bias_oaflx(n) 
    255    
    256   plots, x,y,psym=8,symsize=1.5,color=0 
    257 endfor 
    258  
    259 ;------------------------------------------------------------ 
     279 
     280  plots, x,y,psym=8,symsize=1.5,color=0 
     281endfor 
     282 
    260283 
    261284plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int,/noer, marge=marge1, $ 
    … …  
    281304  y=lat(n) 
    282305  c=bias_olr(n) 
    283    
     306 
    284307  plots, x,y,psym=8,symsize=1.5,color=0 
    285308endfor 
    … …  
    287310erase 
    288311 
    289 ;------------------------------------------------------------ 
    290312 
    291313plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 
    … …  
    311333  y=lat(n) 
    312334  c=rmsd_trop(n) 
    313    
    314   plots, x,y,psym=8,symsize=1.5,color=0 
    315 endfor 
    316 ;------------------------------------------------------------ 
     335 
     336  plots, x,y,psym=8,symsize=1.5,color=0 
     337endfor 
    317338plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 
    318339   title='2) RMSD - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge 
    … …  
    337358  y=lat(n) 
    338359  c=rmsd_era(n) 
    339    
    340   plots, x,y,psym=8,symsize=1.5,color=0 
    341 endfor 
    342 ;------------------------------------------------------------ 
     360 
     361  plots, x,y,psym=8,symsize=1.5,color=0 
     362endfor 
    343363plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 
    344364   title='3) RMSD - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge 
    … …  
    363383  y=lat(n) 
    364384  c=rmsd_oaflx(n) 
    365    
    366   plots, x,y,psym=8,symsize=1.5,color=0 
    367 endfor 
    368  
    369 ;------------------------------------------------------------ 
     385 
     386  plots, x,y,psym=8,symsize=1.5,color=0 
     387endfor 
     388 
    370389plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 
    371390   title='4) RMSD - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1 
    … …  
    390409  y=lat(n) 
    391410  c=rmsd_olr(n) 
    392    
    393   plots, x,y,psym=8,symsize=1.5,color=0 
    394 endfor 
    395 ;------------------------------------------------------------ 
     411 
     412  plots, x,y,psym=8,symsize=1.5,color=0 
     413endfor 
    396414erase 
    397415plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 
    … …  
    420438  plots, x,y,psym=8,symsize=1.5,color=0 
    421439endfor 
    422 ;------------------------------------------------------------ 
    423440plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 
    424441   title='2) STD ratop - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge 
    … …  
    447464endfor 
    448465 
    449 ;------------------------------------------------------------ 
    450466plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 
    451467   title='3) STD rato - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge 
    … …  
    474490endfor 
    475491 
    476 ;------------------------------------------------------------ 
    477492plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 
    478493   title='4) STD rato - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1 
    … …  
    500515  plots, x,y,psym=8,symsize=1.5,color=0 
    501516endfor 
    502 ;------------------------------------------------------------ 
    503517closeps 
    504518fig='swr_statistics_map_2000_2009_v50.ps' 
    505519spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    506520spawn, 'gv '+cpsdir+fig 
    507 ;------------------------------------------------------------ 
    508521return 
    509522end 
    510  

  • trunk/src/paper01/fig13/net_flux_validation_scatter_2000_2007.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro net_flux_validation_scatter_2000_2007,date1,date2 
     1;+ 
     2; .. _net_flux_validation_scatter_2000_2007.pro 
     3; 
     4; ========================================= 
     5; net_flux_validation_scatter_2000_2007.pro 
     6; ========================================= 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> net_flux_validation_scatter_2000_2007, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36pro net_flux_validation_scatter_2000_2007, date1, date2 
    337@common 
    4 ;------------------------------------------------------------ 
    538 
    639reinitplt, /z,/invert 
    … …  
    841 
    942openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1143 
    1244sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    … …  
    2759 
    2860ocean='net_flux' 
    29 nsmooth=1.  
    30 ;------------------------------------------------------------------------------------------------------------------------ 
     61nsmooth=1. 
    3162;;   This program will create the following text files with statistics of respective variables 
    32 ;------------------------------------------------------------------------------------------------------------------------ 
    3363file_trop='/Users/pkb/data/TropFlux/TropFlux_19890101_20091231.nc' 
    3464initncdf, file_trop 
    3565shf_trop=read_ncdf("shf", date1, date2, file=file_trop,/nostr) 
    3666lhf_trop=read_ncdf("lhf", date1, date2, file=file_trop,/nostr) 
    37 lwr_trop=read_ncdf("lwr", date1, date2, file=file_trop,/nostr)  
     67lwr_trop=read_ncdf("lwr", date1, date2, file=file_trop,/nostr) 
    3868 
    3969file_trop='/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20071231_DT_v50.nc' 
    … …  
    71101file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/shf_oafluxgrid_1985_2008.nc' 
    72102initncdf, file_oaf 
    73 shf=read_ncdf("shf", date1, date2, file=file_oaf,/nostr)  
     103shf=read_ncdf("shf", date1, date2, file=file_oaf,/nostr) 
    74104file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/lhf_oafluxgrid_1985_2008.nc' 
    75105initncdf, file_oaf 
    76 lhf=read_ncdf("lhf", date1, date2, file=file_oaf,/nostr)  
     106lhf=read_ncdf("lhf", date1, date2, file=file_oaf,/nostr) 
    77107file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/lwr_oafluxgrid_1985_2007.nc' 
    78108initncdf, file_oaf 
    79 lwr=read_ncdf("lwr", date1, date2, file=file_oaf,/nostr)  
     109lwr=read_ncdf("lwr", date1, date2, file=file_oaf,/nostr) 
    80110file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/swr_oafluxgrid_1985_2007.nc' 
    81111initncdf, file_oaf 
    … …  
    94124help, ncep1 
    95125 
    96 ;------------------------------------------------------------------------------------------------------------------------ 
    97126close,/all 
    98127 
    … …  
    118147printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    119148 
    120 ;------------------------------------------------------------------------------------------------------------------------ 
    121149nn=n_elements(sitelist) 
    122150 
    … …  
    133161    if (y ge 0. and y le 30.) then y=y+360. 
    134162    dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx] 
    135      
     163 
    136164    read_variables_net_flux, csite,date1,date2,nsmooth, $ 
    137165         net 
    138166    tao=net & ind=where(finite(tao)) & valid=n_elements(ind) 
    139167 
    140 ;;     
     168;; 
    141169    if (valid ge 180. ) then begin 
    142170 
    … …  
    155183           extract_flux_tropflux, trop, box, $ 
    156184              tropflux 
    157            tropflux_new=tropflux  
     185           tropflux_new=tropflux 
    158186 
    159187           extract_flux_tropflux, ncep1, box, $ 
    … …  
    165193           tropdf=tropflux 
    166194 
    167            ind=where(finite(tao)) & tao=tao(ind) & oaf_net=oaflux(ind) & nce_net=ncepflux(ind) & trop_net=tropflux_new(ind)  
     195           ind=where(finite(tao)) & tao=tao(ind) & oaf_net=oaflux(ind) & nce_net=ncepflux(ind) & trop_net=tropflux_new(ind) 
    168196           erai_net=eraiflux(ind) & nce1_net=ncepflux1(ind) & tropd_net=tropdf(ind) 
    169197 
    … …  
    383411oplot, [-100,200], [-100,200] 
    384412 
    385 ;----------------------------------------------------------- 
    386413closeps 
    387414 
    … …  
    391418return 
    392419end 
    393 ;-------------------------------------------------------------------------- 
    394420function x_site_location, site 
    395421    n1=strpos(site, 's') 
    … …  
    406432return, float(x) 
    407433end 
    408 ;-------------------------------------------------------------------------- 
    409434function y_site_location, site 
    410435    n1=strpos(site, 'e') 
    … …  
    432457return,float(y) 
    433458end 
    434  
    435 ;-------------------------------------------------------------------------- 
    436  

  • trunk/src/paper01/fig14/as_validation_net_1994_95_v10.pro

    r41 r43  
     1;+ 
     2; .. _as_validation_net_1994_95_v10.pro 
     3; 
     4; ================================= 
     5; as_validation_net_1994_95_v10.pro 
     6; ================================= 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> as_validation_net_1994_95_v10 
     23; 
     24; 
     25; EVOLUTIONS 
     26; ========== 
     27; 
     28; $Id$ 
     29; 
     30; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     31; 
     32;   * minimal header 
     33; 
     34;- 
    135pro as_validation_net_1994_95_v10 
    236@common 
    3 ;------------------------------------------------------------ 
    437reinitplt, /z,/invert 
    538key_portrait = 1 
    … …  
    740 
    841openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1042; partie a changer 
    11 ;------------------------------------------------------------ 
    1243asbox=[61,62,15,16] 
    1344date1=19941016 & date2=19951019 
    … …  
    71102file="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20091231_NRT_v50.nc" 
    72103initncdf, file 
    73 olr1=read_ncdf("sw", date1-1, date2, file=file,/nostr, box=asbox) ;; & olr1=grossemoyenne(olr1, "xy",/nan, box=asbox) &  
     104olr1=read_ncdf("sw", date1-1, date2, file=file,/nostr, box=asbox) ;; & olr1=grossemoyenne(olr1, "xy",/nan, box=asbox) & 
    74105olr1=reform(olr1) 
    75106 
    … …  
    140171nce1=nce1(ind) 
    141172help, nce1 
    142 ;------------------------------------------------------------ 
    143173;ind=where(finite(net)) 
    144174;net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce1=nce1(ind) & nce2=nce2(ind) 
    … …  
    175205cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    176206 
    177 ;------------------------------------------------------------ 
    178207splot, net, era, charsize=1.1, title='Net Flux - AS Vs ERAI', $ 
    179208     xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 
    … …  
    207236 
    208237stop 
    209 ;------------------------------------------------------------ 
    210238closeps 
    211239fig='as_validation_net_1994_95_v10.ps' 
    212240spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    213241spawn, 'gv '+cpsdir+fig 
    214 ;------------------------------------------------------------ 
    215242return 
    216243end 
    217  

  • trunk/src/paper01/fig14/coare_validation_net_1992_93_all_v10.pro

    r41 r43  
     1;+ 
     2; .. _coare_validation_net_1992_93_all_v10.pro: 
     3; 
     4; ======================================== 
     5; coare_validation_net_1992_93_all_v10.pro 
     6; ======================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> coare_validation_net_1992_93_all_v10 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro coare_validation_net_1992_93_all_v10 
    235@common 
    3 ;------------------------------------------------------------ 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242 
    1343box=[155.5, 156.5, -2, -1] 
    … …  
    143173shf4=read_ncdf('shf',st4,en4,file=file,/nostr, box=box) & shf4=grossemoyenne(shf4, "xy",/nan, box=box) & shf4=reform(shf4) 
    144174 
    145 net1=swr1+lwr1+lhf1+shf1 & net2=swr2+lwr2+lhf2+shf2  
     175net1=swr1+lwr1+lhf1+shf1 & net2=swr2+lwr2+lhf2+shf2 
    146176net3=swr3+lwr3+lhf3+shf3 & net4=swr4+lwr4+lhf4+shf4 
    147177 
    … …  
    278308 
    279309help, net, era, trop, oaf, nce1, nce2, olr 
    280 ;------------------------------------------------------------ 
    281310 
    282311statistics, net, era, $ 
    … …  
    310339cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    311340 
    312 ;------------------------------------------------------------ 
    313341 
    314342splot, net, era, charsize=1.1, title='Net Flux - COARE Vs ERAI', $ 
    … …  
    343371 
    344372stop 
    345 ;------------------------------------------------------------ 
    346373closeps 
    347374fig="coare_validation_net_1992_93_all_v11.ps" 
    348375spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    349376spawn, 'gv '+cpsdir+fig 
    350 ;------------------------------------------------------------ 
    351377return 
    352378end 
    353  

  • trunk/src/paper01/fig14/stratus_validation_net_2001_2002_all_v10.pro

    r41 r43  
     1;+ 
     2; .. _stratus_validation_net_2001_2002_all_v10.pro 
     3; 
     4; ============================================ 
     5; stratus_validation_net_2001_2002_all_v10.pro 
     6; ============================================ 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> stratus_validation_net_2001_2002_all_v10 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro stratus_validation_net_2001_2002_all_v10 
    235@common 
    3 ;------------------------------------------------------------ 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242 
    1343;; stratus in the south pacific cold tongue 2000-01 
    … …  
    122152net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce2=nce2(ind) & nce1=nce1(ind) 
    123153olr=olr(ind) 
    124 ;------------------------------------------------------------ 
    125154 
    126155statistics, net, era, $ 
    … …  
    154183cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    155184 
    156 ;------------------------------------------------------------ 
    157185 
    158186splot, net, era, charsize=1.1, title='SWR - Stratus Vs ERAI', $ 
    … …  
    187215 
    188216stop 
    189 ;------------------------------------------------------------ 
    190217closeps 
    191218fig="stratus_validation_net_2001_2002_all_v11.ps" 
    192219spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    193220spawn, 'gv '+cpsdir+fig 
    194 ;------------------------------------------------------------ 
    195221return 
    196222end 
    197  

  • trunk/src/paper01/fig14/subdctn_cs_validation_net_1991_93_all_v10.pro

    r41 r43  
     1;+ 
     2; .. _subdctn_cs_validation_net_1991_93_all_v10.pro 
     3; 
     4; ============================================= 
     5; subdctn_cs_validation_net_1991_93_all_v10.pro 
     6; ============================================= 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> subdctn_cs_validation_net_1991_93_all_v10 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro subdctn_cs_validation_net_1991_93_all_v10 
    235@common 
    3 ;------------------------------------------------------------ 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242 
    1343;; subduction experiment in the atlantic in 1991-93 
    … …  
    136166net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce2=nce2(ind) & nce1=nce1(ind) 
    137167olr=olr(ind) 
    138 ;------------------------------------------------------------ 
    139168 
    140169statistics, net, era, $ 
    … …  
    168197cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    169198 
    170 ;------------------------------------------------------------ 
    171199 
    172200splot, net, era, charsize=1.1, title='Net Flux - Subdn. cs Vs ERAI', $ 
    … …  
    201229 
    202230stop 
    203 ;------------------------------------------------------------ 
    204231closeps 
    205232fig="subdctn_cs_validation_net_1991_93_all_v10.ps" 
    206233spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    207234spawn, 'gv '+cpsdir+fig 
    208 ;------------------------------------------------------------ 
    209235return 
    210236end 
    211  

  • trunk/src/paper01/fig14/subdctn_ses_validation_net_1991_93_all_v10.pro

    r41 r43  
     1;+ 
     2; .. _subdctn_ses_validation_net_1991_93_all_v10.pro 
     3; 
     4; ============================================== 
     5; subdctn_ses_validation_net_1991_93_all_v10.pro 
     6; ============================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> subdctn_ses_validation_net_1991_93_all_v10 
     23; 
     24; 
     25; EVOLUTIONS 
     26; ========== 
     27; 
     28; $Id$ 
     29; 
     30; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     31; 
     32;   * minimal header 
     33; 
     34;- 
    135pro subdctn_ses_validation_net_1991_93_all_v10 
    236@common 
    3 ;------------------------------------------------------------ 
    437reinitplt, /z,/invert 
    538key_portrait = 1 
    … …  
    740 
    841openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1042; partie a changer 
    11 ;------------------------------------------------------------ 
    1243 
    1344;; subduction experiment in the atlantic in 1991-93 
    … …  
    93124 
    94125net1=swr1+lwr1+(shf1+lhf1) & net2=swr2+lwr2+(shf2+lhf2) 
    95 neto_1=olr1+lwr1+lhf1+shf1 & neto_2=olr2+lwr2+lhf2+shf2  
     126neto_1=olr1+lwr1+lhf1+shf1 & neto_2=olr2+lwr2+lhf2+shf2 
    96127trop=[net1,net2] & olr=[neto_1,neto_2] 
    97128help, trop, olr 
    … …  
    163194net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce2=nce2(ind) & nce1=nce1(ind) 
    164195olr=olr(ind) 
    165 ;------------------------------------------------------------ 
    166196 
    167197statistics, net, era, $ 
    … …  
    195225cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    196226 
    197 ;------------------------------------------------------------ 
    198227 
    199228splot, net, era, charsize=1.1, title='Net Flux - Subdn. ses Vs ERAI', $ 
    … …  
    228257 
    229258stop 
    230 ;------------------------------------------------------------ 
    231259closeps 
    232260fig="subdctn_ses_validation_net_1991_93_all_v10.ps" 
    233261spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    234262spawn, 'gv '+cpsdir+fig 
    235 ;------------------------------------------------------------ 
    236263return 
    237264end 
    238  

  • trunk/src/paper01/fig14/subdctn_sws_validation_net_1991_93_all_v10.pro

    r41 r43  
     1;+ 
     2; .. _subdctn_sws_validation_net_1991_93_all_v10.pro 
     3; 
     4; ============================================== 
     5; subdctn_sws_validation_net_1991_93_all_v10.pro 
     6; ============================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> subdctn_sws_validation_net_1991_93_all_v10 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro subdctn_sws_validation_net_1991_93_all_v10 
    235@common 
    3 ;------------------------------------------------------------ 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242 
    1343;; subduction experiment in the atlantic in 1991-93 
    … …  
    90120era=[net1,net2,net3,net4] 
    91121lhf1=0. & lhf2=0. & lhf3=0. & lhf4=0. 
    92 shf1=0. & shf2=0. & shf3=0. & shf4=0.  
    93 lwr1=0. & lwr2=0. & lwr3=0. & lwr4=0. 
    94 swr1=0. & swr2=0. & swr3=0. & swr4=0.  
     122shf1=0. & shf2=0. & shf3=0. & shf4=0. 
     123lwr1=0. & lwr2=0. & lwr3=0. & lwr4=0. 
     124swr1=0. & swr2=0. & swr3=0. & swr4=0. 
    95125net1=0. & net2=0. & net3=0. & net4=0 
    96126help, era 
    … …  
    132162lhf1=0. & lhf2=0. & lhf3=0. & lhf4=0. 
    133163shf1=0. & shf2=0. & shf3=0. & shf4=0. 
    134 lwr1=0. & lwr2=0. & lwr3=0. & lwr4=0.  
     164lwr1=0. & lwr2=0. & lwr3=0. & lwr4=0. 
    135165swr1=0. & swr2=0. & swr3=0. & swr4=0. 
    136166net1=0. & net2=0. & net3=0. & net4=0. 
    … …  
    252282net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce2=nce2(ind) & nce1=nce1(ind) 
    253283olr=olr(ind) 
    254 ;------------------------------------------------------------ 
    255284 
    256285statistics, net, era, $ 
    … …  
    284313cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    285314 
    286 ;------------------------------------------------------------ 
    287315 
    288316splot, net, era, charsize=1.1, title='SWR - Subdctn. sws Vs ERAI', $ 
    … …  
    317345 
    318346stop 
    319 ;------------------------------------------------------------ 
    320347closeps 
    321348fig="subdctn_sws_validation_net_1991_93_all_v10.ps" 
    322349spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    323350spawn, 'gv '+cpsdir+fig 
    324 ;------------------------------------------------------------ 
    325351return 
    326352end 
    327  

  • trunk/src/paper01/fig14/whots_validation_net_2004_2007_all_v10.pro

    r41 r43  
     1;+ 
     2; .. _whots_validation_net_2004_2007_all_v10.pro 
     3; 
     4; ========================================== 
     5; whots_validation_net_2004_2007_all_v10.pro 
     6; ========================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> whots_validation_net_2004_2007_all_v10 
     23; 
     24;  EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro whots_validation_net_2004_2007_all_v10 
    235@common 
    3 ;------------------------------------------------------------ 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242 
    1343;; WHOTS in the NC pacific  2004-09 
    … …  
    5080;; reading other data sets 
    5181 
    52 st1=20040813 & en1=20050725  
     82st1=20040813 & en1=20050725 
    5383st2=20050728 & en2=20060624 
    54 st3=20060627 & en3=20070628  
    55 ;st4=20070626 & en4=20080606  
     84st3=20060627 & en3=20070628 
     85;st4=20070626 & en4=20080606 
    5686;st5=20080605 & en5=20090715 
    5787 
    … …  
    84114lhf1=0. & lhf2=0. & lhf3=0. 
    85115shf1=0. & shf2=0. & shf3=0. 
    86 lwr1=0. & lwr2=0. & lwr3=0.  
    87 swr1=0. & swr2=0. & swr3=0. 
    88 net1=0. & net2=0. & net3=0.  
     116lwr1=0. & lwr2=0. & lwr3=0. 
     117swr1=0. & swr2=0. & swr3=0. 
     118net1=0. & net2=0. & net3=0. 
    89119help, era 
    90120 
    … …  
    115145 
    116146net1=swr1+lwr1+(shf1+lhf1) & net2=swr2+lwr2+(shf2+lhf2) & net3=swr3+lwr3+(shf3+lhf3) 
    117 neto_1=olr1+lwr1+lhf1+shf1 & neto_2=olr2+lwr2+lhf2+shf2 & neto_3=olr3+lwr3+lhf3+shf3  
     147neto_1=olr1+lwr1+lhf1+shf1 & neto_2=olr2+lwr2+lhf2+shf2 & neto_3=olr3+lwr3+lhf3+shf3 
    118148 
    119149trop=[net1,net2,net3] & olr=[neto_1,neto_2,neto_3] 
    120150 
    121151lhf1=0. & lhf2=0. & lhf3=0. 
    122 shf1=0. & shf2=0. & shf3=0.  
    123 lwr1=0. & lwr2=0. & lwr3=0. 
    124 swr1=0. & swr2=0. & swr3=0.  
    125 net1=0. & net2=0. & net3=0.  
     152shf1=0. & shf2=0. & shf3=0. 
     153lwr1=0. & lwr2=0. & lwr3=0. 
     154swr1=0. & swr2=0. & swr3=0. 
     155net1=0. & net2=0. & net3=0. 
    126156neto_1=0 & neto_2=0 & neto_3=0 
    127157 
    … …  
    155185 
    156186oaf=[net1,net2,net3] 
    157 lhf1=0. & lhf2=0. & lhf3=0.  
    158 shf1=0. & shf2=0. & shf3=0.  
     187lhf1=0. & lhf2=0. & lhf3=0. 
     188shf1=0. & shf2=0. & shf3=0. 
    159189lwr1=0. & lwr2=0. & lwr3=0. 
    160190swr1=0. & swr2=0. & swr3=0. 
    … …  
    191221help, nce2 
    192222 
    193 lhf1=0. & lhf2=0. & lhf3=0.  
    194 shf1=0. & shf2=0. & shf3=0. 
    195 lwr1=0. & lwr2=0. & lwr3=0. 
    196 swr1=0. & swr2=0. & swr3=0.  
     223lhf1=0. & lhf2=0. & lhf3=0. 
     224shf1=0. & shf2=0. & shf3=0. 
     225lwr1=0. & lwr2=0. & lwr3=0. 
     226swr1=0. & swr2=0. & swr3=0. 
    197227net1=0. & net2=0. & net3=0. 
    198228 
    … …  
    220250help, nce1 
    221251 
    222 lhf1=0. & lhf2=0. & lhf3=0.  
    223 shf1=0. & shf2=0. & shf3=0. 
    224 lwr1=0. & lwr2=0. & lwr3=0. 
    225 swr1=0. & swr2=0. & swr3=0.  
     252lhf1=0. & lhf2=0. & lhf3=0. 
     253shf1=0. & shf2=0. & shf3=0. 
     254lwr1=0. & lwr2=0. & lwr3=0. 
     255swr1=0. & swr2=0. & swr3=0. 
    226256net1=0. & net2=0. & net3=0. 
    227257 
    228258ind=where(finite(net)) 
    229 net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind)  
     259net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) 
    230260nce2=nce2(ind) & nce1=nce1(ind) & olr=olr(ind) 
    231261 
    232 ;------------------------------------------------------------ 
    233262 
    234263statistics, net, era, $ 
    … …  
    262291cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    263292 
    264 ;------------------------------------------------------------ 
    265293 
    266294splot, net, era, charsize=1.1, title='Net Flux - WHOTS Vs ERAI', $ 
    … …  
    294322xyouts, 0,-200, cstat_olr, charsize=1. 
    295323 
    296 ;------------------------------------------------------------ 
    297324closeps 
    298325fig="WHOTS_validation_net_2004_2007_all_v10.ps" 
    299326spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    300327spawn, 'gv '+cpsdir+fig 
    301 ;------------------------------------------------------------ 
    302328return 
    303329end 
    304  

  • trunk/src/paper01/fig15/fig15_timeseries_mjo.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _fig15_timeseries_mjo.pro: 
     3; 
     4; ======================== 
     5; fig15_timeseries_mjo.pro 
     6; ======================== 
     7; 
     8; 
     9; DESCRIPTION 
     10; =========== 
     11; 
     12; SEE ALSO 
     13; ======== 
     14; 
     15; :ref:`tropflux_profile.sh` 
     16; 
     17; EXAMPLES 
     18; ======== 
     19; 
     20; :: 
     21; 
     22;  IDL> @tropflux_init 
     23;  IDL> fig15_timeseries_mjo 
     24; 
     25; EVOLUTIONS 
     26; ========== 
     27; 
     28; $Id$ 
     29; 
     30; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     31; 
     32;   * minimal header 
     33; 
     34;- 
    235pro fig15_timeseries_mjo 
    336@common 
    4 ;------------------------------------------------------------ 
    537 
    638reinitplt, /z,/invert 
    … …  
    840 
    941openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1142date1=20071001 & date2=20080331 
    1243date1=19990101 & date2=20021231 
    … …  
    79110 
    80111save, olr,trop, nce2, era, nce, filename="netflux_sctr_1999_2002.idl" 
    81 ;------------------------------------------------------------------------------------------------------------------------ 
    82112 
    83113trop_hf=trop*0. & nce2_hf=trop*0. 
    … …  
    130160pltt, nce_hf, "t",/ov1d, color=50 
    131161 
    132 ;---------------------------------------------------------- 
    133162closeps 
    134163fig='fig15_timeseries_mjo.ps' 
    … …  
    137166return 
    138167end 
    139 ;-------------------------------------------------------------------------- 
     168 
    140169function x_site_location, site 
    141170    n1=strpos(site, 's') 
    … …  
    152181return, float(x) 
    153182end 
    154 ;-------------------------------------------------------------------------- 
     183 
    155184function y_site_location, site 
    156185    n1=strpos(site, 'e') 
    … …  
    178207return,float(y) 
    179208end 
    180  
    181 ;-------------------------------------------------------------------------- 
    182  

  • trunk/src/paper01/fig16/fig16_timeseries_nino3.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _fig16_timeseries_nino3.pro: 
     3; 
     4; ========================== 
     5; fig16_timeseries_nino3.pro 
     6; ========================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> fig16_timeseries_nino3 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    234pro fig16_timeseries_nino3 
    335@common 
    4 ;------------------------------------------------------------ 
    536 
    637reinitplt, /z,/invert 
    … …  
    839 
    940openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1141date1=19940101 & date2=20071231 
    1242box=[240,330,-5,5] 
    … …  
    1848restore, filename="/Users/pkb/work/MY_SAXO/netflux_nino3_1989_2007.idl" 
    1949help, trop, oafl, nce2, era, nce, sst 
    20 ;------------------------------------------------------------------------------------------------------------------------ 
    2150;; creating climatology 
    2251ny=(2007-1994)+1 
    … …  
    99128trop_ano=trop_ano-trop_hf & oafl_ano=oafl_ano-oafl_hf 
    100129nce2_ano=nce2_ano-nce2_hf & era_ano=era_ano-era_hf 
    101 nce_ano=nce_ano-nce_hf     
     130nce_ano=nce_ano-nce_hf 
    102131 
    103132nsmooth=30. 
    … …  
    127156pltt, ts_smooth(nce_ano,nsmooth,/nan), "t",/ov1d, color=65 
    128157 
    129 ;---------------------------------------------------------- 
    130158closeps 
    131159fig='fig16_timeseries_nino3.ps' 
    … …  
    134162return 
    135163end 
    136 ;-------------------------------------------------------------------------- 
    137164function x_site_location, site 
    138165    n1=strpos(site, 's') 
    … …  
    149176return, float(x) 
    150177end 
    151 ;-------------------------------------------------------------------------- 
    152178function y_site_location, site 
    153179    n1=strpos(site, 'e') 
    … …  
    175201return,float(y) 
    176202end 
    177  
    178 ;-------------------------------------------------------------------------- 
    179  

  • trunk/src/paper01/fig2/fig2_timeline_diagram.pro

    r41 r43  
     1;+ 
     2; .. _fig2_timeline_diagram.pro: 
     3; 
     4; ========================= 
     5; fig2_timeline_diagram.pro 
     6; ========================= 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> fig2_timeline_diagram 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro fig2_timeline_diagram 
    235@common 
    … …  
    115148    lwr=lw & ind=where(finite(lwr)) & valid=n_elements(ind) 
    116149    if (valid ge 10) then begin 
    117      
     150 
    118151       ind=where(finite(lw)) &  lw(ind)=1. 
    119152       ind=where(finite(lw,/nan)) & lw(ind)=0. 
    … …  
    180213return,float(y) 
    181214end 
    182  
    183  
    184  

  • trunk/src/paper01/fig2/read_lh.pro

    r41 r43  
    1 ;------------------------------------------------------------------------------------------------- 
    2 pro read_lh, csite,date1,date2,nsmooth, $ 
    3                       lh 
    4  
     1;+ 
     2; .. _read_lh.pro 
    53; 
    6 ;------------------------------------------------------------------------------------------------- 
     4; =========== 
     5; read_lh.pro 
     6; =========== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> csite=++ 
     23;  IDL> ++ 
     24;  IDL> read_lh, csite, date1, date2, nsmooth, lh 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36pro read_lh, csite, date1, date2, nsmooth, lh 
     37; 
    738@common 
    839 
    … …  
    4980;; Replace missing values by "NaN" 
    5081; 
    51 tsvars=['lh']  
     82tsvars=['lh'] 
    5283vars=[tsvars] 
    5384nn=n_elements(vars) 
    … …  
    5990 
    6091nsmooth=nsmooth 
    61 lh=smooth(lh,nsmooth,/nan)  
     92lh=smooth(lh,nsmooth,/nan) 
    6293 
    63 ;------------------------------------------------------------------------------------------------- 
    6494end 
    65 ;------------------------------------------------------------------------------------------------- 
    6695 
    6796 
    68 ;------------------------------------------------------------------------------------------------- 
    6997function time_lec, fi 
    7098tt=ncdf_lec(fi,var='time') 
    … …  
    78106return, tt 
    79107end 
    80  
    81 ;------------------------------------------------------------------------------------------------- 

  • trunk/src/paper01/fig2/read_lw.pro

    r41 r43  
    1 ;------------------------------------------------------------------------------------------------- 
    2 pro read_lw, csite,date1,date2,box, $ 
    3                         lw 
     1;+ 
     2; .. _read_lw.pro 
    43; 
    5 ;------------------------------------------------------------------------------------------------- 
     4; =========== 
     5; read_lw.pro 
     6; =========== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> csite=++ 
     23;  IDL> ++ 
     24;  IDL> read_lw, csite, date1, date2, box, lw 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36pro read_lw, csite, date1, date2, box, lw 
     37; 
    638@common 
    739 
    … …  
    4072; 
    4173;tsvars=['at','bf','bp','dyn','emp','evap','heat','iso','lw','lwnet','qlat','qsen','sw', $ 
    42 ;        'swnet','rain','rf','rh','ssd','sss','sst','swnet','tx','ty','wu','wv','ws','lat','lon']  
     74;        'swnet','rain','rf','rh','ssd','sss','sst','swnet','tx','ty','wu','wv','ws','lat','lon'] 
    4375;depvars=['u','v','d','t','s'] 
    4476vars=['lw'] 
    … …  
    5486lw=smooth(lw,nsmooth,/nan) 
    5587 
    56 ;------------------------------------------------------------------------------------------------- 
    5788end 
    58 ;------------------------------------------------------------------------------------------------- 
    5989 
    6090 
    61 ;------------------------------------------------------------------------------------------------- 
    6291function time_lec, fi 
    6392tt=ncdf_lec(fi,var='time') 
    … …  
    71100return, tt 
    72101end 
    73 ;------------------------------------------------------------------------------------------------- 

  • trunk/src/paper01/fig2/read_sw.pro

    r41 r43  
    1 ;------------------------------------------------------------------------------------------------- 
    2 pro read_sw, csite,date1,date2,box, $ 
    3                         sw 
     1;+ 
     2; .. _read_sw.pro 
    43; 
    5 ;------------------------------------------------------------------------------------------------- 
     4; =========== 
     5; read_sw.pro 
     6; =========== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> csite=++ 
     23;  IDL> ++ 
     24;  IDL> read_sw, csite, date1, date2, box, sw 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36pro read_sw, csite, date1, date2, box, sw 
     37; 
    638@common 
    739 
    … …  
    4880; 
    4981;tsvars=['at','bf','bp','dyn','emp','evap','heat','iso','sw','swnet','qlat','qsen','sw', $ 
    50 ;        'swnet','rain','rf','rh','ssd','sss','sst','swnet','tx','ty','wu','wv','ws','lat','lon']  
     82;        'swnet','rain','rf','rh','ssd','sss','sst','swnet','tx','ty','wu','wv','ws','lat','lon'] 
    5183;depvars=['u','v','d','t','s'] 
    5284vars=['sw'] 
    … …  
    6395sw=smooth(sw,nsmooth,/nan) 
    6496 
    65 ;------------------------------------------------------------------------------------------------- 
    6697 
    6798 
    68 ;------------------------------------------------------------------------------------------------- 
    6999function time_lec, fi 
    70100tt=ncdf_lec(fi,var='time') 
    … …  
    78108return, tt 
    79109end 
    80 ;------------------------------------------------------------------------------------------------- 

  • trunk/src/paper01/fig3/air_validation_scatter_2000_2009_v50.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro air_validation_scatter_2000_2009_v50,date1,date2 
     1;+ 
     2; _air_validation_scatter_2000_2009_v50.pro: 
     3; 
     4; ======================================== 
     5; air_validation_scatter_2000_2009_v50.pro 
     6; ======================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> air_validation_scatter_2000_2009_v50, date1, date2 
     25; 
     26; 
     27; EVOLUTIONS 
     28; ========== 
     29; 
     30; $Id$ 
     31; 
     32; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     33; 
     34;   * minimal header 
     35; 
     36;- 
     37 
     38pro air_validation_scatter_2000_2009_v50, date1, date2 
    339@common 
    4 ;------------------------------------------------------------ 
    540 
    641reinitplt, /z,/invert 
    … …  
    843 
    944openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1145;; Before running this program, you have to compile the following subrutines 
    1246;;   1.  .r read_era_total 
    … …  
    3165 
    3266nsmooth=1.    ;; statistics are with 7 day smoothed 
    33 ;------------------------------------------------------------------------------------------------------------------------ 
    3467;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3668close,/all 
    3769 
    … …  
    5486printf,5, 'x     y      cor    bias     std     rmsd    mean_tao' 
    5587 
    56 ;------------------------------------------------------------------------------------------------------------------------ 
    5788;; first reading the whole ERAI uncorrected and corrected data 
    5889 
    5990file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc' 
    6091initncdf, file 
    61 unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15  
     92unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 
    6293help, unc 
    6394 
    … …  
    75106fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 
    76107initncdf, fi 
    77 nce=read_ncdf("air", date1-1, date2, file=fi,/nostr)  
     108nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 
    78109nce=nce-273.15 
    79110help, nce 
    … …  
    81112file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/t2m_ncep1_19890101_20091231.nc' 
    82113initncdf, file 
    83 nce1=read_ncdf("t2m", date1, date2, file=file,/nostr)  
     114nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 
    84115nce1=nce1-273.15 
    85116help, nce1 
    86117 
    87 ;------------------------------------------------------------------------------------------------------------------------ 
    88118nn=n_elements(sitelist) 
    89119date1=date1 
    … …  
    107137                                                ;;        lh  ->  latent heat flux 
    108138                                                ;;        rh  ->  relative humidity 
    109                                                 ;;   wu,wv,ws ->  wind speed     
    110 ;;     
     139                                                ;;   wu,wv,ws ->  wind speed 
     140;; 
    111141 
    112142 
    … …  
    118148    uncr=tropflux 
    119149 
    120      extract_flux_tropflux,cor,box, $   
     150     extract_flux_tropflux,cor,box, $ 
    121151        tropflux 
    122152     corr=tropflux 
    123      
     153 
    124154     extract_flux_tropflux,oaf,box, $ 
    125155        tropflux 
    … …  
    135165 
    136166 
    137     ind=where(finite(air)) & air=air(ind) & uncr_air=uncr(ind)  & corr_air=corr(ind)  
     167    ind=where(finite(air)) & air=air(ind) & uncr_air=uncr(ind)  & corr_air=corr(ind) 
    138168    oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 
    139169    mean_tao=total(air,/nan)/n_elements(ind) 
    140       
     170 
    141171    statistics_3var_v1, air, uncr_air, corr_air, $ 
    142172         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    143   
     173 
    144174    printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    145175    printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    … …  
    155185 
    156186endfor 
    157 ;----------------------------------------------------------------------------------------- 
    158 close,/all   
     187close,/all 
    159188 
    160189fi_air_erai='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_v50.txt' 
    … …  
    294323oplot, mean_tao, yfit, color=250, thick=2 
    295324 
    296 ;----------------------------------------------------------- 
    297325closeps 
    298326 
    … …  
    302330return 
    303331end 
    304 ;-------------------------------------------------------------------------- 
    305332function x_site_location, site 
    306333    n1=strpos(site, 's') 
    … …  
    317344return, float(x) 
    318345end 
    319 ;-------------------------------------------------------------------------- 
    320346function y_site_location, site 
    321347    n1=strpos(site, 'e') 
    … …  
    343369return,float(y) 
    344370end 
    345  
    346 ;-------------------------------------------------------------------------- 
    347  

  • trunk/src/paper01/fig3/q2m_validation_scatter_2000_2009_v50.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro q2m_validation_scatter_2000_2009_v50,date1,date2 
     1;+ 
     2; .. _q2m_validation_scatter_2000_2009_v50.pro 
     3; 
     4; ======================================== 
     5; q2m_validation_scatter_2000_2009_v50.pro 
     6; ======================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> q2m_validation_scatter_2000_2009_v50, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36pro q2m_validation_scatter_2000_2009_v50, date1, date2 
    337@common 
    4 ;------------------------------------------------------------ 
    538 
    639reinitplt, /z,/invert 
    … …  
    841 
    942openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1143;; Before running this program, you have to compile the following subrutines 
    1244;;   1.  .r read_era_total 
    … …  
    3163 
    3264nsmooth=1.    ;; statistics are with 7 day smoothed 
    33 ;------------------------------------------------------------------------------------------------------------------------ 
    3465;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3666close,/all 
    3767 
    … …  
    5484printf,5, 'x     y      cor    bias     std     rmsd    mean_tao' 
    5585 
    56 ;------------------------------------------------------------------------------------------------------------------------ 
    5786;; first reading the whole ERAI uncorrected and corrected data 
    5887 
    … …  
    6493file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc' 
    6594initncdf, file 
    66 cor=read_ncdf('q2m',date1,date2,file=file,/nostr)  
     95cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 
    6796help, cor 
    6897 
    … …  
    74103fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sphum_ncep2_oafluxgrid_19890101_20091231.nc' 
    75104initncdf, fi 
    76 nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr)  
     105nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 
    77106nce=nce*1000 
    78107help, nce 
    … …  
    84113help, nce1 
    85114 
    86 ;------------------------------------------------------------------------------------------------------------------------ 
    87115nn=n_elements(sitelist) 
    88116date1=date1 
    … …  
    106134                                                ;;        lh  ->  latent heat flux 
    107135                                                ;;        rh  ->  relative humidity 
    108                                                 ;;   wu,wv,ws ->  wind speed     
    109 ;;     
     136                                                ;;   wu,wv,ws ->  wind speed 
     137;; 
    110138 
    111139 
    … …  
    117145    uncr=tropflux 
    118146 
    119      extract_flux_tropflux,cor,box, $   
     147     extract_flux_tropflux,cor,box, $ 
    120148        tropflux 
    121149     corr=tropflux 
    122      
     150 
    123151     extract_flux_tropflux,oaf,box, $ 
    124152        tropflux 
    … …  
    133161     ncep1=tropflux 
    134162 
    135     ind=where(finite(q2m)) & q2m=q2m(ind) & uncr_q2m=uncr(ind)  & corr_q2m=corr(ind)  
    136     oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind)  
     163    ind=where(finite(q2m)) & q2m=q2m(ind) & uncr_q2m=uncr(ind)  & corr_q2m=corr(ind) 
     164    oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 
    137165    mean_tao=total(q2m,/nan)/n_elements(ind) 
    138       
     166 
    139167    statistics_3var_v1, q2m, uncr_q2m, corr_q2m, $ 
    140168         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    141   
     169 
    142170    printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    143171    printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    … …  
    153181 
    154182endfor 
    155 ;----------------------------------------------------------------------------------------- 
    156 close,/all   
     183close,/all 
    157184 
    158185fi_q2m_erai='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_erai_v50.txt' 
    … …  
    291318oplot, mean_tao, yfit, color=250, thick=2 
    292319 
    293 ;----------------------------------------------------------- 
    294320closeps 
    295321 
    … …  
    299325return 
    300326end 
    301 ;-------------------------------------------------------------------------- 
    302327function x_site_location, site 
    303328    n1=strpos(site, 's') 
    … …  
    314339return, float(x) 
    315340end 
    316 ;-------------------------------------------------------------------------- 
    317341function y_site_location, site 
    318342    n1=strpos(site, 'e') 
    … …  
    340364return,float(y) 
    341365end 
    342  
    343 ;-------------------------------------------------------------------------- 
    344  

  • trunk/src/paper01/fig3/read_variables_v2.pro

    r41 r43  
    1 ;------------------------------------------------------------------------------------------------- 
     1;+ 
     2; .. _read_variables_v2.pro 
     3; 
     4; ===================== 
     5; read_variables_v2.pro 
     6; ===================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> csite=++ 
     23;  IDL> ++ 
     24;  IDL> read_variables_v2, csite,date1,date2,nsmooth, at, sw,rh,sst,wu,wv,ws, lh 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
    236pro read_variables_v2, csite,date1,date2,nsmooth, $ 
    337                     at, sw,rh,sst,wu,wv,ws, lh 
    438 
    539; 
    6 ;------------------------------------------------------------------------------------------------- 
    740@common 
    841 
    … …  
    75108    lon(tt0(ind))=lon0(ind) 
    76109  endif 
    77 endif  
     110endif 
    78111 
    79112 
    … …  
    122155ind_at=where(at_q ne 1 and at_q ne 2)   &  ind_ws=where(ws_q ne 1 and ws_q ne 2) 
    123156ind_rh=where(rh_q ne 1 and rh_q ne 2)   &  ind_sst=where(sst_q ne 1 and sst_q ne 2) 
    124 ind_lh=where(lh_q ne 1 and lh_q ne 2)  
    125   
     157ind_lh=where(lh_q ne 1 and lh_q ne 2) 
     158 
    126159ind=union(ind_at, union(ind_rh, union(ind_ws, union(ind_lh, ind_sst)))) 
    127160 
    … …  
    135168;; Replace missing values by "NaN" 
    136169; 
    137 tsvars=['at','sw','rh','sst','wu','wv','ws','lat','lon']  
     170tsvars=['at','sw','rh','sst','wu','wv','ws','lat','lon'] 
    138171vars=[tsvars] 
    139172nn=n_elements(vars) 
    … …  
    150183ws=smooth(ws,nsmooth,/nan) & sw=smooth(sw,nsmooth,/nan) 
    151184 
    152 ;------------------------------------------------------------------------------------------------- 
    153185end 
    154 ;------------------------------------------------------------------------------------------------- 
    155186 
    156187 
    157 ;------------------------------------------------------------------------------------------------- 
    158188function time_lec, fi 
    159189tt=ncdf_lec(fi,var='time') 
    … …  
    167197return, tt 
    168198end 
    169  
    170 ;------------------------------------------------------------------------------------------------- 

  • trunk/src/paper01/fig3/sst_validation_scatter_2000_2009_v50.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro sst_validation_scatter_2000_2009_v50,date1,date2 
     1;+ 
     2; .. _sst_validation_scatter_2000_2009_v50.pro 
     3; 
     4; ======================================== 
     5; sst_validation_scatter_2000_2009_v50.pro 
     6; ======================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> sst_validation_scatter_2000_2009_v50, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36pro sst_validation_scatter_2000_2009_v50, date1, date2 
    337@common 
    4 ;------------------------------------------------------------ 
    538 
    639reinitplt, /z,/invert 
    … …  
    841 
    942openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1143;; Before running this program, you have to compile the following subrutines 
    1244;;   1.  .r read_era_total 
    … …  
    3163 
    3264nsmooth=1.    ;; statistics are with 7 day smoothed 
    33 ;------------------------------------------------------------------------------------------------------------------------ 
    3465;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3666close,/all 
    3767 
    … …  
    5888printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    5989 
    60 ;------------------------------------------------------------------------------------------------------------------------ 
    6190;; first reading the whole ERAI uncorrected and corrected data 
    6291 
    6392file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc' 
    6493initncdf, file 
    65 unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15  
     94unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15 
    6695help, unc 
    6796 
    … …  
    78107fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sst_ncep2_oafluxgrid_19890101_20091231.nc' 
    79108initncdf, fi 
    80 nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr)  
     109nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 
    81110nce=nce-273.15 
    82111help, nce 
    … …  
    93122help, nce1 
    94123 
    95 ;------------------------------------------------------------------------------------------------------------------------ 
    96124nn=n_elements(sitelist) 
    97125date1=date1 
    … …  
    122150      uncr=tropflux 
    123151 
    124       extract_flux_tropflux,cor,box, $   
     152      extract_flux_tropflux,cor,box, $ 
    125153        tropflux 
    126154      corr=tropflux 
    127      
     155 
    128156      extract_flux_tropflux,oaf,box, $ 
    129157        tropflux 
    … …  
    143171 
    144172 
    145       ind=where(finite(sst)) & sst=sst(ind) & uncr_sst=uncr(ind)  & corr_sst=corr(ind)  
     173      ind=where(finite(sst)) & sst=sst(ind) & uncr_sst=uncr(ind)  & corr_sst=corr(ind) 
    146174      oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 
    147175      mean_tao=total(sst,/nan)/n_elements(ind) 
    148       
     176 
    149177      statistics_3var_v1, sst, uncr_sst, corr_sst, $ 
    150178         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    151   
     179 
    152180      printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    153181      printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    … …  
    166194 
    167195endfor 
    168 ;----------------------------------------------------------------------------------------- 
    169 close,/all   
     196close,/all 
    170197 
    171198fi_sst_erai='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_v50_erai.txt' 
    … …  
    333360oplot, mean_tao, yfit, color=250, thick=2 
    334361 
    335 ;----------------------------------------------------------- 
    336362closeps 
    337363 
    … …  
    341367return 
    342368end 
    343 ;-------------------------------------------------------------------------- 
    344369function x_site_location, site 
    345370    n1=strpos(site, 's') 
    … …  
    356381return, float(x) 
    357382end 
    358 ;-------------------------------------------------------------------------- 
    359383function y_site_location, site 
    360384    n1=strpos(site, 'e') 
    … …  
    382406return,float(y) 
    383407end 
    384  
    385 ;-------------------------------------------------------------------------- 
    386  

  • trunk/src/paper01/fig3/statistics_3var_v1.pro

    r41 r43  
     1;+ 
     2; .. _statistics_3var_v1.pro 
     3; 
     4; ====================== 
     5; statistics_3var_v1.pro 
     6; ====================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> tao=++ 
     23;  IDL> ++ 
     24;  IDL> statistics_3var_v1, tao, var1, var2, cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36 
     37 
    138pro statistics_3var_v1, tao,var1,var2, $ 
    239    cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    340 
    4   
     41 
    542@common 
    6 ;-------------------------------------------- 
    7 ;-------------------------------------------- 
    843x=tao & y=var1 & z=var2 
    944 
    … …  
    1247 
    1348if (ind ne -1.) then begin 
    14    x(ind)=!Values.f_nan & y(ind)=!Values.f_nan & z(ind)=!Values.f_nan  
     49   x(ind)=!Values.f_nan & y(ind)=!Values.f_nan & z(ind)=!Values.f_nan 
    1550endif 
    1651 
    … …  
    4277 
    4378end 
    44 ;--------------------------------------------    
    45  

  • trunk/src/paper01/fig3/ws_validation_scatter_2000_2009_v50.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _ws_validation_scatter_2000_2009_v50.pro 
     3; 
     4; ======================================= 
     5; ws_validation_scatter_2000_2009_v50.pro 
     6; ======================================= 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> ws_validation_scatter_2000_2009_v50, date1, date2 
     25; 
     26; 
     27; EVOLUTIONS 
     28; ========== 
     29; 
     30; $Id$ 
     31; 
     32; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     33; 
     34;   * minimal header 
     35; 
     36;- 
    237pro ws_validation_scatter_2000_2009_v50,date1,date2 
    338@common 
    4 ;------------------------------------------------------------ 
    539 
    640reinitplt, /z,/invert 
    … …  
    842 
    943openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1144;; Before running this program, you have to compile the following subrutines 
    1245;;   1.  .r read_era_total 
    … …  
    3265da1=10000101 & da2=10081231 
    3366nsmooth=1.    ;; statistics are with 7 day smoothed 
    34 ;------------------------------------------------------------------------------------------------------------------------ 
    3567;;   This program will create the following text files with statistics of respective variables 
    36 ;------------------------------------------------------------------------------------------------------------------------ 
    3768close,/all 
    3869 
    … …  
    5889printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    5990 
    60 ;------------------------------------------------------------------------------------------------------------------------ 
    6191;; first reading the whole ERAI uncorrected and corrected data 
    6292 
    6393file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc' 
    6494initncdf, file 
    65 u=read_ncdf('u10',date1,date2,file=file,/nostr)  
    66 v=read_ncdf('v10',date1,date2,file=file,/nostr)      
    67 unc=sqrt(u*u+v*v)  
     95u=read_ncdf('u10',date1,date2,file=file,/nostr) 
     96v=read_ncdf('v10',date1,date2,file=file,/nostr) 
     97unc=sqrt(u*u+v*v) 
    6898help, unc 
    6999 
    70100file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc' 
    71101initncdf, file 
    72 cor=read_ncdf('ws',date1,date2,file=file,/nostr)  
     102cor=read_ncdf('ws',date1,date2,file=file,/nostr) 
    73103help, cor 
    74104 
    … …  
    80110fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/uwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    81111initncdf, fi 
    82 u=read_ncdf("u", date1-1, date2, file=fi,/nostr)  
     112u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 
    83113fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/vwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    84114initncdf, fi 
    85 v=read_ncdf("v", date1-1, date2, file=fi,/nostr)    
     115v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 
    86116nce=sqrt(u*u+v*v) 
    87117help, nce 
    … …  
    104134help, nce1 
    105135 
    106 ;------------------------------------------------------------------------------------------------------------------------ 
    107136nn=n_elements(sitelist) 
    108137date1=date1 
    … …  
    130159    uncr=tropflux 
    131160 
    132      extract_flux_tropflux,cor,box, $   
     161     extract_flux_tropflux,cor,box, $ 
    133162        tropflux 
    134163     corr=tropflux 
    135      
     164 
    136165     extract_flux_tropflux,oaf,box, $ 
    137166        tropflux 
    … …  
    150179     ncep1=tropflux 
    151180 
    152     ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind)  
     181    ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind) 
    153182    oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 
    154183 
    155184    mean_tao=total(ws)/n_elements(ws) 
    156   
     185 
    157186    statistics_3var_v1, ws, uncr_ws, corr_ws, $ 
    158187         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    159   
     188 
    160189    printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f5.2)' 
    161190    printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f5.2)' 
    … …  
    172201 
    173202endfor 
    174 ;----------------------------------------------------------------------------------------- 
    175 close,/all   
     203close,/all 
    176204fi_ws_erai='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_erai_v50.txt' 
    177205res=read_ascii(fi_ws_erai,data_start=1) 
    … …  
    338366oplot, mean_tao, yfit, color=250, thick=2 
    339367 
    340 ;----------------------------------------------------------- 
    341368closeps 
    342369 
    … …  
    346373return 
    347374end 
    348 ;-------------------------------------------------------------------------- 
    349375function x_site_location, site 
    350376    n1=strpos(site, 's') 
    … …  
    361387return, float(x) 
    362388end 
    363 ;-------------------------------------------------------------------------- 
    364389function y_site_location, site 
    365390    n1=strpos(site, 'e') 
    … …  
    387412return,float(y) 
    388413end 
    389  
    390 ;-------------------------------------------------------------------------- 
    391  

  • trunk/src/paper01/fig4/fig4_met_var_correction_scatter.pro

    r41 r43  
     1;+ 
     2; .. _fig4_met_var_correction_scatter.pro 
     3; 
     4; =================================== 
     5; fig4_met_var_correction_scatter.pro 
     6; =================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> fig4_met_var_correction_scatter 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
     34 
    135pro fig4_met_var_correction_scatter 
    236@common 
    3 ;------------------------------------------------------------ 
    437reinitplt, /z,/invert 
    538key_portrait = 1 
    639 
    740openps, FILENAME = 'idl.ps' 
    8 ;------------------------------------------------------------ 
    941fi='/Users/pkb/work/MY_SAXO/flux_automat/sst_correction_final.txt' 
    1042res=read_ascii(fi,data_start=0) 
    … …  
    3466significance_test_99,mean,bias 
    3567 
    36 ;----------------------- 
    3768fi='/Users/pkb/work/MY_SAXO/flux_automat/t2m_correction_final.txt' 
    3869res=read_ascii(fi,data_start=0) 
    … …  
    6394significance_test_99,mean,bias 
    6495 
    65 ;----------------------- 
    6696fi='/Users/pkb/work/MY_SAXO/flux_automat/q2m_correction_final.txt' 
    6797res=read_ascii(fi,data_start=0) 
    … …  
    90120;; significant test 
    91121significance_test_99,mean,bias 
    92 ;----------------------- 
    93122fi='/Users/pkb/work/MY_SAXO/flux_automat/ws_correction_final.txt' 
    94123res=read_ascii(fi,data_start=0) 
    … …  
    118147;; significant test 
    119148significance_test_99,mean,bias 
    120 ;----------------------------------------------------------------------------------------- 
    121149 
    122150closeps 
    … …  
    127155 
    128156end 
    129 ;-------------------------------------------------------------------------- 
    130  

  • trunk/src/paper01/fig4/significance_test_99.pro

    r41 r43  
    1 pro significance_test_99, x,y, $ 
    2    Sig_value 
     1;+ 
     2; .. _significance_test_99.pro 
     3; 
     4; ======================== 
     5; significance_test_99.pro 
     6; ======================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> x=++ 
     23;  IDL> y=++ 
     24;  IDL> sig_value=++ 
     25;  IDL> significance_test_99, x, y, sig_value 
     26; 
     27; EVOLUTIONS 
     28; ========== 
     29; 
     30; $Id$ 
     31; 
     32; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     33; 
     34;   * minimal header 
     35; 
     36;- 
     37pro significance_test_99, x, y, sig_value 
    338 
    439   X=x & Y=y 
    … …  
    2156   print, 'Statistics are not valid at 99% significance level' 
    2257endelse 
    23   
     58 
    2459end 
    25  

  • trunk/src/paper01/fig5/fig5_gustiness_scatter.pro

    r41 r43  
     1;+ 
     2; .. _fig5_gustiness_scatter.pro 
     3; 
     4; ========================== 
     5; fig5_gustiness_scatter.pro 
     6; ========================== 
     7; 
     8; ; 
     9; DESCRIPTION 
     10; =========== 
     11; 
     12; SEE ALSO 
     13; ======== 
     14; 
     15; :ref:`tropflux_profile.sh` 
     16; 
     17; EXAMPLES 
     18; ======== 
     19; 
     20; :: 
     21; 
     22;  IDL> @tropflux_init 
     23;  IDL> fig5_gustiness_scatter 
     24; 
     25; EVOLUTIONS 
     26; ========== 
     27; 
     28; $Id$ 
     29; 
     30; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     31; 
     32;   * minimal header 
     33; 
     34;- 
    135pro fig5_gustiness_scatter 
    236@common 
    3 ;----------------------------------------------------------------------- 
    437reinitplt, /z,/invert 
    538key_portrait = 1 
    639 
    740openps, FILENAME = 'idl.ps' 
    8 ;----------------------------------------------------------------------- 
    941fi_ws='/Users/pkb/work/MY_SAXO/flux_automat/wind_gustiness_bias.txt' 
    1042res=read_ascii(fi_ws,data_start=1) 
    … …  
    2254print, correlate(erai_sst,gust) 
    2355 
    24 ;----------------------------------------------------------------------- 
    2556splot, erai_wind, gust, small=[2,3,1], title='a) ERAI Wind Vs Gustiness', $ 
    2657      xtitle='ERAI wind', ytitle='Cronin gustiness', charsize=1.2, psym=2 , $ 
    … …  
    4475xyouts, 21,-0.4, 'y = 0.175712 x -3.17303', charsize=1.2 
    4576print, a,b 
    46 ;----------------------------------------------------------------------- 
    4777closeps 
    4878fig='fig5_gustiness_scatter.ps' 
    … …  
    5080spawn, 'gv '+cpsdir+fig 
    5181end 
    52  

  • trunk/src/paper01/fig6/icoads_q2m_stats_paper.pro

    r41 r43  
     1;+ 
     2; .. _icoads_q2m_stats_paper.pro 
     3; 
     4; ========================== 
     5; icoads_q2m_stats_paper.pro 
     6; ========================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> icoads_q2m_stats_paper 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro icoads_q2m_stats_paper 
    235@common 
    3 ;----------------------------------------------- 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242marge=[-2,-2, -4,2] 
    1343st=19890101 & en=20091231 
    … …  
    178208xyouts, 10, 4, abs_ncep1, charsize=1 
    179209 
    180 ;------------------------------------------------------------ 
    181210closeps 
    182211fig="icoads_q2m_stats_paper.ps" 
    183212spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    184213spawn, 'gv '+cpsdir+fig 
    185 ;------------------------------------------------------------ 
    186214return 
    187215end 

  • trunk/src/paper01/fig6/icoads_sst_stats_paper.pro

    r41 r43  
     1;+ 
     2; .. _icoads_sst_stats_paper.pro 
     3; 
     4; ========================== 
     5; icoads_sst_stats_paper.pro 
     6; ========================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> icoads_sst_stats_paper 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro icoads_sst_stats_paper 
    235@common 
    3 ;----------------------------------------------- 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242marge=[-2,-2, -4,2] 
    1343st=19890101 & en=20091231 
    … …  
    95125trop=trop*mask & oaf=oaf*mask & erai=erai*mask 
    96126ncep2=ncep2*mask & ncep1=ncep1*mask 
    97 noc=noc*mask  
     127noc=noc*mask 
    98128 
    99129ind=where(finite(param) and finite(erai) and finite(trop) and finite(oaf) and finite(ncep2) and finite(ncep1)) 
    … …  
    205235xyouts, xp, yp-1, abs_noc, charsize=1 
    206236 
    207 ;------------------------------------------------------------ 
    208237closeps 
    209238fig="icoads_sst_stats_paper.ps" 
    210239spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    211240spawn, 'gv '+cpsdir+fig 
    212 ;------------------------------------------------------------ 
    213241return 
    214 end 

  • trunk/src/paper01/fig6/icoads_t2m_stats_paper.pro

    r41 r43  
     1;+ 
     2; .. _icoads_t2m_stats_paper.pro 
     3; 
     4; ========================== 
     5; icoads_t2m_stats_paper.pro 
     6; ========================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> icoads_t2m_stats_paper 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro icoads_t2m_stats_paper 
    235@common 
    3 ;----------------------------------------------- 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242marge=[-2,-2, -4,2] 
    1343st=19890101 & en=20091231 
    … …  
    186216xyouts, xp, yp-1, abs_ncep1, charsize=1 
    187217 
    188 ;------------------------------------------------------------ 
    189218closeps 
    190219fig="icoads_t2m_stats_paper.ps" 
    191220spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    192221spawn, 'gv '+cpsdir+fig 
    193 ;------------------------------------------------------------ 
    194222return 
    195223end 

  • trunk/src/paper01/fig6/icoads_ws_stats_paper.pro

    r41 r43  
     1;+ 
     2; .. _icoads_ws_stats_paper.pro 
     3; 
     4; ========================= 
     5; icoads_ws_stats_paper.pro 
     6; ========================= 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> icoads_ws_stats_paper 
     23; 
     24; EVOLUTIONS 
     25; ========== 
     26; 
     27; $Id$ 
     28; 
     29; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     30; 
     31;   * minimal header 
     32; 
     33;- 
    134pro icoads_ws_stats_paper 
    235@common 
    3 ;----------------------------------------------- 
    436reinitplt, /z,/invert 
    537key_portrait = 1 
    … …  
    739 
    840openps, FILENAME = 'idl.ps' 
    9 ;------------------------------------------------------------ 
    1041; partie a changer 
    11 ;------------------------------------------------------------ 
    1242marge=[-2,-2, -4,2] 
    1343st=19890101 & en=20091231 
    … …  
    202232xyouts, xp, yp-1, abs_noc, charsize=1 
    203233 
    204 ;------------------------------------------------------------ 
    205234closeps 
    206235fig="icoads_ws_stats_paper.ps" 
    207236spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 
    208237spawn, 'gv '+cpsdir+fig 
    209 ;------------------------------------------------------------ 
    210238return 
    211239end 

  • trunk/src/paper01/fig6/statistics.pro

    r41 r43  
     1;+ 
     2; .. _statistics.pro 
     3; 
     4; ============== 
     5; statistics.pro 
     6; ============== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> moor=++ 
     23;  IDL> ++ 
     24;  IDL> statistics, moor, erai, cor, bias, std, rmsd 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36; 
    137pro statistics, moor, erai, $ 
    238   cor, bias, std, rmsd 
    339@common 
    4 ;-------------------------------------------- 
    540;; mean1= mean of mooring variable      std1=std of mooring variable 
    641;; mean2= mean of erai variable         std2=std of erai variable 
    7 ;-------------------------------------------- 
    842x=moor & y=erai 
    943ind1=where(finite(x,/nan)) & ind2=where(finite(y,/nan)) 
    … …  
    2660rmsd=sqrt(total((x-y)*(x-y))/n_elements(x)) 
    2761end 
    28 ;--------------------------------------------    
    29  

  • trunk/src/paper01/fig7/air_validation_scatter_2000_2009_basin.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _air_validation_scatter_2000_2009.pro 
     3; 
     4; ==================================== 
     5; air_validation_scatter_2000_2009.pro 
     6; ==================================== 
     7; DESCRIPTION 
     8; =========== 
     9; 
     10; SEE ALSO 
     11; ======== 
     12; 
     13; :ref:`tropflux_profile.sh` 
     14; 
     15; EXAMPLES 
     16; ======== 
     17; 
     18; :: 
     19; 
     20;  IDL> @tropflux_init 
     21;  IDL> .r read_era_total 
     22;  IDL> .r read_variables_v2 
     23;  IDL> .r read_era_box 
     24;  IDL> date1=++ 
     25;  IDL> date2=++ 
     26;  IDL> air_validation_scatter_2000_2009, date1, date2 
     27; 
     28; EVOLUTIONS 
     29; ========== 
     30; 
     31; $Id$ 
     32; 
     33; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     34; 
     35;   * minimal header 
     36; 
     37;- 
    238pro air_validation_scatter_2000_2009,date1,date2 
    339@common 
    4 ;------------------------------------------------------------ 
    540 
    641reinitplt, /z,/invert 
    … …  
    843 
    944openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    11 ;; Before running this program, you have to compile the following subrutines 
    12 ;;   1.  .r read_era_total 
    13 ;;   2.  .r read_variables_v2 
    14 ;;   3.  .r read_era_box 
    1545 
    1646;; Give the location of mooring for validation of basic meteorological variables 
    … …  
    6595 
    6696nsmooth=1.    ;; statistics are with 7 day smoothed 
    67 ;------------------------------------------------------------------------------------------------------------------------ 
    6897;;   This program will create the following text files with statistics of respective variables 
    69 ;------------------------------------------------------------------------------------------------------------------------ 
    7098close,/all 
    7199 
    … …  
    88116printf,5, 'x     y      cor    bias     std     rmsd    mean_tao' 
    89117 
    90 ;------------------------------------------------------------------------------------------------------------------------ 
    91118;; first reading the whole ERAI uncorrected and corrected data 
    92119 
    93120file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc' 
    94121initncdf, file 
    95 unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15  
     122unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 
    96123help, unc 
    97124 
    … …  
    109136;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 
    110137;initncdf, fi 
    111 ;nce=read_ncdf("air", date1-1, date2, file=fi,/nostr)  
     138;nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 
    112139;nce=nce-273.15 
    113140;help, nce 
    … …  
    115142;file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/t2m_ncep1_19890101_20091231.nc' 
    116143;initncdf, file 
    117 ;nce1=read_ncdf("t2m", date1, date2, file=file,/nostr)  
     144;nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 
    118145;nce1=nce1-273.15 
    119146;help, nce1 
    … …  
    123150help, unc, cor, oaf, nce, nce1 
    124151 
    125 ;------------------------------------------------------------------------------------------------------------------------ 
    126152nn=n_elements(sitelist) 
    127153date1=date1 
    … …  
    145171                                                ;;        lh  ->  latent heat flux 
    146172                                                ;;        rh  ->  relative humidity 
    147                                                 ;;   wu,wv,ws ->  wind speed     
    148 ;;     
     173                                                ;;   wu,wv,ws ->  wind speed 
     174;; 
    149175 
    150176 
    … …  
    156182    uncr=tropflux 
    157183 
    158      extract_flux_tropflux,cor,box, $   
     184     extract_flux_tropflux,cor,box, $ 
    159185        tropflux 
    160186     corr=tropflux 
    161      
     187 
    162188     extract_flux_tropflux,oaf,box, $ 
    163189        tropflux 
    … …  
    173199 
    174200 
    175     ind=where(finite(air)) & air=air(ind) & uncr_air=uncr(ind)  & corr_air=corr(ind)  
     201    ind=where(finite(air)) & air=air(ind) & uncr_air=uncr(ind)  & corr_air=corr(ind) 
    176202    oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 
    177203    mean_tao=total(air,/nan)/n_elements(ind) 
    178       
     204 
    179205    statistics_3var_v1, air, uncr_air, corr_air, $ 
    180206         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    181   
     207 
    182208    printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    183209    printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    … …  
    193219 
    194220endfor 
    195 ;----------------------------------------------------------------------------------------- 
    196 close,/all   
     221close,/all 
    197222 
    198223fi_air_erai='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_'+ocean+'.txt' 
    … …  
    332357oplot, mean_tao, yfit, color=250, thick=2 
    333358 
    334 ;----------------------------------------------------------- 
    335359closeps 
    336360 
    … …  
    340364return 
    341365end 
    342 ;-------------------------------------------------------------------------- 
    343366function x_site_location, site 
    344367    n1=strpos(site, 's') 
    … …  
    355378return, float(x) 
    356379end 
    357 ;-------------------------------------------------------------------------- 
    358380function y_site_location, site 
    359381    n1=strpos(site, 'e') 
    … …  
    381403return,float(y) 
    382404end 
    383  
    384 ;-------------------------------------------------------------------------- 
    385  

  • trunk/src/paper01/fig7/q2m_validation_scatter_2000_2009_basin.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro q2m_validation_scatter_2000_2009_basin,date1,date2 
     1;+ 
     2; .. _q2m_validation_scatter_2000_2009_basin.pro 
     3; 
     4; ========================================== 
     5; q2m_validation_scatter_2000_2009_basin.pro 
     6; ========================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> q2m_validation_scatter_2000_2009_basin, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36 
     37pro q2m_validation_scatter_2000_2009_basin, date1, date2 
    338@common 
    4 ;------------------------------------------------------------ 
    539 
    640reinitplt, /z,/invert 
    … …  
    842 
    943openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1144;; Before running this program, you have to compile the following subrutines 
    1245;;   1.  .r read_era_total 
    … …  
    6598 
    6699nsmooth=1.    ;; statistics are with 7 day smoothed 
    67 ;------------------------------------------------------------------------------------------------------------------------ 
    68100;;   This program will create the following text files with statistics of respective variables 
    69 ;------------------------------------------------------------------------------------------------------------------------ 
    70101close,/all 
    71102 
    … …  
    88119printf,5, 'x     y      cor    bias     std     rmsd    mean_tao' 
    89120 
    90 ;------------------------------------------------------------------------------------------------------------------------ 
    91121;; first reading the whole ERAI uncorrected and corrected data 
    92122 
    … …  
    98128;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc' 
    99129;initncdf, file 
    100 ;cor=read_ncdf('q2m',date1,date2,file=file,/nostr)  
     130;cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 
    101131;help, cor 
    102132 
    … …  
    108138;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sphum_ncep2_oafluxgrid_19890101_20091231.nc' 
    109139;initncdf, fi 
    110 ;nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr)  
     140;nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 
    111141;nce=nce*1000 
    112142;help, nce 
    … …  
    122152help, unc, cor, oaf, nce, nce1 
    123153 
    124 ;------------------------------------------------------------------------------------------------------------------------ 
    125154nn=n_elements(sitelist) 
    126155date1=date1 
    … …  
    144173                                                ;;        lh  ->  latent heat flux 
    145174                                                ;;        rh  ->  relative humidity 
    146                                                 ;;   wu,wv,ws ->  wind speed     
    147 ;;     
     175                                                ;;   wu,wv,ws ->  wind speed 
     176;; 
    148177 
    149178 
    … …  
    155184    uncr=tropflux 
    156185 
    157      extract_flux_tropflux,cor,box, $   
     186     extract_flux_tropflux,cor,box, $ 
    158187        tropflux 
    159188     corr=tropflux 
    160      
     189 
    161190     extract_flux_tropflux,oaf,box, $ 
    162191        tropflux 
    … …  
    171200     ncep1=tropflux 
    172201 
    173     ind=where(finite(q2m)) & q2m=q2m(ind) & uncr_q2m=uncr(ind)  & corr_q2m=corr(ind)  
    174     oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind)  
     202    ind=where(finite(q2m)) & q2m=q2m(ind) & uncr_q2m=uncr(ind)  & corr_q2m=corr(ind) 
     203    oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 
    175204    mean_tao=total(q2m,/nan)/n_elements(ind) 
    176       
     205 
    177206    statistics_3var_v1, q2m, uncr_q2m, corr_q2m, $ 
    178207         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    179   
     208 
    180209    printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    181210    printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    … …  
    191220 
    192221endfor 
    193 ;----------------------------------------------------------------------------------------- 
    194 close,/all   
     222close,/all 
    195223 
    196224fi_q2m_erai='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_erai_'+ocean+'.txt' 
    … …  
    329357oplot, mean_tao, yfit, color=250, thick=2 
    330358 
    331 ;----------------------------------------------------------- 
    332359closeps 
    333360 
    … …  
    337364return 
    338365end 
    339 ;-------------------------------------------------------------------------- 
    340366function x_site_location, site 
    341367    n1=strpos(site, 's') 
    … …  
    352378return, float(x) 
    353379end 
    354 ;-------------------------------------------------------------------------- 
    355380function y_site_location, site 
    356381    n1=strpos(site, 'e') 
    … …  
    378403return,float(y) 
    379404end 
    380  
    381 ;-------------------------------------------------------------------------- 
    382  

  • trunk/src/paper01/fig7/sst_validation_scatter_2000_2009_basin.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _sst_validation_scatter_2000_2009_basin.pro 
     3; 
     4; ========================================== 
     5; sst_validation_scatter_2000_2009_basin.pro 
     6; ========================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> sst_validation_scatter_2000_2009_basin, date1, date2 
     25; 
     26; 
     27; EVOLUTIONS 
     28; ========== 
     29; 
     30; $Id$ 
     31; 
     32; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     33; 
     34;   * minimal header 
     35; 
     36;- 
    237pro sst_validation_scatter_2000_2009_basin,date1,date2 
    338@common 
    4 ;------------------------------------------------------------ 
    539 
    640reinitplt, /z,/invert 
    … …  
    842 
    943openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1144;; Before running this program, you have to compile the following subrutines 
    1245;;   1.  .r read_era_total 
    … …  
    4376          '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w'] 
    4477 
    45 nsmooth=1.     
    46 ;------------------------------------------------------------------------------------------------------------------------ 
     78nsmooth=1. 
    4779;;   This program will create the following text files with statistics of respective variables 
    48 ;------------------------------------------------------------------------------------------------------------------------ 
    4980close,/all 
    5081 
    … …  
    71102printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    72103 
    73 ;------------------------------------------------------------------------------------------------------------------------ 
    74104;; first reading the whole ERAI uncorrected and corrected data 
    75105 
    76106;file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc' 
    77107;initncdf, file 
    78 ;unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15  
     108;unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15 
    79109;help, unc 
    80110 
    … …  
    91121;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sst_ncep2_oafluxgrid_19890101_20091231.nc' 
    92122;initncdf, fi 
    93 ;nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr)  
     123;nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 
    94124;nce=nce-273.15 
    95125;help, nce 
    … …  
    112142help, unc, cor, oaf, nce, sst_tmi, nce1 
    113143;; data from 20000101, 20090931 
    114 ;------------------------------------------------------------------------------------------------------------------------ 
    115144nn=n_elements(sitelist) 
    116145date1=date1 
    … …  
    141170      uncr=tropflux 
    142171 
    143       extract_flux_tropflux,cor,box, $   
     172      extract_flux_tropflux,cor,box, $ 
    144173        tropflux 
    145174      corr=tropflux 
    146      
     175 
    147176      extract_flux_tropflux,oaf,box, $ 
    148177        tropflux 
    … …  
    162191 
    163192 
    164       ind=where(finite(sst)) & sst=sst(ind) & uncr_sst=uncr(ind)  & corr_sst=corr(ind)  
     193      ind=where(finite(sst)) & sst=sst(ind) & uncr_sst=uncr(ind)  & corr_sst=corr(ind) 
    165194      oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 
    166195      mean_tao=total(sst,/nan)/n_elements(ind) 
    167       
     196 
    168197      statistics_3var_v1, sst, uncr_sst, corr_sst, $ 
    169198         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    170   
     199 
    171200      printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    172201      printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f6.2)' 
    … …  
    185214 
    186215endfor 
    187 ;----------------------------------------------------------------------------------------- 
    188 close,/all   
     216close,/all 
    189217 
    190218fi_sst_erai='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_erai_'+ocean+'.txt' 
    … …  
    352380oplot, mean_tao, yfit, color=250, thick=2 
    353381 
    354 ;----------------------------------------------------------- 
    355382closeps 
    356383 
    … …  
    360387return 
    361388end 
    362 ;-------------------------------------------------------------------------- 
    363389function x_site_location, site 
    364390    n1=strpos(site, 's') 
    … …  
    375401return, float(x) 
    376402end 
    377 ;-------------------------------------------------------------------------- 
    378403function y_site_location, site 
    379404    n1=strpos(site, 'e') 
    … …  
    401426return,float(y) 
    402427end 
    403  
    404 ;-------------------------------------------------------------------------- 
    405  

  • trunk/src/paper01/fig7/ws_validation_scatter_2000_2009_basin.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; 
     3; .. _ws_validation_scatter_2000_2009_basin.pro 
     4; 
     5; ========================================= 
     6; ws_validation_scatter_2000_2009_basin.pro 
     7; ========================================= 
     8; 
     9; DESCRIPTION 
     10; =========== 
     11; 
     12; SEE ALSO 
     13; ======== 
     14; 
     15; :ref:`tropflux_profile.sh` 
     16; 
     17; EXAMPLES 
     18; ======== 
     19; 
     20; :: 
     21; 
     22;  IDL> @tropflux_init 
     23;  IDL> date1=++ 
     24;  IDL> date2=++ 
     25;  IDL> ws_validation_scatter_2000_2009_basin, date1, date2 
     26; 
     27; EVOLUTIONS 
     28; ========== 
     29; 
     30; $Id$ 
     31; 
     32; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     33; 
     34;   * minimal header 
     35; 
     36;- 
    237pro ws_validation_scatter_2000_2009_basin,date1,date2 
    338@common 
    4 ;------------------------------------------------------------ 
    539 
    640reinitplt, /z,/invert 
    … …  
    842 
    943openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1144;; Before running this program, you have to compile the following subrutines 
    1245;;   1.  .r read_era_total 
    … …  
    4982sitelist=['8s125w', '8s110w', '8s95w','5s155w', '5s140w', '5s125w', '5s110w', '5s95w','2s155w', '2s140w', $ 
    5083          '2s125w', '2s110w', '2s95w', '0n155w', '0n140w', '0n125w', '0n110w', '0n95w', '2n155w', '2n140w', $ 
    51           '2n125w', '2n110w', '2n95w',  '5n155w', '5n140w', '5n125w', '5n110w', '5n95w']  
     84          '2n125w', '2n110w', '2n95w',  '5n155w', '5n140w', '5n125w', '5n110w', '5n95w'] 
    5285ocean='Pac_coldtongue' 
    5386 
    … …  
    6598 
    6699nsmooth=1 
    67 ;------------------------------------------------------------------------------------------------------------------------ 
    68100;;   This program will create the following text files with statistics of respective variables 
    69 ;------------------------------------------------------------------------------------------------------------------------ 
    70101close,/all 
    71102 
    … …  
    91122printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    92123 
    93 ;------------------------------------------------------------------------------------------------------------------------ 
    94124;; first reading the whole ERAI uncorrected and corrected data 
    95125;; date1=20000101 & date2=20090931 
    … …  
    97127file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc' 
    98128initncdf, file 
    99 u=read_ncdf('u10',date1,date2,file=file,/nostr)  
    100 ;v=read_ncdf('v10',date1,date2,file=file,/nostr)      
    101 ;unc=sqrt(u*u+v*v)  
     129u=read_ncdf('u10',date1,date2,file=file,/nostr) 
     130;v=read_ncdf('v10',date1,date2,file=file,/nostr) 
     131;unc=sqrt(u*u+v*v) 
    102132;help, unc 
    103133 
    104134;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc' 
    105135;initncdf, file 
    106 ;cor=read_ncdf('ws',date1,date2,file=file,/nostr)  
     136;cor=read_ncdf('ws',date1,date2,file=file,/nostr) 
    107137;help, cor 
    108138 
    … …  
    114144;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/uwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    115145;initncdf, fi 
    116 ;u=read_ncdf("u", date1-1, date2, file=fi,/nostr)  
     146;u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 
    117147;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/vwind_ncep2_oafluxgrid_19890101_20091231.nc' 
    118148;initncdf, fi 
    119 ;v=read_ncdf("v", date1-1, date2, file=fi,/nostr)    
     149;v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 
    120150;nce=sqrt(u*u+v*v) 
    121151;help, nce 
    … …  
    142172help, unc, cor, oaf, nce, ws_tmi, nce1 
    143173 
    144 ;------------------------------------------------------------------------------------------------------------------------ 
    145174nn=n_elements(sitelist) 
    146175date1=date1 
    … …  
    168197    uncr=tropflux 
    169198 
    170      extract_flux_tropflux,cor,box, $   
     199     extract_flux_tropflux,cor,box, $ 
    171200        tropflux 
    172201     corr=tropflux 
    173      
     202 
    174203     extract_flux_tropflux,oaf,box, $ 
    175204        tropflux 
    … …  
    188217     ncep1=tropflux 
    189218 
    190     ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind)  
     219    ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind) 
    191220    oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 
    192221 
    193222    mean_tao=total(ws)/n_elements(ws) 
    194   
     223 
    195224    statistics_3var_v1, ws, uncr_ws, corr_ws, $ 
    196225         cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    197   
     226 
    198227    printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f5.2)' 
    199228    printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f5.2)' 
    … …  
    210239 
    211240endfor 
    212 ;----------------------------------------------------------------------------------------- 
    213 close,/all   
     241close,/all 
    214242fi_ws_erai='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_erai_'+ocean+'.txt' 
    215243res=read_ascii(fi_ws_erai,data_start=1) 
    … …  
    376404oplot, mean_tao, yfit, color=250, thick=2 
    377405 
    378 ;----------------------------------------------------------- 
    379406closeps 
    380407 
    … …  
    384411return 
    385412end 
    386 ;-------------------------------------------------------------------------- 
    387413function x_site_location, site 
    388414    n1=strpos(site, 's') 
    … …  
    399425return, float(x) 
    400426end 
    401 ;-------------------------------------------------------------------------- 
    402427function y_site_location, site 
    403428    n1=strpos(site, 'e') 
    … …  
    425450return,float(y) 
    426451end 
    427  
    428 ;-------------------------------------------------------------------------- 
    429  

  • trunk/src/paper01/fig8/lhf_validation_scatter_2000_2009.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _lhf_validation_scatter_2000_2009.pro 
     3; 
     4; ==================================== 
     5; lhf_validation_scatter_2000_2009.pro 
     6; ==================================== 
     7; 
     8;  DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> lhf_validation_scatter_2000_2009, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
    236pro lhf_validation_scatter_2000_2009,date1,date2 
    337@common 
    4 ;------------------------------------------------------------ 
    538 
    639reinitplt, /z,/invert 
    … …  
    841 
    942openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1143;; Before running this program, you have to compile the following subrutines 
    1244 
    … …  
    3163ocean='global' 
    3264 
    33 ;------------------------------------------------------------------------------------------------------------------------ 
    3465;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3666close,/all 
    3767 
    … …  
    5383printf,4, 'x     y      cor    bias     std     rmsd      mean_tao' 
    5484 
    55 ;------------------------------------------------------------------------------------------------------------------------ 
    5685;; first reading the whole ERAI uncorrected and corrected data 
    5786 
    5887file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_lhf_19890101_20091231_oafluxgrid.nc' 
    5988initncdf, file 
    60 unc=read_ncdf('lhf',date1,date2,file=file,/nostr)  
     89unc=read_ncdf('lhf',date1,date2,file=file,/nostr) 
    6190help, unc 
    6291 
    … …  
    74103fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/lhf_ncep2_oafluxgrid_19890101_20091231.nc' 
    75104initncdf, fi 
    76 nce=read_ncdf("lhf", date1, date2, file=fi,/nostr)  
     105nce=read_ncdf("lhf", date1, date2, file=fi,/nostr) 
    77106help, nce 
    78107 
    79108file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/fluxe_ncep1_19890101_20091231.nc' 
    80109initncdf, file 
    81 nce1=-1*read_ncdf("lhf", date1, date2, file=file,/nostr)  
     110nce1=-1*read_ncdf("lhf", date1, date2, file=file,/nostr) 
    82111help, nce1 
    83112 
    84113nsmooth=1. 
    85 ;------------------------------------------------------------------------------------------------------------------------ 
    86114nn=n_elements(sitelist) 
    87115date1=date1 
    … …  
    101129 
    102130    lhf=lh & ind=where(finite(lhf)) & valid=n_elements(ind) 
    103   
     131 
    104132    if (valid ge 180. ) then begin 
    105133 
    … …  
    108136        uncr=tropflux 
    109137 
    110         extract_flux_tropflux,cor,box, $   
     138        extract_flux_tropflux,cor,box, $ 
    111139                 tropflux 
    112140        corr=tropflux 
    113      
     141 
    114142        extract_flux_tropflux,oaf,box, $ 
    115143                tropflux 
    … …  
    130158        statistics_3var_v1, lhf, uncr_lhf, corr_lhf, $ 
    131159                 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    132   
     160 
    133161        printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    134162        printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    … …  
    170198 
    171199splot, mean_tao, mean_erai, charsize=1.1, title='LHF - TAO Vs ERAI', $ 
    172      xrange=[20,200], yrange=[20,200], small=[2,3,1], psym=2  
     200     xrange=[20,200], yrange=[20,200], small=[2,3,1], psym=2 
    173201xyouts, 50,25, cstat, charsize=1. 
    174202xyouts, 50,10, 'cor    bias    std    rmsd', charsize=1. 
    … …  
    305333oplot, [20,200], [20,200] 
    306334 
    307 ;----------------------------------------------------------- 
    308335closeps 
    309336 
    … …  
    313340return 
    314341end 
    315 ;-------------------------------------------------------------------------- 
    316342function x_site_location, site 
    317343    n1=strpos(site, 's') 
    … …  
    328354return, float(x) 
    329355end 
    330 ;-------------------------------------------------------------------------- 
    331356function y_site_location, site 
    332357    n1=strpos(site, 'e') 
    … …  
    354379return,float(y) 
    355380end 
    356  
    357 ;-------------------------------------------------------------------------- 
    358  

  • trunk/src/paper01/fig8/shf_validation_scatter_2000_2009.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _shf_validation_scatter_2000_2009.pro 
     3; 
     4; ==================================== 
     5; shf_validation_scatter_2000_2009.pro 
     6; ==================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> shf_validation_scatter_2000_2009, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
    236pro shf_validation_scatter_2000_2009,date1,date2 
    337@common 
    4 ;------------------------------------------------------------ 
    538 
    639reinitplt, /z,/invert 
    … …  
    841 
    942openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1143;; Before running this program, you have to compile the following subrutines 
    1244;;   1.  .r read_variables_sh 
    … …  
    3062ocean='global' 
    3163 
    32 nsmooth=1.  
    33 ;------------------------------------------------------------------------------------------------------------------------ 
     64nsmooth=1. 
    3465;;   This program will create the following text files with statistics of respective variables 
    35 ;------------------------------------------------------------------------------------------------------------------------ 
    3666close,/all 
    3767 
    … …  
    5585printf,5, 'x     y      cor    bias     std     rmsd      mean_tao' 
    5686 
    57 ;------------------------------------------------------------------------------------------------------------------------ 
    5887;; first reading the whole ERAI uncorrected and corrected data 
    5988 
    6089file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_shf_19890101_20091231_oafluxgrid.nc' 
    6190initncdf, file 
    62 unc=read_ncdf('shf',date1,date2,file=file,/nostr)  
     91unc=read_ncdf('shf',date1,date2,file=file,/nostr) 
    6392help, unc 
    6493 
    … …  
    76105fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/shf_ncep2_oafluxgrid_19890101_20091231.nc' 
    77106initncdf, fi 
    78 nce=read_ncdf("shf", date1, date2, file=fi,/nostr)  
     107nce=read_ncdf("shf", date1, date2, file=fi,/nostr) 
    79108help, nce 
    80109 
    … …  
    84113help, nce1 
    85114 
    86 ;------------------------------------------------------------------------------------------------------------------------ 
    87115nn=n_elements(sitelist) 
    88116date1=date1 
    … …  
    102130         at, sw,rh,sst,wu,wv,ws,sh 
    103131    shf=sh & ind=where(finite(shf)) & valid=n_elements(ind) 
    104   
     132 
    105133    if (valid ge 180. ) then begin 
    106134 
    … …  
    109137        uncr=tropflux 
    110138 
    111         extract_flux_tropflux,cor,box, $   
     139        extract_flux_tropflux,cor,box, $ 
    112140                 tropflux 
    113141        corr=tropflux 
    114      
     142 
    115143        extract_flux_tropflux,oaf,box, $ 
    116144                tropflux 
    … …  
    130158        statistics_3var_v1, shf, uncr_shf, corr_shf, $ 
    131159                 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    132   
     160 
    133161        printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    134162        printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    … …  
    305333oplot, [-5,25], [-5,25] 
    306334 
    307 ;----------------------------------------------------------- 
    308335closeps 
    309336 
    … …  
    313340return 
    314341end 
    315 ;-------------------------------------------------------------------------- 
    316342function x_site_location, site 
    317343    n1=strpos(site, 's') 
    … …  
    328354return, float(x) 
    329355end 
    330 ;-------------------------------------------------------------------------- 
    331356function y_site_location, site 
    332357    n1=strpos(site, 'e') 
    … …  
    354379return,float(y) 
    355380end 
    356  
    357 ;-------------------------------------------------------------------------- 
    358  

  • trunk/src/paper01/fig9/lwr_validation_scatter_2000_2007.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
     1;+ 
     2; .. _lwr_validation_scatter_2000_2007.pro 
     3; 
     4; ==================================== 
     5; lwr_validation_scatter_2000_2007.pro 
     6; ==================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> lwr_validation_scatter_2000_2007, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
    236pro lwr_validation_scatter_2000_2007,date1,date2 
    337@common 
    4 ;------------------------------------------------------------ 
    538 
    639reinitplt, /z,/invert 
    … …  
    841 
    942openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1143;; Before running this program, you have to compile the following subrutines 
    1244;;   1.  .r read_lw_global_v2 
    … …  
    3466 
    3567da1=20000101 & da2=20071231 
    36 nsmooth=1.     
    37 ;------------------------------------------------------------------------------------------------------------------------ 
     68nsmooth=1. 
    3869;;   This program will create the following text files with statistics of respective variables 
    39 ;------------------------------------------------------------------------------------------------------------------------ 
    4070close,/all 
    4171 
    … …  
    6090printf,6, 'x     y      cor    bias     std     rmsd' 
    6191 
    62 ;------------------------------------------------------------------------------------------------------------------------ 
    6392;; first reading the whole ERAI uncorrected and corrected data 
    6493 
    6594file='/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/erai_fluxes_19930101_20090801_TROP_oafluxgrid.nc' 
    6695initncdf, file 
    67 unc=read_ncdf('lwr',date1,date2,file=file,/nostr)   
     96unc=read_ncdf('lwr',date1,date2,file=file,/nostr) 
    6897help, unc 
    6998 
    … …  
    80109fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/NCEP2_flux_19890101_20090729.nc' 
    81110initncdf, fi 
    82 nce=read_ncdf("lwr", date1, date2, file=fi,/nostr) ;;  & nce=-1*nce  
     111nce=read_ncdf("lwr", date1, date2, file=fi,/nostr) ;;  & nce=-1*nce 
    83112help, nce 
    84113 
    … …  
    93122help, clark 
    94123 
    95 ;------------------------------------------------------------------------------------------------------------------------ 
    96124nn=n_elements(sitelist) 
    97125date1=date1 
    … …  
    111139 
    112140    lwr=lw & ind=where(finite(lwr)) & valid=n_elements(ind) 
    113   
     141 
    114142    if (valid ge 180. ) then begin 
    115143 
    … …  
    118146        uncr=tropflux 
    119147 
    120         extract_flux_tropflux,cor,box, $   
     148        extract_flux_tropflux,cor,box, $ 
    121149                 tropflux 
    122150        corr=tropflux 
    123      
     151 
    124152        extract_flux_tropflux,oaf,box, $ 
    125153                tropflux 
    … …  
    145173        statistics_3var_v1, lwr, uncr_lwr, corr_lwr, $ 
    146174                 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    147   
     175 
    148176        printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f6.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    149177        printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao,format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f6.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    … …  
    161189    endif        
    162190endfor 
    163 ;----------------------------------------------------------------------------------------- 
    164 close,/all   
     191close,/all 
    165192 
    166193fi_lwr_erai='/Users/pkb/work/MY_SAXO/flux_automat/lwr_2000_2007_erai.txt' 
    … …  
    362389significance_test, mean_tao, mean_clark 
    363390 
    364 ;----------------------------------------------------------- 
    365391closeps 
    366392 
    … …  
    370396return 
    371397end 
    372 ;-------------------------------------------------------------------------- 
    373398function x_site_location, site 
    374399    n1=strpos(site, 's') 
    … …  
    385410return, float(x) 
    386411end 
    387 ;-------------------------------------------------------------------------- 
    388412function y_site_location, site 
    389413    n1=strpos(site, 'e') 
    … …  
    411435return,float(y) 
    412436end 
    413  
    414 ;-------------------------------------------------------------------------- 
    415  

  • trunk/src/paper01/fig9/swr_validation_scatter_2000_2007.pro

    r41 r43  
    1 ;------------------------------------------------------------ 
    2 pro swr_validation_scatter_2000_2007,date1,date2 
     1;+ 
     2; .. _swr_validation_scatter_2000_2007.pro 
     3; 
     4; ==================================== 
     5; swr_validation_scatter_2000_2007.pro 
     6; ==================================== 
     7; 
     8; DESCRIPTION 
     9; =========== 
     10; 
     11; SEE ALSO 
     12; ======== 
     13; 
     14; :ref:`tropflux_profile.sh` 
     15; 
     16; EXAMPLES 
     17; ======== 
     18; 
     19; :: 
     20; 
     21;  IDL> @tropflux_init 
     22;  IDL> date1=++ 
     23;  IDL> date2=++ 
     24;  IDL> swr_validation_scatter_2000_2007, date1, date2 
     25; 
     26; EVOLUTIONS 
     27; ========== 
     28; 
     29; $Id$ 
     30; 
     31; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 
     32; 
     33;   * minimal header 
     34; 
     35;- 
     36pro swr_validation_scatter_2000_2007, date1, date2 
    337@common 
    4 ;------------------------------------------------------------ 
    538 
    639reinitplt, /z,/invert 
    … …  
    841 
    942openps, FILENAME = 'idl.ps' 
    10 ;------------------------------------------------------------ 
    1143;; Before running this program, you have to compile the following subrutines 
    1244;;   1.  .r read_global_sw_v50 
    … …  
    3264 
    3365;da1=20000101 & da2=20071231 
    34 nsmooth=1.    
    35 ;------------------------------------------------------------------------------------------------------------------------ 
     66nsmooth=1. 
    3667;;   This program will create the following text files with statistics of respective variables 
    37 ;------------------------------------------------------------------------------------------------------------------------ 
    3868close,/all 
    3969 
    … …  
    5989printf,6, 'x     y      cor    bias     std     rmsd      mean_tao' 
    6090 
    61 ;------------------------------------------------------------------------------------------------------------------------ 
    6291;; first reading the whole ERAI uncorrected and corrected data 
    6392 
    6493file='/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/erai_fluxes_20000101_20090801_TROP_oafluxgrid.nc' 
    6594initncdf, file 
    66 unc=read_ncdf('swr',date1,date2,file=file,/nostr) & unc=-1*unc  
     95unc=read_ncdf('swr',date1,date2,file=file,/nostr) & unc=-1*unc 
    6796help, unc 
    6897 
    … …  
    79108fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/swr_ncep2_oaflxgrid_19890101_20091231.nc' 
    80109initncdf, fi 
    81 nce=read_ncdf("swr", date1, date2, file=fi,/nostr) & nce=0.94*nce  
     110nce=read_ncdf("swr", date1, date2, file=fi,/nostr) & nce=0.94*nce 
    82111help, nce 
    83112 
    … …  
    91120sw_olr=read_ncdf("sw", date1-1, date2, file=file,/nostr) 
    92121help, sw_olr 
    93 ;------------------------------------------------------------------------------------------------------------------------ 
    94122nn=n_elements(sitelist) 
    95123date1=date1 
    … …  
    110138 
    111139    swr=sw & ind=where(finite(swr)) & valid=n_elements(ind) 
    112   
     140 
    113141    if (valid ge 180. ) then begin 
    114142 
    … …  
    117145        uncr=tropflux 
    118146 
    119         extract_flux_tropflux,cor,box, $   
     147        extract_flux_tropflux,cor,box, $ 
    120148                 tropflux 
    121149        corr=tropflux 
    122      
     150 
    123151        extract_flux_tropflux,oaf,box, $ 
    124152                tropflux 
    … …  
    133161        ncep1=tropflux 
    134162 
    135      
     163 
    136164        extract_flux_tropflux,sw_olr,box, $ 
    137165                tropflux 
    … …  
    143171        statistics_3var_v1, swr, uncr_swr, corr_swr, $ 
    144172                 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 
    145   
     173 
    146174        printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    147175        printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f4.2,3x,f7.2,3x,f4.2,3x,f5.2,3x,f6.2)' 
    … …  
    354382 
    355383 
    356 ;----------------------------------------------------------- 
    357384closeps 
    358385 
    … …  
    362389return 
    363390end 
    364 ;-------------------------------------------------------------------------- 
    365391function x_site_location, site 
    366392    n1=strpos(site, 's') 
    … …  
    377403return, float(x) 
    378404end 
    379 ;-------------------------------------------------------------------------- 
    380405function y_site_location, site 
    381406    n1=strpos(site, 'e') 
    … …  
    403428return,float(y) 
    404429end 
    405  
    406 ;-------------------------------------------------------------------------- 
    407  
Note: See TracChangeset for help on using the changeset viewer.