[43] | 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 | ; |
---|
[46] | 11 | ; .. graphviz:: |
---|
| 12 | ; |
---|
| 13 | ; digraph coare_validation_net_1992_93_all_v10 { |
---|
[85] | 14 | ; |
---|
[50] | 15 | ; met_data [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/met_data_as_expt_199*_p*.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 | ; shf_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_shf_19910101_20091231_oafluxgrid.nc"]; |
---|
| 20 | ; tropflux_general [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_19890101_20091231.nc"]; |
---|
| 21 | ; swr_tropflux_nrt [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_swr_19890101_20091231_NRT_v50.nc"]; |
---|
| 22 | ; swr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_oafluxgrid_1985_2007.nc"]; |
---|
| 23 | ; lwr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lwr_oafluxgrid_1985_2007.nc"]; |
---|
| 24 | ; lhf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_oafluxgrid_1985_2007.nc"]; |
---|
| 25 | ; shf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_oafluxgrid_1985_2007.nc"]; |
---|
| 26 | ; swr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_ncep2_oaflxgrid_19890101_20091231.nc"]; |
---|
| 27 | ; lwr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/NCEP2_flux_19890101_20090729.nc"]; |
---|
| 28 | ; lhf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_ncep2_oafluxgrid_19890101_20091231.nc"]; |
---|
| 29 | ; shf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_ncep2_oafluxgrid_19890101_20091231.nc"]; |
---|
[46] | 30 | ; |
---|
[50] | 31 | ; ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/ncep1_flux_19890101_20091231.nc"]; |
---|
[46] | 32 | ; |
---|
[50] | 33 | ; figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/coare_validation_net_1992_93_all_v11.ps"]; |
---|
[46] | 34 | ; |
---|
| 35 | ; coare_validation_net_1992_93_all_v10 [shape=box, |
---|
| 36 | ; fontname=Courier, |
---|
| 37 | ; color=blue, |
---|
| 38 | ; URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig14/coare_validation_net_1992_93_all_v10.pro", |
---|
| 39 | ; label="${TROPFLUX}/src/paper01/fig14/coare_validation_net_1992_93_all_v10.pro"]; |
---|
| 40 | ; |
---|
| 41 | ; {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} -> {coare_validation_net_1992_93_all_v10} -> {figure} |
---|
| 42 | ; } |
---|
| 43 | ; |
---|
[43] | 44 | ; SEE ALSO |
---|
| 45 | ; ======== |
---|
| 46 | ; |
---|
[50] | 47 | ; :ref:`project_profile.sh` |
---|
| 48 | ; :ref:`project_init.pro` |
---|
| 49 | ; :ref:`cm_project.pro` |
---|
[43] | 50 | ; |
---|
[46] | 51 | ; :ref:`statistics.pro` |
---|
| 52 | ; |
---|
[43] | 53 | ; EXAMPLES |
---|
| 54 | ; ======== |
---|
| 55 | ; |
---|
[203] | 56 | ; .. code-block:: idl |
---|
[43] | 57 | ; |
---|
[203] | 58 | ; coare_validation_net_1992_93_all_v10 |
---|
[43] | 59 | ; |
---|
[46] | 60 | ; TODO |
---|
| 61 | ; ==== |
---|
| 62 | ; |
---|
| 63 | ; make it work on cratos : missing data |
---|
| 64 | ; |
---|
| 65 | ; coding rules |
---|
| 66 | ; |
---|
| 67 | ; complete description |
---|
| 68 | ; |
---|
| 69 | ; handle IO error |
---|
| 70 | ; |
---|
[43] | 71 | ; EVOLUTIONS |
---|
| 72 | ; ========== |
---|
| 73 | ; |
---|
| 74 | ; $Id$ |
---|
| 75 | ; |
---|
[85] | 76 | ; $URL$ |
---|
| 77 | ; |
---|
[46] | 78 | ; - fplod 20110420T120449Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 79 | ; |
---|
| 80 | ; * remove hard coding path |
---|
| 81 | ; * add graphviz |
---|
| 82 | ; |
---|
[43] | 83 | ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 84 | ; |
---|
| 85 | ; * minimal header |
---|
| 86 | ; |
---|
| 87 | ;- |
---|
[41] | 88 | pro coare_validation_net_1992_93_all_v10 |
---|
[46] | 89 | @cm_general |
---|
[50] | 90 | @cm_project |
---|
[41] | 91 | reinitplt, /z,/invert |
---|
| 92 | key_portrait = 1 |
---|
| 93 | coefpalit=.9 |
---|
[97] | 94 | ; |
---|
[85] | 95 | openps, FILENAME = project_od_env+"coare_validation_net_1992_93_all_v11.ps" |
---|
[41] | 96 | ; partie a changer |
---|
[97] | 97 | ; |
---|
[41] | 98 | box=[155.5, 156.5, -2, -1] |
---|
[97] | 99 | ; |
---|
| 100 | ; first leg |
---|
| 101 | ; start data - 11/11/92 and end date 03/12/92 |
---|
[50] | 102 | fi=project_id_env+"met_data_coare_expt_1992_p1.txt" |
---|
[41] | 103 | res=read_ascii(fi,data_start=1) |
---|
| 104 | ff=res.field1 |
---|
| 105 | date=reform(ff(0,*)) |
---|
| 106 | sw92_1=reform(ff(7,*))/0.93 |
---|
| 107 | lw92_1=reform(ff(8,*)) |
---|
| 108 | lh92_1=reform(ff(5,*)) |
---|
| 109 | sh92_1=reform(ff(6,*)) |
---|
[97] | 110 | ; |
---|
[41] | 111 | net1=sw92_1+lw92_1+lh92_1+sh92_1 |
---|
[97] | 112 | ; |
---|
| 113 | st1=921111 |
---|
[94] | 114 | en1=921203 |
---|
[97] | 115 | ; |
---|
| 116 | ; second leg |
---|
| 117 | ; |
---|
[50] | 118 | fi=project_id_env+"met_data_coare_expt_1992_p2.txt" |
---|
[41] | 119 | res=read_ascii(fi,data_start=1) |
---|
| 120 | ff=res.field1 |
---|
| 121 | sw92_2=reform(ff(7,*))/0.93 |
---|
| 122 | lw92_2=reform(ff(8,*)) |
---|
| 123 | lh92_2=reform(ff(5,*)) |
---|
| 124 | sh92_2=reform(ff(6,*)) |
---|
[97] | 125 | ; |
---|
[41] | 126 | net2=sw92_2+lw92_2+lh92_2+sh92_2 |
---|
[97] | 127 | ; |
---|
| 128 | ; start data - 17/12/92 and end date 31/12/92 |
---|
| 129 | st2=19921217 |
---|
[94] | 130 | en2=19921231 |
---|
[97] | 131 | ; |
---|
| 132 | ; third leg |
---|
| 133 | ; |
---|
[50] | 134 | fi=project_id_env+"met_data_coare_expt_1993_p1.txt" |
---|
[41] | 135 | res=read_ascii(fi,data_start=1) |
---|
| 136 | ff=res.field1 |
---|
| 137 | sw93_1=reform(ff(7,*))/0.93 |
---|
| 138 | lw93_1=reform(ff(8,*)) |
---|
| 139 | lh93_1=reform(ff(5,*)) |
---|
| 140 | sh93_1=reform(ff(6,*)) |
---|
[97] | 141 | ; |
---|
[41] | 142 | net3=sw93_1+lw93_1+lh93_1+sh93_1 |
---|
[97] | 143 | ; |
---|
| 144 | ; start data - 01/01/93 and end date 11/01/93 |
---|
| 145 | st3=19930101 |
---|
[94] | 146 | en3=19930111 |
---|
[97] | 147 | ; |
---|
| 148 | ; fourth leg |
---|
| 149 | ; |
---|
[50] | 150 | fi=project_id_env+"met_data_coare_expt_1993_p2.txt" |
---|
[41] | 151 | res=read_ascii(fi,data_start=1) |
---|
| 152 | ff=res.field1 |
---|
| 153 | sw93_2=reform(ff(7,*)) |
---|
| 154 | lw93_2=reform(ff(8,*)) |
---|
| 155 | lh93_2=reform(ff(5,*)) |
---|
| 156 | sh93_2=reform(ff(6,*)) |
---|
[97] | 157 | ; |
---|
[41] | 158 | net4=sw93_2+lw93_2+lh93_2+sh93_2 |
---|
[97] | 159 | ; |
---|
| 160 | ; start data - 28/01/93 and end date 16/02/93 |
---|
| 161 | st4=19930128 |
---|
[94] | 162 | en4=19930216 |
---|
[97] | 163 | ; |
---|
[41] | 164 | net=[net1, net2, net3, net4] |
---|
| 165 | help, net |
---|
[97] | 166 | ; |
---|
| 167 | ; reading other data sets |
---|
| 168 | ; |
---|
| 169 | st1=19921111 |
---|
[94] | 170 | en1=19921203 |
---|
[97] | 171 | st2=19921217 |
---|
[94] | 172 | en2=19921231 |
---|
[97] | 173 | st3=19930101 |
---|
[94] | 174 | en3=19930111 |
---|
[97] | 175 | st4=19930128 |
---|
[94] | 176 | en4=19930216 |
---|
[97] | 177 | ; |
---|
[50] | 178 | file=project_id_env+"erai_swr_19910101_20091231_oafluxgrid.nc" |
---|
[41] | 179 | initncdf, file |
---|
[97] | 180 | swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) |
---|
| 181 | swr1=grossemoyenne(swr1, "xy", box=box,/nan) |
---|
[94] | 182 | swr1=reform(swr1) |
---|
[97] | 183 | swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) |
---|
| 184 | swr2=grossemoyenne(swr2, "xy", box=box,/nan) |
---|
[94] | 185 | swr2=reform(swr2) |
---|
[97] | 186 | swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box) |
---|
| 187 | swr3=grossemoyenne(swr3, "xy", box=box,/nan) |
---|
[94] | 188 | swr3=reform(swr3) |
---|
[97] | 189 | swr4=read_ncdf('swr',st4,en4,file=file,/nostr, box=box) |
---|
| 190 | swr4=grossemoyenne(swr4, "xy", box=box,/nan) |
---|
[94] | 191 | swr4=reform(swr4) |
---|
[97] | 192 | ; |
---|
[50] | 193 | file=project_id_env+'erai_lwr_19890101_20091231_oafluxgrid.nc' |
---|
[41] | 194 | initncdf, file |
---|
[97] | 195 | lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box) |
---|
| 196 | lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) |
---|
[94] | 197 | lwr1=reform(lwr1) |
---|
[97] | 198 | lwr2=read_ncdf('lwr',st2-1,en2,file=file,/nostr, box=box) |
---|
| 199 | lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) |
---|
[94] | 200 | lwr2=reform(lwr2) |
---|
[97] | 201 | lwr3=read_ncdf('lwr',st3-1,en3,file=file,/nostr, box=box) |
---|
| 202 | lwr3=grossemoyenne(lwr3, "xy", box=box,/nan) |
---|
[94] | 203 | lwr3=reform(lwr3) |
---|
[97] | 204 | lwr4=read_ncdf('lwr',st4-1,en4,file=file,/nostr, box=box) |
---|
| 205 | lwr4=grossemoyenne(lwr4, "xy", box=box,/nan) |
---|
[94] | 206 | lwr4=reform(lwr4) |
---|
[97] | 207 | ; |
---|
[50] | 208 | file=project_id_env+'erai_lhf_19890101_20091231_oafluxgrid.nc' |
---|
[41] | 209 | initncdf, file |
---|
[97] | 210 | lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) |
---|
| 211 | lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) |
---|
[94] | 212 | lhf1=reform(lhf1) |
---|
[97] | 213 | lhf2=-1*read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) |
---|
| 214 | lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) |
---|
[94] | 215 | lhf2=reform(lhf2) |
---|
[97] | 216 | lhf3=-1*read_ncdf('lhf',st3,en3,file=file,/nostr, box=box) |
---|
| 217 | lhf3=grossemoyenne(lhf3, "xy", box=box,/nan) |
---|
[94] | 218 | lhf3=reform(lhf3) |
---|
[97] | 219 | lhf4=-1*read_ncdf('lhf',st4,en4,file=file,/nostr, box=box) |
---|
| 220 | lhf4=grossemoyenne(lhf4, "xy", box=box,/nan) |
---|
[94] | 221 | lhf4=reform(lhf4) |
---|
[97] | 222 | ; |
---|
[50] | 223 | file=project_id_env+'erai_shf_19890101_20091231_oafluxgrid.nc' |
---|
[41] | 224 | initncdf, file |
---|
[97] | 225 | shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box) |
---|
| 226 | shf1=grossemoyenne(shf1, "xy", box=box,/nan) |
---|
[94] | 227 | shf1=reform(shf1) |
---|
[97] | 228 | shf2=-1*read_ncdf('shf',st2,en2,file=file,/nostr, box=box) |
---|
| 229 | shf2=grossemoyenne(shf2, "xy", box=box,/nan) |
---|
[94] | 230 | shf2=reform(shf2) |
---|
[97] | 231 | shf3=-1*read_ncdf('shf',st3,en3,file=file,/nostr, box=box) |
---|
| 232 | shf3=grossemoyenne(shf3, "xy", box=box,/nan) |
---|
[94] | 233 | shf3=reform(shf3) |
---|
[97] | 234 | shf4=-1*read_ncdf('shf',st4,en4,file=file,/nostr, box=box) |
---|
| 235 | shf4=grossemoyenne(shf4, "xy", box=box,/nan) |
---|
[94] | 236 | shf4=reform(shf4) |
---|
[97] | 237 | ; |
---|
| 238 | net1=swr1+lwr1+lhf1+shf1 |
---|
[94] | 239 | net2=swr2+lwr2+lhf2+shf2 |
---|
[97] | 240 | net3=swr3+lwr3+lhf3+shf3 |
---|
[94] | 241 | net4=swr4+lwr4+lhf4+shf4 |
---|
[41] | 242 | era=[net1,net2,net3,net4] |
---|
[97] | 243 | ; |
---|
[41] | 244 | help, era |
---|
[97] | 245 | swr1=0 |
---|
| 246 | swr2=0 |
---|
| 247 | swr3=0 |
---|
| 248 | swr4=0 |
---|
| 249 | lwr1=0 |
---|
| 250 | lwr2=0 |
---|
| 251 | lwr3=0 |
---|
[94] | 252 | lwr4=0 |
---|
[97] | 253 | lhf1=0 |
---|
| 254 | lhf2=0 |
---|
| 255 | lhf3=0 |
---|
| 256 | lhf4=0 |
---|
| 257 | shf1=0 |
---|
| 258 | shf2=0 |
---|
| 259 | shf3=0 |
---|
[94] | 260 | shf4=0 |
---|
[97] | 261 | net1=0 |
---|
| 262 | net2=0 |
---|
| 263 | net3=0 |
---|
[94] | 264 | net4=0 |
---|
[97] | 265 | ; |
---|
[50] | 266 | file=project_id_env+"TropFlux_19890101_20091231.nc" |
---|
[41] | 267 | initncdf, file |
---|
[97] | 268 | swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) |
---|
| 269 | swr1=grossemoyenne(swr1, "xy",/nan, box=box) |
---|
[94] | 270 | swr1=reform(swr1) |
---|
[97] | 271 | swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) |
---|
| 272 | swr2=grossemoyenne(swr2, "xy",/nan, box=box) |
---|
[94] | 273 | swr2=reform(swr2) |
---|
[97] | 274 | swr3=read_ncdf('swr',st3,en3,file=file,/nostr, box=box) |
---|
| 275 | swr3=grossemoyenne(swr3, "xy",/nan, box=box) |
---|
[94] | 276 | swr3=reform(swr3) |
---|
[97] | 277 | swr4=read_ncdf('swr',st4,en4,file=file,/nostr, box=box) |
---|
| 278 | swr4=grossemoyenne(swr4, "xy",/nan, box=box) |
---|
[94] | 279 | swr4=reform(swr4) |
---|
[97] | 280 | ; |
---|
| 281 | lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) |
---|
| 282 | lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) |
---|
[94] | 283 | lwr1=reform(lwr1) |
---|
[97] | 284 | lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box) |
---|
| 285 | lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) |
---|
[94] | 286 | lwr2=reform(lwr2) |
---|
[97] | 287 | lwr3=read_ncdf('lwr',st3,en3,file=file,/nostr, box=box) |
---|
| 288 | lwr3=grossemoyenne(lwr3, "xy",/nan, box=box) |
---|
[94] | 289 | lwr3=reform(lwr3) |
---|
[97] | 290 | lwr4=read_ncdf('lwr',st4,en4,file=file,/nostr, box=box) |
---|
| 291 | lwr4=grossemoyenne(lwr4, "xy",/nan, box=box) |
---|
[94] | 292 | lwr4=reform(lwr4) |
---|
[97] | 293 | ; |
---|
| 294 | lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) |
---|
| 295 | lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) |
---|
[94] | 296 | lhf1=reform(lhf1) |
---|
[97] | 297 | lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) |
---|
| 298 | lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) |
---|
[94] | 299 | lhf2=reform(lhf2) |
---|
[97] | 300 | lhf3=read_ncdf('lhf',st3,en3,file=file,/nostr, box=box) |
---|
| 301 | lhf3=grossemoyenne(lhf3, "xy",/nan, box=box) |
---|
[94] | 302 | lhf3=reform(lhf3) |
---|
[97] | 303 | lhf4=read_ncdf('lhf',st4,en4,file=file,/nostr, box=box) |
---|
| 304 | lhf4=grossemoyenne(lhf4, "xy",/nan, box=box) |
---|
[94] | 305 | lhf4=reform(lhf4) |
---|
[97] | 306 | ; |
---|
| 307 | shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) |
---|
| 308 | shf1=grossemoyenne(shf1, "xy",/nan, box=box) |
---|
[94] | 309 | shf1=reform(shf1) |
---|
[97] | 310 | shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box) |
---|
| 311 | shf2=grossemoyenne(shf2, "xy",/nan, box=box) |
---|
[94] | 312 | shf2=reform(shf2) |
---|
[97] | 313 | shf3=read_ncdf('shf',st3,en3,file=file,/nostr, box=box) |
---|
| 314 | shf3=grossemoyenne(shf3, "xy",/nan, box=box) |
---|
[94] | 315 | shf3=reform(shf3) |
---|
[97] | 316 | shf4=read_ncdf('shf',st4,en4,file=file,/nostr, box=box) |
---|
| 317 | shf4=grossemoyenne(shf4, "xy",/nan, box=box) |
---|
[94] | 318 | shf4=reform(shf4) |
---|
[97] | 319 | ; |
---|
| 320 | net1=swr1+lwr1+lhf1+shf1 |
---|
[94] | 321 | net2=swr2+lwr2+lhf2+shf2 |
---|
[97] | 322 | net3=swr3+lwr3+lhf3+shf3 |
---|
[94] | 323 | net4=swr4+lwr4+lhf4+shf4 |
---|
[97] | 324 | ; |
---|
[50] | 325 | file=project_id_env+"TropFlux_swr_19890101_20091231_NRT_v50.nc" |
---|
[41] | 326 | initncdf, file |
---|
[97] | 327 | olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box) |
---|
| 328 | olr1=grossemoyenne(olr1, "xy",/nan, box=box) |
---|
[94] | 329 | olr1=reform(olr1) |
---|
[97] | 330 | olr2=read_ncdf("sw", st2-1, en2, file=file,/nostr, box=box) |
---|
| 331 | olr2=grossemoyenne(olr2, "xy",/nan, box=box) |
---|
[94] | 332 | olr2=reform(olr2) |
---|
[97] | 333 | olr3=read_ncdf("sw", st3-1, en3, file=file,/nostr, box=box) |
---|
| 334 | olr3=grossemoyenne(olr3, "xy",/nan, box=box) |
---|
[94] | 335 | olr3=reform(olr3) |
---|
[97] | 336 | olr4=read_ncdf("sw", st4-1, en4, file=file,/nostr, box=box) |
---|
| 337 | olr4=grossemoyenne(olr4, "xy",/nan, box=box) |
---|
[94] | 338 | olr4=reform(olr4) |
---|
[97] | 339 | ; |
---|
| 340 | neto_1=olr1+lwr1+lhf1+shf1 |
---|
[94] | 341 | neto_2=olr2+lwr2+lhf2+shf2 |
---|
[97] | 342 | neto_3=olr3+lwr3+lhf3+shf3 |
---|
[94] | 343 | neto_4=olr4+lwr4+lhf4+shf4 |
---|
[97] | 344 | ; |
---|
[41] | 345 | trop=[net1,net2,net3,net4] |
---|
| 346 | olr=[neto_1,neto_2,neto_3,neto_4] |
---|
| 347 | help, trop, olr |
---|
[97] | 348 | ; |
---|
| 349 | swr1=0 |
---|
| 350 | swr2=0 |
---|
| 351 | swr3=0 |
---|
| 352 | swr4=0 |
---|
| 353 | lwr1=0 |
---|
| 354 | lwr2=0 |
---|
| 355 | lwr3=0 |
---|
[94] | 356 | lwr4=0 |
---|
[97] | 357 | lhf1=0 |
---|
| 358 | lhf2=0 |
---|
| 359 | lhf3=0 |
---|
| 360 | lhf4=0 |
---|
| 361 | shf1=0 |
---|
| 362 | shf2=0 |
---|
| 363 | shf3=0 |
---|
[94] | 364 | shf4=0 |
---|
[97] | 365 | net1=0 |
---|
| 366 | net2=0 |
---|
| 367 | net3=0 |
---|
| 368 | net4=0 |
---|
| 369 | neto_1=0 |
---|
| 370 | neto_2=0 |
---|
| 371 | neto_3=0 |
---|
[94] | 372 | neto_4=0 |
---|
[97] | 373 | ; |
---|
[50] | 374 | file=project_id_env+'swr_oafluxgrid_1985_2007.nc' |
---|
[41] | 375 | initncdf, file |
---|
[97] | 376 | swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box) |
---|
| 377 | swr1=grossemoyenne(swr1, "xy",/nan, box=box) |
---|
[94] | 378 | swr1=reform(swr1) |
---|
[97] | 379 | swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box) |
---|
| 380 | swr2=grossemoyenne(swr2, "xy",/nan, box=box) |
---|
[94] | 381 | swr2=reform(swr2) |
---|
[97] | 382 | swr3=read_ncdf("swr", st3, en3, file=file,/nostr, box=box) |
---|
| 383 | swr3=grossemoyenne(swr3, "xy",/nan, box=box) |
---|
[94] | 384 | swr3=reform(swr3) |
---|
[97] | 385 | swr4=read_ncdf("swr", st4, en4, file=file,/nostr, box=box) |
---|
| 386 | swr4=grossemoyenne(swr4, "xy",/nan, box=box) |
---|
[94] | 387 | swr4=reform(swr4) |
---|
[97] | 388 | ; |
---|
[50] | 389 | file=project_id_env+'lwr_oafluxgrid_1985_2007.nc' |
---|
[41] | 390 | initncdf, file |
---|
[97] | 391 | lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box) |
---|
| 392 | lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) |
---|
[94] | 393 | lwr1=reform(lwr1) |
---|
[97] | 394 | lwr2=-1*read_ncdf("lwr", st2, en2, file=file,/nostr, box=box) |
---|
| 395 | lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) |
---|
[94] | 396 | lwr2=reform(lwr2) |
---|
[97] | 397 | lwr3=-1*read_ncdf("lwr", st3, en3, file=file,/nostr, box=box) |
---|
| 398 | lwr3=grossemoyenne(lwr3, "xy",/nan, box=box) |
---|
[94] | 399 | lwr3=reform(lwr3) |
---|
[97] | 400 | lwr4=-1*read_ncdf("lwr", st4, en4, file=file,/nostr, box=box) |
---|
| 401 | lwr4=grossemoyenne(lwr4, "xy",/nan, box=box) |
---|
[94] | 402 | lwr4=reform(lwr4) |
---|
[97] | 403 | ; |
---|
[50] | 404 | file=project_id_env+'lhf_oafluxgrid_1985_2009.nc' |
---|
[41] | 405 | initncdf, file |
---|
[97] | 406 | lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box) |
---|
| 407 | lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) |
---|
[94] | 408 | lhf1=reform(lhf1) |
---|
[97] | 409 | lhf2=-1*read_ncdf("lhf", st2, en2, file=file,/nostr, box=box) |
---|
| 410 | lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) |
---|
[94] | 411 | lhf2=reform(lhf2) |
---|
[97] | 412 | lhf3=-1*read_ncdf("lhf", st3, en3, file=file,/nostr, box=box) |
---|
[94] | 413 | lhf3=grossemoyenne(lhf3, "xy",/nan, box=box) |
---|
| 414 | lhf3=reform(lhf3) |
---|
[97] | 415 | lhf4=-1*read_ncdf("lhf", st4, en4, file=file,/nostr, box=box) |
---|
| 416 | lhf4=grossemoyenne(lhf4, "xy",/nan, box=box) |
---|
[94] | 417 | lhf4=reform(lhf4) |
---|
[97] | 418 | ; |
---|
[50] | 419 | file=project_id_env+'shf_oafluxgrid_1985_2009.nc' |
---|
[41] | 420 | initncdf, file |
---|
[97] | 421 | shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box) |
---|
| 422 | shf1=grossemoyenne(shf1, "xy",/nan, box=box) |
---|
[94] | 423 | shf1=reform(shf1) |
---|
[97] | 424 | shf2=-1*read_ncdf("shf", st2, en2, file=file,/nostr, box=box) |
---|
| 425 | shf2=grossemoyenne(shf2, "xy",/nan, box=box) |
---|
[94] | 426 | shf2=reform(shf2) |
---|
[97] | 427 | shf3=-1*read_ncdf("shf", st3, en3, file=file,/nostr, box=box) |
---|
| 428 | shf3=grossemoyenne(shf3, "xy",/nan, box=box) |
---|
[94] | 429 | shf3=reform(shf3) |
---|
[97] | 430 | shf4=-1*read_ncdf("shf", st4, en4, file=file,/nostr, box=box) |
---|
| 431 | shf4=grossemoyenne(shf4, "xy",/nan, box=box) |
---|
[94] | 432 | shf4=reform(shf4) |
---|
[97] | 433 | ; |
---|
| 434 | net1=swr1+lwr1+lhf1+shf1 |
---|
[94] | 435 | net2=swr2+lwr2+lhf2+shf2 |
---|
[97] | 436 | net3=swr3+lwr3+lhf3+shf3 |
---|
[94] | 437 | net4=swr4+lwr4+lhf4+shf4 |
---|
[97] | 438 | ; |
---|
[41] | 439 | oaf=[net1,net2,net3,net4] |
---|
[97] | 440 | ; |
---|
[41] | 441 | help, oaf |
---|
[97] | 442 | swr1=0 |
---|
| 443 | swr2=0 |
---|
| 444 | swr3=0 |
---|
| 445 | swr4=0 |
---|
| 446 | lwr1=0 |
---|
| 447 | lwr2=0 |
---|
| 448 | lwr3=0 |
---|
[94] | 449 | lwr4=0 |
---|
[97] | 450 | lhf1=0 |
---|
| 451 | lhf2=0 |
---|
| 452 | lhf3=0 |
---|
| 453 | lhf4=0 |
---|
| 454 | shf1=0 |
---|
| 455 | shf2=0 |
---|
[94] | 456 | shf3=0 |
---|
| 457 | shf4=0 |
---|
[97] | 458 | net1=0 |
---|
| 459 | net2=0 |
---|
| 460 | net3=0 |
---|
[94] | 461 | net4=0 |
---|
[97] | 462 | ; |
---|
[50] | 463 | fi=project_id_env+'swr_ncep2_oaflxgrid_19890101_20091231.nc' |
---|
[41] | 464 | initncdf, fi |
---|
[97] | 465 | swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box) |
---|
| 466 | swr1=grossemoyenne(swr1, "xy",/nan, box=box) |
---|
[94] | 467 | swr1=reform(swr1) |
---|
[97] | 468 | swr2=read_ncdf("swr", st2, en2, file=fi,/nostr, box=box) |
---|
| 469 | swr2=grossemoyenne(swr2, "xy",/nan, box=box) |
---|
[94] | 470 | swr2=reform(swr2) |
---|
[97] | 471 | swr3=read_ncdf("swr", st3, en3, file=fi,/nostr, box=box) |
---|
| 472 | swr3=grossemoyenne(swr3, "xy",/nan, box=box) |
---|
[94] | 473 | swr3=reform(swr3) |
---|
[97] | 474 | swr4=read_ncdf("swr", st4, en4, file=fi,/nostr, box=box) |
---|
| 475 | swr4=grossemoyenne(swr4, "xy",/nan, box=box) |
---|
[94] | 476 | swr4=reform(swr4) |
---|
[97] | 477 | ; |
---|
[50] | 478 | fi=project_id_env+"NCEP2_flux_19890101_20090729.nc" |
---|
[41] | 479 | initncdf, fi |
---|
[97] | 480 | lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box) |
---|
| 481 | lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) |
---|
[94] | 482 | lwr1=reform(lwr1) |
---|
[97] | 483 | lwr2=-1*read_ncdf("lwr", st2, en2, file=fi,/nostr, box=box) |
---|
| 484 | lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) |
---|
[94] | 485 | lwr2=reform(lwr2) |
---|
[97] | 486 | lwr3=-1*read_ncdf("lwr", st3, en3, file=fi,/nostr, box=box) |
---|
| 487 | lwr3=grossemoyenne(lwr3, "xy",/nan, box=box) |
---|
[94] | 488 | lwr3=reform(lwr3) |
---|
[97] | 489 | lwr4=-1*read_ncdf("lwr", st4, en4, file=fi,/nostr, box=box) |
---|
| 490 | lwr4=grossemoyenne(lwr4, "xy",/nan, box=box) |
---|
[94] | 491 | lwr4=reform(lwr4) |
---|
[97] | 492 | ; |
---|
[50] | 493 | fi=project_id_env+'lhf_ncep2_oafluxgrid_19890101_20091231.nc' |
---|
[41] | 494 | initncdf, fi |
---|
[97] | 495 | lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box) |
---|
| 496 | lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) |
---|
[94] | 497 | lhf1=reform(lhf1) |
---|
[97] | 498 | lhf2=-1*read_ncdf("lhf", st2, en2, file=fi,/nostr, box=box) |
---|
| 499 | lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) |
---|
[94] | 500 | lhf2=reform(lhf2) |
---|
[97] | 501 | lhf3=-1*read_ncdf("lhf", st3, en3, file=fi,/nostr, box=box) |
---|
| 502 | lhf3=grossemoyenne(lhf3, "xy",/nan, box=box) |
---|
[94] | 503 | lhf3=reform(lhf3) |
---|
[97] | 504 | lhf4=-1*read_ncdf("lhf", st4, en4, file=fi,/nostr, box=box) |
---|
| 505 | lhf4=grossemoyenne(lhf4, "xy",/nan, box=box) |
---|
[94] | 506 | lhf4=reform(lhf4) |
---|
[97] | 507 | ; |
---|
[50] | 508 | fi=project_id_env+'shf_ncep2_oafluxgrid_19890101_20091231.nc' |
---|
[41] | 509 | initncdf, fi |
---|
[97] | 510 | shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box) |
---|
| 511 | shf1=grossemoyenne(shf1, "xy",/nan, box=box) |
---|
[94] | 512 | shf1=reform(shf1) |
---|
[97] | 513 | shf2=-1*read_ncdf("shf", st2, en2, file=fi,/nostr, box=box) |
---|
| 514 | shf2=grossemoyenne(shf2, "xy",/nan, box=box) |
---|
[94] | 515 | shf2=reform(shf2) |
---|
[97] | 516 | shf3=-1*read_ncdf("shf", st3, en3, file=fi,/nostr, box=box) |
---|
| 517 | shf3=grossemoyenne(shf3, "xy",/nan, box=box) |
---|
[94] | 518 | shf3=reform(shf3) |
---|
[97] | 519 | shf4=-1*read_ncdf("shf", st4, en4, file=fi,/nostr, box=box) |
---|
| 520 | shf4=grossemoyenne(shf4, "xy",/nan, box=box) |
---|
[94] | 521 | shf4=reform(shf4) |
---|
[97] | 522 | ; |
---|
| 523 | net1=swr1+lwr1+lhf1+shf1 |
---|
[94] | 524 | net2=swr2+lwr2+lhf2+shf2 |
---|
[97] | 525 | net3=swr3+lwr3+lhf3+shf3 |
---|
[94] | 526 | net4=swr4+lwr4+lhf4+shf4 |
---|
[97] | 527 | ; |
---|
[41] | 528 | nce2=[net1,net2,net3,net4] |
---|
[97] | 529 | ; |
---|
[41] | 530 | help, nce2 |
---|
[97] | 531 | swr1=0 |
---|
| 532 | swr2=0 |
---|
| 533 | swr3=0 |
---|
| 534 | swr4=0 |
---|
| 535 | lwr1=0 |
---|
| 536 | lwr2=0 |
---|
| 537 | lwr3=0 |
---|
[94] | 538 | lwr4=0 |
---|
[97] | 539 | lhf1=0 |
---|
| 540 | lhf2=0 |
---|
| 541 | lhf3=0 |
---|
| 542 | lhf4=0 |
---|
| 543 | shf1=0 |
---|
| 544 | shf2=0 |
---|
| 545 | shf3=0 |
---|
[94] | 546 | shf4=0 |
---|
[97] | 547 | net1=0 |
---|
| 548 | net2=0 |
---|
| 549 | net3=0 |
---|
[94] | 550 | net4=0 |
---|
[97] | 551 | ; |
---|
[50] | 552 | file=project_id_env+"ncep1_flux_19890101_20091231.nc" |
---|
[41] | 553 | initncdf, file |
---|
[97] | 554 | swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box) |
---|
| 555 | swr1=grossemoyenne(swr1, "xy",/nan, box=box) |
---|
[94] | 556 | swr1=reform(swr1) |
---|
[97] | 557 | swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box) |
---|
| 558 | swr2=grossemoyenne(swr2, "xy",/nan, box=box) |
---|
[94] | 559 | swr2=reform(swr2) |
---|
[97] | 560 | swr3=read_ncdf("swr", st3, en3, file=file,/nostr, box=box) |
---|
| 561 | swr3=grossemoyenne(swr3, "xy",/nan, box=box) |
---|
[94] | 562 | swr3=reform(swr3) |
---|
[97] | 563 | swr4=read_ncdf("swr", st4, en4, file=file,/nostr, box=box) |
---|
| 564 | swr4=grossemoyenne(swr4, "xy",/nan, box=box) |
---|
[94] | 565 | swr4=reform(swr4) |
---|
[97] | 566 | ; |
---|
| 567 | lwr1=read_ncdf("lwr", st1, en1, file=file,/nostr, box=box) |
---|
| 568 | lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) |
---|
[94] | 569 | lwr1=reform(lwr1) |
---|
[97] | 570 | lwr2=read_ncdf("lwr", st2, en2, file=file,/nostr, box=box) |
---|
| 571 | lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) |
---|
[94] | 572 | lwr2=reform(lwr2) |
---|
[97] | 573 | lwr3=read_ncdf("lwr", st3, en3, file=file,/nostr, box=box) |
---|
| 574 | lwr3=grossemoyenne(lwr3, "xy",/nan, box=box) |
---|
[94] | 575 | lwr3=reform(lwr3) |
---|
[97] | 576 | lwr4=read_ncdf("lwr", st4, en4, file=file,/nostr, box=box) |
---|
| 577 | lwr4=grossemoyenne(lwr4, "xy",/nan, box=box) |
---|
[94] | 578 | lwr4=reform(lwr4) |
---|
[97] | 579 | ; |
---|
| 580 | lhf1=read_ncdf("lhf", st1, en1, file=file,/nostr, box=box) |
---|
| 581 | lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) |
---|
[94] | 582 | lhf1=reform(lhf1) |
---|
[97] | 583 | lhf2=read_ncdf("lhf", st2, en2, file=file,/nostr, box=box) |
---|
| 584 | lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) |
---|
[94] | 585 | lhf2=reform(lhf2) |
---|
[97] | 586 | lhf3=read_ncdf("lhf", st3, en3, file=file,/nostr, box=box) |
---|
| 587 | lhf3=grossemoyenne(lhf3, "xy",/nan, box=box) |
---|
[94] | 588 | lhf3=reform(lhf3) |
---|
[97] | 589 | lhf4=read_ncdf("lhf", st4, en4, file=file,/nostr, box=box) |
---|
| 590 | lhf4=grossemoyenne(lhf4, "xy",/nan, box=box) |
---|
[94] | 591 | lhf4=reform(lhf4) |
---|
[97] | 592 | ; |
---|
| 593 | shf1=read_ncdf("shf", st1, en1, file=file,/nostr, box=box) |
---|
| 594 | shf1=grossemoyenne(shf1, "xy",/nan, box=box) |
---|
[94] | 595 | shf1=reform(shf1) |
---|
[97] | 596 | shf2=read_ncdf("shf", st2, en2, file=file,/nostr, box=box) |
---|
| 597 | shf2=grossemoyenne(shf2, "xy",/nan, box=box) |
---|
[94] | 598 | shf2=reform(shf2) |
---|
[97] | 599 | shf3=read_ncdf("shf", st3, en3, file=file,/nostr, box=box) |
---|
| 600 | shf3=grossemoyenne(shf3, "xy",/nan, box=box) |
---|
[94] | 601 | shf3=reform(shf3) |
---|
[97] | 602 | shf4=read_ncdf("shf", st4, en4, file=file,/nostr, box=box) |
---|
| 603 | shf4=grossemoyenne(shf4, "xy",/nan, box=box) |
---|
[94] | 604 | shf4=reform(shf4) |
---|
[97] | 605 | ; |
---|
| 606 | net1=swr1+lwr1+lhf1+shf1 |
---|
[94] | 607 | net2=swr2+lwr2+lhf2+shf2 |
---|
[97] | 608 | net3=swr3+lwr3+lhf3+shf3 |
---|
[94] | 609 | net4=swr4+lwr4+lhf4+shf4 |
---|
[97] | 610 | ; |
---|
| 611 | nce1=[net1,net2,net3,net4] |
---|
[94] | 612 | nce1=nce1 |
---|
[97] | 613 | ; |
---|
[41] | 614 | help, nce1 |
---|
[97] | 615 | swr1=0 |
---|
| 616 | swr2=0 |
---|
| 617 | swr3=0 |
---|
| 618 | swr4=0 |
---|
| 619 | lwr1=0 |
---|
| 620 | lwr2=0 |
---|
| 621 | lwr3=0 |
---|
[94] | 622 | lwr4=0 |
---|
[97] | 623 | lhf1=0 |
---|
| 624 | lhf2=0 |
---|
| 625 | lhf3=0 |
---|
| 626 | lhf4=0 |
---|
| 627 | shf1=0 |
---|
| 628 | shf2=0 |
---|
| 629 | shf3=0 |
---|
[94] | 630 | shf4=0 |
---|
[97] | 631 | net1=0 |
---|
| 632 | net2=0 |
---|
| 633 | net3=0 |
---|
[94] | 634 | net4=0 |
---|
[97] | 635 | ; |
---|
[41] | 636 | ind=where(finite(net)) |
---|
[97] | 637 | net=net(ind) |
---|
[94] | 638 | era=era(ind) |
---|
[97] | 639 | trop=trop(ind) |
---|
| 640 | oaf=oaf(ind) |
---|
| 641 | nce1=nce1(ind) |
---|
[94] | 642 | nce2=nce2(ind) |
---|
[41] | 643 | olr=olr(ind) |
---|
[97] | 644 | ; |
---|
[41] | 645 | help, net, era, trop, oaf, nce1, nce2, olr |
---|
[97] | 646 | ; |
---|
[41] | 647 | statistics, net, era, $ |
---|
| 648 | cor, bias, std, rmsd |
---|
| 649 | print, cor, bias, std, rmsd |
---|
| 650 | cstat_era=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 651 | ; |
---|
[41] | 652 | statistics, net, trop, $ |
---|
| 653 | cor, bias, std, rmsd |
---|
| 654 | print, cor, bias, std, rmsd |
---|
| 655 | cstat_trop=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 656 | ; |
---|
[41] | 657 | statistics, net, oaf, $ |
---|
| 658 | cor, bias, std, rmsd |
---|
| 659 | print, cor, bias, std, rmsd |
---|
| 660 | cstat_oaf=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 661 | ; |
---|
[41] | 662 | statistics, net, nce1, $ |
---|
| 663 | cor, bias, std, rmsd |
---|
| 664 | print, cor, bias, std, rmsd |
---|
| 665 | cstat_nce1=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 666 | ; |
---|
[41] | 667 | statistics, net, nce2, $ |
---|
| 668 | cor, bias, std, rmsd |
---|
| 669 | print, cor, bias, std, rmsd |
---|
| 670 | cstat_nce2=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 671 | ; |
---|
[41] | 672 | statistics, net, olr, $ |
---|
| 673 | cor, bias, std, rmsd |
---|
| 674 | print, cor, bias, std, rmsd |
---|
| 675 | cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)') |
---|
[97] | 676 | ; |
---|
[41] | 677 | splot, net, era, charsize=1.1, title='Net Flux - COARE Vs ERAI', $ |
---|
| 678 | xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1 |
---|
| 679 | oplot, [-200,200], [-200,200] |
---|
| 680 | xyouts, -150,-175, cstat_era, charsize=1. |
---|
[97] | 681 | ; |
---|
[41] | 682 | splot, net, trop, charsize=1.1, title='Net Flux - COARE Vs TropFlux',/noer, $ |
---|
| 683 | xrange=[-200,200], yrange=[-200,200], small=[2,3,2], psym=2, xmin=1, ymin=1 |
---|
| 684 | oplot, [-200,200], [-200,200] |
---|
| 685 | xyouts, -150,-175, cstat_trop, charsize=1. |
---|
[97] | 686 | ; |
---|
[41] | 687 | splot, net, oaf, charsize=1.1, title='Net Flux - COARE Vs OAFlux',/noer, $ |
---|
| 688 | xrange=[-200,200], yrange=[-200,200], small=[2,3,3], psym=2, xmin=1, ymin=1 |
---|
| 689 | oplot, [-200,200], [-200,200] |
---|
| 690 | xyouts, -150,-175, cstat_oaf, charsize=1. |
---|
[97] | 691 | ; |
---|
[41] | 692 | splot, net, nce1, charsize=1.1, title='Net Flux - COARE Vs NCEP1',/noer, $ |
---|
| 693 | xrange=[-200,200], yrange=[-200,200], small=[2,3,4], psym=2, xmin=1, ymin=1 |
---|
| 694 | oplot, [-200,200], [-200,200] |
---|
| 695 | xyouts, -150,-175, cstat_nce1, charsize=1. |
---|
[97] | 696 | ; |
---|
[41] | 697 | splot, net, nce2, charsize=1.1, title='Net Flux - COARE Vs NCEP2',/noer, $ |
---|
| 698 | xrange=[-200,200], yrange=[-200,200], small=[2,3,5], psym=2, xmin=1, ymin=1 |
---|
| 699 | oplot, [-200,200], [-200,200] |
---|
| 700 | xyouts, -150,-175, cstat_nce2, charsize=1. |
---|
[97] | 701 | ; |
---|
[41] | 702 | splot, net, olr, charsize=1.1, title='Net Flux - COARE Vs TropFlux_NRT',/noer, $ |
---|
| 703 | xrange=[-200,200], yrange=[-200,200], small=[2,3,6], psym=2, xmin=1, ymin=1 |
---|
| 704 | oplot, [-200,200], [-200,200] |
---|
| 705 | xyouts, -150,-175, cstat_olr, charsize=1. |
---|
[97] | 706 | ; |
---|
[41] | 707 | closeps |
---|
[97] | 708 | ; |
---|
[41] | 709 | end |
---|