Changeset 43 for trunk/src/paper01/fig14
- Timestamp:
- 04/11/11 18:05:51 (13 years ago)
- Location:
- trunk/src/paper01/fig14
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
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
Note: See TracChangeset
for help on using the changeset viewer.