[43] | 1 | ;+ |
---|
[46] | 2 | ; .. _as_validation_net_1994_95_v10.pro: |
---|
[43] | 3 | ; |
---|
| 4 | ; ================================= |
---|
| 5 | ; as_validation_net_1994_95_v10.pro |
---|
| 6 | ; ================================= |
---|
| 7 | ; |
---|
| 8 | ; DESCRIPTION |
---|
| 9 | ; =========== |
---|
| 10 | ; |
---|
[46] | 11 | ; .. graphviz:: |
---|
| 12 | ; |
---|
| 13 | ; digraph as_validation_net_1994_95_v10 { |
---|
[85] | 14 | ; |
---|
[50] | 15 | ; met_data [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/met_data_as_expt_19*.txt"]; |
---|
| 16 | ; swr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_swr_19910101_20091231_oafluxgrid.nc"]; |
---|
| 17 | ; lwr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_lwr_19910101_20091231_oafluxgrid.nc"]; |
---|
| 18 | ; lhf_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_lhf_19910101_20091231_oafluxgrid.nc"]; |
---|
| 19 | ; tropflux_general [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_19890101_20091231.nc"]; |
---|
| 20 | ; swr_tropflux_nrt [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_swr_19890101_20091231_NRT_v50.nc"]; |
---|
| 21 | ; swr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_oafluxgrid_1985_2007.nc"]; |
---|
| 22 | ; lwr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lwr_oafluxgrid_1985_2007.nc"]; |
---|
| 23 | ; lhf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_oafluxgrid_1985_2007.nc"]; |
---|
| 24 | ; shf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_oafluxgrid_1985_2007.nc"]; |
---|
| 25 | ; swr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_ncep2_oaflxgrid_19890101_20091231.nc"]; |
---|
| 26 | ; lwr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/NCEP2_flux_19890101_20090729.nc"]; |
---|
| 27 | ; lhf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_ncep2_oafluxgrid_19890101_20091231.nc"]; |
---|
| 28 | ; shf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_ncep2_oafluxgrid_19890101_20091231.nc"]; |
---|
[46] | 29 | ; |
---|
[50] | 30 | ; ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/ncep1_flux_19890101_20091231.nc"]; |
---|
[46] | 31 | ; |
---|
[50] | 32 | ; figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/as_validation_net_1994_95_v10.ps"]; |
---|
[46] | 33 | ; |
---|
| 34 | ; as_validation_net_1994_95_v10 [shape=box, |
---|
| 35 | ; fontname=Courier, |
---|
| 36 | ; color=blue, |
---|
| 37 | ; URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig14/as_validation_net_1994_95_v10.pro", |
---|
| 38 | ; label="${TROPFLUX}/src/paper01/fig14/as_validation_net_1994_95_v10.pro"]; |
---|
[50] | 39 | ; {met_data swr_erai lwr_erai lhf_erai tropflux_general swr_tropflux_nrt swr_oafluxgrid lwr_oafluxgrid lhf_oafluxgrid shf_oafluxgrid swr_ncep2 lwr_ncep2 lhf_ncep2 shf_ncep2 ncep1} -> {as_validation_net_1994_95_v10} -> {figure} |
---|
[46] | 40 | ; } |
---|
| 41 | ; |
---|
[43] | 42 | ; SEE ALSO |
---|
| 43 | ; ======== |
---|
| 44 | ; |
---|
[50] | 45 | ; :ref:`project_profile.sh` |
---|
| 46 | ; :ref:`project_init.pro` |
---|
| 47 | ; :ref:`cm_project.pro` |
---|
[43] | 48 | ; |
---|
[46] | 49 | ; :ref:`statistics.pro` |
---|
| 50 | ; |
---|
[43] | 51 | ; EXAMPLES |
---|
| 52 | ; ======== |
---|
| 53 | ; |
---|
| 54 | ; :: |
---|
| 55 | ; |
---|
[180] | 56 | ; as_validation_net_1994_95_v10 |
---|
[43] | 57 | ; |
---|
[46] | 58 | ; TODO |
---|
| 59 | ; ==== |
---|
| 60 | ; |
---|
[50] | 61 | ; make it work on cratos : missing data erai_swr_19910101_20091231_oafluxgrid.nc, etc |
---|
[47] | 62 | ; |
---|
[46] | 63 | ; coding rules |
---|
| 64 | ; |
---|
| 65 | ; complete description |
---|
| 66 | ; |
---|
| 67 | ; handle IO error |
---|
| 68 | ; |
---|
[43] | 69 | ; EVOLUTIONS |
---|
| 70 | ; ========== |
---|
| 71 | ; |
---|
| 72 | ; $Id$ |
---|
| 73 | ; |
---|
[85] | 74 | ; $URL$ |
---|
| 75 | ; |
---|
[46] | 76 | ; - fplod 20110420T115926Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 77 | ; |
---|
| 78 | ; * remove hard coding path |
---|
| 79 | ; * add graphviz |
---|
| 80 | ; |
---|
[43] | 81 | ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 82 | ; |
---|
| 83 | ; * minimal header |
---|
| 84 | ; |
---|
| 85 | ;- |
---|
[41] | 86 | pro as_validation_net_1994_95_v10 |
---|
[46] | 87 | @cm_general |
---|
[50] | 88 | @cm_project |
---|
[41] | 89 | reinitplt, /z,/invert |
---|
| 90 | key_portrait = 1 |
---|
| 91 | coefpalit=.9 |
---|
[97] | 92 | ; |
---|
[85] | 93 | openps, FILENAME = project_od_env+'as_validation_net_1994_95_v10.ps' |
---|
[41] | 94 | ; partie a changer |
---|
| 95 | asbox=[61,62,15,16] |
---|
[97] | 96 | date1=19941016 |
---|
[94] | 97 | date2=19951019 |
---|
[97] | 98 | ; |
---|
[50] | 99 | fi=project_id_env+"met_data_as_expt_1994.txt" |
---|
[41] | 100 | res=read_ascii(fi,data_start=1) |
---|
| 101 | ff=res.field1 |
---|
| 102 | swr94=reform(ff(7,*))*0.93 |
---|
| 103 | lwr94=reform(ff(8,*)) |
---|
| 104 | lhf94=reform(ff(5,*)) |
---|
| 105 | shf94=reform(ff(6,*)) |
---|
| 106 | net94=swr94+lwr94+lhf94+shf94 |
---|
[97] | 107 | ; |
---|
[50] | 108 | fi=project_id_env+"met_data_as_expt_1995.txt" |
---|
[41] | 109 | res=read_ascii(fi,data_start=1) |
---|
| 110 | ff=res.field1 |
---|
| 111 | swr95=reform(ff(7,*))*0.93 |
---|
| 112 | lwr95=reform(ff(8,*)) |
---|
| 113 | lhf95=reform(ff(5,*)) |
---|
| 114 | shf95=reform(ff(6,*)) |
---|
| 115 | net95=swr95+lwr95+lhf95+shf95 |
---|
[97] | 116 | ; |
---|
[41] | 117 | net=[net94, net95] |
---|
| 118 | ind=where(finite(net)) |
---|
| 119 | net=net(ind) |
---|
[97] | 120 | ; |
---|
| 121 | ; start data - 16/10/94 and end date 31/12/94 |
---|
| 122 | ; start data - 01/01/95 and end date 19/10/95 |
---|
| 123 | ; |
---|
[41] | 124 | help, net |
---|
[97] | 125 | ; reading other data sets |
---|
[50] | 126 | file=project_id_env+"erai_swr_19910101_20091231_oafluxgrid.nc" |
---|
[41] | 127 | initncdf, file |
---|
[97] | 128 | swr=read_ncdf('swr',date1,date2,file=file,/nostr, box=asbox) |
---|
[94] | 129 | swr=reform(swr) |
---|
[50] | 130 | file=project_id_env+'erai_lwr_19890101_20091231_oafluxgrid.nc' |
---|
[41] | 131 | initncdf, file |
---|
[97] | 132 | lwr=read_ncdf('lwr',date1-1,date2,file=file,/nostr, box=asbox) |
---|
[94] | 133 | lwr=reform(lwr) |
---|
[50] | 134 | file=project_id_env+'erai_shf_19890101_20091231_oafluxgrid.nc' |
---|
[41] | 135 | initncdf, file |
---|
[97] | 136 | shf=-1*read_ncdf('shf',date1,date2,file=file,/nostr, box=asbox) |
---|
[94] | 137 | shf=reform(shf) |
---|
[50] | 138 | file=project_id_env+'erai_lhf_19890101_20091231_oafluxgrid.nc' |
---|
[41] | 139 | initncdf, file |
---|
[97] | 140 | lhf=-1*read_ncdf('lhf',date1,date2,file=file,/nostr, box=asbox) |
---|
[94] | 141 | lhf=reform(lhf) |
---|
[97] | 142 | ; |
---|
[41] | 143 | era=swr+lwr+lhf+shf |
---|
| 144 | help, era |
---|
| 145 | era=era(ind) |
---|
| 146 | help, era |
---|
[97] | 147 | ; |
---|
| 148 | swr=0. |
---|
| 149 | lwr=0. |
---|
| 150 | shf=0. |
---|
[94] | 151 | lhf=0. |
---|
[97] | 152 | ; |
---|
[50] | 153 | file=project_id_env+"TropFlux_19890101_20091231.nc" |
---|
[41] | 154 | initncdf, file |
---|
[97] | 155 | swr=read_ncdf('swr',date1,date2,file=file,/nostr, box=asbox) |
---|
[94] | 156 | swr=reform(swr) |
---|
[97] | 157 | lwr=read_ncdf('lwr',date1,date2,file=file,/nostr, box=asbox) |
---|
[94] | 158 | lwr=reform(lwr) |
---|
[97] | 159 | shf=read_ncdf('shf',date1,date2,file=file,/nostr, box=asbox) |
---|
[94] | 160 | shf=reform(shf) |
---|
[97] | 161 | lhf=read_ncdf('lhf',date1,date2,file=file,/nostr, box=asbox) |
---|
[94] | 162 | lhf=reform(lhf) |
---|
[41] | 163 | trop=swr+lwr+shf+lhf |
---|
| 164 | help, trop |
---|
[97] | 165 | ; |
---|
[50] | 166 | file=project_id_env+"TropFlux_swr_19890101_20091231_NRT_v50.nc" |
---|
[41] | 167 | initncdf, file |
---|
[97] | 168 | olr1=read_ncdf("sw", date1-1, date2, file=file,/nostr, box=asbox) |
---|
[94] | 169 | olr1=grossemoyenne(olr1, "xy",/nan, box=asbox) |
---|
[41] | 170 | olr1=reform(olr1) |
---|
[97] | 171 | ; |
---|
[41] | 172 | olr=olr1+lwr+shf+lhf |
---|
[97] | 173 | ; |
---|
| 174 | swr=0. |
---|
| 175 | lwr=0. |
---|
| 176 | shf=0. |
---|
| 177 | lhf=0. |
---|
[94] | 178 | olr1=0. |
---|
[97] | 179 | ; |
---|
| 180 | trop=trop(ind) |
---|
[94] | 181 | olr=olr(ind) |
---|
[41] | 182 | help, trop, olr |
---|
[97] | 183 | ; |
---|
[50] | 184 | file=project_id_env+'swr_oafluxgrid_1985_2007.nc' |
---|
[41] | 185 | initncdf, file |
---|
[97] | 186 | swr=read_ncdf("swr", date1, date2, file=file,/nostr, box=asbox) |
---|
[94] | 187 | swr=reform(swr) |
---|
[97] | 188 | ; |
---|
[50] | 189 | file=project_id_env+'lwr_oafluxgrid_1985_2007.nc' |
---|
[41] | 190 | initncdf, file |
---|
[97] | 191 | lwr=read_ncdf("lwr", date1, date2, file=file,/nostr, box=asbox) |
---|
[94] | 192 | lwr=-1*reform(lwr) |
---|
[97] | 193 | ; |
---|
[50] | 194 | file=project_id_env+'lhf_oafluxgrid_1985_2009.nc' |
---|
[41] | 195 | initncdf, file |
---|
[97] | 196 | lhf=read_ncdf("lhf", date1, date2, file=file,/nostr, box=asbox) |
---|
[94] | 197 | lhf=-1*reform(lhf) |
---|
[97] | 198 | ; |
---|
[50] | 199 | file=project_id_env+'shf_oafluxgrid_1985_2009.nc' |
---|
[41] | 200 | initncdf, file |
---|
[97] | 201 | shf=read_ncdf("shf", date1, date2, file=file,/nostr, box=asbox) |
---|
[94] | 202 | shf=-1*reform(shf) |
---|
[97] | 203 | ; |
---|
[41] | 204 | oaf=swr+lwr+shf+lhf |
---|
| 205 | oaf=oaf(ind) |
---|
| 206 | help, oaf |
---|
[97] | 207 | ; |
---|
| 208 | swr=0. |
---|
| 209 | lwr=0. |
---|
| 210 | shf=0. |
---|
[94] | 211 | lhf=0. |
---|
[97] | 212 | ; |
---|
[50] | 213 | fi=project_id_env+'swr_ncep2_oaflxgrid_19890101_20091231.nc' |
---|
[41] | 214 | initncdf, fi |
---|
[97] | 215 | swr=read_ncdf("swr", date1, date2, file=fi,/nostr, box=asbox) |
---|
[94] | 216 | swr=reform(swr) |
---|
[97] | 217 | ; |
---|
[50] | 218 | fi=project_id_env+"NCEP2_flux_19890101_20090729.nc" |
---|
[41] | 219 | initncdf, fi |
---|
[97] | 220 | lwr=read_ncdf("lwr", date1, date2, file=fi,/nostr, box=asbox) |
---|
| 221 | lwr=grossemoyenne(lwr, "xy",/nan, box=asbox) |
---|
[94] | 222 | lwr=reform(lwr) |
---|
[97] | 223 | ; |
---|
[41] | 224 | lwr=-1*lwr |
---|
[97] | 225 | ; |
---|
[50] | 226 | fi=project_id_env+'lhf_ncep2_oafluxgrid_19890101_20091231.nc' |
---|
[41] | 227 | initncdf, fi |
---|
[97] | 228 | lhf=read_ncdf("lhf", date1, date2, file=fi,/nostr, box=asbox) |
---|
[94] | 229 | lhf=-1*reform(lhf) |
---|
[97] | 230 | ; |
---|
[50] | 231 | fi=project_id_env+'shf_ncep2_oafluxgrid_19890101_20091231.nc' |
---|
[41] | 232 | initncdf, fi |
---|
[97] | 233 | shf=read_ncdf("shf", date1, date2, file=fi,/nostr, box=asbox) |
---|
[94] | 234 | shf=-1*reform(shf) |
---|
[97] | 235 | ; |
---|
[41] | 236 | nce2=swr+lwr+shf+lhf |
---|
| 237 | help, nce2 |
---|
[97] | 238 | ; |
---|
[41] | 239 | nce2=nce2(ind) |
---|
| 240 | help, nce2 |
---|
[97] | 241 | swr=0. |
---|
| 242 | lwr=0. |
---|
| 243 | shf=0. |
---|
[94] | 244 | lhf=0. |
---|
[97] | 245 | ; |
---|
[50] | 246 | file=project_id_env+"ncep1_flux_19890101_20091231.nc" |
---|
[41] | 247 | initncdf, file |
---|
[94] | 248 | swr=read_ncdf("swr", date1, date2, file=file,/nostr, box=asbox) |
---|
| 249 | swr=reform(swr) |
---|
[97] | 250 | lwr=read_ncdf("lwr", date1, date2, file=file,/nostr, box=asbox) |
---|
[94] | 251 | lwr=reform(lwr) |
---|
[97] | 252 | lhf=read_ncdf("lhf", date1, date2, file=file,/nostr, box=asbox) |
---|
[94] | 253 | lhf=reform(lhf) |
---|
[97] | 254 | shf=read_ncdf("shf", date1, date2, file=file,/nostr, box=asbox) |
---|
[94] | 255 | shf=reform(shf) |
---|
[41] | 256 | nce1=swr+lwr+shf+lhf |
---|
[97] | 257 | ; |
---|
[41] | 258 | help, nce1 |
---|
[97] | 259 | ; |
---|
[41] | 260 | nce1=nce1(ind) |
---|
| 261 | help, nce1 |
---|
| 262 | ;ind=where(finite(net)) |
---|
[97] | 263 | ;net=net(ind) |
---|
| 264 | ;era=era(ind) |
---|
| 265 | ;trop=trop(ind) |
---|
| 266 | ;oaf=oaf(ind) |
---|
| 267 | ;nce1=nce1(ind) |
---|
[94] | 268 | ;nce2=nce2(ind) |
---|
[41] | 269 | ;olr=olr(ind) |
---|
[97] | 270 | ; |
---|
[41] | 271 | statistics, net, era, $ |
---|
| 272 | cor, bias, std, rmsd |
---|
| 273 | print, cor, bias, std, rmsd |
---|
| 274 | cstat_era=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 275 | ; |
---|
[41] | 276 | statistics, net, trop, $ |
---|
| 277 | cor, bias, std, rmsd |
---|
| 278 | print, cor, bias, std, rmsd |
---|
| 279 | cstat_trop=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 280 | ; |
---|
[41] | 281 | statistics, net, oaf, $ |
---|
| 282 | cor, bias, std, rmsd |
---|
| 283 | print, cor, bias, std, rmsd |
---|
| 284 | cstat_oaf=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 285 | ; |
---|
[41] | 286 | statistics, net, nce1, $ |
---|
| 287 | cor, bias, std, rmsd |
---|
| 288 | print, cor, bias, std, rmsd |
---|
| 289 | cstat_nce1=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 290 | ; |
---|
[41] | 291 | statistics, net, nce2, $ |
---|
| 292 | cor, bias, std, rmsd |
---|
| 293 | print, cor, bias, std, rmsd |
---|
| 294 | cstat_nce2=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 295 | ; |
---|
[41] | 296 | statistics, net, olr, $ |
---|
| 297 | cor, bias, std, rmsd |
---|
| 298 | print, cor, bias, std, rmsd |
---|
| 299 | cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 300 | ; |
---|
[41] | 301 | splot, net, era, charsize=1.1, title='Net Flux - AS Vs ERAI', $ |
---|
| 302 | xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 |
---|
| 303 | oplot, [-200,200], [-200,200] |
---|
| 304 | xyouts, -150,-175, cstat_era, charsize=1. |
---|
[97] | 305 | ; |
---|
[41] | 306 | splot, net, trop, charsize=1.1, title='Net Flux - AS Vs TropFlux',/noer, $ |
---|
| 307 | xrange=[-200,200], yrange=[-200,200], small=[2,3,2], psym=2, xmin=1, ymin=1 |
---|
| 308 | oplot, [-200,200], [-200,200] |
---|
| 309 | xyouts, -150,-175, cstat_trop, charsize=1. |
---|
[97] | 310 | ; |
---|
[41] | 311 | splot, net, oaf, charsize=1.1, title='Net Flux - AS Vs OAFlux',/noer, $ |
---|
| 312 | xrange=[-200,200], yrange=[-200,200], small=[2,3,3], psym=2, xmin=1, ymin=1 |
---|
| 313 | oplot, [-200,200], [-200,200] |
---|
| 314 | xyouts, -150,-175, cstat_oaf, charsize=1. |
---|
[97] | 315 | ; |
---|
[41] | 316 | splot, net, nce1, charsize=1.1, title='Net Flux - AS Vs NCEP1',/noer, $ |
---|
| 317 | xrange=[-200,200], yrange=[-200,200], small=[2,3,4], psym=2, xmin=1, ymin=1 |
---|
| 318 | oplot, [-200,200], [-200,200] |
---|
| 319 | xyouts, -150,-175, cstat_nce1, charsize=1. |
---|
[97] | 320 | ; |
---|
[41] | 321 | splot, net, nce2, charsize=1.1, title='Net Flux - AS Vs NCEP2',/noer, $ |
---|
| 322 | xrange=[-200,200], yrange=[-200,200], small=[2,3,5], psym=2, xmin=1, ymin=1 |
---|
| 323 | oplot, [-200,200], [-200,200] |
---|
| 324 | xyouts, -150,-175, cstat_nce2, charsize=1. |
---|
[97] | 325 | ; |
---|
[41] | 326 | splot, net, olr, charsize=1.1, title='Net Flux - AS Vs TropFlux_NRT',/noer, $ |
---|
| 327 | xrange=[-200,200], yrange=[-200,200], small=[2,3,6], psym=2, xmin=1, ymin=1 |
---|
| 328 | oplot, [-200,200], [-200,200] |
---|
| 329 | xyouts, -150,-175, cstat_olr, charsize=1. |
---|
[97] | 330 | ; |
---|
[41] | 331 | closeps |
---|
[97] | 332 | ; |
---|
[41] | 333 | end |
---|