Changeset 94 for trunk/src/paper01/fig6


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

get rid of multistatement lines

Location:
trunk/src/paper01/fig6
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/paper01/fig6/icoads_q2m_stats_paper.pro

    r85 r94  
    8484; partie a changer 
    8585marge=[-2,-2, -4,2] 
    86 st=19890101 & en=20091231 
     86st=19890101  
     87en=20091231 
    8788domdef, box 
    8889 
    … …  
    118119 
    119120 
    120 si=size(trop) & nt=si(3) 
     121si=size(trop)  
     122nt=si(3) 
    121123 
    122124lat_rama=[67, 55, 55, 80.5, 80.5, 80.5, 80.5, 90, 90, 90, 90, 90, 90, 90, 95] 
    … …  
    136138lon=[lon_rama, lon_pirata, lon_tao] 
    137139 
    138 lat_tao=lat & lon_tao=lon 
     140lat_tao=lat  
     141lon_tao=lon 
    139142nn=n_elements(lat) 
    140143 
    … …  
    143146  x=lat_tao(n) 
    144147  y=lon_tao(n) 
    145   dx=abs(reform(glamt-x)) & dy=abs(reform(gphit)-y) 
     148  dx=abs(reform(glamt-x))  
     149  dy=abs(reform(gphit)-y) 
    146150  ind=where((dx le 0.5) and (dy le 0.5)) 
    147151  if (total(ind) ge 0.) then begin 
    … …  
    160164 
    161165param=icoads*mask 
    162 trop=trop*mask & oaf=oaf*mask & erai=erai*mask 
    163 ncep2=ncep2*mask & ncep1=ncep1*mask 
     166trop=trop*mask  
     167oaf=oaf*mask  
     168erai=erai*mask 
     169ncep2=ncep2*mask  
     170ncep1=ncep1*mask 
    164171 
    165172 
    166173ind=where(finite(param) and finite(erai) and finite(trop) and finite(oaf) and finite(ncep2) and finite(ncep1)) 
    167 x=param(ind) & y=trop(ind) 
    168 param=param(ind) & trop=trop(ind) & erai=erai(ind) & oaf=oaf(ind) & ncep2=ncep2(ind) & ncep1=ncep1(ind) 
     174x=param(ind)  
     175y=trop(ind) 
     176param=param(ind)  
     177trop=trop(ind)  
     178erai=erai(ind)  
     179oaf=oaf(ind)  
     180ncep2=ncep2(ind)  
     181ncep1=ncep1(ind) 
    169182 
    170183; 
    … …  
    208221 
    209222; 
    210 mio=4 & mao=24 & inx=1 & nx=(mao-mio)/inx+1l & xx=mio+indgen(nx)*inx 
    211 mis=4 & mas=24 & iny=1 & ny=(mas-mis)/iny+1l & yy=mis+indgen(ny)*iny 
     223mio=4  
     224mao=24  
     225inx=1  
     226nx=(mao-mio)/inx+1l  
     227xx=mio+indgen(nx)*inx 
     228mis=4  
     229mas=24  
     230iny=1  
     231ny=(mas-mis)/iny+1l  
     232yy=mis+indgen(ny)*iny 
    212233; 
    213234pdf, param,trop,mio,mao,inx,mis,mas,iny, $ 

  • trunk/src/paper01/fig6/icoads_sst_stats_paper.pro

    r85 r94  
    8484; partie a changer 
    8585marge=[-2,-2, -4,2] 
    86 st=19890101 & en=20091231 
     86st=19890101  
     87en=20091231 
    8788domdef, box 
    8889 
    … …  
    123124 
    124125 
    125 si=size(trop) & nt=si(3) 
     126si=size(trop)  
     127nt=si(3) 
    126128 
    127129lat_rama=[67, 55, 55, 80.5, 80.5, 80.5, 80.5, 90, 90, 90, 90, 90, 90, 90, 95] 
    … …  
    141143lon=[lon_rama, lon_pirata, lon_tao] 
    142144 
    143 lat_tao=lat & lon_tao=lon 
     145lat_tao=lat  
     146lon_tao=lon 
    144147nn=n_elements(lat) 
    145148 
    … …  
    148151  x=lat_tao(n) 
    149152  y=lon_tao(n) 
    150   dx=abs(reform(glamt-x)) & dy=abs(reform(gphit)-y) 
     153  dx=abs(reform(glamt-x))  
     154  dy=abs(reform(gphit)-y) 
    151155  ind=where((dx le 0.5) and (dy le 0.5)) 
    152156  if (total(ind) ge 0.) then begin 
    … …  
    165169 
    166170param=icoads*mask 
    167 trop=trop*mask & oaf=oaf*mask & erai=erai*mask 
    168 ncep2=ncep2*mask & ncep1=ncep1*mask 
     171trop=trop*mask  
     172oaf=oaf*mask  
     173erai=erai*mask 
     174ncep2=ncep2*mask  
     175ncep1=ncep1*mask 
    169176noc=noc*mask 
    170177 
    171178ind=where(finite(param) and finite(erai) and finite(trop) and finite(oaf) and finite(ncep2) and finite(ncep1)) 
    172 x=param(ind) & y=trop(ind) 
    173 param=param(ind) & trop=trop(ind) & erai=erai(ind) & oaf=oaf(ind) & ncep2=ncep2(ind) & ncep1=ncep1(ind) 
     179x=param(ind)  
     180y=trop(ind) 
     181param=param(ind)  
     182trop=trop(ind)  
     183erai=erai(ind)  
     184oaf=oaf(ind)  
     185ncep2=ncep2(ind)  
     186ncep1=ncep1(ind) 
    174187noc=noc(ind) 
    175188 
    … …  
    220233 
    221234; 
    222 mio=15 & mao=32 & inx=1 & nx=(mao-mio)/inx+1l & xx=mio+indgen(nx)*inx 
    223 mis=15 & mas=32 & iny=1 & ny=(mas-mis)/iny+1l & yy=mis+indgen(ny)*iny 
    224 xp=18 & yp=16.5 
    225 x1=15 & x2=32 & y1=x1 & y2=x2 
     235mio=15  
     236mao=32  
     237inx=1 
     238nx=(mao-mio)/inx+1l  
     239xx=mio+indgen(nx)*inx 
     240mis=15  
     241mas=32  
     242iny=1  
     243ny=(mas-mis)/iny+1l  
     244yy=mis+indgen(ny)*iny 
     245xp=18  
     246yp=16.5 
     247x1=15  
     248x2=32  
     249y1=x1  
     250y2=x2 
    226251; 
    227252pdf, param,trop,mio,mao,inx,mis,mas,iny, $ 

  • trunk/src/paper01/fig6/icoads_t2m_stats_paper.pro

    r85 r94  
    8585; partie a changer 
    8686marge=[-2,-2, -4,2] 
    87 st=19890101 & en=20091231 
     87st=19890101  
     88en=20091231 
    8889domdef, box 
    8990 
    … …  
    119120 
    120121 
    121 si=size(trop) & nt=si(3) 
     122si=size(trop)  
     123nt=si(3) 
    122124 
    123125lat_rama=[67, 55, 55, 80.5, 80.5, 80.5, 80.5, 90, 90, 90, 90, 90, 90, 90, 95] 
    … …  
    137139lon=[lon_rama, lon_pirata, lon_tao] 
    138140 
    139 lat_tao=lat & lon_tao=lon 
     141lat_tao=lat  
     142lon_tao=lon 
    140143nn=n_elements(lat) 
    141144 
    … …  
    144147  x=lat_tao(n) 
    145148  y=lon_tao(n) 
    146   dx=abs(reform(glamt-x)) & dy=abs(reform(gphit)-y) 
     149  dx=abs(reform(glamt-x))  
     150  dy=abs(reform(gphit)-y) 
    147151  ind=where((dx le 0.5) and (dy le 0.5)) 
    148152  if (total(ind) ge 0.) then begin 
    … …  
    161165 
    162166param=icoads*mask 
    163 trop=trop*mask & oaf=oaf*mask & erai=erai*mask 
    164 ncep2=ncep2*mask & ncep1=ncep1*mask 
     167trop=trop*mask  
     168oaf=oaf*mask  
     169erai=erai*mask 
     170ncep2=ncep2*mask  
     171ncep1=ncep1*mask 
    165172 
    166173 
    167174ind=where(finite(param) and finite(erai) and finite(trop) and finite(oaf) and finite(ncep2) and finite(ncep1)) 
    168 x=param(ind) & y=trop(ind) 
    169 param=param(ind) & trop=trop(ind) & erai=erai(ind) & oaf=oaf(ind) & ncep2=ncep2(ind) & ncep1=ncep1(ind) 
     175x=param(ind)  
     176y=trop(ind) 
     177param=param(ind)  
     178trop=trop(ind)  
     179erai=erai(ind)  
     180oaf=oaf(ind)  
     181ncep2=ncep2(ind)  
     182ncep1=ncep1(ind) 
    170183 
    171184; 
    … …  
    209222 
    210223; 
    211 mio=15 & mao=32 & inx=1 & nx=(mao-mio)/inx+1l & xx=mio+indgen(nx)*inx 
    212 mis=15 & mas=32 & iny=1 & ny=(mas-mis)/iny+1l & yy=mis+indgen(ny)*iny 
    213 xp=16 & yp=14 
    214  
    215 x1=13 & x2=32 & y1=x1 & y2=x2 
     224mio=15  
     225mao=32  
     226inx=1  
     227nx=(mao-mio)/inx+1l  
     228xx=mio+indgen(nx)*inx 
     229mis=15  
     230mas=32  
     231iny=1  
     232ny=(mas-mis)/iny+1l  
     233yy=mis+indgen(ny)*iny 
     234xp=16  
     235yp=14 
     236 
     237x1=13  
     238x2=32  
     239y1=x1  
     240y2=x2 
    216241; 
    217242pdf, param,trop,mio,mao,inx,mis,mas,iny, $ 

  • 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, $ 

  • trunk/src/paper01/fig6/statistics.pro

    r85 r94  
    6262;; mean1= mean of mooring variable      std1=std of mooring variable 
    6363;; mean2= mean of erai variable         std2=std of erai variable 
    64 x=moor & y=erai 
    65 ind1=where(finite(x,/nan)) & ind2=where(finite(y,/nan)) 
     64x=moor  
     65y=erai 
     66ind1=where(finite(x,/nan))  
     67ind2=where(finite(y,/nan)) 
    6668ind=inter(ind1,ind2) 
    6769if (ind ne -1 ) then begin 
    68    x(ind)=!Values.f_nan &y(ind)=!Values.f_nan 
     70   x(ind)=!Values.f_nan  
     71   y(ind)=!Values.f_nan 
    6972endif 
    7073 
    7174ind=where(finite(x)) 
    72 x=x(ind) & y=y(ind) 
     75x=x(ind)  
     76y=y(ind) 
    7377 
    7478nn=n_elements(ind) 
    … …  
    7680cor=correlate(x,y) 
    7781xxx=moment(x) 
    78 mean1=xxx(0) & std1=sqrt(xxx(1)) 
     82mean1=xxx(0)  
     83std1=sqrt(xxx(1)) 
    7984xxx=moment(y) 
    80 mean2=xxx(0) & std2=sqrt(xxx(1)) 
    81 bias=mean2-mean1 & std=std2/std1 
     85mean2=xxx(0)  
     86std2=sqrt(xxx(1)) 
     87bias=mean2-mean1  
     88std=std2/std1 
    8289rmsd=sqrt(total((x-y)*(x-y))/n_elements(x)) 
    8390end 
Note: See TracChangeset for help on using the changeset viewer.