- Timestamp:
- 08/30/11 14:38:29 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/paper01/fig14/subdctn_cs_validation_net_1991_93_all_v10.pro
r94 r97 92 92 key_portrait = 1 93 93 coefpalit=.9 94 94 ; 95 95 openps, FILENAME = project_od_env+"subdctn_cs_validation_net_1991_93_all_v10.ps" 96 96 ; partie a changer 97 98 ; ;subduction experiment in the atlantic in 1991-9397 ; 98 ; subduction experiment in the atlantic in 1991-93 99 99 box=[330.5, 331.5, 25., 26.] 100 101 ; ;first leg102 100 ; 101 ; first leg 102 ; 103 103 fi=project_id_env+"met_data_subduct_cs_91.txt" 104 104 res=read_ascii(fi,data_start=1) … … 109 109 lwr_91=reform(ff(8,*)) 110 110 net_91=[swr_91+lwr_91+shf_91+lhf_91] 111 112 ; ;start data - 06/02/92 and end date 02/06/92 then from 11/10/92 to 31/12/92113 ; ;total no of daily records -114 ; ;second leg115 111 ; 112 ; start data - 06/02/92 and end date 02/06/92 then from 11/10/92 to 31/12/92 113 ; total no of daily records - 114 ; second leg 115 ; 116 116 fi=project_id_env+"met_data_subduct_cs_92.txt" 117 117 res=read_ascii(fi,data_start=1) … … 122 122 lwr_92=reform(ff(8,*)) 123 123 net_92=[swr_92+lwr_92+shf_92+lhf_92] 124 125 ; ;second leg126 124 ; 125 ; second leg 126 ; 127 127 fi=project_id_env+"met_data_subduct_cs_93.txt" 128 128 res=read_ascii(fi,data_start=1) … … 133 133 lwr_93=reform(ff(8,*)) 134 134 net_93=[swr_93+lwr_93+shf_93+lhf_93] 135 135 ; 136 136 net=[net_91, net_92, net_93] 137 137 help, net 138 139 ; ;reading other data sets140 141 st1=19910623 138 ; 139 ; reading other data sets 140 ; 141 st1=19910623 142 142 en1=19930616 143 143 ; 144 144 file=project_id_env+"erai_swr_19910101_20091231_oafluxgrid.nc" 145 145 initncdf, file 146 swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 147 swr1=grossemoyenne(swr1, "xy", box=box,/nan) 146 swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 147 swr1=grossemoyenne(swr1, "xy", box=box,/nan) 148 148 swr1=reform(swr1) 149 149 file=project_id_env+'erai_lhf_19890101_20091231_oafluxgrid.nc' 150 150 initncdf, file 151 lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 152 lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 151 lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 152 lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 153 153 lhf1=reform(lhf1) 154 154 file=project_id_env+'erai_shf_19890101_20091231_oafluxgrid.nc' 155 155 initncdf, file 156 shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 157 shf1=grossemoyenne(shf1, "xy", box=box,/nan) 156 shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 157 shf1=grossemoyenne(shf1, "xy", box=box,/nan) 158 158 shf1=reform(shf1) 159 159 file=project_id_env+'erai_lwr_19890101_20091231_oafluxgrid.nc' 160 160 initncdf, file 161 lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box) 162 lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 163 lwr1=reform(lwr1) 164 161 lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box) 162 lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 163 lwr1=reform(lwr1) 164 ; 165 165 era=swr1+lwr1+lhf1+shf1 166 swr1=0 167 lwr1=0 168 shf1=0 166 swr1=0 167 lwr1=0 168 shf1=0 169 169 lhf1=0 170 170 help, era 171 171 ; 172 172 file=project_id_env+"TropFlux_19890101_20091231.nc" 173 173 initncdf, file 174 lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 175 lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 176 lhf1=reform(lhf1) 177 shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 178 shf1=grossemoyenne(shf1, "xy", box=box,/nan) 179 shf1=reform(shf1) 180 lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) 181 lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 182 lwr1=reform(lwr1) 183 swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 184 swr1=grossemoyenne(swr1, "xy", box=box,/nan) 174 lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 175 lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 176 lhf1=reform(lhf1) 177 shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 178 shf1=grossemoyenne(shf1, "xy", box=box,/nan) 179 shf1=reform(shf1) 180 lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) 181 lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 182 lwr1=reform(lwr1) 183 swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 184 swr1=grossemoyenne(swr1, "xy", box=box,/nan) 185 185 swr1=reform(swr1) 186 186 trop=swr1+lwr1+lhf1+shf1 187 187 ; 188 188 file=project_id_env+"TropFlux_swr_19890101_20091231_NRT_v50.nc" 189 189 initncdf, file 190 olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box) 191 olr1=grossemoyenne(olr1, "xy",/nan, box=box) 190 olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box) 191 olr1=grossemoyenne(olr1, "xy",/nan, box=box) 192 192 olr1=reform(olr1) 193 193 ; 194 194 olr=olr1+lwr1+shf1+lhf1 195 195 help, trop, olr 196 swr1=0 197 lwr1=0 198 shf1=0 199 lhf1=0 196 swr1=0 197 lwr1=0 198 shf1=0 199 lhf1=0 200 200 olr1=0 201 201 ; 202 202 file=project_id_env+'swr_oafluxgrid_1985_2007.nc' 203 203 initncdf, file 204 swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box) 205 swr1=grossemoyenne(swr1, "xy",/nan, box=box) 204 swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box) 205 swr1=grossemoyenne(swr1, "xy",/nan, box=box) 206 206 swr1=reform(swr1) 207 207 file=project_id_env+'lwr_oafluxgrid_1985_2007.nc' 208 208 initncdf, file 209 lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box) 210 lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) 209 lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box) 210 lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) 211 211 lwr1=reform(lwr1) 212 212 file=project_id_env+'lhf_oafluxgrid_1985_2009.nc' 213 213 initncdf, file 214 lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box) 215 lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) 214 lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box) 215 lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) 216 216 lhf1=reform(lhf1) 217 217 file=project_id_env+'shf_oafluxgrid_1985_2009.nc' 218 218 initncdf, file 219 shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box) 220 shf1=grossemoyenne(shf1, "xy",/nan, box=box) 221 shf1=reform(shf1) 222 219 shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box) 220 shf1=grossemoyenne(shf1, "xy",/nan, box=box) 221 shf1=reform(shf1) 222 ; 223 223 oaf=swr1+lwr1+lhf1+shf1 224 224 help, oaf 225 swr1=0 226 lwr1=0 227 shf1=0 228 lhf1=0 225 swr1=0 226 lwr1=0 227 shf1=0 228 lhf1=0 229 229 olr1=0 230 230 ; 231 231 fi=project_id_env+'swr_ncep2_oaflxgrid_19890101_20091231.nc' 232 232 initncdf, fi 233 swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box) 234 swr1=grossemoyenne(swr1, "xy",/nan, box=box) 233 swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box) 234 swr1=grossemoyenne(swr1, "xy",/nan, box=box) 235 235 swr1=reform(swr1) 236 236 fi=project_id_env+"NCEP2_flux_19890101_20090729.nc" 237 237 initncdf, fi 238 lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box) 239 lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) 238 lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box) 239 lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) 240 240 lwr1=reform(lwr1) 241 241 fi=project_id_env+'lhf_ncep2_oafluxgrid_19890101_20091231.nc' 242 242 initncdf, fi 243 lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box) 244 lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) 243 lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box) 244 lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) 245 245 lhf1=reform(lhf1) 246 246 fi=project_id_env+'shf_ncep2_oafluxgrid_19890101_20091231.nc' 247 247 initncdf, fi 248 shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box) 249 shf1=grossemoyenne(shf1, "xy",/nan, box=box) 250 shf1=reform(shf1) 251 248 shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box) 249 shf1=grossemoyenne(shf1, "xy",/nan, box=box) 250 shf1=reform(shf1) 251 ; 252 252 nce2=swr1+lwr1+lhf1+shf1 253 253 help, nce2 254 swr1=0 255 lwr1=0 256 shf1=0 257 lhf1=0 254 swr1=0 255 lwr1=0 256 shf1=0 257 lhf1=0 258 258 olr1=0 259 259 ; 260 260 file=project_id_env+"ncep1_flux_19890101_20091231.nc" 261 261 initncdf, file 262 263 lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 264 lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 265 lhf1=reform(lhf1) 266 shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 267 shf1=grossemoyenne(shf1, "xy", box=box,/nan) 268 shf1=reform(shf1) 269 lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) 270 lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 271 lwr1=reform(lwr1) 272 swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 273 swr1=grossemoyenne(swr1, "xy", box=box,/nan) 262 ; 263 lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 264 lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 265 lhf1=reform(lhf1) 266 shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 267 shf1=grossemoyenne(shf1, "xy", box=box,/nan) 268 shf1=reform(shf1) 269 lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) 270 lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 271 lwr1=reform(lwr1) 272 swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 273 swr1=grossemoyenne(swr1, "xy", box=box,/nan) 274 274 swr1=reform(swr1) 275 275 nce1=swr1+(lwr1+shf1+lhf1) 276 276 help, nce1 277 277 ; 278 278 ind=where(finite(net)) 279 net=net(ind) 280 era=era(ind) 281 trop=trop(ind) 282 oaf=oaf(ind) 283 nce2=nce2(ind) 279 net=net(ind) 280 era=era(ind) 281 trop=trop(ind) 282 oaf=oaf(ind) 283 nce2=nce2(ind) 284 284 nce1=nce1(ind) 285 285 olr=olr(ind) 286 286 ; 287 287 statistics, net, era, $ 288 288 cor, bias, std, rmsd 289 289 print, cor, bias, std, rmsd 290 290 cstat_era=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 291 291 ; 292 292 statistics, net, trop, $ 293 293 cor, bias, std, rmsd 294 294 print, cor, bias, std, rmsd 295 295 cstat_trop=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 296 296 ; 297 297 statistics, net, oaf, $ 298 298 cor, bias, std, rmsd 299 299 print, cor, bias, std, rmsd 300 300 cstat_oaf=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 301 301 ; 302 302 statistics, net, nce1, $ 303 303 cor, bias, std, rmsd 304 304 print, cor, bias, std, rmsd 305 305 cstat_nce1=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 306 306 ; 307 307 statistics, net, nce2, $ 308 308 cor, bias, std, rmsd 309 309 print, cor, bias, std, rmsd 310 310 cstat_nce2=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 311 311 ; 312 312 statistics, net, olr, $ 313 313 cor, bias, std, rmsd 314 314 print, cor, bias, std, rmsd 315 315 cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') 316 317 316 ; 317 ; 318 318 splot, net, era, charsize=1.1, title='Net Flux - Subdn. cs Vs ERAI', $ 319 319 xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 320 320 oplot, [-200,200], [-200,200] 321 321 xyouts, 0,-175, cstat_era, charsize=1. 322 322 ; 323 323 splot, net, trop, charsize=1.1, title='Net Flux - Subdn. cs Vs TropFlux',/noer, $ 324 324 xrange=[-200,200], yrange=[-200,200], small=[2,3,2], psym=2, xmin=1, ymin=1 325 325 oplot, [-200,200], [-200,200] 326 326 xyouts, 0,-175, cstat_trop, charsize=1. 327 327 ; 328 328 splot, net, oaf, charsize=1.1, title='Net Flux - Subdn. cs Vs OAFlux',/noer, $ 329 329 xrange=[-200,200], yrange=[-200,200], small=[2,3,3], psym=2, xmin=1, ymin=1 330 330 oplot, [-200,200], [-200,200] 331 331 xyouts, 0,-175, cstat_oaf, charsize=1. 332 332 ; 333 333 splot, net, nce1, charsize=1.1, title='Net Flux - Subdn. cs Vs NCEP1',/noer, $ 334 334 xrange=[-200,200], yrange=[-200,200], small=[2,3,4], psym=2, xmin=1, ymin=1 335 335 oplot, [-200,200], [-200,200] 336 336 xyouts, 0,-175, cstat_nce1, charsize=1. 337 337 ; 338 338 splot, net, nce2, charsize=1.1, title='Net Flux - Subdn. cs Vs NCEP2',/noer, $ 339 339 xrange=[-200,200], yrange=[-200,200], small=[2,3,5], psym=2, xmin=1, ymin=1 340 340 oplot, [-200,200], [-200,200] 341 341 xyouts, 0,-175, cstat_nce2, charsize=1. 342 342 ; 343 343 splot, net, olr, charsize=1.1, title='Net Flux - Subdn. cs Vs TropFlux_NRT',/noer, $ 344 344 xrange=[-200,200], yrange=[-200,200], small=[2,3,6], psym=2, xmin=1, ymin=1 345 345 oplot, [-200,200], [-200,200] 346 346 xyouts, 0,-175, cstat_olr, charsize=1. 347 347 ; 348 348 closeps 349 ; 349 350 end
Note: See TracChangeset
for help on using the changeset viewer.