[43] | 1 | ;+ |
---|
[46] | 2 | ; .. _whots_validation_net_2004_2007_all_v10.pro: |
---|
[43] | 3 | ; |
---|
| 4 | ; ========================================== |
---|
| 5 | ; whots_validation_net_2004_2007_all_v10.pro |
---|
| 6 | ; ========================================== |
---|
| 7 | ; |
---|
| 8 | ; DESCRIPTION |
---|
| 9 | ; =========== |
---|
| 10 | ; |
---|
[46] | 11 | ; .. graphviz:: |
---|
| 12 | ; |
---|
| 13 | ; digraph whots_validation_net_2004_2007_all_v10 { |
---|
| 14 | ; graph [ |
---|
| 15 | ; rankdir="LR", |
---|
| 16 | ; ] |
---|
| 17 | ; met_data [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/met_data_WHOTS_200?_200?.txt"]; |
---|
| 18 | ; swr_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_swr_19910101_20091231_oafluxgrid.nc"]; |
---|
| 19 | ; lwr_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_lwr_19910101_20091231_oafluxgrid.nc"]; |
---|
| 20 | ; lhf_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_lhf_19910101_20091231_oafluxgrid.nc"]; |
---|
| 21 | ; shf_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_shf_19910101_20091231_oafluxgrid.nc"]; |
---|
| 22 | ; tropflux_general [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_19890101_20091231.nc"]; |
---|
| 23 | ; swr_tropflux_nrt [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_swr_19890101_20091231_NRT_v50.nc"]; |
---|
| 24 | ; swr_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/swr_oafluxgrid_1985_2007.nc"]; |
---|
| 25 | ; lwr_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lwr_oafluxgrid_1985_2007.nc"]; |
---|
| 26 | ; lhf_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lhf_oafluxgrid_1985_2007.nc"]; |
---|
| 27 | ; shf_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/shf_oafluxgrid_1985_2007.nc"]; |
---|
| 28 | ; swr_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/swr_ncep2_oaflxgrid_19890101_20091231.nc"]; |
---|
| 29 | ; lwr_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/NCEP2_flux_19890101_20090729.nc"]; |
---|
| 30 | ; lhf_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lhf_ncep2_oafluxgrid_19890101_20091231.nc"]; |
---|
| 31 | ; shf_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/shf_ncep2_oafluxgrid_19890101_20091231.nc"]; |
---|
| 32 | ; |
---|
| 33 | ; ncep1 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/ncep1_flux_19890101_20091231.nc"]; |
---|
| 34 | ; |
---|
| 35 | ; figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/WHOTS_validation_net_2004_2007_all_v10.ps"]; |
---|
| 36 | ; |
---|
| 37 | ; whots_validation_net_2004_2007_all_v10 [shape=box, |
---|
| 38 | ; fontname=Courier, |
---|
| 39 | ; color=blue, |
---|
| 40 | ; URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig14/whots_validation_net_2004_2007_all_v10.pro", |
---|
| 41 | ; label="${TROPFLUX}/src/paper01/fig14/whots_validation_net_2004_2007_all_v10.pro"]; |
---|
| 42 | ; |
---|
| 43 | ; {met_data swr_erai lwr_erai lhf_erai shf_erai tropflux_general swr_tropflux_nrt swr_oafluxgrid lwr_oafluxgrid lhf_oafluxgrid shf_oafluxgrid swr_ncep2 lwr_ncep2 lhf_ncep2 shf_ncep2 ncep1} -> {whots_validation_net_2004_2007_all_v10} -> {figure} |
---|
| 44 | ; } |
---|
| 45 | ; |
---|
[43] | 46 | ; SEE ALSO |
---|
| 47 | ; ======== |
---|
| 48 | ; |
---|
| 49 | ; :ref:`tropflux_profile.sh` |
---|
[46] | 50 | ; :ref:`tropflux_init.pro` |
---|
| 51 | ; :ref:`cm_tropflux.pro` |
---|
[43] | 52 | ; |
---|
[46] | 53 | ; :ref:`statistics.pro` |
---|
| 54 | ; |
---|
[43] | 55 | ; EXAMPLES |
---|
| 56 | ; ======== |
---|
| 57 | ; |
---|
| 58 | ; :: |
---|
| 59 | ; |
---|
| 60 | ; IDL> @tropflux_init |
---|
| 61 | ; IDL> whots_validation_net_2004_2007_all_v10 |
---|
| 62 | ; |
---|
[46] | 63 | ; TODO |
---|
| 64 | ; ==== |
---|
| 65 | ; |
---|
| 66 | ; make it work on cratos : missing data |
---|
| 67 | ; |
---|
| 68 | ; coding rules |
---|
| 69 | ; |
---|
| 70 | ; get rid of spwan |
---|
| 71 | ; |
---|
| 72 | ; complete description |
---|
| 73 | ; |
---|
| 74 | ; handle IO error |
---|
| 75 | ; |
---|
| 76 | ; EVOLUTIONS |
---|
[43] | 77 | ; ========== |
---|
| 78 | ; |
---|
| 79 | ; $Id$ |
---|
| 80 | ; |
---|
[46] | 81 | ; - fplod 20110420T122057Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 82 | ; |
---|
| 83 | ; * remove hard coding path |
---|
| 84 | ; * add graphviz |
---|
| 85 | ; |
---|
[43] | 86 | ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 87 | ; |
---|
| 88 | ; * minimal header |
---|
| 89 | ; |
---|
| 90 | ;- |
---|
[41] | 91 | pro whots_validation_net_2004_2007_all_v10 |
---|
[46] | 92 | @cm_general |
---|
| 93 | @cm_tropflux |
---|
[41] | 94 | reinitplt, /z,/invert |
---|
| 95 | key_portrait = 1 |
---|
| 96 | coefpalit=.9 |
---|
| 97 | |
---|
| 98 | openps, FILENAME = 'idl.ps' |
---|
| 99 | ; partie a changer |
---|
| 100 | |
---|
| 101 | ;; WHOTS in the NC pacific 2004-09 |
---|
| 102 | ;22.45N, 157.54W (202.46E) |
---|
| 103 | box=[202, 203, 22., 23.] |
---|
| 104 | |
---|
| 105 | ;; first leg |
---|
| 106 | |
---|
[46] | 107 | fi=tropflux_id_env+"met_data_WHOTS_2004_2005.txt" |
---|
[41] | 108 | res=read_ascii(fi,data_start=1) |
---|
| 109 | ff=res.field1 |
---|
| 110 | lhf_04=reform(ff(5,*)) |
---|
| 111 | shf_04=reform(ff(6,*)) |
---|
| 112 | swr_04=reform(ff(7,*))/0.93 |
---|
| 113 | lwr_04=reform(ff(8,*)) |
---|
| 114 | net_04=[lhf_04+shf_04+swr_04+lwr_04] |
---|
| 115 | |
---|
[46] | 116 | fi=tropflux_id_env+"met_data_WHOTS_2005_2006.txt" |
---|
[41] | 117 | res=read_ascii(fi,data_start=1) |
---|
| 118 | ff=res.field1 |
---|
| 119 | lhf_05=reform(ff(5,*)) |
---|
| 120 | shf_05=reform(ff(6,*)) |
---|
| 121 | swr_05=reform(ff(7,*))/0.93 |
---|
| 122 | lwr_05=reform(ff(8,*)) |
---|
| 123 | net_05=[lhf_05+shf_05+swr_05+lwr_05] |
---|
| 124 | |
---|
[46] | 125 | fi=tropflux_id_env+"met_data_WHOTS_2006_2007.txt" |
---|
[41] | 126 | res=read_ascii(fi,data_start=1) |
---|
| 127 | ff=res.field1 |
---|
| 128 | lhf_06=reform(ff(5,*)) |
---|
| 129 | shf_06=reform(ff(6,*)) |
---|
| 130 | swr_06=reform(ff(7,*))/0.93 |
---|
| 131 | lwr_06=reform(ff(8,*)) |
---|
| 132 | net_06=[lhf_06+shf_06+swr_06+lwr_06] |
---|
| 133 | |
---|
| 134 | net=[net_04, net_05, net_06] |
---|
| 135 | |
---|
| 136 | help, net |
---|
| 137 | |
---|
| 138 | ;; reading other data sets |
---|
| 139 | |
---|
[43] | 140 | st1=20040813 & en1=20050725 |
---|
[41] | 141 | st2=20050728 & en2=20060624 |
---|
[43] | 142 | st3=20060627 & en3=20070628 |
---|
| 143 | ;st4=20070626 & en4=20080606 |
---|
[41] | 144 | ;st5=20080605 & en5=20090715 |
---|
| 145 | |
---|
[46] | 146 | file=tropflux_id_env+'erai_lhf_19890101_20091231_oafluxgrid.nc' |
---|
[41] | 147 | initncdf, file |
---|
| 148 | lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) & lhf1=reform(lhf1) |
---|
| 149 | lhf2=-1*read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) & lhf2=reform(lhf2) |
---|
| 150 | lhf3=-1*read_ncdf('lhf',st3,en3,file=file,/nostr, box=box) & lhf3=grossemoyenne(lhf3, "xy", box=box,/nan) & lhf3=reform(lhf3) |
---|
| 151 | |
---|
[46] | 152 | file=tropflux_id_env+'erai_shf_19890101_20091231_oafluxgrid.nc' |
---|
[41] | 153 | initncdf, file |
---|
| 154 | shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy", box=box,/nan) & shf1=reform(shf1) |
---|
| 155 | shf2=-1*read_ncdf('shf',st2,en2,file=file,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy", box=box,/nan) & shf2=reform(shf2) |
---|
| 156 | shf3=-1*read_ncdf('shf',st3,en3,file=file,/nostr, box=box) & shf3=grossemoyenne(shf3, "xy", box=box,/nan) & shf3=reform(shf3) |
---|
| 157 | |
---|
[46] | 158 | file=tropflux_id_env+'erai_lwr_19890101_20091231_oafluxgrid.nc' |
---|
[41] | 159 | initncdf, file |
---|
| 160 | lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) & lwr1=reform(lwr1) |
---|
| 161 | lwr2=read_ncdf('lwr',st2-1,en2,file=file,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) & lwr2=reform(lwr2) |
---|
| 162 | lwr3=read_ncdf('lwr',st3-1,en3,file=file,/nostr, box=box) & lwr3=grossemoyenne(lwr3, "xy", box=box,/nan) & lwr3=reform(lwr3) |
---|
| 163 | |
---|
[46] | 164 | file=tropflux_id_env+"erai_swr_19910101_20091231_oafluxgrid.nc" |
---|
[41] | 165 | initncdf, file |
---|
| 166 | swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy", box=box,/nan) & swr1=reform(swr1) |
---|
| 167 | swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy", box=box,/nan) & swr2=reform(swr2) |
---|
| 168 | swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box) & swr3=grossemoyenne(swr3, "xy", box=box,/nan) & swr3=reform(swr3) |
---|
| 169 | |
---|
| 170 | net1=swr1+(lwr1+shf1+lhf1) & net2=swr2+(lwr2+shf2+lhf2) & net3=swr3+(lwr3+shf3+lhf3) |
---|
| 171 | era=[net1,net2,net3] |
---|
| 172 | lhf1=0. & lhf2=0. & lhf3=0. |
---|
| 173 | shf1=0. & shf2=0. & shf3=0. |
---|
[43] | 174 | lwr1=0. & lwr2=0. & lwr3=0. |
---|
[41] | 175 | swr1=0. & swr2=0. & swr3=0. |
---|
[43] | 176 | net1=0. & net2=0. & net3=0. |
---|
[41] | 177 | help, era |
---|
| 178 | |
---|
[46] | 179 | file=tropflux_id_env+"TropFlux_19890101_20091231.nc" |
---|
[41] | 180 | initncdf, file |
---|
| 181 | |
---|
| 182 | lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) & lhf1=reform(lhf1) |
---|
| 183 | lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) & lhf2=reform(lhf2) |
---|
| 184 | lhf3=read_ncdf('lhf',st3,en3,file=file,/nostr, box=box) & lhf3=grossemoyenne(lhf3, "xy", box=box,/nan) & lhf3=reform(lhf3) |
---|
| 185 | |
---|
| 186 | shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy", box=box,/nan) & shf1=reform(shf1) |
---|
| 187 | shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy", box=box,/nan) & shf2=reform(shf2) |
---|
| 188 | shf3=read_ncdf('shf',st3,en3,file=file,/nostr, box=box) & shf3=grossemoyenne(shf3, "xy", box=box,/nan) & shf3=reform(shf3) |
---|
| 189 | |
---|
| 190 | lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) & lwr1=reform(lwr1) |
---|
| 191 | lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) & lwr2=reform(lwr2) |
---|
| 192 | lwr3=read_ncdf('lwr',st3,en3,file=file,/nostr, box=box) & lwr3=grossemoyenne(lwr3, "xy", box=box,/nan) & lwr3=reform(lwr3) |
---|
| 193 | |
---|
| 194 | swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy", box=box,/nan) & swr1=reform(swr1) |
---|
| 195 | swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy", box=box,/nan) & swr2=reform(swr2) |
---|
| 196 | swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box) & swr3=grossemoyenne(swr3, "xy", box=box,/nan) & swr3=reform(swr3) |
---|
| 197 | |
---|
[46] | 198 | file=tropflux_id_env+"TropFlux_swr_19890101_20091231_NRT_v50.nc" |
---|
[41] | 199 | initncdf, file |
---|
| 200 | olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box) & olr1=grossemoyenne(olr1, "xy",/nan, box=box) & olr1=reform(olr1) |
---|
| 201 | olr2=read_ncdf("sw", st2-1, en2, file=file,/nostr, box=box) & olr2=grossemoyenne(olr2, "xy",/nan, box=box) & olr2=reform(olr2) |
---|
| 202 | olr3=read_ncdf("sw", st3-1, en3, file=file,/nostr, box=box) & olr3=grossemoyenne(olr3, "xy",/nan, box=box) & olr3=reform(olr3) |
---|
| 203 | |
---|
| 204 | net1=swr1+lwr1+(shf1+lhf1) & net2=swr2+lwr2+(shf2+lhf2) & net3=swr3+lwr3+(shf3+lhf3) |
---|
[43] | 205 | neto_1=olr1+lwr1+lhf1+shf1 & neto_2=olr2+lwr2+lhf2+shf2 & neto_3=olr3+lwr3+lhf3+shf3 |
---|
[41] | 206 | |
---|
| 207 | trop=[net1,net2,net3] & olr=[neto_1,neto_2,neto_3] |
---|
| 208 | |
---|
| 209 | lhf1=0. & lhf2=0. & lhf3=0. |
---|
[43] | 210 | shf1=0. & shf2=0. & shf3=0. |
---|
[41] | 211 | lwr1=0. & lwr2=0. & lwr3=0. |
---|
[43] | 212 | swr1=0. & swr2=0. & swr3=0. |
---|
| 213 | net1=0. & net2=0. & net3=0. |
---|
[41] | 214 | neto_1=0 & neto_2=0 & neto_3=0 |
---|
| 215 | |
---|
| 216 | help, trop, olr |
---|
| 217 | |
---|
[46] | 218 | file=tropflux_id_env+'swr_oafluxgrid_1985_2007.nc' |
---|
[41] | 219 | initncdf, file |
---|
| 220 | swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy",/nan, box=box) & swr1=reform(swr1) |
---|
| 221 | swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy",/nan, box=box) & swr2=reform(swr2) |
---|
| 222 | swr3=read_ncdf("swr", st3, en3, file=file,/nostr, box=box) & swr3=grossemoyenne(swr3, "xy",/nan, box=box) & swr3=reform(swr3) |
---|
| 223 | |
---|
[46] | 224 | file=tropflux_id_env+'lwr_oafluxgrid_1985_2007.nc' |
---|
[41] | 225 | initncdf, file |
---|
| 226 | lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) & lwr1=reform(lwr1) |
---|
| 227 | lwr2=-1*read_ncdf("lwr", st2, en2, file=file,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) & lwr2=reform(lwr2) |
---|
| 228 | lwr3=-1*read_ncdf("lwr", st3, en3, file=file,/nostr, box=box) & lwr3=grossemoyenne(lwr3, "xy",/nan, box=box) & lwr3=reform(lwr3) |
---|
| 229 | |
---|
[46] | 230 | file=tropflux_id_env+'lhf_oafluxgrid_1985_2009.nc' |
---|
[41] | 231 | initncdf, file |
---|
| 232 | lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) & lhf1=reform(lhf1) |
---|
| 233 | lhf2=-1*read_ncdf("lhf", st2, en2, file=file,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) & lhf2=reform(lhf2) |
---|
| 234 | lhf3=-1*read_ncdf("lhf", st3, en3, file=file,/nostr, box=box) & lhf3=grossemoyenne(lhf3, "xy",/nan, box=box) & lhf3=reform(lhf3) |
---|
| 235 | |
---|
[46] | 236 | file=tropflux_id_env+'shf_oafluxgrid_1985_2009.nc' |
---|
[41] | 237 | initncdf, file |
---|
| 238 | shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy",/nan, box=box) & shf1=reform(shf1) |
---|
| 239 | shf2=-1*read_ncdf("shf", st2, en2, file=file,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy",/nan, box=box) & shf2=reform(shf2) |
---|
| 240 | shf3=-1*read_ncdf("shf", st3, en3, file=file,/nostr, box=box) & shf3=grossemoyenne(shf3, "xy",/nan, box=box) & shf3=reform(shf3) |
---|
| 241 | |
---|
| 242 | net1=swr1+lwr1+lhf1+shf1 & net2=swr2+lwr2+lhf2+shf2 & net3=swr3+lwr3+lhf3+shf3 |
---|
| 243 | |
---|
| 244 | oaf=[net1,net2,net3] |
---|
[43] | 245 | lhf1=0. & lhf2=0. & lhf3=0. |
---|
| 246 | shf1=0. & shf2=0. & shf3=0. |
---|
[41] | 247 | lwr1=0. & lwr2=0. & lwr3=0. |
---|
| 248 | swr1=0. & swr2=0. & swr3=0. |
---|
| 249 | net1=0. & net2=0. & net3=0. |
---|
| 250 | help, oaf |
---|
| 251 | |
---|
[46] | 252 | fi=tropflux_id_env+'swr_ncep2_oaflxgrid_19890101_20091231.nc' |
---|
[41] | 253 | initncdf, fi |
---|
| 254 | swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy",/nan, box=box) & swr1=reform(swr1) |
---|
| 255 | swr2=read_ncdf("swr", st2, en2, file=fi,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy",/nan, box=box) & swr2=reform(swr2) |
---|
| 256 | swr3=read_ncdf("swr", st3, en3, file=fi,/nostr, box=box) & swr3=grossemoyenne(swr3, "xy",/nan, box=box) & swr3=reform(swr3) |
---|
| 257 | |
---|
[46] | 258 | fi=tropflux_id_env+'NCEP2_flux_19890101_20090729.nc" |
---|
[41] | 259 | initncdf, fi |
---|
| 260 | lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) & lwr1=reform(lwr1) |
---|
| 261 | lwr2=-1*read_ncdf("lwr", st2, en2, file=fi,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) & lwr2=reform(lwr2) |
---|
| 262 | lwr3=-1*read_ncdf("lwr", st3, en3, file=fi,/nostr, box=box) & lwr3=grossemoyenne(lwr3, "xy",/nan, box=box) & lwr3=reform(lwr3) |
---|
| 263 | |
---|
[46] | 264 | fi=tropflux_id_env+'lhf_ncep2_oafluxgrid_19890101_20091231.nc' |
---|
[41] | 265 | initncdf, fi |
---|
| 266 | lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) & lhf1=reform(lhf1) |
---|
| 267 | lhf2=-1*read_ncdf("lhf", st2, en2, file=fi,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) & lhf2=reform(lhf2) |
---|
| 268 | lhf3=-1*read_ncdf("lhf", st3, en3, file=fi,/nostr, box=box) & lhf3=grossemoyenne(lhf3, "xy",/nan, box=box) & lhf3=reform(lhf3) |
---|
| 269 | |
---|
[46] | 270 | fi=tropflux_id_env+'shf_ncep2_oafluxgrid_19890101_20091231.nc' |
---|
[41] | 271 | initncdf, fi |
---|
| 272 | shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy",/nan, box=box) & shf1=reform(shf1) |
---|
| 273 | shf2=-1*read_ncdf("shf", st2, en2, file=fi,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy",/nan, box=box) & shf2=reform(shf2) |
---|
| 274 | shf3=-1*read_ncdf("shf", st3, en3, file=fi,/nostr, box=box) & shf3=grossemoyenne(shf3, "xy",/nan, box=box) & shf3=reform(shf3) |
---|
| 275 | |
---|
| 276 | net1=swr1+lwr1+lhf1+shf1 & net2=swr2+lwr2+lhf2+shf2 & net3=swr3+lwr3+lhf3+shf3 |
---|
| 277 | |
---|
| 278 | nce2=[net1,net2,net3] |
---|
| 279 | help, nce2 |
---|
| 280 | |
---|
[43] | 281 | lhf1=0. & lhf2=0. & lhf3=0. |
---|
[41] | 282 | shf1=0. & shf2=0. & shf3=0. |
---|
| 283 | lwr1=0. & lwr2=0. & lwr3=0. |
---|
[43] | 284 | swr1=0. & swr2=0. & swr3=0. |
---|
[41] | 285 | net1=0. & net2=0. & net3=0. |
---|
| 286 | |
---|
[46] | 287 | file=tropflux_id_env+"ncep1_flux_19890101_20091231.nc" |
---|
[41] | 288 | initncdf, file |
---|
| 289 | |
---|
| 290 | lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) & lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) & lhf1=reform(lhf1) |
---|
| 291 | lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) & lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) & lhf2=reform(lhf2) |
---|
| 292 | lhf3=read_ncdf('lhf',st3,en3,file=file,/nostr, box=box) & lhf3=grossemoyenne(lhf3, "xy", box=box,/nan) & lhf3=reform(lhf3) |
---|
| 293 | |
---|
| 294 | shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) & shf1=grossemoyenne(shf1, "xy", box=box,/nan) & shf1=reform(shf1) |
---|
| 295 | shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box) & shf2=grossemoyenne(shf2, "xy", box=box,/nan) & shf2=reform(shf2) |
---|
| 296 | shf3=read_ncdf('shf',st3,en3,file=file,/nostr, box=box) & shf3=grossemoyenne(shf3, "xy", box=box,/nan) & shf3=reform(shf3) |
---|
| 297 | |
---|
| 298 | lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) & lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) & lwr1=reform(lwr1) |
---|
| 299 | lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box) & lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) & lwr2=reform(lwr2) |
---|
| 300 | lwr3=read_ncdf('lwr',st3,en3,file=file,/nostr, box=box) & lwr3=grossemoyenne(lwr3, "xy", box=box,/nan) & lwr3=reform(lwr3) |
---|
| 301 | |
---|
| 302 | swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) & swr1=grossemoyenne(swr1, "xy", box=box,/nan) & swr1=reform(swr1) |
---|
| 303 | swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) & swr2=grossemoyenne(swr2, "xy", box=box,/nan) & swr2=reform(swr2) |
---|
| 304 | swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box) & swr3=grossemoyenne(swr3, "xy", box=box,/nan) & swr3=reform(swr3) |
---|
| 305 | |
---|
| 306 | net1=swr1+(lwr1+shf1+lhf1) & net2=swr2+(lwr2+shf2+lhf2) & net3=swr3+(lwr3+shf3+lhf3) |
---|
| 307 | nce1=[net1,net2,net3] |
---|
| 308 | help, nce1 |
---|
| 309 | |
---|
[43] | 310 | lhf1=0. & lhf2=0. & lhf3=0. |
---|
[41] | 311 | shf1=0. & shf2=0. & shf3=0. |
---|
| 312 | lwr1=0. & lwr2=0. & lwr3=0. |
---|
[43] | 313 | swr1=0. & swr2=0. & swr3=0. |
---|
[41] | 314 | net1=0. & net2=0. & net3=0. |
---|
| 315 | |
---|
| 316 | ind=where(finite(net)) |
---|
[43] | 317 | net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) |
---|
[41] | 318 | nce2=nce2(ind) & nce1=nce1(ind) & olr=olr(ind) |
---|
| 319 | |
---|
| 320 | |
---|
| 321 | statistics, net, era, $ |
---|
| 322 | cor, bias, std, rmsd |
---|
| 323 | print, cor, bias, std, rmsd |
---|
| 324 | cstat_era=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
| 325 | |
---|
| 326 | statistics, net, trop, $ |
---|
| 327 | cor, bias, std, rmsd |
---|
| 328 | print, cor, bias, std, rmsd |
---|
| 329 | cstat_trop=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
| 330 | |
---|
| 331 | statistics, net, oaf, $ |
---|
| 332 | cor, bias, std, rmsd |
---|
| 333 | print, cor, bias, std, rmsd |
---|
| 334 | cstat_oaf=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
| 335 | |
---|
| 336 | statistics, net, nce1, $ |
---|
| 337 | cor, bias, std, rmsd |
---|
| 338 | print, cor, bias, std, rmsd |
---|
| 339 | cstat_nce1=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
| 340 | |
---|
| 341 | statistics, net, nce2, $ |
---|
| 342 | cor, bias, std, rmsd |
---|
| 343 | print, cor, bias, std, rmsd |
---|
| 344 | cstat_nce2=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
| 345 | |
---|
| 346 | statistics, net, olr, $ |
---|
| 347 | cor, bias, std, rmsd |
---|
| 348 | print, cor, bias, std, rmsd |
---|
| 349 | cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
| 350 | |
---|
| 351 | |
---|
| 352 | splot, net, era, charsize=1.1, title='Net Flux - WHOTS Vs ERAI', $ |
---|
| 353 | xrange=[-250,250], yrange=[-250,250], small=[2,3,1], psym=2, xmin=1, ymin=1 |
---|
| 354 | oplot, [-250,250], [-250,250] |
---|
| 355 | xyouts, 0,-200, cstat_era, charsize=1. |
---|
| 356 | |
---|
| 357 | splot, net, trop, charsize=1.1, title='Net Flux - WHOTS Vs TropFlux',/noer, $ |
---|
| 358 | xrange=[-250,250], yrange=[-250,250], small=[2,3,2], psym=2, xmin=1, ymin=1 |
---|
| 359 | oplot, [-250,250], [-250,250] |
---|
| 360 | xyouts, 0,-200, cstat_trop, charsize=1. |
---|
| 361 | |
---|
| 362 | splot, net, oaf, charsize=1.1, title='Net Flux - WHOTS Vs OAFlux',/noer, $ |
---|
| 363 | xrange=[-250,250], yrange=[-250,250], small=[2,3,3], psym=2, xmin=1, ymin=1 |
---|
| 364 | oplot, [-250,250], [-250,250] |
---|
| 365 | xyouts, 0,-200, cstat_oaf, charsize=1. |
---|
| 366 | |
---|
| 367 | splot, net, nce1, charsize=1.1, title='Net Flux - WHOTS Vs NCEP1',/noer, $ |
---|
| 368 | xrange=[-250,250], yrange=[-250,250], small=[2,3,4], psym=2, xmin=1, ymin=1 |
---|
| 369 | oplot, [-250,250], [-250,250] |
---|
| 370 | xyouts, 0,-200, cstat_nce1, charsize=1. |
---|
| 371 | |
---|
| 372 | splot, net, nce2, charsize=1.1, title='Net Flux - WHOTS Vs NCEP2',/noer, $ |
---|
| 373 | xrange=[-250,250], yrange=[-250,250], small=[2,3,5], psym=2, xmin=1, ymin=1 |
---|
| 374 | oplot, [-250,250], [-250,250] |
---|
| 375 | xyouts, 0,-200, cstat_nce2, charsize=1. |
---|
| 376 | |
---|
| 377 | splot, net, olr, charsize=1.1, title='Net Flux - WHOTS Vs TropFlux_NRT',/noer, $ |
---|
| 378 | xrange=[-250,250], yrange=[-250,250], small=[2,3,6], psym=2, xmin=1, ymin=1 |
---|
| 379 | oplot, [-250,250], [-250,250] |
---|
| 380 | xyouts, 0,-200, cstat_olr, charsize=1. |
---|
| 381 | |
---|
| 382 | closeps |
---|
[46] | 383 | fig=tropflux_od_env+"WHOTS_validation_net_2004_2007_all_v10.ps" |
---|
| 384 | spawn, 'mv '+psdir+'idl.ps '+fig |
---|
| 385 | spawn, 'gv '+fig |
---|
[41] | 386 | end |
---|