[43] | 1 | ;+ |
---|
[46] | 2 | ; .. _swr_statistics_map_2000_2009_v50.pro: |
---|
[43] | 3 | ; |
---|
| 4 | ; ==================================== |
---|
| 5 | ; swr_statistics_map_2000_2009_v50.pro |
---|
| 6 | ; ==================================== |
---|
| 7 | ; |
---|
| 8 | ; DESCRIPTION |
---|
| 9 | ; =========== |
---|
| 10 | ; |
---|
[46] | 11 | ; .. graphviz:: |
---|
| 12 | ; |
---|
| 13 | ; digraph swr_statistics_map_2000_2009_v50 { |
---|
[85] | 14 | ; |
---|
[50] | 15 | ; swr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_2000_2009_erai_v50.txt"]; |
---|
| 16 | ; swr_trop [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_2000_2009_trop_v50.txt"]; |
---|
| 17 | ; swr_oaflx [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_2000_2009_oaflx_v50.txt"]; |
---|
| 18 | ; swr_olr [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_2000_2009_olr_v50.txt"]; |
---|
| 19 | ; mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/longwave_IO_mask.nc"]; |
---|
[46] | 20 | ; |
---|
[50] | 21 | ; figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/swr_statistics_map_2000_2009_v50.ps"]; |
---|
[46] | 22 | ; |
---|
| 23 | ; swr_statistics_map_2000_2009_v50 [shape=box, |
---|
| 24 | ; fontname=Courier, |
---|
| 25 | ; color=blue, |
---|
| 26 | ; URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig12/swr_statistics_map_2000_2009_v50.pro", |
---|
| 27 | ; label="${TROPFLUX}/src/paper01/fig12/swr_statistics_map_2000_2009_v50.pro"]; |
---|
| 28 | ; |
---|
| 29 | ; {swr_erai swr_trop swr_oaflx swr_olr mask} -> {swr_statistics_map_2000_2009_v50} -> {figure} |
---|
| 30 | ; } |
---|
| 31 | ; |
---|
[43] | 32 | ; SEE ALSO |
---|
| 33 | ; ======== |
---|
| 34 | ; |
---|
[50] | 35 | ; :ref:`project_profile.sh` |
---|
| 36 | ; :ref:`project_init.pro` |
---|
| 37 | ; :ref:`cm_project.pro` |
---|
[43] | 38 | ; |
---|
| 39 | ; EXAMPLES |
---|
| 40 | ; ======== |
---|
| 41 | ; |
---|
| 42 | ; :: |
---|
| 43 | ; |
---|
[180] | 44 | ; swr_statistics_map_2000_2009_v50 |
---|
[43] | 45 | ; |
---|
[46] | 46 | ; TODO |
---|
| 47 | ; ==== |
---|
| 48 | ; |
---|
| 49 | ; make it work on cratos : missing data |
---|
| 50 | ; |
---|
| 51 | ; coding rules |
---|
| 52 | ; |
---|
| 53 | ; complete description |
---|
| 54 | ; |
---|
| 55 | ; handle IO error |
---|
| 56 | ; |
---|
[43] | 57 | ; EVOLUTIONS |
---|
| 58 | ; ========== |
---|
| 59 | ; |
---|
| 60 | ; $Id$ |
---|
| 61 | ; |
---|
[85] | 62 | ; $URL$ |
---|
| 63 | ; |
---|
[46] | 64 | ; - fplod 20110420T113846Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 65 | ; |
---|
| 66 | ; * remove stop |
---|
| 67 | ; * remove hard coding path |
---|
| 68 | ; * add graphviz |
---|
| 69 | ; |
---|
[43] | 70 | ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 71 | ; |
---|
| 72 | ; * minimal header |
---|
| 73 | ; |
---|
| 74 | ;- |
---|
[41] | 75 | pro swr_statistics_map_2000_2009_v50 |
---|
[46] | 76 | @cm_general |
---|
[50] | 77 | @cm_project |
---|
[41] | 78 | reinitplt, /z,/invert |
---|
| 79 | key_portrait = 1 |
---|
| 80 | coefpalit=.9 |
---|
[97] | 81 | ; |
---|
[85] | 82 | openps, FILENAME = project_od_env+'swr_statistics_map_2000_2009_v50.ps' |
---|
[41] | 83 | ; partie a changer |
---|
[97] | 84 | bias_mi=-20 |
---|
| 85 | bias_ma=20 |
---|
[94] | 86 | bias_int=2 |
---|
[97] | 87 | std_mi=0.7 |
---|
| 88 | std_ma=1.31 |
---|
[94] | 89 | std_int=0.05 |
---|
[97] | 90 | rmsd_mi=10 |
---|
| 91 | rmsd_ma=30 |
---|
[94] | 92 | rmsd_int=1.25 |
---|
[97] | 93 | cor_mi=0.5 |
---|
| 94 | cor_ma=1. |
---|
[94] | 95 | cor_int=0.02 |
---|
[50] | 96 | fi_swr_erai=project_id_env+'swr_2000_2009_erai_v50.txt' |
---|
| 97 | fi_swr_trop=project_id_env+'swr_2000_2009_trop_v50.txt' |
---|
| 98 | fi_swr_oaflx=project_id_env+'swr_2000_2009_oaflx_v50.txt' |
---|
| 99 | fi_swr_olr=project_id_env+'swr_2000_2009_olr_v50.txt' |
---|
[97] | 100 | ; |
---|
[41] | 101 | res=read_ascii(fi_swr_erai,data_start=1) |
---|
| 102 | ff=res.field1 |
---|
| 103 | lat=reform(ff(0,*)) |
---|
| 104 | lon=reform(ff(1,*)) |
---|
| 105 | cor_era=reform(ff(2,*)) |
---|
| 106 | bias_era=reform(ff(3,*)) |
---|
| 107 | std_era=reform(ff(4,*)) |
---|
| 108 | rmsd_era=reform(ff(5,*)) |
---|
[97] | 109 | ; |
---|
| 110 | ind=where(bias_era ge bias_ma) |
---|
[94] | 111 | bias_era(ind)=bias_ma-0.5 |
---|
[97] | 112 | ind=where(rmsd_era ge rmsd_ma) |
---|
[94] | 113 | rmsd_era(ind)=rmsd_ma-0.5 |
---|
[97] | 114 | ; |
---|
[41] | 115 | res=read_ascii(fi_swr_trop,data_start=1) |
---|
| 116 | ff=res.field1 |
---|
| 117 | lat=reform(ff(0,*)) |
---|
| 118 | lon=reform(ff(1,*)) |
---|
| 119 | cor_trop=reform(ff(2,*)) |
---|
| 120 | bias_trop=reform(ff(3,*)) |
---|
| 121 | std_trop=reform(ff(4,*)) |
---|
| 122 | rmsd_trop=reform(ff(5,*)) |
---|
[97] | 123 | ind=where(rmsd_trop ge rmsd_ma) |
---|
[94] | 124 | rmsd_trop(ind)=rmsd_ma-0.5 |
---|
[97] | 125 | ; |
---|
[41] | 126 | res=read_ascii(fi_swr_oaflx,data_start=1) |
---|
| 127 | ff=res.field1 |
---|
| 128 | lat=reform(ff(0,*)) |
---|
| 129 | lon=reform(ff(1,*)) |
---|
| 130 | cor_oaflx=reform(ff(2,*)) |
---|
| 131 | bias_oaflx=reform(ff(3,*)) |
---|
| 132 | std_oaflx=reform(ff(4,*)) |
---|
| 133 | rmsd_oaflx=reform(ff(5,*)) |
---|
[97] | 134 | ind=where(rmsd_oaflx ge rmsd_ma) |
---|
[94] | 135 | rmsd_oaflx(ind)=rmsd_ma-0.5 |
---|
[97] | 136 | ; |
---|
[41] | 137 | res=read_ascii(fi_swr_olr,data_start=1) |
---|
| 138 | ff=res.field1 |
---|
| 139 | lat=reform(ff(0,*)) |
---|
| 140 | lon=reform(ff(1,*)) |
---|
| 141 | cor_olr=reform(ff(2,*)) |
---|
| 142 | bias_olr=reform(ff(3,*)) |
---|
| 143 | std_olr=reform(ff(4,*)) |
---|
| 144 | rmsd_olr=reform(ff(5,*)) |
---|
[97] | 145 | ind=where(rmsd_olr ge rmsd_ma) |
---|
[94] | 146 | rmsd_olr(ind)=rmsd_ma-0.5 |
---|
[97] | 147 | ; |
---|
| 148 | ;ind=where(std_olr ge std_ma) |
---|
[94] | 149 | ;std_olr(ind)=std_ma-0.01 |
---|
[50] | 150 | file=project_id_env+'longwave_IO_mask.nc' |
---|
[41] | 151 | initncdf, file |
---|
| 152 | domdef, 30,390,-30,30 |
---|
| 153 | msk=ncdf_lec(file,var='msk') |
---|
| 154 | ;box=[0,290,-30,30] |
---|
| 155 | marge=[0,0,-0.8,0] |
---|
| 156 | marge1=[0,0,-5,0] |
---|
| 157 | plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,box=box,lct=64,cor_mi, cor_ma,int=cor_int, $ |
---|
| 158 | title='1) SWR Correlation - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge |
---|
[97] | 159 | ; |
---|
[41] | 160 | NN=n_elements(lat) |
---|
[97] | 161 | ; |
---|
[41] | 162 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 163 | ; |
---|
[41] | 164 | for n=0,NN-1 do begin |
---|
| 165 | x=lon(n) |
---|
| 166 | y=lat(n) |
---|
| 167 | c=cor_trop(n) |
---|
[97] | 168 | cmi=cor_mi |
---|
| 169 | cma=cor_ma |
---|
[94] | 170 | dc=cma-cmi |
---|
[41] | 171 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 172 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 173 | endfor |
---|
[97] | 174 | ; |
---|
[41] | 175 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 176 | ; |
---|
[41] | 177 | for n=0,NN-1 do begin |
---|
| 178 | x=lon(n) |
---|
[43] | 179 | y=lat(n) |
---|
| 180 | c=cor_trop(n) |
---|
[97] | 181 | ; |
---|
[41] | 182 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 183 | endfor |
---|
| 184 | plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $ |
---|
| 185 | title='2) SWR Correlation - ERAI', subtitle='', small=[1,4,2],/noer,/rempl, marge=marge |
---|
[97] | 186 | ; |
---|
[41] | 187 | NN=n_elements(lat) |
---|
[97] | 188 | ; |
---|
[41] | 189 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 190 | ; |
---|
[41] | 191 | for n=0,NN-1 do begin |
---|
| 192 | x=lon(n) |
---|
| 193 | y=lat(n) |
---|
| 194 | c=cor_era(n) |
---|
[97] | 195 | cmi=cor_mi |
---|
| 196 | cma=cor_ma |
---|
[94] | 197 | dc=cma-cmi |
---|
[41] | 198 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 199 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 200 | endfor |
---|
[97] | 201 | ; |
---|
[41] | 202 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 203 | ; |
---|
[41] | 204 | for n=0,NN-1 do begin |
---|
| 205 | x=lon(n) |
---|
[43] | 206 | y=lat(n) |
---|
| 207 | c=cor_era(n) |
---|
[97] | 208 | ; |
---|
[41] | 209 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 210 | endfor |
---|
| 211 | plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $ |
---|
| 212 | title='3) SWR Correlation - OAFlux', subtitle='', small=[1,4,3],/noer,/rempl, marge=marge |
---|
[97] | 213 | ; |
---|
[41] | 214 | NN=n_elements(lat) |
---|
[97] | 215 | ; |
---|
[41] | 216 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 217 | ; |
---|
[41] | 218 | for n=0,NN-1 do begin |
---|
| 219 | x=lon(n) |
---|
| 220 | y=lat(n) |
---|
| 221 | c=cor_oaflx(n) |
---|
[97] | 222 | cmi=cor_mi |
---|
| 223 | cma=cor_ma |
---|
[94] | 224 | dc=cma-cmi |
---|
[41] | 225 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 226 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 227 | endfor |
---|
[97] | 228 | ; |
---|
[41] | 229 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 230 | ; |
---|
[41] | 231 | for n=0,NN-1 do begin |
---|
| 232 | x=lon(n) |
---|
| 233 | y=lat(n) |
---|
| 234 | c=cor_oaflx(n) |
---|
[97] | 235 | ; |
---|
[41] | 236 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 237 | endfor |
---|
| 238 | plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int, $ |
---|
| 239 | title='4) SWR Correlation - SWR_olr', subtitle='', small=[1,4,4],/noer,/rempl, marge=marge1 |
---|
[97] | 240 | ; |
---|
[41] | 241 | NN=n_elements(lat) |
---|
[97] | 242 | ; |
---|
[41] | 243 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 244 | ; |
---|
[41] | 245 | for n=0,NN-1 do begin |
---|
| 246 | x=lon(n) |
---|
| 247 | y=lat(n) |
---|
| 248 | c=cor_olr(n) |
---|
[97] | 249 | cmi=cor_mi |
---|
| 250 | cma=cor_ma |
---|
[94] | 251 | dc=cma-cmi |
---|
[41] | 252 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 253 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 254 | endfor |
---|
[97] | 255 | ; |
---|
[41] | 256 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 257 | ; |
---|
[41] | 258 | for n=0,NN-1 do begin |
---|
| 259 | x=lon(n) |
---|
| 260 | y=lat(n) |
---|
| 261 | c=cor_olr(n) |
---|
[97] | 262 | ; |
---|
[41] | 263 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 264 | endfor |
---|
[97] | 265 | ; |
---|
[41] | 266 | erase |
---|
[97] | 267 | ; |
---|
[41] | 268 | plt, -5+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int, marge=marge,/nocolorb, $ |
---|
| 269 | title='1) SWR Mean bias - TropFlux', subtitle='', small=[1,4,1],/rempl,format='(i3)' |
---|
[97] | 270 | ; |
---|
[41] | 271 | NN=n_elements(lat) |
---|
[97] | 272 | ; |
---|
[41] | 273 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 274 | ; |
---|
[41] | 275 | for n=0,NN-1 do begin |
---|
| 276 | x=lon(n) |
---|
| 277 | y=lat(n) |
---|
| 278 | c=bias_trop(n) |
---|
[97] | 279 | cmi=bias_mi |
---|
| 280 | cma=bias_ma |
---|
[94] | 281 | dc=cma-cmi |
---|
[41] | 282 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 283 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 284 | endfor |
---|
[97] | 285 | ; |
---|
[41] | 286 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 287 | ; |
---|
[41] | 288 | for n=0,NN-1 do begin |
---|
| 289 | x=lon(n) |
---|
[43] | 290 | y=lat(n) |
---|
| 291 | c=bias_trop(n) |
---|
[97] | 292 | ; |
---|
[41] | 293 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 294 | endfor |
---|
| 295 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int,/noer, marge=marge, $ |
---|
| 296 | title='2) SWR Mean bias - ERAI', subtitle='', small=[1,4,2],/rempl,format='(i3)',/nocolorb |
---|
[97] | 297 | ; |
---|
[41] | 298 | NN=n_elements(lat) |
---|
[97] | 299 | ; |
---|
[41] | 300 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 301 | ; |
---|
[41] | 302 | for n=0,NN-1 do begin |
---|
| 303 | x=lon(n) |
---|
| 304 | y=lat(n) |
---|
| 305 | c=bias_era(n) |
---|
[97] | 306 | cmi=bias_mi |
---|
| 307 | cma=bias_ma |
---|
[94] | 308 | dc=cma-cmi |
---|
[41] | 309 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 310 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 311 | endfor |
---|
[97] | 312 | ; |
---|
[41] | 313 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 314 | ; |
---|
[41] | 315 | for n=0,NN-1 do begin |
---|
| 316 | x=lon(n) |
---|
| 317 | y=lat(n) |
---|
| 318 | c=bias_era(n) |
---|
[97] | 319 | ; |
---|
[41] | 320 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 321 | endfor |
---|
[97] | 322 | ; |
---|
[41] | 323 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int,/noer, marge=marge, $ |
---|
| 324 | title='3) SWR Mean bias - OAFlux', subtitle='', small=[1,4,3],/rempl,format='(i3)',/nocolorb |
---|
[97] | 325 | ; |
---|
[41] | 326 | NN=n_elements(lat) |
---|
[97] | 327 | ; |
---|
[41] | 328 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 329 | ; |
---|
[41] | 330 | for n=0,NN-1 do begin |
---|
| 331 | x=lon(n) |
---|
| 332 | y=lat(n) |
---|
| 333 | c=bias_oaflx(n) |
---|
[97] | 334 | cmi=bias_mi |
---|
| 335 | cma=bias_ma |
---|
[94] | 336 | dc=cma-cmi |
---|
[41] | 337 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 338 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 339 | endfor |
---|
[97] | 340 | ; |
---|
[41] | 341 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 342 | ; |
---|
[41] | 343 | for n=0,NN-1 do begin |
---|
| 344 | x=lon(n) |
---|
| 345 | y=lat(n) |
---|
| 346 | c=bias_oaflx(n) |
---|
[97] | 347 | ; |
---|
[41] | 348 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 349 | endfor |
---|
[97] | 350 | ; |
---|
[41] | 351 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,bias_mi, bias_ma, int=bias_int,/noer, marge=marge1, $ |
---|
[97] | 352 | title='4) SWR Mean bias - SWR_olr', subtitle='', small=[1,4,4],/rempl ;,format='(i4)' |
---|
| 353 | ; |
---|
[41] | 354 | NN=n_elements(lat) |
---|
[97] | 355 | ; |
---|
[41] | 356 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 357 | ; |
---|
[41] | 358 | for n=0,NN-1 do begin |
---|
| 359 | x=lon(n) |
---|
| 360 | y=lat(n) |
---|
| 361 | c=bias_olr(n) |
---|
[97] | 362 | cmi=bias_mi |
---|
| 363 | cma=bias_ma |
---|
[94] | 364 | dc=cma-cmi |
---|
[41] | 365 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 366 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 367 | endfor |
---|
[97] | 368 | ; |
---|
[41] | 369 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 370 | ; |
---|
[41] | 371 | for n=0,NN-1 do begin |
---|
| 372 | x=lon(n) |
---|
| 373 | y=lat(n) |
---|
| 374 | c=bias_olr(n) |
---|
[97] | 375 | ; |
---|
[41] | 376 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 377 | endfor |
---|
[97] | 378 | ; |
---|
[41] | 379 | erase |
---|
[97] | 380 | ; |
---|
[41] | 381 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ |
---|
| 382 | title='1) RMSD - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge |
---|
[97] | 383 | ; |
---|
[41] | 384 | NN=n_elements(lat) |
---|
[97] | 385 | ; |
---|
[41] | 386 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 387 | ; |
---|
[41] | 388 | for n=0,NN-1 do begin |
---|
| 389 | x=lon(n) |
---|
| 390 | y=lat(n) |
---|
| 391 | c=rmsd_trop(n) |
---|
[97] | 392 | cmi=rmsd_mi |
---|
| 393 | cma=rmsd_ma |
---|
[94] | 394 | dc=cma-cmi |
---|
[41] | 395 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 396 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 397 | endfor |
---|
[97] | 398 | ; |
---|
[41] | 399 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 400 | ; |
---|
[41] | 401 | for n=0,NN-1 do begin |
---|
| 402 | x=lon(n) |
---|
| 403 | y=lat(n) |
---|
| 404 | c=rmsd_trop(n) |
---|
[97] | 405 | ; |
---|
[41] | 406 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 407 | endfor |
---|
| 408 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ |
---|
| 409 | title='2) RMSD - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge |
---|
[97] | 410 | ; |
---|
[41] | 411 | NN=n_elements(lat) |
---|
[97] | 412 | ; |
---|
[41] | 413 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 414 | ; |
---|
[41] | 415 | for n=0,NN-1 do begin |
---|
| 416 | x=lon(n) |
---|
| 417 | y=lat(n) |
---|
| 418 | c=rmsd_era(n) |
---|
[97] | 419 | cmi=rmsd_mi |
---|
| 420 | cma=rmsd_ma |
---|
[94] | 421 | dc=cma-cmi |
---|
[41] | 422 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 423 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 424 | endfor |
---|
[97] | 425 | ; |
---|
[41] | 426 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 427 | ; |
---|
[41] | 428 | for n=0,NN-1 do begin |
---|
| 429 | x=lon(n) |
---|
| 430 | y=lat(n) |
---|
| 431 | c=rmsd_era(n) |
---|
[97] | 432 | ; |
---|
[41] | 433 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 434 | endfor |
---|
| 435 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ |
---|
| 436 | title='3) RMSD - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge |
---|
[97] | 437 | ; |
---|
[41] | 438 | NN=n_elements(lat) |
---|
[97] | 439 | ; |
---|
[41] | 440 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 441 | ; |
---|
[41] | 442 | for n=0,NN-1 do begin |
---|
| 443 | x=lon(n) |
---|
| 444 | y=lat(n) |
---|
| 445 | c=rmsd_oaflx(n) |
---|
[97] | 446 | cmi=rmsd_mi |
---|
| 447 | cma=rmsd_ma |
---|
[94] | 448 | dc=cma-cmi |
---|
[41] | 449 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 450 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 451 | endfor |
---|
[97] | 452 | ; |
---|
[41] | 453 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 454 | ; |
---|
[41] | 455 | for n=0,NN-1 do begin |
---|
| 456 | x=lon(n) |
---|
| 457 | y=lat(n) |
---|
| 458 | c=rmsd_oaflx(n) |
---|
[97] | 459 | ; |
---|
[41] | 460 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 461 | endfor |
---|
[97] | 462 | ; |
---|
[41] | 463 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ |
---|
| 464 | title='4) RMSD - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1 |
---|
[97] | 465 | ; |
---|
[41] | 466 | NN=n_elements(lat) |
---|
[97] | 467 | ; |
---|
[41] | 468 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 469 | ; |
---|
[41] | 470 | for n=0,NN-1 do begin |
---|
| 471 | x=lon(n) |
---|
| 472 | y=lat(n) |
---|
| 473 | c=rmsd_olr(n) |
---|
[97] | 474 | cmi=rmsd_mi |
---|
| 475 | cma=rmsd_ma |
---|
[94] | 476 | dc=cma-cmi |
---|
[41] | 477 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 478 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 479 | endfor |
---|
[97] | 480 | ; |
---|
[41] | 481 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 482 | ; |
---|
[41] | 483 | for n=0,NN-1 do begin |
---|
| 484 | x=lon(n) |
---|
| 485 | y=lat(n) |
---|
| 486 | c=rmsd_olr(n) |
---|
[97] | 487 | ; |
---|
[41] | 488 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 489 | endfor |
---|
| 490 | erase |
---|
| 491 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ |
---|
| 492 | title='1) STD ratio - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge |
---|
[97] | 493 | ; |
---|
[41] | 494 | NN=n_elements(lat) |
---|
[97] | 495 | ; |
---|
[41] | 496 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 497 | ; |
---|
[41] | 498 | for n=0,NN-1 do begin |
---|
| 499 | x=lon(n) |
---|
| 500 | y=lat(n) |
---|
| 501 | c=std_trop(n) |
---|
[97] | 502 | cmi=std_mi |
---|
| 503 | cma=std_ma |
---|
[94] | 504 | dc=cma-cmi |
---|
[41] | 505 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 506 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 507 | endfor |
---|
[97] | 508 | ; |
---|
[41] | 509 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 510 | ; |
---|
[41] | 511 | for n=0,NN-1 do begin |
---|
| 512 | x=lon(n) |
---|
| 513 | y=lat(n) |
---|
| 514 | c=std_trop(n) |
---|
[97] | 515 | ; |
---|
[41] | 516 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 517 | endfor |
---|
| 518 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ |
---|
| 519 | title='2) STD ratop - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge |
---|
[97] | 520 | ; |
---|
[41] | 521 | NN=n_elements(lat) |
---|
[97] | 522 | ; |
---|
[41] | 523 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 524 | ; |
---|
[41] | 525 | for n=0,NN-1 do begin |
---|
| 526 | x=lon(n) |
---|
| 527 | y=lat(n) |
---|
| 528 | c=std_era(n) |
---|
[97] | 529 | cmi=std_mi |
---|
| 530 | cma=std_ma |
---|
[94] | 531 | dc=cma-cmi |
---|
[41] | 532 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 533 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 534 | endfor |
---|
[97] | 535 | ; |
---|
[41] | 536 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 537 | ; |
---|
[41] | 538 | for n=0,NN-1 do begin |
---|
| 539 | x=lon(n) |
---|
| 540 | y=lat(n) |
---|
| 541 | c=std_era(n) |
---|
[97] | 542 | ; |
---|
[41] | 543 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 544 | endfor |
---|
[97] | 545 | ; |
---|
[41] | 546 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ |
---|
| 547 | title='3) STD rato - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge |
---|
[97] | 548 | ; |
---|
[41] | 549 | NN=n_elements(lat) |
---|
[97] | 550 | ; |
---|
[41] | 551 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 552 | ; |
---|
[41] | 553 | for n=0,NN-1 do begin |
---|
| 554 | x=lon(n) |
---|
| 555 | y=lat(n) |
---|
| 556 | c=std_oaflx(n) |
---|
[97] | 557 | cmi=std_mi |
---|
| 558 | cma=std_ma |
---|
[94] | 559 | dc=cma-cmi |
---|
[41] | 560 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 561 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 562 | endfor |
---|
[97] | 563 | ; |
---|
[41] | 564 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 565 | ; |
---|
[41] | 566 | for n=0,NN-1 do begin |
---|
| 567 | x=lon(n) |
---|
| 568 | y=lat(n) |
---|
| 569 | c=std_oaflx(n) |
---|
[97] | 570 | ; |
---|
[41] | 571 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 572 | endfor |
---|
[97] | 573 | ; |
---|
[41] | 574 | plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ |
---|
| 575 | title='4) STD rato - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1 |
---|
[97] | 576 | ; |
---|
[41] | 577 | NN=n_elements(lat) |
---|
[97] | 578 | ; |
---|
[41] | 579 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill |
---|
[97] | 580 | ; |
---|
[41] | 581 | for n=0,NN-1 do begin |
---|
| 582 | x=lon(n) |
---|
| 583 | y=lat(n) |
---|
| 584 | c=std_olr(n) |
---|
[97] | 585 | cmi=std_mi |
---|
| 586 | cma=std_ma |
---|
[94] | 587 | dc=cma-cmi |
---|
[41] | 588 | col=((10+244*(c-cmi)/dc) > 10) < 254 |
---|
| 589 | plots, x,y,psym=8,symsize=1.5,color=col |
---|
| 590 | endfor |
---|
[97] | 591 | ; |
---|
[41] | 592 | usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] |
---|
[97] | 593 | ; |
---|
[41] | 594 | for n=0,NN-1 do begin |
---|
| 595 | x=lon(n) |
---|
| 596 | y=lat(n) |
---|
| 597 | c=std_olr(n) |
---|
[97] | 598 | ; |
---|
[41] | 599 | plots, x,y,psym=8,symsize=1.5,color=0 |
---|
| 600 | endfor |
---|
[97] | 601 | ; |
---|
[41] | 602 | closeps |
---|
[97] | 603 | ; |
---|
[41] | 604 | end |
---|