Changeset 97 for trunk/src/paper01/fig1


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/fig1/fig1_no_obs_map.pro

    r94 r97  
    1313; has been produced by :ref:`oaflux_mask_30N30S.pro`. 
    1414; 
    15 ; Figure of  
     15; Figure of 
    1616; ++ 
    1717; is saved in 
     
    107107reinitplt, /z,/invert 
    108108key_portrait = 1 
    109  
     109; 
    110110openps, FILENAME = project_od_env + 'fig1_no_obs_map.ps' 
    111111close,/all 
    112  
     112; 
    113113yr=365. 
    114 lh_mi=365/yr  
    115 lh_ma=3285/yr  
     114lh_mi=365/yr 
     115lh_ma=3285/yr 
    116116lh_int=0.25 
    117 sw_mi=365/yr  
    118 sw_ma=3285/yr  
     117sw_mi=365/yr 
     118sw_ma=3285/yr 
    119119sw_int=0.25 
    120 lw_mi=365/yr  
    121 lw_ma=1460/yr  
     120lw_mi=365/yr 
     121lw_ma=1460/yr 
    122122lw_int=0.125 
    123123lct=65 
    124  
     124; 
    125125file=project_id_env + 'mask_oaflux_30N30S.nc' 
    126126initncdf, file 
    127127domdef, 30,390,-30,30 
    128128msk=ncdf_lec(file,var='msk') 
    129  
     129; 
    130130fi=project_id_env + 'lhf_no_obs.txt' 
    131131res=read_ascii(fi,data_start=1) 
     
    135135no=reform(ff(2,*)) 
    136136no=no/yr 
    137 ;ind=where(no ge lh_ma)  
     137;ind=where(no ge lh_ma) 
    138138; no(ind)=no-0.5 
    139  
     139; 
    140140marge1=[0,0,-5,0] 
    141141plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,box=box,lct=lct, lh_mi, lh_ma, int=lh_int, $ 
    142142   title='a) Turbulent Flux validation sites', subtitle='', small=[1,4,1],/rempl, marge=marge1, format='(i5)' 
    143143NN=n_elements(lat) 
    144  
     144; 
    145145usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    146  
    147 for n=0,NN-1 do begin 
    148   x=lon[n] 
    149   y=lat[n] 
    150   c=no[n] 
    151   cmi=lh_mi  
    152   cma=lh_ma  
     146; 
     147for n=0,NN-1 do begin 
     148  x=lon[n] 
     149  y=lat[n] 
     150  c=no[n] 
     151  cmi=lh_mi 
     152  cma=lh_ma 
    153153  dc=cma-cmi 
    154154  col=((10+244*(c-cmi)/dc) > 10) < 254 
    155155  plots, x,y,psym=8,symsize=1.5,color=col 
    156156endfor 
    157  
     157; 
    158158usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    159  
    160 for n=0,NN-1 do begin 
    161   x=lon[n] 
    162   y=lat[n] 
    163   c=no[n] 
    164  
     159; 
     160for n=0,NN-1 do begin 
     161  x=lon[n] 
     162  y=lat[n] 
     163  c=no[n] 
     164; 
    165165  plots, x,y,psym=8,symsize=1.5,color=0 
    166166endfor 
    167  
     167; 
    168168fi=project_id_env + 'sw_no_obs.txt' 
    169169res=read_ascii(fi,data_start=1) 
     
    173173no=reform(ff(2,*)) 
    174174no=no/yr 
    175 ;ind=where(no ge sw_ma)  
     175;ind=where(no ge sw_ma) 
    176176;no(ind)=sw_ma-0.5 
    177177marge1=[0,0,-5,0] 
     
    179179   title='b) SWR validation sites', subtitle='', small=[1,4,2],/noer,/rempl, marge=marge1, format='(i5)' 
    180180NN=n_elements(lat) 
    181  
     181; 
    182182usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    183  
    184 for n=0,NN-1 do begin 
    185   x=lon[n] 
    186   y=lat[n] 
    187   c=no[n] 
    188   cmi=sw_mi  
    189   cma=sw_ma  
     183; 
     184for n=0,NN-1 do begin 
     185  x=lon[n] 
     186  y=lat[n] 
     187  c=no[n] 
     188  cmi=sw_mi 
     189  cma=sw_ma 
    190190  dc=cma-cmi 
    191191  col=((10+244*(c-cmi)/dc) > 10) < 254 
    192192  plots, x,y,psym=8,symsize=1.5,color=col 
    193193endfor 
    194  
     194; 
    195195usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    196  
    197 for n=0,NN-1 do begin 
    198   x=lon[n] 
    199   y=lat[n] 
    200   c=no[n] 
    201  
     196; 
     197for n=0,NN-1 do begin 
     198  x=lon[n] 
     199  y=lat[n] 
     200  c=no[n] 
     201; 
    202202  plots, x,y,psym=8,symsize=1.5,color=0 
    203203endfor 
    204  
     204; 
    205205fi=project_id_env + 'lw_no_obs.txt' 
    206206res=read_ascii(fi,data_start=1) 
     
    210210no=reform(ff(2,*)) 
    211211no=no/yr 
    212 ;ind=where(no ge lw_ma)  
     212;ind=where(no ge lw_ma) 
    213213;no(ind)=lw_ma-0.25 
    214  
     214; 
    215215marge1=[0,0,-5,0] 
    216216plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,/noer,xminor=1,yminor=1,box=box,lct=lct,lw_mi, lw_ma, int=lw_int, $ 
    217217   title='c) LWR validation sites', subtitle='', small=[1,4,3],/rempl, marge=marge1, format='(i5)' 
    218218NN=n_elements(lat) 
    219  
     219; 
    220220usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 
    221  
    222 for n=0,NN-1 do begin 
    223   x=lon[n] 
    224   y=lat[n] 
    225   c=no[n] 
    226   cmi=lw_mi  
    227   cma=lw_ma  
     221; 
     222for n=0,NN-1 do begin 
     223  x=lon[n] 
     224  y=lat[n] 
     225  c=no[n] 
     226  cmi=lw_mi 
     227  cma=lw_ma 
    228228  dc=cma-cmi 
    229229  col=((10+244*(c-cmi)/dc) > 10) < 254 
    230230  plots, x,y,psym=8,symsize=1.5,color=col 
    231231endfor 
    232  
     232; 
    233233usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 
    234  
    235 for n=0,NN-1 do begin 
    236   x=lon[n] 
    237   y=lat[n] 
    238   c=no[n] 
    239  
     234; 
     235for n=0,NN-1 do begin 
     236  x=lon[n] 
     237  y=lat[n] 
     238  c=no[n] 
     239; 
    240240  plots, x,y,psym=8,symsize=1.5,color=0 
    241241endfor 
    242  
    243 ;; independent validation sites 
    244  
    245 plots, 61.5, 15.5,psym=2      ;; as location 
    246 plots, 331, 25.5,psym=2       ;; subduction (central) 
    247 plots, 326, 18,psym=2         ;; subduction (sw) 
    248 plots, 338, 18,psym=2         ;; subduction (se) 
    249 plots, 156, -1.7 , psym=2     ;; coare 
    250 plots, 274.2, -21, psym=2     ;; stratus 
    251 plots, 202.46, 22.45, psym=2  ;; WHOTS 
    252  
     242; 
     243; independent validation sites 
     244; 
     245plots, 61.5, 15.5,psym=2      ; as location 
     246plots, 331, 25.5,psym=2       ; subduction (central) 
     247plots, 326, 18,psym=2         ; subduction (sw) 
     248plots, 338, 18,psym=2         ; subduction (se) 
     249plots, 156, -1.7 , psym=2     ; coare 
     250plots, 274.2, -21, psym=2     ; stratus 
     251plots, 202.46, 22.45, psym=2  ; WHOTS 
     252; 
    253253closeps 
    254  
     254; 
    255255end 
Note: See TracChangeset for help on using the changeset viewer.