Changeset 99 for trunk/src


Ignore:
Timestamp:
08/30/11 16:21:34 (13 years ago)
Author:
pinsard
Message:

suppress successive semi-column lines

Location:
trunk/src
Files:
48 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/TropFlux_19890101_20091231.pro

    r98 r99  
    444444time=tt 
    445445jpt=n_elements(time) 
    446 ; 
    447446; 
    448447; Constants for flux computation 
    … …  
    506505;stop 
    507506  y=cor30a(wn,us,ts,t,Qs,Q,Rs,Rl,rain,zi,P,zu,zt,zq,ylat,jcool,jwave,twave,hwave) 
    508 ; 
    509 ; 
    510507; A few punctual missing values (coare does not converge): filled by spatial extrapolation 
    511508  tab=fltarr(jpi,jpj)+!values.f_nan 

  • trunk/src/TropFlux_NRT_ncdf.pro

    r97 r99  
    107107@cm_4data 
    108108@cm_project 
    109 ; 
    110109; 
    111110; test if ${PROJECT_OD} defined 

  • trunk/src/calc_cloud_vlat.pro

    r97 r99  
    1818;       CLD:  cloud fraction - value between 0-1. 
    1919;       Qclsky: clearsky Qsw = Q0 in Reeds notation 
    20 ; 
    2120; 
    2221; SET COEFFICIENTS FOR THE Q0 CALCULATION 

  • trunk/src/dew_temp_airt_rel_humidity_jv.pro

    r97 r99  
    99;     :param t: air temperature in degC 
    1010;     :param td: dew point temperature in degC 
    11 ; 
    1211; 
    1312; SEE ALSO 
    … …  
    2221;    IDL> td=+todo+ 
    2322;    IDL> result=dew_temp_airt_rel_humidity_jv(t, td) 
    24 ; 
    2523; 
    2624; TODO 

  • trunk/src/interp_erai_dewt_1989_2009.pro

    r97 r99  
    267267endfor 
    268268; 
    269 timein=timein 
    270269jptin=jpt 
    271270; 

  • trunk/src/interp_erai_msl_1989_2009.pro

    r97 r99  
    269269endfor 
    270270; 
    271 timein=timein 
    272271jptin=jpt 
    273272; 

  • trunk/src/interp_erai_sst_1989_2009.pro

    r98 r99  
    263263endfor 
    264264; 
    265 timein=timein 
    266265jptin=jpt 
    267266; 

  • trunk/src/interp_erai_t2m_1989_2009.pro

    r97 r99  
    183183;       Cannot allocate memory 
    184184; 
    185 ; 
    186185;     When timestep keyword is not used we can see these messages:: 
    187186; 
    … …  
    192191; 
    193192;        % Array used to subscript array contains out of range subscript: DATAIN 
    194 ; 
    195193; 
    196194; - fplod 20101222T163216Z aedon.locean-ipsl.upmc.fr (Darwin) 

  • trunk/src/interp_olr_30n30s_1989_2009.pro

    r98 r99  
    233233olrin=read_ncdf("olr",19881231, 20091231, file=fullfilename,/nostr) 
    234234; 
    235 ; 
    236235timein=time 
    237236jptin=jpt 
    … …  
    262261endfor 
    263262; 
    264 ; 
    265 timein=timein 
    266263jptin=jpt 
    267 ; 
    268264timein=timein-julday(1,1,1950,00,00) 
    269265; 

  • trunk/src/oaflux_mask_30n30s.pro

    r94 r99  
    136136ENDIF 
    137137; 
    138 ; 
    139138fullfilename_out=project_od_env+'mask_oaflux_30N30S.nc' 
    140139; in order to avoid unexpected overwritten 

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

    r97 r99  
    9595mean=reform(ff(4,*)) 
    9696; 
    97 ; 
    9897splot, mean, bias, small=[2,3,1], xtitle='Mean ERA-I SWR', ytitle='ERA-I SWR bias', $ 
    9998       title='', psym=2, charsize=1.2,  $ 

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

    r97 r99  
    321321endfor 
    322322; 
    323 ; 
    324323plt,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, $ 
    325324   title='3) SWR Mean bias - OAFlux', subtitle='', small=[1,4,3],/rempl,format='(i3)',/nocolorb 
    … …  
    350349endfor 
    351350; 
    352 ; 
    353351plt,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, $ 
    354352   title='4) SWR Mean bias - SWR_olr', subtitle='', small=[1,4,4],/rempl ;,format='(i4)' 
    … …  
    381379erase 
    382380; 
    383 ; 
    384381plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 
    385382   title='1) RMSD - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge 

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

    r98 r99  
    218218openw,6,fi_net_NRT_tropd 
    219219; 
    220 ; 
    221220printf,1, 'x     y      cor    bias     std     rmsd      mean_tao' 
    222221printf,2, 'x     y      cor    bias     std     rmsd      mean_tao' 
    … …  
    344343oplot, [-100,200], [-100,200] 
    345344; 
    346 ; 
    347345fi_net_NRT_trop=project_od_env+'net_2000_2007_trop.txt' 
    348346res=read_ascii(fi_net_NRT_trop,data_start=1) 
    … …  
    383381oplot, [-100,200], [-100,200] 
    384382; 
    385 ; 
    386383fi_net_NRT_oaflx=project_od_env+'net_2000_2007_oaflx.txt' 
    387384res=read_ascii(fi_net_NRT_oaflx,data_start=1) 
    … …  
    423420oplot, [-100,200], [-100,200] 
    424421; 
    425 ; 
    426422fi_net_NRT_ncep=project_od_env+'net_2000_2007_ncep.txt' 
    427423res=read_ascii(fi_net_NRT_ncep,data_start=1) 
    … …  
    462458oplot, [-100,200], [-100,200] 
    463459; 
    464 ; 
    465460fi_net_NRT_ncep1=project_od_env+'net_2000_2007_ncep1.txt' 
    466461res=read_ascii(fi_net_NRT_ncep1,data_start=1) 
    … …  
    501496oplot, [-100,200], [-100,200] 
    502497; 
    503 ; 
    504498fi_net_NRT_tropd=project_od_env+'net_2000_2007_tropd.txt' 
    505499res=read_ascii(fi_net_NRT_tropd,data_start=1) 

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

    r97 r99  
    5656;  IDL> as_validation_net_1994_95_v10 
    5757; 
    58 ; 
    5958; TODO 
    6059; ==== 

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

    r97 r99  
    675675cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    676676; 
    677 ; 
    678677splot, net, era, charsize=1.1, title='Net Flux - COARE Vs ERAI', $ 
    679678     xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 

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

    r97 r99  
    292292cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    293293; 
    294 ; 
    295294splot, net, era, charsize=1.1, title='SWR - Stratus Vs ERAI', $ 
    296295     xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 

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

    r97 r99  
    315315cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    316316; 
    317 ; 
    318317splot, net, era, charsize=1.1, title='Net Flux - Subdn. cs Vs ERAI', $ 
    319318     xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 

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

    r97 r99  
    384384cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    385385; 
    386 ; 
    387386splot, net, era, charsize=1.1, title='Net Flux - Subdn. ses Vs ERAI', $ 
    388387     xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 

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

    r97 r99  
    642642cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    643643; 
    644 ; 
    645644splot, net, era, charsize=1.1, title='SWR - Subdctn. sws Vs ERAI', $ 
    646645     xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 

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

    r97 r99  
    516516olr=olr(ind) 
    517517; 
    518 ; 
    519518statistics, net, era, $ 
    520519                 cor, bias, std, rmsd 
    … …  
    547546cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 
    548547; 
    549 ; 
    550548splot, net, era, charsize=1.1, title='Net Flux - WHOTS Vs ERAI', $ 
    551549     xrange=[-250,250], yrange=[-250,250], small=[2,3,1], psym=2, xmin=1, ymin=1 

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

    r97 r99  
    193193    nce_hf=float(fft(ff_ts_f,/inverse)) 
    194194; 
    195 ; 
    196195help, trop_hf, nce2_hf, era_hf, nce_hf 
    197196nsmooth=15. 

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

    r97 r99  
    6666; 
    6767; coding rules 
    68 ; 
    6968; 
    7069; complete description 

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

    r97 r99  
    8484time=jda1+dindgen(jpt) 
    8585; 
    86 ; 
    8786; FIRST READ ALL DATA FROM SITE 
    8887;   - if file exists, extract correct time axis 
    8988;   - if not, fill variable with missing values 
    90 ; 
    91 ; 
    9289; 
    9390;LHF 
    … …  
    112109ind=where(lh_q ne 1 and lh_q ne 2) 
    113110; 
    114 ; 
    115111if (ind(0) ne -1) then begin 
    116112  lh(ind)=!Values.f_nan 

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

    r97 r99  
    2424;       {lwnet} -> {read_lw} -> {lw} 
    2525;    } 
    26 ; 
    2726; 
    2827; SEE ALSO 
    … …  
    8483time=jda1+dindgen(jpt) 
    8584; 
    86 ; 
    8785; FIRST READ ALL DATA FROM SITE 
    8886;   - if file exists, extract correct time axis 
    8987;   - if not, fill variable with missing values 
    90 ; 
    9188; 
    9289; LWR 
    … …  
    105102  print, fi, ' not found' 
    106103endelse 
    107 ; 
    108 ; 
    109104; 
    110105; Replace missing values by "NaN" 

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

    r97 r99  
    8585time=jda1+dindgen(jpt) 
    8686; 
    87 ; 
    8887; FIRST READ ALL DATA FROM SITE 
    8988;   - if file exists, extract correct time axis 
    9089;   - if not, fill variable with missing values 
    91 ; 
    9290; 
    9391; SWR 
    … …  
    114112sw(ind1)=!Values.f_nan 
    115113; 
    116 ; 
    117 ; 
    118114; Replace missing values by "NaN" 
    119115; 

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

    r97 r99  
    106106; 
    107107; Give the location of mooring for validation of basic meteorological variables 
    108 ; 
    109108; 
    110109sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    … …  
    135134fi_air_ncep1=project_id_env+'air_2000_2009_ncep1_v50.txt' 
    136135openw,5,fi_air_ncep1 
    137 ; 
    138136; 
    139137printf,1, 'x     y      cor    bias     std     rmsd    mean_tao' 
    … …  
    191189         at, sw,rh,sst,wu,wv,ws,lh 
    192190; 
    193 ; 
    194191;    air=rh_to_sphum(rh,at,1008)               ;        at  ->  air temperature 
    195192    air=at                                               ;        air ->  sea surface temperature 
    … …  
    220217        tropflux 
    221218     ncep1=tropflux 
    222 ; 
    223219; 
    224220    ind=where(finite(air)) 

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

    r97 r99  
    108108; 
    109109; Give the location of mooring for validation of basic meteorological variables 
    110 ; 
    111110; 
    112111sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    … …  
    137136fi_q2m_ncep1=project_id_env+'q2m_2000_2009_ncep1_v50.txt' 
    138137openw,5,fi_q2m_ncep1 
    139 ; 
    140138; 
    141139printf,1, 'x     y      cor    bias     std     rmsd    mean_tao' 
    … …  
    190188    read_variables_v2, csite,date1,date2,nsmooth, $ 
    191189         at, sw,rh,sst,wu,wv,ws,lh 
    192 ; 
    193190; 
    194191    q2m=rh_to_sphum(rh,at,1008)               ;        at  ->  q2m temperature 
    … …  
    346343oplot, mean_tao, yfit, color=250, thick=2 
    347344; 
    348 ; 
    349345fi_q2m_ncep=project_id_env+'q2m_2000_2009_ncep_v50.txt' 
    350346res=read_ascii(fi_q2m_ncep,data_start=1) 

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

    r98 r99  
    9999time=jda1+dindgen(jpt) 
    100100; 
    101 ; 
    102101; FIRST READ ALL DATA FROM SITE 
    103102;   - if file exists, extract correct time axis 
    … …  
    163162  print, fi, ' not found' 
    164163endelse 
    165 ; 
    166164; 
    167165;MET PARAMETERS 

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

    r97 r99  
    100100; Give the location of mooring for validation of basic meteorological variables 
    101101; 
    102 ; 
    103102sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    104103           '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $ 
    … …  
    131130openw,6,fi_sst_ncep1 
    132131; 
    133 ; 
    134 ; 
    135132printf,1, 'x     y      cor    bias     std     rmsd      mean_tao' 
    136133printf,2, 'x     y      cor    bias     std     rmsd      mean_tao' 
    … …  
    193190         at, sw,rh,sst,wu,wv,ws,lh 
    194191; 
    195 ; 
    196192    ind=where(finite(sst)) 
    197193    valid=n_elements(ind) 
    … …  
    225221        tropflux 
    226222      ncep1=tropflux 
    227 ; 
    228223; 
    229224      ind=where(finite(sst)) 
    … …  
    356351oplot, mean_tao, yfit, color=250, thick=2 
    357352; 
    358 ; 
    359 ; 
    360353fi_sst_ncep=project_id_env+'sst_2000_2009_v50_ncep.txt' 
    361354res=read_ascii(fi_sst_ncep,data_start=1) 
    … …  
    391384oplot, mean_tao, yfit, color=250, thick=2 
    392385; 
    393 ; 
    394386fi_sst_tmi=project_id_env+'sst_2000_2009_v50_tmi.txt' 
    395387res=read_ascii(fi_sst_tmi,data_start=1) 

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

    r97 r99  
    9898; 
    9999; Give the location of mooring for validation of basic meteorological variables 
    100 ; 
    101100; 
    102101sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    … …  
    132131openw,6,fi_ws_ncep1 
    133132; 
    134 ; 
    135133printf,1, 'x     y      cor    bias     std     rmsd      mean_tao' 
    136134printf,2, 'x     y      cor    bias     std     rmsd      mean_tao' 
    … …  
    357355oplot, mean_tao, yfit, color=250, thick=2 
    358356; 
    359 ; 
    360 ; 
    361357fi_ws_ncep=project_id_env+'ws_2000_2009_ncep_v50.txt' 
    362358res=read_ascii(fi_ws_ncep,data_start=1) 
    … …  
    425421oplot, mean_tao, yfit, color=250, thick=2 
    426422; 
    427 ; 
    428423fi_ws_ncep1=project_id_env+'ws_2000_2009_ncep1_v50.txt' 
    429424res=read_ascii(fi_ws_ncep1,data_start=1) 

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

    r97 r99  
    182182print, xx, yy 
    183183; 
    184 ; 
    185184splot, mean, bias, small=[2,3,4],/noer, xtitle='Mean ERAI WS', ytitle='ERAI WS bias', $ 
    186185        psym=2, charsize=1.,  $ 

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

    r97 r99  
    117117ncep1=1000*read_ncdf("q2m", st, en , file=file,/nostr, box=box) 
    118118help, ncep1 
    119 ; 
    120119; 
    121120si=size(trop) 
    … …  
    170169ncep1=ncep1*mask 
    171170; 
    172 ; 
    173171ind=where(finite(param) and finite(erai) and finite(trop) and finite(oaf) and finite(ncep2) and finite(ncep1)) 
    174172x=param(ind) 
    … …  
    181179ncep1=ncep1(ind) 
    182180; 
    183 ; 
    184181statistics, param, trop, $ 
    185182                 cor, bias, std, rmsd 
    … …  
    196193abs_oaf=string(abs_diff, format='(f5.2)') 
    197194; 
    198 ; 
    199195statistics, param, erai, $ 
    200196                 cor, bias, std, rmsd 
    … …  
    204200abs_erai=string(abs_diff, format='(f5.2)') 
    205201; 
    206 ; 
    207202statistics, param, ncep2, $ 
    208203                 cor, bias, std, rmsd 
    … …  
    212207abs_ncep2=string(abs_diff, format='(f5.2)') 
    213208; 
    214 ; 
    215209statistics, param, ncep1, $ 
    216210                 cor, bias, std, rmsd 
    … …  
    219213abs_diff=total(abs(ncep1-param))/n_elements(ncep1) 
    220214abs_ncep1=string(abs_diff, format='(f5.2)') 
    221 ; 
    222215; 
    223216mio=4 
    … …  
    239232xyouts, 10, 4, abs_trop, charsize=1 
    240233; 
    241 ; 
    242234pdf, param,oaf,mio,mao,inx,mis,mas,iny, $ 
    243235    pdf,xx,yy 
    … …  
    247239xyouts, 10, 4, abs_oaf, charsize=1 
    248240; 
    249 ; 
    250241pdf, param,erai,mio,mao,inx,mis,mas,iny, $ 
    251242    pdf,xx,yy 
    … …  
    255246xyouts, 10, 4, abs_erai, charsize=1 
    256247; 
    257 ; 
    258248pdf, param,ncep2,mio,mao,inx,mis,mas,iny, $ 
    259249    pdf,xx,yy 
    … …  
    263253xyouts, 10, 4, abs_ncep2, charsize=1 
    264254; 
    265 ; 
    266255pdf, param,ncep1,mio,mao,inx,mis,mas,iny, $ 
    267256    pdf,xx,yy 

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

    r97 r99  
    122122ncep1=read_ncdf("sst", st, en , file=file,/nostr, box=box)-273.15 
    123123help, ncep1 
    124 ; 
    125124; 
    126125si=size(trop) 
    … …  
    187186noc=noc(ind) 
    188187; 
    189 ; 
    190188statistics, param, trop, $ 
    191189                 cor, bias, std, rmsd 
    … …  
    202200abs_oaf=string(abs_diff, format='(f5.2)') 
    203201; 
    204 ; 
    205202statistics, param, erai, $ 
    206203                 cor, bias, std, rmsd 
    … …  
    210207abs_erai=string(abs_diff, format='(f5.2)') 
    211208; 
    212 ; 
    213209statistics, param, ncep2, $ 
    214210                 cor, bias, std, rmsd 
    … …  
    231227abs_diff=total(abs(noc-param))/n_elements(noc) 
    232228abs_noc=string(abs_diff, format='(f5.2)') 
    233 ; 
    234229; 
    235230mio=15 
    … …  
    258253xyouts, xp, yp-1, abs_trop, charsize=1 
    259254; 
    260 ; 
    261255pdf, param,oaf,mio,mao,inx,mis,mas,iny, $ 
    262256    pdf,xx,yy 
    … …  
    267261xyouts, xp, yp-1, abs_oaf, charsize=1 
    268262; 
    269 ; 
    270263pdf, param,erai,mio,mao,inx,mis,mas,iny, $ 
    271264    pdf,xx,yy 
    … …  
    276269xyouts, xp, yp-1, abs_erai, charsize=1 
    277270; 
    278 ; 
    279271pdf, param,ncep2,mio,mao,inx,mis,mas,iny, $ 
    280272    pdf,xx,yy 
    … …  
    285277xyouts, xp, yp-1, abs_ncep2, charsize=1 
    286278; 
    287 ; 
    288279pdf, param,ncep1,mio,mao,inx,mis,mas,iny, $ 
    289280    pdf,xx,yy 

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

    r97 r99  
    4545; 
    4646;  IDL> icoads_t2m_stats_paper 
    47 ; 
    4847; 
    4948; TODO 
    … …  
    118117ncep1=read_ncdf("t2m", st, en , file=file,/nostr, box=box)-273.15 
    119118help, ncep1 
    120 ; 
    121119; 
    122120si=size(trop) 
    … …  
    171169ncep1=ncep1*mask 
    172170; 
    173 ; 
    174171ind=where(finite(param) and finite(erai) and finite(trop) and finite(oaf) and finite(ncep2) and finite(ncep1)) 
    175172x=param(ind) 
    … …  
    182179ncep1=ncep1(ind) 
    183180; 
    184 ; 
    185181statistics, param, trop, $ 
    186182                 cor, bias, std, rmsd 
    … …  
    197193abs_oaf=string(abs_diff, format='(f5.2)') 
    198194; 
    199 ; 
    200195statistics, param, erai, $ 
    201196                 cor, bias, std, rmsd 
    … …  
    205200abs_erai=string(abs_diff, format='(f5.2)') 
    206201; 
    207 ; 
    208202statistics, param, ncep2, $ 
    209203                 cor, bias, std, rmsd 
    … …  
    213207abs_ncep2=string(abs_diff, format='(f5.2)') 
    214208; 
    215 ; 
    216209statistics, param, ncep1, $ 
    217210                 cor, bias, std, rmsd 
    … …  
    220213abs_diff=total(abs(ncep1-param))/n_elements(ncep1) 
    221214abs_ncep1=string(abs_diff, format='(f5.2)') 
    222 ; 
    223215; 
    224216mio=15 
    … …  
    248240xyouts, xp, yp-1, abs_trop, charsize=1 
    249241; 
    250 ; 
    251242pdf, param,oaf,mio,mao,inx,mis,mas,iny, $ 
    252243    pdf,xx,yy 
    … …  
    257248xyouts, xp, yp-1, abs_oaf, charsize=1 
    258249; 
    259 ; 
    260250pdf, param,erai,mio,mao,inx,mis,mas,iny, $ 
    261251    pdf,xx,yy 
    … …  
    266256xyouts, xp, yp-1, abs_erai, charsize=1 
    267257; 
    268 ; 
    269258pdf, param,ncep2,mio,mao,inx,mis,mas,iny, $ 
    270259    pdf,xx,yy 
    … …  
    275264xyouts, xp, yp-1, abs_ncep2, charsize=1 
    276265; 
    277 ; 
    278266pdf, param,ncep1,mio,mao,inx,mis,mas,iny, $ 
    279267    pdf,xx,yy 

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

    r97 r99  
    4646; 
    4747;  IDL> icoads_ws_stats_paper 
    48 ; 
    4948; 
    5049; TODO 
    … …  
    188187noc=noc(ind) 
    189188; 
    190 ; 
    191189statistics, param, trop, $ 
    192190                 cor, bias, std, rmsd 
    … …  
    209207abs_diff=total(abs(erai-param))/n_elements(erai) 
    210208abs_erai=string(abs_diff, format='(f5.2)') 
    211 ; 
    212209; 
    213210statistics, param, ncep2, $ 
    … …  
    257254xyouts, xp, yp-1, abs_trop, charsize=1 
    258255; 
    259 ; 
    260256pdf, param,oaf,mio,mao,inx,mis,mas,iny, $ 
    261257    pdf,xx,yy 
    … …  
    266262xyouts, xp, yp-1, abs_oaf, charsize=1 
    267263; 
    268 ; 
    269264pdf, param,erai,mio,mao,inx,mis,mas,iny, $ 
    270265    pdf,xx,yy 
    … …  
    275270xyouts, xp, yp-1, abs_erai, charsize=1 
    276271; 
    277 ; 
    278272pdf, param,ncep2,mio,mao,inx,mis,mas,iny, $ 
    279273    pdf,xx,yy 
    … …  
    284278xyouts, xp, yp-1, abs_ncep2, charsize=1 
    285279; 
    286 ; 
    287280pdf, param,ncep1,mio,mao,inx,mis,mas,iny, $ 
    288281    pdf,xx,yy 

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

    r97 r99  
    9494; 
    9595; Give the location of mooring for validation of basic meteorological variables 
    96 ; 
    9796; 
    9897sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    … …  
    157156fi_air_ncep1=project_id_env+'air_2000_2009_ncep1_'+ocean+'.txt' 
    158157openw,5,fi_air_ncep1 
    159 ; 
    160158; 
    161159printf,1, 'x     y      cor    bias     std     rmsd    mean_tao' 
    … …  
    193191    read_variables_v2, csite,date1,date2,nsmooth, $ 
    194192         at, sw,rh,sst,wu,wv,ws,lh 
    195 ; 
    196 ; 
    197193;    air=rh_to_sphum(rh,at,1008)               ;        at  ->  air temperature 
    198194    air=at                                               ;        air ->  sea surface temperature 
    … …  
    201197                                                ;   wu,wv,ws ->  wind speed 
    202198; 
    203 ; 
    204 ; 
    205199;  extracting the corrected and uncorrected ERAI data at the locations 
    206200    nsmooth=1. 
    … …  
    225219        tropflux 
    226220     ncep1=tropflux 
    227 ; 
    228221; 
    229222    ind=where(finite(air)) 
    … …  
    352345oplot, mean_tao, yfit, color=250, thick=2 
    353346; 
    354 ; 
    355 ; 
    356347fi_air_ncep=project_id_env+'air_2000_2009_ncep2_'+ocean+'.txt' 
    357348res=read_ascii(fi_air_ncep,data_start=1) 

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

    r97 r99  
    156156openw,5,fi_q2m_ncep1 
    157157; 
    158 ; 
    159158printf,1, 'x     y      cor    bias     std     rmsd    mean_tao' 
    160159printf,2, 'x     y      cor    bias     std     rmsd    mean_tao' 
    … …  
    190189    read_variables_v2, csite,date1,date2,nsmooth, $ 
    191190         at, sw,rh,sst,wu,wv,ws,lh 
    192 ; 
    193191; 
    194192    q2m=rh_to_sphum(rh,at,1008)               ;        at  ->  q2m temperature 
    … …  
    346344oplot, mean_tao, yfit, color=250, thick=2 
    347345; 
    348 ; 
    349346fi_q2m_ncep=project_id_env+'q2m_2000_2009_ncep_'+ocean+'.txt' 
    350347res=read_ascii(fi_q2m_ncep,data_start=1) 

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

    r97 r99  
    136136openw,6,fi_sst_ncep1 
    137137; 
    138 ; 
    139 ; 
    140138printf,1, 'x     y      cor    bias     std     rmsd      mean_tao' 
    141139printf,2, 'x     y      cor    bias     std     rmsd      mean_tao' 
    … …  
    170168         at, sw,rh,sst,wu,wv,ws,lh 
    171169; 
    172 ; 
    173170    ind=where(finite(sst)) 
    174171    valid=n_elements(ind) 
    … …  
    202199        tropflux 
    203200      ncep1=tropflux 
    204 ; 
    205201; 
    206202      ind=where(finite(sst)) 
    … …  
    333329oplot, mean_tao, yfit, color=250, thick=2 
    334330; 
    335 ; 
    336 ; 
    337331fi_sst_ncep=project_id_env+'sst_2000_2009_ncep_'+ocean+'.txt' 
    338332res=read_ascii(fi_sst_ncep,data_start=1) 
    … …  
    368362oplot, mean_tao, yfit, color=250, thick=2 
    369363; 
    370 ; 
    371364fi_sst_tmi=project_id_env+'sst_2000_2009_tmi_'+ocean+'.txt' 
    372365res=read_ascii(fi_sst_tmi,data_start=1) 

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

    r97 r99  
    122122          '2n125w', '2n110w', '2n95w',  '5n155w', '5n140w', '5n125w', '5n110w', '5n95w'] 
    123123ocean='Pac_coldtongue' 
    124 ; 
    125124; 
    126125;sitelist=['0n0e', '0n10w', '0n23w', '0n35w', '10s10w', '12n23w', '12n38w', '14s32w', '15n38w', '19s34w', '20n38w', $ 
    … …  
    152151openw,6,fi_ws_ncep1 
    153152; 
    154 ; 
    155153printf,1, 'x     y      cor    bias     std     rmsd      mean_tao' 
    156154printf,2, 'x     y      cor    bias     std     rmsd      mean_tao' 
    … …  
    348346oplot, mean_tao, yfit, color=250, thick=2 
    349347; 
    350 ; 
    351 ; 
    352348fi_ws_ncep=project_id_env+'ws_2000_2009_ncep_'+ocean+'.txt' 
    353349res=read_ascii(fi_ws_ncep,data_start=1) 
    … …  
    416412oplot, mean_tao, yfit, color=250, thick=2 
    417413; 
    418 ; 
    419414fi_ws_ncep1=project_id_env+'ws_2000_2009_ncep1_'+ocean+'.txt' 
    420415res=read_ascii(fi_ws_ncep1,data_start=1) 

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

    r97 r99  
    9999; location of moorings 
    100100; 
    101 ; 
    102101sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    103102           '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $ 
    … …  
    126125fi_lhf_ncep1=project_id_env+'lhf_2000_2009_ncep1_v52.txt' 
    127126openw,5,fi_lhf_ncep1 
    128 ; 
    129127; 
    130128printf,1, 'x     y      cor    bias     std     rmsd      mean_tao' 
    … …  
    229227        printf, 5, 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)' 
    230228; 
    231 ; 
    232229    endif 
    233230endfor 
    … …  
    311308oplot, [20,200], [20,200] 
    312309; 
    313 ; 
    314310fi_lhf_oaflx=project_id_env+'lhf_2000_2009_oaflx_v52.txt' 
    315311res=read_ascii(fi_lhf_oaflx,data_start=1) 

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

    r97 r99  
    129129openw,5,fi_shf_ncep1 
    130130; 
    131 ; 
    132131printf,1, 'x     y      cor    bias     std     rmsd      mean_tao' 
    133132printf,2, 'x     y      cor    bias     std     rmsd      mean_tao' 
    … …  
    274273oplot, [-5,25], [-5,25] 
    275274; 
    276 ; 
    277275fi_shf_trop=project_id_env+'shf_2000_2009_trop_v52.txt' 
    278276res=read_ascii(fi_shf_trop,data_start=1) 
    … …  
    313311oplot, [-5,25], [-5,25] 
    314312; 
    315 ; 
    316313fi_shf_oaflx=project_id_env+'shf_2000_2009_oaflx_v52.txt' 
    317314res=read_ascii(fi_shf_oaflx,data_start=1) 

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

    r97 r99  
    102102; location of mooring 
    103103; 
    104 ; 
    105104sitelist=['5n165e','8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 
    106105           '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $ 
    … …  
    135134openw,6,fi_swr_olr 
    136135; 
    137 ; 
    138136printf,1, 'x     y      cor    bias     std     rmsd      mean_tao' 
    139137printf,2, 'x     y      cor    bias     std     rmsd      mean_tao' 
    … …  
    293291oplot, [100,300], [100,300] 
    294292; 
    295 ; 
    296293fi_swr_trop=project_id_env+'swr_2000_2007_trop_v50.txt' 
    297294res=read_ascii(fi_swr_trop,data_start=1) 
    … …  
    332329oplot, [100,300], [100,300] 
    333330; 
    334 ; 
    335331fi_swr_oaflx=project_id_env+'swr_2000_2007_oaflx_v50.txt' 
    336332res=read_ascii(fi_swr_oaflx,data_start=1) 
    … …  
    371367oplot, [100,300], [100,300] 
    372368; 
    373 ; 
    374369fi_swr_ncep=project_id_env+'swr_2000_2007_ncep_v50.txt' 
    375370res=read_ascii(fi_swr_ncep,data_start=1) 
    … …  
    410405oplot, [100,300], [100,300] 
    411406; 
    412 ; 
    413407fi_swr_ncep1=project_id_env+'swr_2000_2007_ncep1_v50.txt' 
    414408res=read_ascii(fi_swr_ncep1,data_start=1) 
    … …  
    449443oplot, [100,300], [100,300] 
    450444; 
    451 ; 
    452445fi_swr_olr=project_id_env+'swr_2000_2007_olr_v50.txt' 
    453446res=read_ascii(fi_swr_olr,data_start=1) 

  • trunk/src/plot_tropflux.pro

    r91 r99  
    1111; 
    1212; Produce +todo+ figures with PROJECT data 
    13 ; 
    1413; 
    1514; .. graphviz:: 

  • trunk/src/project_init.pro

    r97 r99  
    2323;  $ IDL 
    2424;  IDL> @project_init 
    25 ; 
    2625; 
    2726; SEE ALSO 

  • trunk/src/sst_correction_ncdf.pro

    r98 r99  
    5252; 
    5353;        } 
    54 ; 
    5554; 
    5655; SEE ALSO 

  • trunk/src/time_axe_modif_TropFlux.pro

    r97 r99  
    2222; 
    2323;        } 
    24 ; 
    2524; 
    2625; SEE ALSO 
    … …  
    8786print, ' ' 
    8887; 
    89 ; 
    9088; opening the file to modify the time axis 
    9189; 
    … …  
    10199; 
    102100; 
    103 ; 
    104101ncdf_control, cdfid ,/endef 
    105102ntt=n_elements(tt) 

  • trunk/src/time_axe_modif_tropflux_1d_1989_2010.pro

    r97 r99  
    2222; 
    2323;        } 
    24 ; 
    2524; 
    2625; SEE ALSO 
    … …  
    7574print, ' ' 
    7675; 
    77 ; 
    7876; opening the file to modify the time axis 
    7977; ouverture fichier dont on modifie le timeaxe 
    … …  
    9290;varid=ncdf_varid(cdfid,'tt') 
    9391; 
    94 ; 
    95 ; 
    9692ncdf_control, cdfid ,/endef 
    9793ntt=n_elements(tt) 

  • trunk/src/ws_tropflux_1d_to_1m_ncdf.pro

    r98 r99  
    172172endfor 
    173173; 
    174 ; 
    175174;----------------------------------------------------------------------------- 
    176175cda0=string(jul2date(time(0)),format='(i8.8)') 
    … …  
    197196    , Website:'http://www.locean-ipsl.upmc.fr/tropflux/'} 
    198197; 
    199 ; 
    200198ncfields = 'ws[longitude,latitude,time]=wsm:ws_attr; ' $ 
    201199                      + 'longitude[]=xlon:lon_attr; ' $ 
Note: See TracChangeset for help on using the changeset viewer.