- Timestamp:
- 08/25/11 13:02:20 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/paper01/fig3/ws_validation_scatter_2000_2009_v50.pro
r85 r94 113 113 ocean='global' 114 114 115 da1=10000101 & da2=10081231 115 da1=10000101 116 da2=10081231 116 117 nsmooth=1. ;; statistics are with 7 day smoothed 117 118 ;; This program will create the following text files with statistics of respective variables … … 188 189 ;; reading data from mooring 189 190 190 site=sitelist(n) & csite=site 191 site=sitelist(n) 192 csite=site 191 193 print, csite 192 194 x=x_site_location(site) 193 195 y=y_site_location(site) 194 196 if (y ge 0. and y le 30.) then y=y+360. 195 dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx] 197 dx=0.5 198 dy=0.5 199 box=[y-dy, y+dy, x-dx, x+dx] 196 200 read_variables_v2, csite,date1,date2,nsmooth, $ 197 201 at, sw,rh,sst,wu,wv,ws, lh … … 226 230 ncep1=tropflux 227 231 228 ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind) & corr_ws=corr(ind) 229 oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind) 232 ind=where(finite(ws)) 233 ws=ws(ind) 234 uncr_ws=uncr(ind) 235 corr_ws=corr(ind) 236 oafl=oafl(ind) 237 ncep=ncep(ind) 238 tmi=tmi(ind) 239 ncep1=ncep1(ind) 230 240 231 241 mean_tao=total(ws)/n_elements(ws) … … 254 264 lat=reform(ff(0,*)) 255 265 lon=reform(ff(1,*)) 256 cor_era=reform(ff(2,*)) & cor_erai=total(cor_era)/n_elements(cor_era) 257 bias_era=reform(ff(3,*)) & bias_erai=total(bias_era)/n_elements(bias_era) 258 std_era=reform(ff(4,*)) & std_erai=total(std_era)/n_elements(std_era) 259 rmsd_era=reform(ff(5,*)) & rmsd_erai=total(rmsd_era)/n_elements(rmsd_era) 260 mean_tao=reform(ff(6,*)) & mean_erai=bias_era+mean_tao 266 cor_era=reform(ff(2,*)) 267 cor_erai=total(cor_era)/n_elements(cor_era) 268 bias_era=reform(ff(3,*)) 269 bias_erai=total(bias_era)/n_elements(bias_era) 270 std_era=reform(ff(4,*)) 271 std_erai=total(std_era)/n_elements(std_era) 272 rmsd_era=reform(ff(5,*)) 273 rmsd_erai=total(rmsd_era)/n_elements(rmsd_era) 274 mean_tao=reform(ff(6,*)) 275 mean_erai=bias_era+mean_tao 261 276 262 277 print, '' … … 272 287 oplot, [2,10], [2,10] 273 288 ab=linfit(mean_tao, mean_erai,yfit=yfit) 274 a=float(ab(0)) & b=float(ab(1)) 289 a=float(ab(0)) 290 b=float(ab(1)) 275 291 oplot, mean_tao, yfit, color=250, thick=2 276 292 … … 280 296 lat=reform(ff(0,*)) 281 297 lon=reform(ff(1,*)) 282 cor_tro=reform(ff(2,*)) & cor_trop=total(cor_tro)/n_elements(cor_tro) 283 bias_tro=reform(ff(3,*)) & bias_trop=total(bias_tro)/n_elements(bias_tro) 284 std_tro=reform(ff(4,*)) & std_trop=total(std_tro)/n_elements(std_tro) 285 rmsd_tro=reform(ff(5,*)) & rmsd_trop=total(rmsd_tro)/n_elements(rmsd_tro) 286 mean_tao=reform(ff(6,*)) & mean_trop=bias_tro+mean_tao 298 cor_tro=reform(ff(2,*)) 299 cor_trop=total(cor_tro)/n_elements(cor_tro) 300 bias_tro=reform(ff(3,*)) 301 bias_trop=total(bias_tro)/n_elements(bias_tro) 302 std_tro=reform(ff(4,*)) 303 std_trop=total(std_tro)/n_elements(std_tro) 304 rmsd_tro=reform(ff(5,*)) 305 rmsd_trop=total(rmsd_tro)/n_elements(rmsd_tro) 306 mean_tao=reform(ff(6,*)) 307 mean_trop=bias_tro+mean_tao 287 308 288 309 print, '' … … 299 320 300 321 ab=linfit(mean_tao, mean_trop,yfit=yfit) 301 a=float(ab(0)) & b=float(ab(1)) 322 a=float(ab(0)) 323 b=float(ab(1)) 302 324 oplot, mean_tao, yfit, color=250, thick=2 303 325 … … 307 329 lat=reform(ff(0,*)) 308 330 lon=reform(ff(1,*)) 309 cor_oaf=reform(ff(2,*)) & cor_oafl=total(cor_oaf)/n_elements(cor_oaf) 310 bias_oaf=reform(ff(3,*)) & bias_oafl=total(bias_oaf)/n_elements(bias_oaf) 311 std_oaf=reform(ff(4,*)) & std_oafl=total(std_oaf)/n_elements(std_oaf) 312 rmsd_oaf=reform(ff(5,*)) & rmsd_oafl=total(rmsd_oaf)/n_elements(rmsd_oaf) 313 mean_tao=reform(ff(6,*)) & mean_oafl=bias_oaf+mean_tao 331 cor_oaf=reform(ff(2,*)) 332 cor_oafl=total(cor_oaf)/n_elements(cor_oaf) 333 bias_oaf=reform(ff(3,*)) 334 bias_oafl=total(bias_oaf)/n_elements(bias_oaf) 335 std_oaf=reform(ff(4,*)) 336 std_oafl=total(std_oaf)/n_elements(std_oaf) 337 rmsd_oaf=reform(ff(5,*)) 338 rmsd_oafl=total(rmsd_oaf)/n_elements(rmsd_oaf) 339 mean_tao=reform(ff(6,*)) 340 mean_oafl=bias_oaf+mean_tao 314 341 315 342 print, '' … … 326 353 327 354 ab=linfit(mean_tao, mean_oafl,yfit=yfit) 328 a=float(ab(0)) & b=float(ab(1)) 355 a=float(ab(0)) 356 b=float(ab(1)) 329 357 oplot, mean_tao, yfit, color=250, thick=2 330 358 … … 336 364 lat=reform(ff(0,*)) 337 365 lon=reform(ff(1,*)) 338 cor_nce=reform(ff(2,*)) & cor_ncep=total(cor_nce)/n_elements(cor_nce) 339 bias_nce=reform(ff(3,*)) & bias_ncep=total(bias_nce)/n_elements(bias_nce) 340 std_nce=reform(ff(4,*)) & std_ncep=total(std_nce)/n_elements(std_nce) 341 rmsd_nce=reform(ff(5,*)) & rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce) 342 mean_tao=reform(ff(6,*)) & mean_ncep=bias_nce+mean_tao 366 cor_nce=reform(ff(2,*)) 367 cor_ncep=total(cor_nce)/n_elements(cor_nce) 368 bias_nce=reform(ff(3,*)) 369 bias_ncep=total(bias_nce)/n_elements(bias_nce) 370 std_nce=reform(ff(4,*)) 371 std_ncep=total(std_nce)/n_elements(std_nce) 372 rmsd_nce=reform(ff(5,*)) 373 rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce) 374 mean_tao=reform(ff(6,*)) 375 mean_ncep=bias_nce+mean_tao 343 376 344 377 print, '' … … 355 388 356 389 ab=linfit(mean_tao, mean_ncep,yfit=yfit) 357 a=float(ab(0)) & b=float(ab(1)) 390 a=float(ab(0)) 391 b=float(ab(1)) 358 392 oplot, mean_tao, yfit, color=250, thick=2 359 393 … … 363 397 lat=reform(ff(0,*)) 364 398 lon=reform(ff(1,*)) 365 cor_tm=reform(ff(2,*)) & cor_tmi=total(cor_tm)/n_elements(cor_tm) 366 bias_tm=reform(ff(3,*)) & bias_tmi=total(bias_tm)/n_elements(bias_tm) 367 std_tm=reform(ff(4,*)) & std_tmi=total(std_tm)/n_elements(std_tm) 368 rmsd_tm=reform(ff(5,*)) & rmsd_tmi=total(rmsd_tm)/n_elements(rmsd_tm) 369 mean_tao=reform(ff(6,*)) & mean_tmi=bias_tm+mean_tao 399 cor_tm=reform(ff(2,*)) 400 cor_tmi=total(cor_tm)/n_elements(cor_tm) 401 bias_tm=reform(ff(3,*)) 402 bias_tmi=total(bias_tm)/n_elements(bias_tm) 403 std_tm=reform(ff(4,*)) 404 std_tmi=total(std_tm)/n_elements(std_tm) 405 rmsd_tm=reform(ff(5,*)) 406 rmsd_tmi=total(rmsd_tm)/n_elements(rmsd_tm) 407 mean_tao=reform(ff(6,*)) 408 mean_tmi=bias_tm+mean_tao 370 409 371 410 print, '' … … 382 421 383 422 ab=linfit(mean_tao, mean_tmi,yfit=yfit) 384 a=float(ab(0)) & b=float(ab(1)) 423 a=float(ab(0)) 424 b=float(ab(1)) 385 425 oplot, mean_tao, yfit, color=250, thick=2 386 426 … … 391 431 lat=reform(ff(0,*)) 392 432 lon=reform(ff(1,*)) 393 cor_nce=reform(ff(2,*)) & cor_ncep=total(cor_nce)/n_elements(cor_nce) 394 bias_nce=reform(ff(3,*)) & bias_ncep=total(bias_nce)/n_elements(bias_nce) 395 std_nce=reform(ff(4,*)) & std_ncep=total(std_nce)/n_elements(std_nce) 396 rmsd_nce=reform(ff(5,*)) & rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce) 397 mean_tao=reform(ff(6,*)) & mean_ncep=bias_nce+mean_tao 433 cor_nce=reform(ff(2,*)) 434 cor_ncep=total(cor_nce)/n_elements(cor_nce) 435 bias_nce=reform(ff(3,*)) 436 bias_ncep=total(bias_nce)/n_elements(bias_nce) 437 std_nce=reform(ff(4,*)) 438 std_ncep=total(std_nce)/n_elements(std_nce) 439 rmsd_nce=reform(ff(5,*)) 440 rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce) 441 mean_tao=reform(ff(6,*)) 442 mean_ncep=bias_nce+mean_tao 398 443 399 444 print, '' … … 410 455 411 456 ab=linfit(mean_tao, mean_ncep,yfit=yfit) 412 a=float(ab(0)) & b=float(ab(1)) 457 a=float(ab(0)) 458 b=float(ab(1)) 413 459 oplot, mean_tao, yfit, color=250, thick=2 414 460
Note: See TracChangeset
for help on using the changeset viewer.