- Timestamp:
- 04/11/11 18:05:51 (13 years ago)
- 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 ;- 2 34 pro fig1_no_obs_map 3 35 @common 4 ;------------------------------------------------------------5 36 6 37 reinitplt, /z,/invert … … 8 39 9 40 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 41 close,/all 12 42 … … 136 166 plots, 202.46, 22.45, psym=2 ;; WHOTS 137 167 138 ;-----------------------------------------------------------139 168 closeps 140 169 … … 144 173 return 145 174 end 146 ;-------------------------------------------------------------------------- 175 147 176 function x_site_location, site 148 177 n1=strpos(site, 's') … … 159 188 return, float(x) 160 189 end 161 ;-------------------------------------------------------------------------- 190 162 191 function y_site_location, site 163 192 n1=strpos(site, 'e') … … 185 214 return,float(y) 186 215 end 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 1 35 pro fig10_swr_correction 2 36 @common 3 ;------------------------------------------------------------4 37 reinitplt, /z,/invert 5 38 key_portrait = 1 6 39 7 40 openps, FILENAME = 'idl.ps' 8 ;------------------------------------------------------------9 41 close,/all 10 42 … … 48 80 significance_test,mean,std 49 81 50 ;;------------------------------------------------------------------------51 82 closeps 52 83 -
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 ;- 1 34 pro swr_isccp_tropflux_correlation_fig10 2 35 @common 3 ;--------------------------------------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 6 38 7 39 openps, FILENAME = 'idl.ps' 8 ;--------------------------------------------------------------------------9 40 st=19890101 & en=20071231 10 41 … … 24 55 contour=sst-273.15, contmin=24, contmax=29, contint=1.5 25 56 26 ;--------------------------------------------------------------------------27 57 closeps 28 58 fig='correlation_rmsd_olra_swra_smooth.ps' … … 32 62 stop 33 63 end 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 ;- 1 35 pro swr_isccp_tropflux_new_v1 2 36 @common 3 ;--------------------------------------------------------------------------4 37 reinitplt, /z,/invert 5 38 key_portrait = 1 6 39 7 40 openps, FILENAME = 'idl.ps' 8 ;--------------------------------------------------------------------------9 41 st=19890101 & en=20071231 10 42 … … 50 82 help, swi, olri 51 83 52 ;--------------------------------------------------------------------------53 84 caldat, time,mon,day,yea 54 85 swm=olr*0. & olrm=olr*0. … … 61 92 swa=sw_isccp-swm ; swr anomaly 62 93 olra=olr-olrm ; olr anomaly 63 ;--------------------------------------------------------------------------64 94 help, swa, olra 65 95 m=reform(sw_trop(*,*,0)) … … 81 111 endfor 82 112 endfor 83 113 84 114 cor=fltarr(jpi,jpj)+!values.f_nan & rms=fltarr(jpi,jpj)+!values.f_nan 85 115 for jj=1,jpj-2 do begin … … 92 122 93 123 save, cor, cor_olr, filename='correlation_isccp_olr.idl' 94 124 95 125 lct=64 & cs=0.9 96 126 marge1=[0,0,-5,0] … … 118 148 subt='',title='b) RMS Difference (OLRA, SWRA)',xminor=1,yminor=1 119 149 120 ;--------------------------------------------------------------------------121 150 closeps 122 151 fig='correlation_rmsd_olra_swra_smooth.ps' … … 127 156 128 157 end 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 ;- 1 34 pro swr_statistics_map_2000_2009_v50 2 35 @common 3 ;------------------------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 bias_mi=-20 & bias_ma=20 & bias_int=2 13 43 std_mi=0.7 & std_ma=1.31 & std_int=0.05 … … 90 120 for n=0,NN-1 do begin 91 121 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 126 endfor 98 127 plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $ 99 128 title='2) SWR Correlation - ERAI', subtitle='', small=[1,4,2],/noer,/rempl, marge=marge … … 116 145 for n=0,NN-1 do begin 117 146 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 151 endfor 124 152 plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $ 125 153 title='3) SWR Correlation - OAFlux', subtitle='', small=[1,4,3],/noer,/rempl, marge=marge … … 144 172 y=lat(n) 145 173 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 176 endfor 150 177 plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int, $ 151 178 title='4) SWR Correlation - SWR_olr', subtitle='', small=[1,4,4],/noer,/rempl, marge=marge1 … … 170 197 y=lat(n) 171 198 c=cor_olr(n) 172 199 173 200 plots, x,y,psym=8,symsize=1.5,color=0 174 201 endfor 175 202 176 203 erase 177 ;------------------------------------------------------------178 204 179 205 plt, -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, $ … … 197 223 for n=0,NN-1 do begin 198 224 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 229 endfor 205 230 plt,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, $ 206 231 title='2) SWR Mean bias - ERAI', subtitle='', small=[1,4,2],/rempl,format='(i3)',/nocolorb … … 225 250 y=lat(n) 226 251 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 254 endfor 255 232 256 233 257 plt,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, $ … … 253 277 y=lat(n) 254 278 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 281 endfor 282 260 283 261 284 plt,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, $ … … 281 304 y=lat(n) 282 305 c=bias_olr(n) 283 306 284 307 plots, x,y,psym=8,symsize=1.5,color=0 285 308 endfor … … 287 310 erase 288 311 289 ;------------------------------------------------------------290 312 291 313 plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ … … 311 333 y=lat(n) 312 334 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 337 endfor 317 338 plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 318 339 title='2) RMSD - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge … … 337 358 y=lat(n) 338 359 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 362 endfor 343 363 plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 344 364 title='3) RMSD - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge … … 363 383 y=lat(n) 364 384 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 387 endfor 388 370 389 plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 371 390 title='4) RMSD - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1 … … 390 409 y=lat(n) 391 410 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 413 endfor 396 414 erase 397 415 plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ … … 420 438 plots, x,y,psym=8,symsize=1.5,color=0 421 439 endfor 422 ;------------------------------------------------------------423 440 plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 424 441 title='2) STD ratop - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge … … 447 464 endfor 448 465 449 ;------------------------------------------------------------450 466 plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 451 467 title='3) STD rato - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge … … 474 490 endfor 475 491 476 ;------------------------------------------------------------477 492 plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 478 493 title='4) STD rato - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1 … … 500 515 plots, x,y,psym=8,symsize=1.5,color=0 501 516 endfor 502 ;------------------------------------------------------------503 517 closeps 504 518 fig='swr_statistics_map_2000_2009_v50.ps' 505 519 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 506 520 spawn, 'gv '+cpsdir+fig 507 ;------------------------------------------------------------508 521 return 509 522 end 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 ;- 36 pro net_flux_validation_scatter_2000_2007, date1, date2 3 37 @common 4 ;------------------------------------------------------------5 38 6 39 reinitplt, /z,/invert … … 8 41 9 42 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 43 12 44 sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ … … 27 59 28 60 ocean='net_flux' 29 nsmooth=1. 30 ;------------------------------------------------------------------------------------------------------------------------ 61 nsmooth=1. 31 62 ;; This program will create the following text files with statistics of respective variables 32 ;------------------------------------------------------------------------------------------------------------------------33 63 file_trop='/Users/pkb/data/TropFlux/TropFlux_19890101_20091231.nc' 34 64 initncdf, file_trop 35 65 shf_trop=read_ncdf("shf", date1, date2, file=file_trop,/nostr) 36 66 lhf_trop=read_ncdf("lhf", date1, date2, file=file_trop,/nostr) 37 lwr_trop=read_ncdf("lwr", date1, date2, file=file_trop,/nostr) 67 lwr_trop=read_ncdf("lwr", date1, date2, file=file_trop,/nostr) 38 68 39 69 file_trop='/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20071231_DT_v50.nc' … … 71 101 file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/shf_oafluxgrid_1985_2008.nc' 72 102 initncdf, file_oaf 73 shf=read_ncdf("shf", date1, date2, file=file_oaf,/nostr) 103 shf=read_ncdf("shf", date1, date2, file=file_oaf,/nostr) 74 104 file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/lhf_oafluxgrid_1985_2008.nc' 75 105 initncdf, file_oaf 76 lhf=read_ncdf("lhf", date1, date2, file=file_oaf,/nostr) 106 lhf=read_ncdf("lhf", date1, date2, file=file_oaf,/nostr) 77 107 file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/lwr_oafluxgrid_1985_2007.nc' 78 108 initncdf, file_oaf 79 lwr=read_ncdf("lwr", date1, date2, file=file_oaf,/nostr) 109 lwr=read_ncdf("lwr", date1, date2, file=file_oaf,/nostr) 80 110 file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/swr_oafluxgrid_1985_2007.nc' 81 111 initncdf, file_oaf … … 94 124 help, ncep1 95 125 96 ;------------------------------------------------------------------------------------------------------------------------97 126 close,/all 98 127 … … 118 147 printf,6, 'x y cor bias std rmsd mean_tao' 119 148 120 ;------------------------------------------------------------------------------------------------------------------------121 149 nn=n_elements(sitelist) 122 150 … … 133 161 if (y ge 0. and y le 30.) then y=y+360. 134 162 dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx] 135 163 136 164 read_variables_net_flux, csite,date1,date2,nsmooth, $ 137 165 net 138 166 tao=net & ind=where(finite(tao)) & valid=n_elements(ind) 139 167 140 ;; 168 ;; 141 169 if (valid ge 180. ) then begin 142 170 … … 155 183 extract_flux_tropflux, trop, box, $ 156 184 tropflux 157 tropflux_new=tropflux 185 tropflux_new=tropflux 158 186 159 187 extract_flux_tropflux, ncep1, box, $ … … 165 193 tropdf=tropflux 166 194 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) 168 196 erai_net=eraiflux(ind) & nce1_net=ncepflux1(ind) & tropd_net=tropdf(ind) 169 197 … … 383 411 oplot, [-100,200], [-100,200] 384 412 385 ;-----------------------------------------------------------386 413 closeps 387 414 … … 391 418 return 392 419 end 393 ;--------------------------------------------------------------------------394 420 function x_site_location, site 395 421 n1=strpos(site, 's') … … 406 432 return, float(x) 407 433 end 408 ;--------------------------------------------------------------------------409 434 function y_site_location, site 410 435 n1=strpos(site, 'e') … … 432 457 return,float(y) 433 458 end 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 ;- 1 35 pro as_validation_net_1994_95_v10 2 36 @common 3 ;------------------------------------------------------------4 37 reinitplt, /z,/invert 5 38 key_portrait = 1 … … 7 40 8 41 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 42 ; partie a changer 11 ;------------------------------------------------------------12 43 asbox=[61,62,15,16] 13 44 date1=19941016 & date2=19951019 … … 71 102 file="/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20091231_NRT_v50.nc" 72 103 initncdf, file 73 olr1=read_ncdf("sw", date1-1, date2, file=file,/nostr, box=asbox) ;; & olr1=grossemoyenne(olr1, "xy",/nan, box=asbox) & 104 olr1=read_ncdf("sw", date1-1, date2, file=file,/nostr, box=asbox) ;; & olr1=grossemoyenne(olr1, "xy",/nan, box=asbox) & 74 105 olr1=reform(olr1) 75 106 … … 140 171 nce1=nce1(ind) 141 172 help, nce1 142 ;------------------------------------------------------------143 173 ;ind=where(finite(net)) 144 174 ;net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce1=nce1(ind) & nce2=nce2(ind) … … 175 205 cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 176 206 177 ;------------------------------------------------------------178 207 splot, net, era, charsize=1.1, title='Net Flux - AS Vs ERAI', $ 179 208 xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 … … 207 236 208 237 stop 209 ;------------------------------------------------------------210 238 closeps 211 239 fig='as_validation_net_1994_95_v10.ps' 212 240 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 213 241 spawn, 'gv '+cpsdir+fig 214 ;------------------------------------------------------------215 242 return 216 243 end 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 ;- 1 34 pro coare_validation_net_1992_93_all_v10 2 35 @common 3 ;------------------------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 13 43 box=[155.5, 156.5, -2, -1] … … 143 173 shf4=read_ncdf('shf',st4,en4,file=file,/nostr, box=box) & shf4=grossemoyenne(shf4, "xy",/nan, box=box) & shf4=reform(shf4) 144 174 145 net1=swr1+lwr1+lhf1+shf1 & net2=swr2+lwr2+lhf2+shf2 175 net1=swr1+lwr1+lhf1+shf1 & net2=swr2+lwr2+lhf2+shf2 146 176 net3=swr3+lwr3+lhf3+shf3 & net4=swr4+lwr4+lhf4+shf4 147 177 … … 278 308 279 309 help, net, era, trop, oaf, nce1, nce2, olr 280 ;------------------------------------------------------------281 310 282 311 statistics, net, era, $ … … 310 339 cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 311 340 312 ;------------------------------------------------------------313 341 314 342 splot, net, era, charsize=1.1, title='Net Flux - COARE Vs ERAI', $ … … 343 371 344 372 stop 345 ;------------------------------------------------------------346 373 closeps 347 374 fig="coare_validation_net_1992_93_all_v11.ps" 348 375 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 349 376 spawn, 'gv '+cpsdir+fig 350 ;------------------------------------------------------------351 377 return 352 378 end 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 ;- 1 34 pro stratus_validation_net_2001_2002_all_v10 2 35 @common 3 ;------------------------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 13 43 ;; stratus in the south pacific cold tongue 2000-01 … … 122 152 net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce2=nce2(ind) & nce1=nce1(ind) 123 153 olr=olr(ind) 124 ;------------------------------------------------------------125 154 126 155 statistics, net, era, $ … … 154 183 cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 155 184 156 ;------------------------------------------------------------157 185 158 186 splot, net, era, charsize=1.1, title='SWR - Stratus Vs ERAI', $ … … 187 215 188 216 stop 189 ;------------------------------------------------------------190 217 closeps 191 218 fig="stratus_validation_net_2001_2002_all_v11.ps" 192 219 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 193 220 spawn, 'gv '+cpsdir+fig 194 ;------------------------------------------------------------195 221 return 196 222 end 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 ;- 1 34 pro subdctn_cs_validation_net_1991_93_all_v10 2 35 @common 3 ;------------------------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 13 43 ;; subduction experiment in the atlantic in 1991-93 … … 136 166 net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce2=nce2(ind) & nce1=nce1(ind) 137 167 olr=olr(ind) 138 ;------------------------------------------------------------139 168 140 169 statistics, net, era, $ … … 168 197 cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 169 198 170 ;------------------------------------------------------------171 199 172 200 splot, net, era, charsize=1.1, title='Net Flux - Subdn. cs Vs ERAI', $ … … 201 229 202 230 stop 203 ;------------------------------------------------------------204 231 closeps 205 232 fig="subdctn_cs_validation_net_1991_93_all_v10.ps" 206 233 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 207 234 spawn, 'gv '+cpsdir+fig 208 ;------------------------------------------------------------209 235 return 210 236 end 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 ;- 1 35 pro subdctn_ses_validation_net_1991_93_all_v10 2 36 @common 3 ;------------------------------------------------------------4 37 reinitplt, /z,/invert 5 38 key_portrait = 1 … … 7 40 8 41 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 42 ; partie a changer 11 ;------------------------------------------------------------12 43 13 44 ;; subduction experiment in the atlantic in 1991-93 … … 93 124 94 125 net1=swr1+lwr1+(shf1+lhf1) & net2=swr2+lwr2+(shf2+lhf2) 95 neto_1=olr1+lwr1+lhf1+shf1 & neto_2=olr2+lwr2+lhf2+shf2 126 neto_1=olr1+lwr1+lhf1+shf1 & neto_2=olr2+lwr2+lhf2+shf2 96 127 trop=[net1,net2] & olr=[neto_1,neto_2] 97 128 help, trop, olr … … 163 194 net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce2=nce2(ind) & nce1=nce1(ind) 164 195 olr=olr(ind) 165 ;------------------------------------------------------------166 196 167 197 statistics, net, era, $ … … 195 225 cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 196 226 197 ;------------------------------------------------------------198 227 199 228 splot, net, era, charsize=1.1, title='Net Flux - Subdn. ses Vs ERAI', $ … … 228 257 229 258 stop 230 ;------------------------------------------------------------231 259 closeps 232 260 fig="subdctn_ses_validation_net_1991_93_all_v10.ps" 233 261 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 234 262 spawn, 'gv '+cpsdir+fig 235 ;------------------------------------------------------------236 263 return 237 264 end 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 ;- 1 34 pro subdctn_sws_validation_net_1991_93_all_v10 2 35 @common 3 ;------------------------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 13 43 ;; subduction experiment in the atlantic in 1991-93 … … 90 120 era=[net1,net2,net3,net4] 91 121 lhf1=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. 122 shf1=0. & shf2=0. & shf3=0. & shf4=0. 123 lwr1=0. & lwr2=0. & lwr3=0. & lwr4=0. 124 swr1=0. & swr2=0. & swr3=0. & swr4=0. 95 125 net1=0. & net2=0. & net3=0. & net4=0 96 126 help, era … … 132 162 lhf1=0. & lhf2=0. & lhf3=0. & lhf4=0. 133 163 shf1=0. & shf2=0. & shf3=0. & shf4=0. 134 lwr1=0. & lwr2=0. & lwr3=0. & lwr4=0. 164 lwr1=0. & lwr2=0. & lwr3=0. & lwr4=0. 135 165 swr1=0. & swr2=0. & swr3=0. & swr4=0. 136 166 net1=0. & net2=0. & net3=0. & net4=0. … … 252 282 net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce2=nce2(ind) & nce1=nce1(ind) 253 283 olr=olr(ind) 254 ;------------------------------------------------------------255 284 256 285 statistics, net, era, $ … … 284 313 cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 285 314 286 ;------------------------------------------------------------287 315 288 316 splot, net, era, charsize=1.1, title='SWR - Subdctn. sws Vs ERAI', $ … … 317 345 318 346 stop 319 ;------------------------------------------------------------320 347 closeps 321 348 fig="subdctn_sws_validation_net_1991_93_all_v10.ps" 322 349 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 323 350 spawn, 'gv '+cpsdir+fig 324 ;------------------------------------------------------------325 351 return 326 352 end 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 ;- 1 34 pro whots_validation_net_2004_2007_all_v10 2 35 @common 3 ;------------------------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 13 43 ;; WHOTS in the NC pacific 2004-09 … … 50 80 ;; reading other data sets 51 81 52 st1=20040813 & en1=20050725 82 st1=20040813 & en1=20050725 53 83 st2=20050728 & en2=20060624 54 st3=20060627 & en3=20070628 55 ;st4=20070626 & en4=20080606 84 st3=20060627 & en3=20070628 85 ;st4=20070626 & en4=20080606 56 86 ;st5=20080605 & en5=20090715 57 87 … … 84 114 lhf1=0. & lhf2=0. & lhf3=0. 85 115 shf1=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. 116 lwr1=0. & lwr2=0. & lwr3=0. 117 swr1=0. & swr2=0. & swr3=0. 118 net1=0. & net2=0. & net3=0. 89 119 help, era 90 120 … … 115 145 116 146 net1=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 147 neto_1=olr1+lwr1+lhf1+shf1 & neto_2=olr2+lwr2+lhf2+shf2 & neto_3=olr3+lwr3+lhf3+shf3 118 148 119 149 trop=[net1,net2,net3] & olr=[neto_1,neto_2,neto_3] 120 150 121 151 lhf1=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. 152 shf1=0. & shf2=0. & shf3=0. 153 lwr1=0. & lwr2=0. & lwr3=0. 154 swr1=0. & swr2=0. & swr3=0. 155 net1=0. & net2=0. & net3=0. 126 156 neto_1=0 & neto_2=0 & neto_3=0 127 157 … … 155 185 156 186 oaf=[net1,net2,net3] 157 lhf1=0. & lhf2=0. & lhf3=0. 158 shf1=0. & shf2=0. & shf3=0. 187 lhf1=0. & lhf2=0. & lhf3=0. 188 shf1=0. & shf2=0. & shf3=0. 159 189 lwr1=0. & lwr2=0. & lwr3=0. 160 190 swr1=0. & swr2=0. & swr3=0. … … 191 221 help, nce2 192 222 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. 223 lhf1=0. & lhf2=0. & lhf3=0. 224 shf1=0. & shf2=0. & shf3=0. 225 lwr1=0. & lwr2=0. & lwr3=0. 226 swr1=0. & swr2=0. & swr3=0. 197 227 net1=0. & net2=0. & net3=0. 198 228 … … 220 250 help, nce1 221 251 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. 252 lhf1=0. & lhf2=0. & lhf3=0. 253 shf1=0. & shf2=0. & shf3=0. 254 lwr1=0. & lwr2=0. & lwr3=0. 255 swr1=0. & swr2=0. & swr3=0. 226 256 net1=0. & net2=0. & net3=0. 227 257 228 258 ind=where(finite(net)) 229 net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) 259 net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) 230 260 nce2=nce2(ind) & nce1=nce1(ind) & olr=olr(ind) 231 261 232 ;------------------------------------------------------------233 262 234 263 statistics, net, era, $ … … 262 291 cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 263 292 264 ;------------------------------------------------------------265 293 266 294 splot, net, era, charsize=1.1, title='Net Flux - WHOTS Vs ERAI', $ … … 294 322 xyouts, 0,-200, cstat_olr, charsize=1. 295 323 296 ;------------------------------------------------------------297 324 closeps 298 325 fig="WHOTS_validation_net_2004_2007_all_v10.ps" 299 326 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 300 327 spawn, 'gv '+cpsdir+fig 301 ;------------------------------------------------------------302 328 return 303 329 end 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 ;- 2 35 pro fig15_timeseries_mjo 3 36 @common 4 ;------------------------------------------------------------5 37 6 38 reinitplt, /z,/invert … … 8 40 9 41 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 42 date1=20071001 & date2=20080331 12 43 date1=19990101 & date2=20021231 … … 79 110 80 111 save, olr,trop, nce2, era, nce, filename="netflux_sctr_1999_2002.idl" 81 ;------------------------------------------------------------------------------------------------------------------------82 112 83 113 trop_hf=trop*0. & nce2_hf=trop*0. … … 130 160 pltt, nce_hf, "t",/ov1d, color=50 131 161 132 ;----------------------------------------------------------133 162 closeps 134 163 fig='fig15_timeseries_mjo.ps' … … 137 166 return 138 167 end 139 ;-------------------------------------------------------------------------- 168 140 169 function x_site_location, site 141 170 n1=strpos(site, 's') … … 152 181 return, float(x) 153 182 end 154 ;-------------------------------------------------------------------------- 183 155 184 function y_site_location, site 156 185 n1=strpos(site, 'e') … … 178 207 return,float(y) 179 208 end 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 ;- 2 34 pro fig16_timeseries_nino3 3 35 @common 4 ;------------------------------------------------------------5 36 6 37 reinitplt, /z,/invert … … 8 39 9 40 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 41 date1=19940101 & date2=20071231 12 42 box=[240,330,-5,5] … … 18 48 restore, filename="/Users/pkb/work/MY_SAXO/netflux_nino3_1989_2007.idl" 19 49 help, trop, oafl, nce2, era, nce, sst 20 ;------------------------------------------------------------------------------------------------------------------------21 50 ;; creating climatology 22 51 ny=(2007-1994)+1 … … 99 128 trop_ano=trop_ano-trop_hf & oafl_ano=oafl_ano-oafl_hf 100 129 nce2_ano=nce2_ano-nce2_hf & era_ano=era_ano-era_hf 101 nce_ano=nce_ano-nce_hf 130 nce_ano=nce_ano-nce_hf 102 131 103 132 nsmooth=30. … … 127 156 pltt, ts_smooth(nce_ano,nsmooth,/nan), "t",/ov1d, color=65 128 157 129 ;----------------------------------------------------------130 158 closeps 131 159 fig='fig16_timeseries_nino3.ps' … … 134 162 return 135 163 end 136 ;--------------------------------------------------------------------------137 164 function x_site_location, site 138 165 n1=strpos(site, 's') … … 149 176 return, float(x) 150 177 end 151 ;--------------------------------------------------------------------------152 178 function y_site_location, site 153 179 n1=strpos(site, 'e') … … 175 201 return,float(y) 176 202 end 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 ;- 1 34 pro fig2_timeline_diagram 2 35 @common … … 115 148 lwr=lw & ind=where(finite(lwr)) & valid=n_elements(ind) 116 149 if (valid ge 10) then begin 117 150 118 151 ind=where(finite(lw)) & lw(ind)=1. 119 152 ind=where(finite(lw,/nan)) & lw(ind)=0. … … 180 213 return,float(y) 181 214 end 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 5 3 ; 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 ;- 36 pro read_lh, csite, date1, date2, nsmooth, lh 37 ; 7 38 @common 8 39 … … 49 80 ;; Replace missing values by "NaN" 50 81 ; 51 tsvars=['lh'] 82 tsvars=['lh'] 52 83 vars=[tsvars] 53 84 nn=n_elements(vars) … … 59 90 60 91 nsmooth=nsmooth 61 lh=smooth(lh,nsmooth,/nan) 92 lh=smooth(lh,nsmooth,/nan) 62 93 63 ;-------------------------------------------------------------------------------------------------64 94 end 65 ;-------------------------------------------------------------------------------------------------66 95 67 96 68 ;-------------------------------------------------------------------------------------------------69 97 function time_lec, fi 70 98 tt=ncdf_lec(fi,var='time') … … 78 106 return, tt 79 107 end 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 4 3 ; 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 ;- 36 pro read_lw, csite, date1, date2, box, lw 37 ; 6 38 @common 7 39 … … 40 72 ; 41 73 ;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'] 43 75 ;depvars=['u','v','d','t','s'] 44 76 vars=['lw'] … … 54 86 lw=smooth(lw,nsmooth,/nan) 55 87 56 ;-------------------------------------------------------------------------------------------------57 88 end 58 ;-------------------------------------------------------------------------------------------------59 89 60 90 61 ;-------------------------------------------------------------------------------------------------62 91 function time_lec, fi 63 92 tt=ncdf_lec(fi,var='time') … … 71 100 return, tt 72 101 end 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 4 3 ; 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 ;- 36 pro read_sw, csite, date1, date2, box, sw 37 ; 6 38 @common 7 39 … … 48 80 ; 49 81 ;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'] 51 83 ;depvars=['u','v','d','t','s'] 52 84 vars=['sw'] … … 63 95 sw=smooth(sw,nsmooth,/nan) 64 96 65 ;-------------------------------------------------------------------------------------------------66 97 67 98 68 ;-------------------------------------------------------------------------------------------------69 99 function time_lec, fi 70 100 tt=ncdf_lec(fi,var='time') … … 78 108 return, tt 79 109 end 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 38 pro air_validation_scatter_2000_2009_v50, date1, date2 3 39 @common 4 ;------------------------------------------------------------5 40 6 41 reinitplt, /z,/invert … … 8 43 9 44 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 45 ;; Before running this program, you have to compile the following subrutines 12 46 ;; 1. .r read_era_total … … 31 65 32 66 nsmooth=1. ;; statistics are with 7 day smoothed 33 ;------------------------------------------------------------------------------------------------------------------------34 67 ;; This program will create the following text files with statistics of respective variables 35 ;------------------------------------------------------------------------------------------------------------------------36 68 close,/all 37 69 … … 54 86 printf,5, 'x y cor bias std rmsd mean_tao' 55 87 56 ;------------------------------------------------------------------------------------------------------------------------57 88 ;; first reading the whole ERAI uncorrected and corrected data 58 89 59 90 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc' 60 91 initncdf, file 61 unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 92 unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 62 93 help, unc 63 94 … … 75 106 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 76 107 initncdf, fi 77 nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 108 nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 78 109 nce=nce-273.15 79 110 help, nce … … 81 112 file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/t2m_ncep1_19890101_20091231.nc' 82 113 initncdf, file 83 nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 114 nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 84 115 nce1=nce1-273.15 85 116 help, nce1 86 117 87 ;------------------------------------------------------------------------------------------------------------------------88 118 nn=n_elements(sitelist) 89 119 date1=date1 … … 107 137 ;; lh -> latent heat flux 108 138 ;; rh -> relative humidity 109 ;; wu,wv,ws -> wind speed 110 ;; 139 ;; wu,wv,ws -> wind speed 140 ;; 111 141 112 142 … … 118 148 uncr=tropflux 119 149 120 extract_flux_tropflux,cor,box, $ 150 extract_flux_tropflux,cor,box, $ 121 151 tropflux 122 152 corr=tropflux 123 153 124 154 extract_flux_tropflux,oaf,box, $ 125 155 tropflux … … 135 165 136 166 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) 138 168 oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 139 169 mean_tao=total(air,/nan)/n_elements(ind) 140 170 141 171 statistics_3var_v1, air, uncr_air, corr_air, $ 142 172 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 143 173 144 174 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)' 145 175 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)' … … 155 185 156 186 endfor 157 ;----------------------------------------------------------------------------------------- 158 close,/all 187 close,/all 159 188 160 189 fi_air_erai='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_v50.txt' … … 294 323 oplot, mean_tao, yfit, color=250, thick=2 295 324 296 ;-----------------------------------------------------------297 325 closeps 298 326 … … 302 330 return 303 331 end 304 ;--------------------------------------------------------------------------305 332 function x_site_location, site 306 333 n1=strpos(site, 's') … … 317 344 return, float(x) 318 345 end 319 ;--------------------------------------------------------------------------320 346 function y_site_location, site 321 347 n1=strpos(site, 'e') … … 343 369 return,float(y) 344 370 end 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 ;- 36 pro q2m_validation_scatter_2000_2009_v50, date1, date2 3 37 @common 4 ;------------------------------------------------------------5 38 6 39 reinitplt, /z,/invert … … 8 41 9 42 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 43 ;; Before running this program, you have to compile the following subrutines 12 44 ;; 1. .r read_era_total … … 31 63 32 64 nsmooth=1. ;; statistics are with 7 day smoothed 33 ;------------------------------------------------------------------------------------------------------------------------34 65 ;; This program will create the following text files with statistics of respective variables 35 ;------------------------------------------------------------------------------------------------------------------------36 66 close,/all 37 67 … … 54 84 printf,5, 'x y cor bias std rmsd mean_tao' 55 85 56 ;------------------------------------------------------------------------------------------------------------------------57 86 ;; first reading the whole ERAI uncorrected and corrected data 58 87 … … 64 93 file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc' 65 94 initncdf, file 66 cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 95 cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 67 96 help, cor 68 97 … … 74 103 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sphum_ncep2_oafluxgrid_19890101_20091231.nc' 75 104 initncdf, fi 76 nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 105 nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 77 106 nce=nce*1000 78 107 help, nce … … 84 113 help, nce1 85 114 86 ;------------------------------------------------------------------------------------------------------------------------87 115 nn=n_elements(sitelist) 88 116 date1=date1 … … 106 134 ;; lh -> latent heat flux 107 135 ;; rh -> relative humidity 108 ;; wu,wv,ws -> wind speed 109 ;; 136 ;; wu,wv,ws -> wind speed 137 ;; 110 138 111 139 … … 117 145 uncr=tropflux 118 146 119 extract_flux_tropflux,cor,box, $ 147 extract_flux_tropflux,cor,box, $ 120 148 tropflux 121 149 corr=tropflux 122 150 123 151 extract_flux_tropflux,oaf,box, $ 124 152 tropflux … … 133 161 ncep1=tropflux 134 162 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) 137 165 mean_tao=total(q2m,/nan)/n_elements(ind) 138 166 139 167 statistics_3var_v1, q2m, uncr_q2m, corr_q2m, $ 140 168 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 141 169 142 170 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)' 143 171 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)' … … 153 181 154 182 endfor 155 ;----------------------------------------------------------------------------------------- 156 close,/all 183 close,/all 157 184 158 185 fi_q2m_erai='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_erai_v50.txt' … … 291 318 oplot, mean_tao, yfit, color=250, thick=2 292 319 293 ;-----------------------------------------------------------294 320 closeps 295 321 … … 299 325 return 300 326 end 301 ;--------------------------------------------------------------------------302 327 function x_site_location, site 303 328 n1=strpos(site, 's') … … 314 339 return, float(x) 315 340 end 316 ;--------------------------------------------------------------------------317 341 function y_site_location, site 318 342 n1=strpos(site, 'e') … … 340 364 return,float(y) 341 365 end 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 ;- 2 36 pro read_variables_v2, csite,date1,date2,nsmooth, $ 3 37 at, sw,rh,sst,wu,wv,ws, lh 4 38 5 39 ; 6 ;-------------------------------------------------------------------------------------------------7 40 @common 8 41 … … 75 108 lon(tt0(ind))=lon0(ind) 76 109 endif 77 endif 110 endif 78 111 79 112 … … 122 155 ind_at=where(at_q ne 1 and at_q ne 2) & ind_ws=where(ws_q ne 1 and ws_q ne 2) 123 156 ind_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 157 ind_lh=where(lh_q ne 1 and lh_q ne 2) 158 126 159 ind=union(ind_at, union(ind_rh, union(ind_ws, union(ind_lh, ind_sst)))) 127 160 … … 135 168 ;; Replace missing values by "NaN" 136 169 ; 137 tsvars=['at','sw','rh','sst','wu','wv','ws','lat','lon'] 170 tsvars=['at','sw','rh','sst','wu','wv','ws','lat','lon'] 138 171 vars=[tsvars] 139 172 nn=n_elements(vars) … … 150 183 ws=smooth(ws,nsmooth,/nan) & sw=smooth(sw,nsmooth,/nan) 151 184 152 ;-------------------------------------------------------------------------------------------------153 185 end 154 ;-------------------------------------------------------------------------------------------------155 186 156 187 157 ;-------------------------------------------------------------------------------------------------158 188 function time_lec, fi 159 189 tt=ncdf_lec(fi,var='time') … … 167 197 return, tt 168 198 end 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 ;- 36 pro sst_validation_scatter_2000_2009_v50, date1, date2 3 37 @common 4 ;------------------------------------------------------------5 38 6 39 reinitplt, /z,/invert … … 8 41 9 42 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 43 ;; Before running this program, you have to compile the following subrutines 12 44 ;; 1. .r read_era_total … … 31 63 32 64 nsmooth=1. ;; statistics are with 7 day smoothed 33 ;------------------------------------------------------------------------------------------------------------------------34 65 ;; This program will create the following text files with statistics of respective variables 35 ;------------------------------------------------------------------------------------------------------------------------36 66 close,/all 37 67 … … 58 88 printf,6, 'x y cor bias std rmsd mean_tao' 59 89 60 ;------------------------------------------------------------------------------------------------------------------------61 90 ;; first reading the whole ERAI uncorrected and corrected data 62 91 63 92 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc' 64 93 initncdf, file 65 unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15 94 unc=read_ncdf('sst',date1,date2,file=file,/nostr) & unc=unc-273.15 66 95 help, unc 67 96 … … 78 107 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sst_ncep2_oafluxgrid_19890101_20091231.nc' 79 108 initncdf, fi 80 nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 109 nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 81 110 nce=nce-273.15 82 111 help, nce … … 93 122 help, nce1 94 123 95 ;------------------------------------------------------------------------------------------------------------------------96 124 nn=n_elements(sitelist) 97 125 date1=date1 … … 122 150 uncr=tropflux 123 151 124 extract_flux_tropflux,cor,box, $ 152 extract_flux_tropflux,cor,box, $ 125 153 tropflux 126 154 corr=tropflux 127 155 128 156 extract_flux_tropflux,oaf,box, $ 129 157 tropflux … … 143 171 144 172 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) 146 174 oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 147 175 mean_tao=total(sst,/nan)/n_elements(ind) 148 176 149 177 statistics_3var_v1, sst, uncr_sst, corr_sst, $ 150 178 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 151 179 152 180 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)' 153 181 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)' … … 166 194 167 195 endfor 168 ;----------------------------------------------------------------------------------------- 169 close,/all 196 close,/all 170 197 171 198 fi_sst_erai='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_v50_erai.txt' … … 333 360 oplot, mean_tao, yfit, color=250, thick=2 334 361 335 ;-----------------------------------------------------------336 362 closeps 337 363 … … 341 367 return 342 368 end 343 ;--------------------------------------------------------------------------344 369 function x_site_location, site 345 370 n1=strpos(site, 's') … … 356 381 return, float(x) 357 382 end 358 ;--------------------------------------------------------------------------359 383 function y_site_location, site 360 384 n1=strpos(site, 'e') … … 382 406 return,float(y) 383 407 end 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 1 38 pro statistics_3var_v1, tao,var1,var2, $ 2 39 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 3 40 4 41 5 42 @common 6 ;--------------------------------------------7 ;--------------------------------------------8 43 x=tao & y=var1 & z=var2 9 44 … … 12 47 13 48 if (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 15 50 endif 16 51 … … 42 77 43 78 end 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 ;- 2 37 pro ws_validation_scatter_2000_2009_v50,date1,date2 3 38 @common 4 ;------------------------------------------------------------5 39 6 40 reinitplt, /z,/invert … … 8 42 9 43 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 44 ;; Before running this program, you have to compile the following subrutines 12 45 ;; 1. .r read_era_total … … 32 65 da1=10000101 & da2=10081231 33 66 nsmooth=1. ;; statistics are with 7 day smoothed 34 ;------------------------------------------------------------------------------------------------------------------------35 67 ;; This program will create the following text files with statistics of respective variables 36 ;------------------------------------------------------------------------------------------------------------------------37 68 close,/all 38 69 … … 58 89 printf,6, 'x y cor bias std rmsd mean_tao' 59 90 60 ;------------------------------------------------------------------------------------------------------------------------61 91 ;; first reading the whole ERAI uncorrected and corrected data 62 92 63 93 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc' 64 94 initncdf, 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) 95 u=read_ncdf('u10',date1,date2,file=file,/nostr) 96 v=read_ncdf('v10',date1,date2,file=file,/nostr) 97 unc=sqrt(u*u+v*v) 68 98 help, unc 69 99 70 100 file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc' 71 101 initncdf, file 72 cor=read_ncdf('ws',date1,date2,file=file,/nostr) 102 cor=read_ncdf('ws',date1,date2,file=file,/nostr) 73 103 help, cor 74 104 … … 80 110 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/uwind_ncep2_oafluxgrid_19890101_20091231.nc' 81 111 initncdf, fi 82 u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 112 u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 83 113 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/vwind_ncep2_oafluxgrid_19890101_20091231.nc' 84 114 initncdf, fi 85 v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 115 v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 86 116 nce=sqrt(u*u+v*v) 87 117 help, nce … … 104 134 help, nce1 105 135 106 ;------------------------------------------------------------------------------------------------------------------------107 136 nn=n_elements(sitelist) 108 137 date1=date1 … … 130 159 uncr=tropflux 131 160 132 extract_flux_tropflux,cor,box, $ 161 extract_flux_tropflux,cor,box, $ 133 162 tropflux 134 163 corr=tropflux 135 164 136 165 extract_flux_tropflux,oaf,box, $ 137 166 tropflux … … 150 179 ncep1=tropflux 151 180 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) 153 182 oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 154 183 155 184 mean_tao=total(ws)/n_elements(ws) 156 185 157 186 statistics_3var_v1, ws, uncr_ws, corr_ws, $ 158 187 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 159 188 160 189 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)' 161 190 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)' … … 172 201 173 202 endfor 174 ;----------------------------------------------------------------------------------------- 175 close,/all 203 close,/all 176 204 fi_ws_erai='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_erai_v50.txt' 177 205 res=read_ascii(fi_ws_erai,data_start=1) … … 338 366 oplot, mean_tao, yfit, color=250, thick=2 339 367 340 ;-----------------------------------------------------------341 368 closeps 342 369 … … 346 373 return 347 374 end 348 ;--------------------------------------------------------------------------349 375 function x_site_location, site 350 376 n1=strpos(site, 's') … … 361 387 return, float(x) 362 388 end 363 ;--------------------------------------------------------------------------364 389 function y_site_location, site 365 390 n1=strpos(site, 'e') … … 387 412 return,float(y) 388 413 end 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 1 35 pro fig4_met_var_correction_scatter 2 36 @common 3 ;------------------------------------------------------------4 37 reinitplt, /z,/invert 5 38 key_portrait = 1 6 39 7 40 openps, FILENAME = 'idl.ps' 8 ;------------------------------------------------------------9 41 fi='/Users/pkb/work/MY_SAXO/flux_automat/sst_correction_final.txt' 10 42 res=read_ascii(fi,data_start=0) … … 34 66 significance_test_99,mean,bias 35 67 36 ;-----------------------37 68 fi='/Users/pkb/work/MY_SAXO/flux_automat/t2m_correction_final.txt' 38 69 res=read_ascii(fi,data_start=0) … … 63 94 significance_test_99,mean,bias 64 95 65 ;-----------------------66 96 fi='/Users/pkb/work/MY_SAXO/flux_automat/q2m_correction_final.txt' 67 97 res=read_ascii(fi,data_start=0) … … 90 120 ;; significant test 91 121 significance_test_99,mean,bias 92 ;-----------------------93 122 fi='/Users/pkb/work/MY_SAXO/flux_automat/ws_correction_final.txt' 94 123 res=read_ascii(fi,data_start=0) … … 118 147 ;; significant test 119 148 significance_test_99,mean,bias 120 ;-----------------------------------------------------------------------------------------121 149 122 150 closeps … … 127 155 128 156 end 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 ;- 37 pro significance_test_99, x, y, sig_value 3 38 4 39 X=x & Y=y … … 21 56 print, 'Statistics are not valid at 99% significance level' 22 57 endelse 23 58 24 59 end 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 ;- 1 35 pro fig5_gustiness_scatter 2 36 @common 3 ;-----------------------------------------------------------------------4 37 reinitplt, /z,/invert 5 38 key_portrait = 1 6 39 7 40 openps, FILENAME = 'idl.ps' 8 ;-----------------------------------------------------------------------9 41 fi_ws='/Users/pkb/work/MY_SAXO/flux_automat/wind_gustiness_bias.txt' 10 42 res=read_ascii(fi_ws,data_start=1) … … 22 54 print, correlate(erai_sst,gust) 23 55 24 ;-----------------------------------------------------------------------25 56 splot, erai_wind, gust, small=[2,3,1], title='a) ERAI Wind Vs Gustiness', $ 26 57 xtitle='ERAI wind', ytitle='Cronin gustiness', charsize=1.2, psym=2 , $ … … 44 75 xyouts, 21,-0.4, 'y = 0.175712 x -3.17303', charsize=1.2 45 76 print, a,b 46 ;-----------------------------------------------------------------------47 77 closeps 48 78 fig='fig5_gustiness_scatter.ps' … … 50 80 spawn, 'gv '+cpsdir+fig 51 81 end 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 ;- 1 34 pro icoads_q2m_stats_paper 2 35 @common 3 ;-----------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 marge=[-2,-2, -4,2] 13 43 st=19890101 & en=20091231 … … 178 208 xyouts, 10, 4, abs_ncep1, charsize=1 179 209 180 ;------------------------------------------------------------181 210 closeps 182 211 fig="icoads_q2m_stats_paper.ps" 183 212 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 184 213 spawn, 'gv '+cpsdir+fig 185 ;------------------------------------------------------------186 214 return 187 215 end -
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 ;- 1 34 pro icoads_sst_stats_paper 2 35 @common 3 ;-----------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 marge=[-2,-2, -4,2] 13 43 st=19890101 & en=20091231 … … 95 125 trop=trop*mask & oaf=oaf*mask & erai=erai*mask 96 126 ncep2=ncep2*mask & ncep1=ncep1*mask 97 noc=noc*mask 127 noc=noc*mask 98 128 99 129 ind=where(finite(param) and finite(erai) and finite(trop) and finite(oaf) and finite(ncep2) and finite(ncep1)) … … 205 235 xyouts, xp, yp-1, abs_noc, charsize=1 206 236 207 ;------------------------------------------------------------208 237 closeps 209 238 fig="icoads_sst_stats_paper.ps" 210 239 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 211 240 spawn, 'gv '+cpsdir+fig 212 ;------------------------------------------------------------213 241 return 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 ;- 1 34 pro icoads_t2m_stats_paper 2 35 @common 3 ;-----------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 marge=[-2,-2, -4,2] 13 43 st=19890101 & en=20091231 … … 186 216 xyouts, xp, yp-1, abs_ncep1, charsize=1 187 217 188 ;------------------------------------------------------------189 218 closeps 190 219 fig="icoads_t2m_stats_paper.ps" 191 220 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 192 221 spawn, 'gv '+cpsdir+fig 193 ;------------------------------------------------------------194 222 return 195 223 end -
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 ;- 1 34 pro icoads_ws_stats_paper 2 35 @common 3 ;-----------------------------------------------4 36 reinitplt, /z,/invert 5 37 key_portrait = 1 … … 7 39 8 40 openps, FILENAME = 'idl.ps' 9 ;------------------------------------------------------------10 41 ; partie a changer 11 ;------------------------------------------------------------12 42 marge=[-2,-2, -4,2] 13 43 st=19890101 & en=20091231 … … 202 232 xyouts, xp, yp-1, abs_noc, charsize=1 203 233 204 ;------------------------------------------------------------205 234 closeps 206 235 fig="icoads_ws_stats_paper.ps" 207 236 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 208 237 spawn, 'gv '+cpsdir+fig 209 ;------------------------------------------------------------210 238 return 211 239 end -
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 ; 1 37 pro statistics, moor, erai, $ 2 38 cor, bias, std, rmsd 3 39 @common 4 ;--------------------------------------------5 40 ;; mean1= mean of mooring variable std1=std of mooring variable 6 41 ;; mean2= mean of erai variable std2=std of erai variable 7 ;--------------------------------------------8 42 x=moor & y=erai 9 43 ind1=where(finite(x,/nan)) & ind2=where(finite(y,/nan)) … … 26 60 rmsd=sqrt(total((x-y)*(x-y))/n_elements(x)) 27 61 end 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 ;- 2 38 pro air_validation_scatter_2000_2009,date1,date2 3 39 @common 4 ;------------------------------------------------------------5 40 6 41 reinitplt, /z,/invert … … 8 43 9 44 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 ;; Before running this program, you have to compile the following subrutines12 ;; 1. .r read_era_total13 ;; 2. .r read_variables_v214 ;; 3. .r read_era_box15 45 16 46 ;; Give the location of mooring for validation of basic meteorological variables … … 65 95 66 96 nsmooth=1. ;; statistics are with 7 day smoothed 67 ;------------------------------------------------------------------------------------------------------------------------68 97 ;; This program will create the following text files with statistics of respective variables 69 ;------------------------------------------------------------------------------------------------------------------------70 98 close,/all 71 99 … … 88 116 printf,5, 'x y cor bias std rmsd mean_tao' 89 117 90 ;------------------------------------------------------------------------------------------------------------------------91 118 ;; first reading the whole ERAI uncorrected and corrected data 92 119 93 120 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc' 94 121 initncdf, file 95 unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 122 unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 96 123 help, unc 97 124 … … 109 136 ;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 110 137 ;initncdf, fi 111 ;nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 138 ;nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) 112 139 ;nce=nce-273.15 113 140 ;help, nce … … 115 142 ;file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/t2m_ncep1_19890101_20091231.nc' 116 143 ;initncdf, file 117 ;nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 144 ;nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) 118 145 ;nce1=nce1-273.15 119 146 ;help, nce1 … … 123 150 help, unc, cor, oaf, nce, nce1 124 151 125 ;------------------------------------------------------------------------------------------------------------------------126 152 nn=n_elements(sitelist) 127 153 date1=date1 … … 145 171 ;; lh -> latent heat flux 146 172 ;; rh -> relative humidity 147 ;; wu,wv,ws -> wind speed 148 ;; 173 ;; wu,wv,ws -> wind speed 174 ;; 149 175 150 176 … … 156 182 uncr=tropflux 157 183 158 extract_flux_tropflux,cor,box, $ 184 extract_flux_tropflux,cor,box, $ 159 185 tropflux 160 186 corr=tropflux 161 187 162 188 extract_flux_tropflux,oaf,box, $ 163 189 tropflux … … 173 199 174 200 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) 176 202 oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind) 177 203 mean_tao=total(air,/nan)/n_elements(ind) 178 204 179 205 statistics_3var_v1, air, uncr_air, corr_air, $ 180 206 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 181 207 182 208 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)' 183 209 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)' … … 193 219 194 220 endfor 195 ;----------------------------------------------------------------------------------------- 196 close,/all 221 close,/all 197 222 198 223 fi_air_erai='/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_'+ocean+'.txt' … … 332 357 oplot, mean_tao, yfit, color=250, thick=2 333 358 334 ;-----------------------------------------------------------335 359 closeps 336 360 … … 340 364 return 341 365 end 342 ;--------------------------------------------------------------------------343 366 function x_site_location, site 344 367 n1=strpos(site, 's') … … 355 378 return, float(x) 356 379 end 357 ;--------------------------------------------------------------------------358 380 function y_site_location, site 359 381 n1=strpos(site, 'e') … … 381 403 return,float(y) 382 404 end 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 37 pro q2m_validation_scatter_2000_2009_basin, date1, date2 3 38 @common 4 ;------------------------------------------------------------5 39 6 40 reinitplt, /z,/invert … … 8 42 9 43 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 44 ;; Before running this program, you have to compile the following subrutines 12 45 ;; 1. .r read_era_total … … 65 98 66 99 nsmooth=1. ;; statistics are with 7 day smoothed 67 ;------------------------------------------------------------------------------------------------------------------------68 100 ;; This program will create the following text files with statistics of respective variables 69 ;------------------------------------------------------------------------------------------------------------------------70 101 close,/all 71 102 … … 88 119 printf,5, 'x y cor bias std rmsd mean_tao' 89 120 90 ;------------------------------------------------------------------------------------------------------------------------91 121 ;; first reading the whole ERAI uncorrected and corrected data 92 122 … … 98 128 ;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc' 99 129 ;initncdf, file 100 ;cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 130 ;cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 101 131 ;help, cor 102 132 … … 108 138 ;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sphum_ncep2_oafluxgrid_19890101_20091231.nc' 109 139 ;initncdf, fi 110 ;nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 140 ;nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) 111 141 ;nce=nce*1000 112 142 ;help, nce … … 122 152 help, unc, cor, oaf, nce, nce1 123 153 124 ;------------------------------------------------------------------------------------------------------------------------125 154 nn=n_elements(sitelist) 126 155 date1=date1 … … 144 173 ;; lh -> latent heat flux 145 174 ;; rh -> relative humidity 146 ;; wu,wv,ws -> wind speed 147 ;; 175 ;; wu,wv,ws -> wind speed 176 ;; 148 177 149 178 … … 155 184 uncr=tropflux 156 185 157 extract_flux_tropflux,cor,box, $ 186 extract_flux_tropflux,cor,box, $ 158 187 tropflux 159 188 corr=tropflux 160 189 161 190 extract_flux_tropflux,oaf,box, $ 162 191 tropflux … … 171 200 ncep1=tropflux 172 201 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) 175 204 mean_tao=total(q2m,/nan)/n_elements(ind) 176 205 177 206 statistics_3var_v1, q2m, uncr_q2m, corr_q2m, $ 178 207 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 179 208 180 209 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)' 181 210 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)' … … 191 220 192 221 endfor 193 ;----------------------------------------------------------------------------------------- 194 close,/all 222 close,/all 195 223 196 224 fi_q2m_erai='/Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_erai_'+ocean+'.txt' … … 329 357 oplot, mean_tao, yfit, color=250, thick=2 330 358 331 ;-----------------------------------------------------------332 359 closeps 333 360 … … 337 364 return 338 365 end 339 ;--------------------------------------------------------------------------340 366 function x_site_location, site 341 367 n1=strpos(site, 's') … … 352 378 return, float(x) 353 379 end 354 ;--------------------------------------------------------------------------355 380 function y_site_location, site 356 381 n1=strpos(site, 'e') … … 378 403 return,float(y) 379 404 end 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 ;- 2 37 pro sst_validation_scatter_2000_2009_basin,date1,date2 3 38 @common 4 ;------------------------------------------------------------5 39 6 40 reinitplt, /z,/invert … … 8 42 9 43 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 44 ;; Before running this program, you have to compile the following subrutines 12 45 ;; 1. .r read_era_total … … 43 76 '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w'] 44 77 45 nsmooth=1. 46 ;------------------------------------------------------------------------------------------------------------------------ 78 nsmooth=1. 47 79 ;; This program will create the following text files with statistics of respective variables 48 ;------------------------------------------------------------------------------------------------------------------------49 80 close,/all 50 81 … … 71 102 printf,6, 'x y cor bias std rmsd mean_tao' 72 103 73 ;------------------------------------------------------------------------------------------------------------------------74 104 ;; first reading the whole ERAI uncorrected and corrected data 75 105 76 106 ;file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc' 77 107 ;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 79 109 ;help, unc 80 110 … … 91 121 ;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sst_ncep2_oafluxgrid_19890101_20091231.nc' 92 122 ;initncdf, fi 93 ;nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 123 ;nce=read_ncdf("sst", date1-1, date2, file=fi,/nostr) 94 124 ;nce=nce-273.15 95 125 ;help, nce … … 112 142 help, unc, cor, oaf, nce, sst_tmi, nce1 113 143 ;; data from 20000101, 20090931 114 ;------------------------------------------------------------------------------------------------------------------------115 144 nn=n_elements(sitelist) 116 145 date1=date1 … … 141 170 uncr=tropflux 142 171 143 extract_flux_tropflux,cor,box, $ 172 extract_flux_tropflux,cor,box, $ 144 173 tropflux 145 174 corr=tropflux 146 175 147 176 extract_flux_tropflux,oaf,box, $ 148 177 tropflux … … 162 191 163 192 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) 165 194 oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 166 195 mean_tao=total(sst,/nan)/n_elements(ind) 167 196 168 197 statistics_3var_v1, sst, uncr_sst, corr_sst, $ 169 198 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 170 199 171 200 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)' 172 201 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)' … … 185 214 186 215 endfor 187 ;----------------------------------------------------------------------------------------- 188 close,/all 216 close,/all 189 217 190 218 fi_sst_erai='/Users/pkb/work/MY_SAXO/flux_automat/sst_2000_2009_erai_'+ocean+'.txt' … … 352 380 oplot, mean_tao, yfit, color=250, thick=2 353 381 354 ;-----------------------------------------------------------355 382 closeps 356 383 … … 360 387 return 361 388 end 362 ;--------------------------------------------------------------------------363 389 function x_site_location, site 364 390 n1=strpos(site, 's') … … 375 401 return, float(x) 376 402 end 377 ;--------------------------------------------------------------------------378 403 function y_site_location, site 379 404 n1=strpos(site, 'e') … … 401 426 return,float(y) 402 427 end 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 ;- 2 37 pro ws_validation_scatter_2000_2009_basin,date1,date2 3 38 @common 4 ;------------------------------------------------------------5 39 6 40 reinitplt, /z,/invert … … 8 42 9 43 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 44 ;; Before running this program, you have to compile the following subrutines 12 45 ;; 1. .r read_era_total … … 49 82 sitelist=['8s125w', '8s110w', '8s95w','5s155w', '5s140w', '5s125w', '5s110w', '5s95w','2s155w', '2s140w', $ 50 83 '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'] 52 85 ocean='Pac_coldtongue' 53 86 … … 65 98 66 99 nsmooth=1 67 ;------------------------------------------------------------------------------------------------------------------------68 100 ;; This program will create the following text files with statistics of respective variables 69 ;------------------------------------------------------------------------------------------------------------------------70 101 close,/all 71 102 … … 91 122 printf,6, 'x y cor bias std rmsd mean_tao' 92 123 93 ;------------------------------------------------------------------------------------------------------------------------94 124 ;; first reading the whole ERAI uncorrected and corrected data 95 125 ;; date1=20000101 & date2=20090931 … … 97 127 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc' 98 128 initncdf, 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) 129 u=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) 102 132 ;help, unc 103 133 104 134 ;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc' 105 135 ;initncdf, file 106 ;cor=read_ncdf('ws',date1,date2,file=file,/nostr) 136 ;cor=read_ncdf('ws',date1,date2,file=file,/nostr) 107 137 ;help, cor 108 138 … … 114 144 ;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/uwind_ncep2_oafluxgrid_19890101_20091231.nc' 115 145 ;initncdf, fi 116 ;u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 146 ;u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 117 147 ;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/vwind_ncep2_oafluxgrid_19890101_20091231.nc' 118 148 ;initncdf, fi 119 ;v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 149 ;v=read_ncdf("v", date1-1, date2, file=fi,/nostr) 120 150 ;nce=sqrt(u*u+v*v) 121 151 ;help, nce … … 142 172 help, unc, cor, oaf, nce, ws_tmi, nce1 143 173 144 ;------------------------------------------------------------------------------------------------------------------------145 174 nn=n_elements(sitelist) 146 175 date1=date1 … … 168 197 uncr=tropflux 169 198 170 extract_flux_tropflux,cor,box, $ 199 extract_flux_tropflux,cor,box, $ 171 200 tropflux 172 201 corr=tropflux 173 202 174 203 extract_flux_tropflux,oaf,box, $ 175 204 tropflux … … 188 217 ncep1=tropflux 189 218 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) 191 220 oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 192 221 193 222 mean_tao=total(ws)/n_elements(ws) 194 223 195 224 statistics_3var_v1, ws, uncr_ws, corr_ws, $ 196 225 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 197 226 198 227 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)' 199 228 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)' … … 210 239 211 240 endfor 212 ;----------------------------------------------------------------------------------------- 213 close,/all 241 close,/all 214 242 fi_ws_erai='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_erai_'+ocean+'.txt' 215 243 res=read_ascii(fi_ws_erai,data_start=1) … … 376 404 oplot, mean_tao, yfit, color=250, thick=2 377 405 378 ;-----------------------------------------------------------379 406 closeps 380 407 … … 384 411 return 385 412 end 386 ;--------------------------------------------------------------------------387 413 function x_site_location, site 388 414 n1=strpos(site, 's') … … 399 425 return, float(x) 400 426 end 401 ;--------------------------------------------------------------------------402 427 function y_site_location, site 403 428 n1=strpos(site, 'e') … … 425 450 return,float(y) 426 451 end 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 ;- 2 36 pro lhf_validation_scatter_2000_2009,date1,date2 3 37 @common 4 ;------------------------------------------------------------5 38 6 39 reinitplt, /z,/invert … … 8 41 9 42 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 43 ;; Before running this program, you have to compile the following subrutines 12 44 … … 31 63 ocean='global' 32 64 33 ;------------------------------------------------------------------------------------------------------------------------34 65 ;; This program will create the following text files with statistics of respective variables 35 ;------------------------------------------------------------------------------------------------------------------------36 66 close,/all 37 67 … … 53 83 printf,4, 'x y cor bias std rmsd mean_tao' 54 84 55 ;------------------------------------------------------------------------------------------------------------------------56 85 ;; first reading the whole ERAI uncorrected and corrected data 57 86 58 87 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_lhf_19890101_20091231_oafluxgrid.nc' 59 88 initncdf, file 60 unc=read_ncdf('lhf',date1,date2,file=file,/nostr) 89 unc=read_ncdf('lhf',date1,date2,file=file,/nostr) 61 90 help, unc 62 91 … … 74 103 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/lhf_ncep2_oafluxgrid_19890101_20091231.nc' 75 104 initncdf, fi 76 nce=read_ncdf("lhf", date1, date2, file=fi,/nostr) 105 nce=read_ncdf("lhf", date1, date2, file=fi,/nostr) 77 106 help, nce 78 107 79 108 file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/fluxe_ncep1_19890101_20091231.nc' 80 109 initncdf, file 81 nce1=-1*read_ncdf("lhf", date1, date2, file=file,/nostr) 110 nce1=-1*read_ncdf("lhf", date1, date2, file=file,/nostr) 82 111 help, nce1 83 112 84 113 nsmooth=1. 85 ;------------------------------------------------------------------------------------------------------------------------86 114 nn=n_elements(sitelist) 87 115 date1=date1 … … 101 129 102 130 lhf=lh & ind=where(finite(lhf)) & valid=n_elements(ind) 103 131 104 132 if (valid ge 180. ) then begin 105 133 … … 108 136 uncr=tropflux 109 137 110 extract_flux_tropflux,cor,box, $ 138 extract_flux_tropflux,cor,box, $ 111 139 tropflux 112 140 corr=tropflux 113 141 114 142 extract_flux_tropflux,oaf,box, $ 115 143 tropflux … … 130 158 statistics_3var_v1, lhf, uncr_lhf, corr_lhf, $ 131 159 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 132 160 133 161 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)' 134 162 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)' … … 170 198 171 199 splot, 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 173 201 xyouts, 50,25, cstat, charsize=1. 174 202 xyouts, 50,10, 'cor bias std rmsd', charsize=1. … … 305 333 oplot, [20,200], [20,200] 306 334 307 ;-----------------------------------------------------------308 335 closeps 309 336 … … 313 340 return 314 341 end 315 ;--------------------------------------------------------------------------316 342 function x_site_location, site 317 343 n1=strpos(site, 's') … … 328 354 return, float(x) 329 355 end 330 ;--------------------------------------------------------------------------331 356 function y_site_location, site 332 357 n1=strpos(site, 'e') … … 354 379 return,float(y) 355 380 end 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 ;- 2 36 pro shf_validation_scatter_2000_2009,date1,date2 3 37 @common 4 ;------------------------------------------------------------5 38 6 39 reinitplt, /z,/invert … … 8 41 9 42 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 43 ;; Before running this program, you have to compile the following subrutines 12 44 ;; 1. .r read_variables_sh … … 30 62 ocean='global' 31 63 32 nsmooth=1. 33 ;------------------------------------------------------------------------------------------------------------------------ 64 nsmooth=1. 34 65 ;; This program will create the following text files with statistics of respective variables 35 ;------------------------------------------------------------------------------------------------------------------------36 66 close,/all 37 67 … … 55 85 printf,5, 'x y cor bias std rmsd mean_tao' 56 86 57 ;------------------------------------------------------------------------------------------------------------------------58 87 ;; first reading the whole ERAI uncorrected and corrected data 59 88 60 89 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_shf_19890101_20091231_oafluxgrid.nc' 61 90 initncdf, file 62 unc=read_ncdf('shf',date1,date2,file=file,/nostr) 91 unc=read_ncdf('shf',date1,date2,file=file,/nostr) 63 92 help, unc 64 93 … … 76 105 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/shf_ncep2_oafluxgrid_19890101_20091231.nc' 77 106 initncdf, fi 78 nce=read_ncdf("shf", date1, date2, file=fi,/nostr) 107 nce=read_ncdf("shf", date1, date2, file=fi,/nostr) 79 108 help, nce 80 109 … … 84 113 help, nce1 85 114 86 ;------------------------------------------------------------------------------------------------------------------------87 115 nn=n_elements(sitelist) 88 116 date1=date1 … … 102 130 at, sw,rh,sst,wu,wv,ws,sh 103 131 shf=sh & ind=where(finite(shf)) & valid=n_elements(ind) 104 132 105 133 if (valid ge 180. ) then begin 106 134 … … 109 137 uncr=tropflux 110 138 111 extract_flux_tropflux,cor,box, $ 139 extract_flux_tropflux,cor,box, $ 112 140 tropflux 113 141 corr=tropflux 114 142 115 143 extract_flux_tropflux,oaf,box, $ 116 144 tropflux … … 130 158 statistics_3var_v1, shf, uncr_shf, corr_shf, $ 131 159 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 132 160 133 161 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)' 134 162 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)' … … 305 333 oplot, [-5,25], [-5,25] 306 334 307 ;-----------------------------------------------------------308 335 closeps 309 336 … … 313 340 return 314 341 end 315 ;--------------------------------------------------------------------------316 342 function x_site_location, site 317 343 n1=strpos(site, 's') … … 328 354 return, float(x) 329 355 end 330 ;--------------------------------------------------------------------------331 356 function y_site_location, site 332 357 n1=strpos(site, 'e') … … 354 379 return,float(y) 355 380 end 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 ;- 2 36 pro lwr_validation_scatter_2000_2007,date1,date2 3 37 @common 4 ;------------------------------------------------------------5 38 6 39 reinitplt, /z,/invert … … 8 41 9 42 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 43 ;; Before running this program, you have to compile the following subrutines 12 44 ;; 1. .r read_lw_global_v2 … … 34 66 35 67 da1=20000101 & da2=20071231 36 nsmooth=1. 37 ;------------------------------------------------------------------------------------------------------------------------ 68 nsmooth=1. 38 69 ;; This program will create the following text files with statistics of respective variables 39 ;------------------------------------------------------------------------------------------------------------------------40 70 close,/all 41 71 … … 60 90 printf,6, 'x y cor bias std rmsd' 61 91 62 ;------------------------------------------------------------------------------------------------------------------------63 92 ;; first reading the whole ERAI uncorrected and corrected data 64 93 65 94 file='/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/erai_fluxes_19930101_20090801_TROP_oafluxgrid.nc' 66 95 initncdf, file 67 unc=read_ncdf('lwr',date1,date2,file=file,/nostr) 96 unc=read_ncdf('lwr',date1,date2,file=file,/nostr) 68 97 help, unc 69 98 … … 80 109 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/NCEP2_flux_19890101_20090729.nc' 81 110 initncdf, fi 82 nce=read_ncdf("lwr", date1, date2, file=fi,/nostr) ;; & nce=-1*nce 111 nce=read_ncdf("lwr", date1, date2, file=fi,/nostr) ;; & nce=-1*nce 83 112 help, nce 84 113 … … 93 122 help, clark 94 123 95 ;------------------------------------------------------------------------------------------------------------------------96 124 nn=n_elements(sitelist) 97 125 date1=date1 … … 111 139 112 140 lwr=lw & ind=where(finite(lwr)) & valid=n_elements(ind) 113 141 114 142 if (valid ge 180. ) then begin 115 143 … … 118 146 uncr=tropflux 119 147 120 extract_flux_tropflux,cor,box, $ 148 extract_flux_tropflux,cor,box, $ 121 149 tropflux 122 150 corr=tropflux 123 151 124 152 extract_flux_tropflux,oaf,box, $ 125 153 tropflux … … 145 173 statistics_3var_v1, lwr, uncr_lwr, corr_lwr, $ 146 174 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 147 175 148 176 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)' 149 177 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)' … … 161 189 endif 162 190 endfor 163 ;----------------------------------------------------------------------------------------- 164 close,/all 191 close,/all 165 192 166 193 fi_lwr_erai='/Users/pkb/work/MY_SAXO/flux_automat/lwr_2000_2007_erai.txt' … … 362 389 significance_test, mean_tao, mean_clark 363 390 364 ;-----------------------------------------------------------365 391 closeps 366 392 … … 370 396 return 371 397 end 372 ;--------------------------------------------------------------------------373 398 function x_site_location, site 374 399 n1=strpos(site, 's') … … 385 410 return, float(x) 386 411 end 387 ;--------------------------------------------------------------------------388 412 function y_site_location, site 389 413 n1=strpos(site, 'e') … … 411 435 return,float(y) 412 436 end 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 ;- 36 pro swr_validation_scatter_2000_2007, date1, date2 3 37 @common 4 ;------------------------------------------------------------5 38 6 39 reinitplt, /z,/invert … … 8 41 9 42 openps, FILENAME = 'idl.ps' 10 ;------------------------------------------------------------11 43 ;; Before running this program, you have to compile the following subrutines 12 44 ;; 1. .r read_global_sw_v50 … … 32 64 33 65 ;da1=20000101 & da2=20071231 34 nsmooth=1. 35 ;------------------------------------------------------------------------------------------------------------------------ 66 nsmooth=1. 36 67 ;; This program will create the following text files with statistics of respective variables 37 ;------------------------------------------------------------------------------------------------------------------------38 68 close,/all 39 69 … … 59 89 printf,6, 'x y cor bias std rmsd mean_tao' 60 90 61 ;------------------------------------------------------------------------------------------------------------------------62 91 ;; first reading the whole ERAI uncorrected and corrected data 63 92 64 93 file='/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/erai_fluxes_20000101_20090801_TROP_oafluxgrid.nc' 65 94 initncdf, file 66 unc=read_ncdf('swr',date1,date2,file=file,/nostr) & unc=-1*unc 95 unc=read_ncdf('swr',date1,date2,file=file,/nostr) & unc=-1*unc 67 96 help, unc 68 97 … … 79 108 fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/swr_ncep2_oaflxgrid_19890101_20091231.nc' 80 109 initncdf, fi 81 nce=read_ncdf("swr", date1, date2, file=fi,/nostr) & nce=0.94*nce 110 nce=read_ncdf("swr", date1, date2, file=fi,/nostr) & nce=0.94*nce 82 111 help, nce 83 112 … … 91 120 sw_olr=read_ncdf("sw", date1-1, date2, file=file,/nostr) 92 121 help, sw_olr 93 ;------------------------------------------------------------------------------------------------------------------------94 122 nn=n_elements(sitelist) 95 123 date1=date1 … … 110 138 111 139 swr=sw & ind=where(finite(swr)) & valid=n_elements(ind) 112 140 113 141 if (valid ge 180. ) then begin 114 142 … … 117 145 uncr=tropflux 118 146 119 extract_flux_tropflux,cor,box, $ 147 extract_flux_tropflux,cor,box, $ 120 148 tropflux 121 149 corr=tropflux 122 150 123 151 extract_flux_tropflux,oaf,box, $ 124 152 tropflux … … 133 161 ncep1=tropflux 134 162 135 163 136 164 extract_flux_tropflux,sw_olr,box, $ 137 165 tropflux … … 143 171 statistics_3var_v1, swr, uncr_swr, corr_swr, $ 144 172 cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2 145 173 146 174 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)' 147 175 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)' … … 354 382 355 383 356 ;-----------------------------------------------------------357 384 closeps 358 385 … … 362 389 return 363 390 end 364 ;--------------------------------------------------------------------------365 391 function x_site_location, site 366 392 n1=strpos(site, 's') … … 377 403 return, float(x) 378 404 end 379 ;--------------------------------------------------------------------------380 405 function y_site_location, site 381 406 n1=strpos(site, 'e') … … 403 428 return,float(y) 404 429 end 405 406 ;--------------------------------------------------------------------------407
Note: See TracChangeset
for help on using the changeset viewer.