Ignore:
Timestamp:
08/30/11 14:38:29 (13 years ago)
Author:
pinsard
Message:

suppress blank lines trailing blank

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/paper01/fig12/swr_statistics_map_2000_2009_v50.pro

    r94 r97  
    7979key_portrait = 1 
    8080coefpalit=.9 
    81  
     81; 
    8282openps, FILENAME = project_od_env+'swr_statistics_map_2000_2009_v50.ps' 
    8383; partie a changer 
    84 bias_mi=-20  
    85 bias_ma=20  
     84bias_mi=-20 
     85bias_ma=20 
    8686bias_int=2 
    87 std_mi=0.7    
    88 std_ma=1.31   
     87std_mi=0.7 
     88std_ma=1.31 
    8989std_int=0.05 
    90 rmsd_mi=10   
    91 rmsd_ma=30  
     90rmsd_mi=10 
     91rmsd_ma=30 
    9292rmsd_int=1.25 
    93 cor_mi=0.5    
    94 cor_ma=1.   
     93cor_mi=0.5 
     94cor_ma=1. 
    9595cor_int=0.02 
    9696fi_swr_erai=project_id_env+'swr_2000_2009_erai_v50.txt' 
    … …  
    9898fi_swr_oaflx=project_id_env+'swr_2000_2009_oaflx_v50.txt' 
    9999fi_swr_olr=project_id_env+'swr_2000_2009_olr_v50.txt' 
    100  
     100; 
    101101res=read_ascii(fi_swr_erai,data_start=1) 
    102102ff=res.field1 
    … …  
    107107std_era=reform(ff(4,*)) 
    108108rmsd_era=reform(ff(5,*)) 
    109  
    110 ind=where(bias_era ge bias_ma)  
     109; 
     110ind=where(bias_era ge bias_ma) 
    111111bias_era(ind)=bias_ma-0.5 
    112 ind=where(rmsd_era ge rmsd_ma)  
     112ind=where(rmsd_era ge rmsd_ma) 
    113113rmsd_era(ind)=rmsd_ma-0.5 
    114  
     114; 
    115115res=read_ascii(fi_swr_trop,data_start=1) 
    116116ff=res.field1 
    … …  
    121121std_trop=reform(ff(4,*)) 
    122122rmsd_trop=reform(ff(5,*)) 
    123 ind=where(rmsd_trop ge rmsd_ma)  
     123ind=where(rmsd_trop ge rmsd_ma) 
    124124rmsd_trop(ind)=rmsd_ma-0.5 
    125  
     125; 
    126126res=read_ascii(fi_swr_oaflx,data_start=1) 
    127127ff=res.field1 
    … …  
    132132std_oaflx=reform(ff(4,*)) 
    133133rmsd_oaflx=reform(ff(5,*)) 
    134 ind=where(rmsd_oaflx ge rmsd_ma)  
     134ind=where(rmsd_oaflx ge rmsd_ma) 
    135135rmsd_oaflx(ind)=rmsd_ma-0.5 
    136  
     136; 
    137137res=read_ascii(fi_swr_olr,data_start=1) 
    138138ff=res.field1 
    … …  
    143143std_olr=reform(ff(4,*)) 
    144144rmsd_olr=reform(ff(5,*)) 
    145 ind=where(rmsd_olr ge rmsd_ma)  
     145ind=where(rmsd_olr ge rmsd_ma) 
    146146rmsd_olr(ind)=rmsd_ma-0.5 
    147  
    148 ;ind=where(std_olr ge std_ma)  
     147; 
     148;ind=where(std_olr ge std_ma) 
    149149;std_olr(ind)=std_ma-0.01 
    150150file=project_id_env+'longwave_IO_mask.nc' 
    … …  
    157157plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,box=box,lct=64,cor_mi, cor_ma,int=cor_int, $ 
    158158   title='1) SWR Correlation - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge 
    159  
    160 NN=n_elements(lat) 
    161  
    162 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    163  
     159; 
     160NN=n_elements(lat) 
     161; 
     162usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     163; 
    164164for n=0,NN-1 do begin 
    165165  x=lon(n) 
    166166  y=lat(n) 
    167167  c=cor_trop(n) 
    168   cmi=cor_mi  
    169   cma=cor_ma  
    170   dc=cma-cmi 
    171   col=((10+244*(c-cmi)/dc) > 10) < 254 
    172   plots, x,y,psym=8,symsize=1.5,color=col 
    173 endfor 
    174  
    175 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    176  
     168  cmi=cor_mi 
     169  cma=cor_ma 
     170  dc=cma-cmi 
     171  col=((10+244*(c-cmi)/dc) > 10) < 254 
     172  plots, x,y,psym=8,symsize=1.5,color=col 
     173endfor 
     174; 
     175usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     176; 
    177177for n=0,NN-1 do begin 
    178178  x=lon(n) 
    179179  y=lat(n) 
    180180  c=cor_trop(n) 
    181  
     181; 
    182182  plots, x,y,psym=8,symsize=1.5,color=0 
    183183endfor 
    184184plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $ 
    185185   title='2) SWR Correlation -  ERAI', subtitle='', small=[1,4,2],/noer,/rempl, marge=marge 
    186  
    187 NN=n_elements(lat) 
    188  
    189 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    190  
     186; 
     187NN=n_elements(lat) 
     188; 
     189usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     190; 
    191191for n=0,NN-1 do begin 
    192192  x=lon(n) 
    193193  y=lat(n) 
    194194  c=cor_era(n) 
    195   cmi=cor_mi  
    196   cma=cor_ma  
    197   dc=cma-cmi 
    198   col=((10+244*(c-cmi)/dc) > 10) < 254 
    199   plots, x,y,psym=8,symsize=1.5,color=col 
    200 endfor 
    201  
    202 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    203  
     195  cmi=cor_mi 
     196  cma=cor_ma 
     197  dc=cma-cmi 
     198  col=((10+244*(c-cmi)/dc) > 10) < 254 
     199  plots, x,y,psym=8,symsize=1.5,color=col 
     200endfor 
     201; 
     202usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     203; 
    204204for n=0,NN-1 do begin 
    205205  x=lon(n) 
    206206  y=lat(n) 
    207207  c=cor_era(n) 
    208  
     208; 
    209209  plots, x,y,psym=8,symsize=1.5,color=0 
    210210endfor 
    211211plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int,/nocolorb, $ 
    212212   title='3) SWR Correlation - OAFlux', subtitle='', small=[1,4,3],/noer,/rempl, marge=marge 
    213  
    214 NN=n_elements(lat) 
    215  
    216 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    217  
     213; 
     214NN=n_elements(lat) 
     215; 
     216usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     217; 
    218218for n=0,NN-1 do begin 
    219219  x=lon(n) 
    220220  y=lat(n) 
    221221  c=cor_oaflx(n) 
    222   cmi=cor_mi  
    223   cma=cor_ma  
    224   dc=cma-cmi 
    225   col=((10+244*(c-cmi)/dc) > 10) < 254 
    226   plots, x,y,psym=8,symsize=1.5,color=col 
    227 endfor 
    228  
    229 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    230  
     222  cmi=cor_mi 
     223  cma=cor_ma 
     224  dc=cma-cmi 
     225  col=((10+244*(c-cmi)/dc) > 10) < 254 
     226  plots, x,y,psym=8,symsize=1.5,color=col 
     227endfor 
     228; 
     229usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     230; 
    231231for n=0,NN-1 do begin 
    232232  x=lon(n) 
    233233  y=lat(n) 
    234234  c=cor_oaflx(n) 
    235  
     235; 
    236236  plots, x,y,psym=8,symsize=1.5,color=0 
    237237endfor 
    238238plt, 12+msk,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,cor_mi, cor_ma,int=cor_int, $ 
    239239   title='4) SWR Correlation - SWR_olr', subtitle='', small=[1,4,4],/noer,/rempl, marge=marge1 
    240  
    241 NN=n_elements(lat) 
    242  
    243 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    244  
     240; 
     241NN=n_elements(lat) 
     242; 
     243usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     244; 
    245245for n=0,NN-1 do begin 
    246246  x=lon(n) 
    247247  y=lat(n) 
    248248  c=cor_olr(n) 
    249   cmi=cor_mi  
    250   cma=cor_ma  
    251   dc=cma-cmi 
    252   col=((10+244*(c-cmi)/dc) > 10) < 254 
    253   plots, x,y,psym=8,symsize=1.5,color=col 
    254 endfor 
    255  
    256 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    257  
     249  cmi=cor_mi 
     250  cma=cor_ma 
     251  dc=cma-cmi 
     252  col=((10+244*(c-cmi)/dc) > 10) < 254 
     253  plots, x,y,psym=8,symsize=1.5,color=col 
     254endfor 
     255; 
     256usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     257; 
    258258for n=0,NN-1 do begin 
    259259  x=lon(n) 
    260260  y=lat(n) 
    261261  c=cor_olr(n) 
    262  
    263   plots, x,y,psym=8,symsize=1.5,color=0 
    264 endfor 
    265  
     262; 
     263  plots, x,y,psym=8,symsize=1.5,color=0 
     264endfor 
     265; 
    266266erase 
    267  
     267; 
    268268plt, -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, $ 
    269269   title='1) SWR Mean bias - TropFlux', subtitle='', small=[1,4,1],/rempl,format='(i3)' 
    270  
    271 NN=n_elements(lat) 
    272  
    273 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    274  
     270; 
     271NN=n_elements(lat) 
     272; 
     273usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     274; 
    275275for n=0,NN-1 do begin 
    276276  x=lon(n) 
    277277  y=lat(n) 
    278278  c=bias_trop(n) 
    279   cmi=bias_mi  
    280   cma=bias_ma  
    281   dc=cma-cmi 
    282   col=((10+244*(c-cmi)/dc) > 10) < 254 
    283   plots, x,y,psym=8,symsize=1.5,color=col 
    284 endfor 
    285  
    286 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    287  
     279  cmi=bias_mi 
     280  cma=bias_ma 
     281  dc=cma-cmi 
     282  col=((10+244*(c-cmi)/dc) > 10) < 254 
     283  plots, x,y,psym=8,symsize=1.5,color=col 
     284endfor 
     285; 
     286usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     287; 
    288288for n=0,NN-1 do begin 
    289289  x=lon(n) 
    290290  y=lat(n) 
    291291  c=bias_trop(n) 
    292  
     292; 
    293293  plots, x,y,psym=8,symsize=1.5,color=0 
    294294endfor 
    295295plt,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, $ 
    296296   title='2) SWR Mean bias - ERAI', subtitle='', small=[1,4,2],/rempl,format='(i3)',/nocolorb 
    297  
    298 NN=n_elements(lat) 
    299  
    300 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    301  
     297; 
     298NN=n_elements(lat) 
     299; 
     300usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     301; 
    302302for n=0,NN-1 do begin 
    303303  x=lon(n) 
    304304  y=lat(n) 
    305305  c=bias_era(n) 
    306   cmi=bias_mi  
    307   cma=bias_ma  
    308   dc=cma-cmi 
    309   col=((10+244*(c-cmi)/dc) > 10) < 254 
    310   plots, x,y,psym=8,symsize=1.5,color=col 
    311 endfor 
    312  
    313 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    314  
     306  cmi=bias_mi 
     307  cma=bias_ma 
     308  dc=cma-cmi 
     309  col=((10+244*(c-cmi)/dc) > 10) < 254 
     310  plots, x,y,psym=8,symsize=1.5,color=col 
     311endfor 
     312; 
     313usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     314; 
    315315for n=0,NN-1 do begin 
    316316  x=lon(n) 
    317317  y=lat(n) 
    318318  c=bias_era(n) 
    319  
    320   plots, x,y,psym=8,symsize=1.5,color=0 
    321 endfor 
    322  
    323  
     319; 
     320  plots, x,y,psym=8,symsize=1.5,color=0 
     321endfor 
     322; 
     323; 
    324324plt,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, $ 
    325325   title='3) SWR Mean bias - OAFlux', subtitle='', small=[1,4,3],/rempl,format='(i3)',/nocolorb 
    326  
    327 NN=n_elements(lat) 
    328  
    329 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    330  
     326; 
     327NN=n_elements(lat) 
     328; 
     329usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     330; 
    331331for n=0,NN-1 do begin 
    332332  x=lon(n) 
    333333  y=lat(n) 
    334334  c=bias_oaflx(n) 
    335   cmi=bias_mi  
    336   cma=bias_ma  
    337   dc=cma-cmi 
    338   col=((10+244*(c-cmi)/dc) > 10) < 254 
    339   plots, x,y,psym=8,symsize=1.5,color=col 
    340 endfor 
    341  
    342 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    343  
     335  cmi=bias_mi 
     336  cma=bias_ma 
     337  dc=cma-cmi 
     338  col=((10+244*(c-cmi)/dc) > 10) < 254 
     339  plots, x,y,psym=8,symsize=1.5,color=col 
     340endfor 
     341; 
     342usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     343; 
    344344for n=0,NN-1 do begin 
    345345  x=lon(n) 
    346346  y=lat(n) 
    347347  c=bias_oaflx(n) 
    348  
    349   plots, x,y,psym=8,symsize=1.5,color=0 
    350 endfor 
    351  
    352  
     348; 
     349  plots, x,y,psym=8,symsize=1.5,color=0 
     350endfor 
     351; 
     352; 
    353353plt,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, $ 
    354    title='4) SWR Mean bias - SWR_olr', subtitle='', small=[1,4,4],/rempl ;;,format='(i4)' 
    355  
    356 NN=n_elements(lat) 
    357  
    358 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    359  
     354   title='4) SWR Mean bias - SWR_olr', subtitle='', small=[1,4,4],/rempl ;,format='(i4)' 
     355; 
     356NN=n_elements(lat) 
     357; 
     358usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     359; 
    360360for n=0,NN-1 do begin 
    361361  x=lon(n) 
    362362  y=lat(n) 
    363363  c=bias_olr(n) 
    364   cmi=bias_mi  
    365   cma=bias_ma  
    366   dc=cma-cmi 
    367   col=((10+244*(c-cmi)/dc) > 10) < 254 
    368   plots, x,y,psym=8,symsize=1.5,color=col 
    369 endfor 
    370  
    371 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    372  
     364  cmi=bias_mi 
     365  cma=bias_ma 
     366  dc=cma-cmi 
     367  col=((10+244*(c-cmi)/dc) > 10) < 254 
     368  plots, x,y,psym=8,symsize=1.5,color=col 
     369endfor 
     370; 
     371usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     372; 
    373373for n=0,NN-1 do begin 
    374374  x=lon(n) 
    375375  y=lat(n) 
    376376  c=bias_olr(n) 
    377  
    378   plots, x,y,psym=8,symsize=1.5,color=0 
    379 endfor 
    380  
     377; 
     378  plots, x,y,psym=8,symsize=1.5,color=0 
     379endfor 
     380; 
    381381erase 
    382  
    383  
     382; 
     383; 
    384384plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 
    385385   title='1) RMSD - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge 
    386  
    387 NN=n_elements(lat) 
    388  
    389 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    390  
     386; 
     387NN=n_elements(lat) 
     388; 
     389usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     390; 
    391391for n=0,NN-1 do begin 
    392392  x=lon(n) 
    393393  y=lat(n) 
    394394  c=rmsd_trop(n) 
    395   cmi=rmsd_mi  
    396   cma=rmsd_ma  
    397   dc=cma-cmi 
    398   col=((10+244*(c-cmi)/dc) > 10) < 254 
    399   plots, x,y,psym=8,symsize=1.5,color=col 
    400 endfor 
    401  
    402 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    403  
     395  cmi=rmsd_mi 
     396  cma=rmsd_ma 
     397  dc=cma-cmi 
     398  col=((10+244*(c-cmi)/dc) > 10) < 254 
     399  plots, x,y,psym=8,symsize=1.5,color=col 
     400endfor 
     401; 
     402usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     403; 
    404404for n=0,NN-1 do begin 
    405405  x=lon(n) 
    406406  y=lat(n) 
    407407  c=rmsd_trop(n) 
    408  
     408; 
    409409  plots, x,y,psym=8,symsize=1.5,color=0 
    410410endfor 
    411411plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 
    412412   title='2) RMSD - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge 
    413  
    414 NN=n_elements(lat) 
    415  
    416 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    417  
     413; 
     414NN=n_elements(lat) 
     415; 
     416usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     417; 
    418418for n=0,NN-1 do begin 
    419419  x=lon(n) 
    420420  y=lat(n) 
    421421  c=rmsd_era(n) 
    422   cmi=rmsd_mi  
    423   cma=rmsd_ma  
    424   dc=cma-cmi 
    425   col=((10+244*(c-cmi)/dc) > 10) < 254 
    426   plots, x,y,psym=8,symsize=1.5,color=col 
    427 endfor 
    428  
    429 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    430  
     422  cmi=rmsd_mi 
     423  cma=rmsd_ma 
     424  dc=cma-cmi 
     425  col=((10+244*(c-cmi)/dc) > 10) < 254 
     426  plots, x,y,psym=8,symsize=1.5,color=col 
     427endfor 
     428; 
     429usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     430; 
    431431for n=0,NN-1 do begin 
    432432  x=lon(n) 
    433433  y=lat(n) 
    434434  c=rmsd_era(n) 
    435  
     435; 
    436436  plots, x,y,psym=8,symsize=1.5,color=0 
    437437endfor 
    438438plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 
    439439   title='3) RMSD - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge 
    440  
    441 NN=n_elements(lat) 
    442  
    443 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    444  
     440; 
     441NN=n_elements(lat) 
     442; 
     443usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     444; 
    445445for n=0,NN-1 do begin 
    446446  x=lon(n) 
    447447  y=lat(n) 
    448448  c=rmsd_oaflx(n) 
    449   cmi=rmsd_mi  
    450   cma=rmsd_ma  
    451   dc=cma-cmi 
    452   col=((10+244*(c-cmi)/dc) > 10) < 254 
    453   plots, x,y,psym=8,symsize=1.5,color=col 
    454 endfor 
    455  
    456 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    457  
     449  cmi=rmsd_mi 
     450  cma=rmsd_ma 
     451  dc=cma-cmi 
     452  col=((10+244*(c-cmi)/dc) > 10) < 254 
     453  plots, x,y,psym=8,symsize=1.5,color=col 
     454endfor 
     455; 
     456usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     457; 
    458458for n=0,NN-1 do begin 
    459459  x=lon(n) 
    460460  y=lat(n) 
    461461  c=rmsd_oaflx(n) 
    462  
    463   plots, x,y,psym=8,symsize=1.5,color=0 
    464 endfor 
    465  
     462; 
     463  plots, x,y,psym=8,symsize=1.5,color=0 
     464endfor 
     465; 
    466466plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,rmsd_mi, rmsd_ma, int=rmsd_int,/noer, $ 
    467467   title='4) RMSD - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1 
    468  
    469 NN=n_elements(lat) 
    470  
    471 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    472  
     468; 
     469NN=n_elements(lat) 
     470; 
     471usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     472; 
    473473for n=0,NN-1 do begin 
    474474  x=lon(n) 
    475475  y=lat(n) 
    476476  c=rmsd_olr(n) 
    477   cmi=rmsd_mi  
    478   cma=rmsd_ma  
    479   dc=cma-cmi 
    480   col=((10+244*(c-cmi)/dc) > 10) < 254 
    481   plots, x,y,psym=8,symsize=1.5,color=col 
    482 endfor 
    483  
    484 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    485  
     477  cmi=rmsd_mi 
     478  cma=rmsd_ma 
     479  dc=cma-cmi 
     480  col=((10+244*(c-cmi)/dc) > 10) < 254 
     481  plots, x,y,psym=8,symsize=1.5,color=col 
     482endfor 
     483; 
     484usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     485; 
    486486for n=0,NN-1 do begin 
    487487  x=lon(n) 
    488488  y=lat(n) 
    489489  c=rmsd_olr(n) 
    490  
     490; 
    491491  plots, x,y,psym=8,symsize=1.5,color=0 
    492492endfor 
    … …  
    494494plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 
    495495   title='1) STD ratio - TropFlux', subtitle='', small=[1,4,1],/rempl,/nocolorb, marge=marge 
    496  
    497 NN=n_elements(lat) 
    498  
    499 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    500  
     496; 
     497NN=n_elements(lat) 
     498; 
     499usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     500; 
    501501for n=0,NN-1 do begin 
    502502  x=lon(n) 
    503503  y=lat(n) 
    504504  c=std_trop(n) 
    505   cmi=std_mi  
    506   cma=std_ma  
    507   dc=cma-cmi 
    508   col=((10+244*(c-cmi)/dc) > 10) < 254 
    509   plots, x,y,psym=8,symsize=1.5,color=col 
    510 endfor 
    511  
    512 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    513  
     505  cmi=std_mi 
     506  cma=std_ma 
     507  dc=cma-cmi 
     508  col=((10+244*(c-cmi)/dc) > 10) < 254 
     509  plots, x,y,psym=8,symsize=1.5,color=col 
     510endfor 
     511; 
     512usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     513; 
    514514for n=0,NN-1 do begin 
    515515  x=lon(n) 
    516516  y=lat(n) 
    517517  c=std_trop(n) 
    518  
     518; 
    519519  plots, x,y,psym=8,symsize=1.5,color=0 
    520520endfor 
    521521plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 
    522522   title='2) STD ratop - ERAI', subtitle='', small=[1,4,2],/rempl,/nocolorb, marge=marge 
    523  
    524 NN=n_elements(lat) 
    525  
    526 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    527  
     523; 
     524NN=n_elements(lat) 
     525; 
     526usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     527; 
    528528for n=0,NN-1 do begin 
    529529  x=lon(n) 
    530530  y=lat(n) 
    531531  c=std_era(n) 
    532   cmi=std_mi  
    533   cma=std_ma  
    534   dc=cma-cmi 
    535   col=((10+244*(c-cmi)/dc) > 10) < 254 
    536   plots, x,y,psym=8,symsize=1.5,color=col 
    537 endfor 
    538  
    539 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    540  
     532  cmi=std_mi 
     533  cma=std_ma 
     534  dc=cma-cmi 
     535  col=((10+244*(c-cmi)/dc) > 10) < 254 
     536  plots, x,y,psym=8,symsize=1.5,color=col 
     537endfor 
     538; 
     539usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     540; 
    541541for n=0,NN-1 do begin 
    542542  x=lon(n) 
    543543  y=lat(n) 
    544544  c=std_era(n) 
    545  
    546   plots, x,y,psym=8,symsize=1.5,color=0 
    547 endfor 
    548  
     545; 
     546  plots, x,y,psym=8,symsize=1.5,color=0 
     547endfor 
     548; 
    549549plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 
    550550   title='3) STD rato - OAFlux', subtitle='', small=[1,4,3],/rempl,/nocolorb, marge=marge 
    551  
    552 NN=n_elements(lat) 
    553  
    554 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    555  
     551; 
     552NN=n_elements(lat) 
     553; 
     554usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     555; 
    556556for n=0,NN-1 do begin 
    557557  x=lon(n) 
    558558  y=lat(n) 
    559559  c=std_oaflx(n) 
    560   cmi=std_mi  
    561   cma=std_ma  
    562   dc=cma-cmi 
    563   col=((10+244*(c-cmi)/dc) > 10) < 254 
    564   plots, x,y,psym=8,symsize=1.5,color=col 
    565 endfor 
    566  
    567 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    568  
     560  cmi=std_mi 
     561  cma=std_ma 
     562  dc=cma-cmi 
     563  col=((10+244*(c-cmi)/dc) > 10) < 254 
     564  plots, x,y,psym=8,symsize=1.5,color=col 
     565endfor 
     566; 
     567usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     568; 
    569569for n=0,NN-1 do begin 
    570570  x=lon(n) 
    571571  y=lat(n) 
    572572  c=std_oaflx(n) 
    573  
    574   plots, x,y,psym=8,symsize=1.5,color=0 
    575 endfor 
    576  
     573; 
     574  plots, x,y,psym=8,symsize=1.5,color=0 
     575endfor 
     576; 
    577577plt,0.8+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,lct=64,std_mi, std_ma, int=std_int,/noer, $ 
    578578   title='4) STD rato - SWR_olr', subtitle='', small=[1,4,4],/rempl, marge=marge1 
    579  
    580 NN=n_elements(lat) 
    581  
    582 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    583  
     579; 
     580NN=n_elements(lat) 
     581; 
     582usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
     583; 
    584584for n=0,NN-1 do begin 
    585585  x=lon(n) 
    586586  y=lat(n) 
    587587  c=std_olr(n) 
    588   cmi=std_mi  
    589   cma=std_ma  
    590   dc=cma-cmi 
    591   col=((10+244*(c-cmi)/dc) > 10) < 254 
    592   plots, x,y,psym=8,symsize=1.5,color=col 
    593 endfor 
    594  
    595 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    596  
     588  cmi=std_mi 
     589  cma=std_ma 
     590  dc=cma-cmi 
     591  col=((10+244*(c-cmi)/dc) > 10) < 254 
     592  plots, x,y,psym=8,symsize=1.5,color=col 
     593endfor 
     594; 
     595usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
     596; 
    597597for n=0,NN-1 do begin 
    598598  x=lon(n) 
    599599  y=lat(n) 
    600600  c=std_olr(n) 
    601  
    602   plots, x,y,psym=8,symsize=1.5,color=0 
    603 endfor 
     601; 
     602  plots, x,y,psym=8,symsize=1.5,color=0 
     603endfor 
     604; 
    604605closeps 
     606; 
    605607end 
Note: See TracChangeset for help on using the changeset viewer.