Ignore:
Timestamp:
08/25/11 13:02:20 (13 years ago)
Author:
pinsard
Message:

get rid of multistatement lines

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/paper01/fig6/icoads_ws_stats_paper.pro

    r85 r94  
    8686; partie a changer 
    8787marge=[-2,-2, -4,2] 
    88 st=19890101 & en=20091231 
     88st=19890101  
     89en=20091231 
    8990domdef, box 
    9091 
     
    124125help, noc 
    125126 
    126 si=size(trop) & nt=si(3) 
     127si=size(trop)  
     128nt=si(3) 
    127129 
    128130lat_rama=[67, 55, 55, 80.5, 80.5, 80.5, 80.5, 90, 90, 90, 90, 90, 90, 90, 95] 
     
    142144lon=[lon_rama, lon_pirata, lon_tao] 
    143145 
    144 lat_tao=lat & lon_tao=lon 
     146lat_tao=lat  
     147lon_tao=lon 
    145148nn=n_elements(lat) 
    146149 
     
    149152  x=lat_tao(n) 
    150153  y=lon_tao(n) 
    151   dx=abs(reform(glamt-x)) & dy=abs(reform(gphit)-y) 
     154  dx=abs(reform(glamt-x))  
     155  dy=abs(reform(gphit)-y) 
    152156  ind=where((dx le 0.5) and (dy le 0.5)) 
    153157  if (total(ind) ge 0.) then begin 
     
    166170 
    167171param=icoads*mask 
    168 trop=trop*mask & oaf=oaf*mask & erai=erai*mask 
    169 ncep2=ncep2*mask & ncep1=ncep1*mask 
     172trop=trop*mask  
     173oaf=oaf*mask  
     174erai=erai*mask 
     175ncep2=ncep2*mask  
     176ncep1=ncep1*mask 
    170177noc=noc*mask 
    171178 
    172179ind=where(finite(param) and finite(noc) and finite(erai) and finite(trop) and finite(oaf) and finite(ncep2) and finite(ncep1)) 
    173 x=param(ind) & y=trop(ind) 
    174 param=param(ind) & trop=trop(ind) & erai=erai(ind) & oaf=oaf(ind) & ncep2=ncep2(ind) & ncep1=ncep1(ind) 
     180x=param(ind)  
     181y=trop(ind) 
     182param=param(ind)  
     183trop=trop(ind)  
     184erai=erai(ind)  
     185oaf=oaf(ind)  
     186ncep2=ncep2(ind)  
     187ncep1=ncep1(ind) 
    175188noc=noc(ind) 
    176189 
     
    219232abs_noc=string(abs_diff, format='(f5.2)') 
    220233; 
    221 mio=0 & mao=15 & inx=1 & nx=(mao-mio)/inx+1l & xx=mio+indgen(nx)*inx 
    222 mis=0 & mas=15 & iny=1 & ny=(mas-mis)/iny+1l & yy=mis+indgen(ny)*iny 
    223 xp=4 & yp=13 
    224 x1=0 & x2=15 & y1=x1 & y2=x2 
     234mio=0  
     235mao=15  
     236inx=1  
     237nx=(mao-mio)/inx+1l  
     238xx=mio+indgen(nx)*inx 
     239mis=0  
     240mas=15  
     241iny=1  
     242ny=(mas-mis)/iny+1l  
     243yy=mis+indgen(ny)*iny 
     244xp=4  
     245yp=13 
     246x1=0  
     247x2=15  
     248y1=x1  
     249y2=x2 
    225250; 
    226251pdf, param,trop,mio,mao,inx,mis,mas,iny, $ 
Note: See TracChangeset for help on using the changeset viewer.