sukye

north idrop

October 3, 2007 on 11:08 am | In Uncategorized | Comments Off

myb=NB_ACS_D10Z_ASCII
myi=NI_ACS_D10Z_ASCII
myv=NV_ACS_D10Z_ASCII
myz=NZ_ACS_D10Z_ASCII

b=myb.field034
i=myi.field034
v=myv.field034
z=myz.field034

sz=myz.field036
sv=myv.field036
si=myi.field036
sb=myb.field036

nb=myb.field037
ni=myi.field037
nv=myv.field037
nz=myz.field037

flagz=myz.field030
flagv=myv.field030
flagi=myi.field030
flagb=myb.field030

radb=myb.field024
radi=myi.field024
radv=myv.field024
radz=myz.field024

clsz=myz.field028
clsv=myv.field028
clsi=myi.field028
clsb=myb.field028

.run
openw,lun,”ni_drop.txt”,/get_lun,width=85
For j=0u,39431u Do Begin
if i[j]-z[j] gt 1.30 and sz[j]/nz[j] gt 5.0 and (sv[j]/nv[j] lt 2.0 or v[j]-z[j] gt 2.80) $
and radz[j] lt 3.40 $
and flagz[j] le 8.0 $
and clsz[j] ge 0.80 $
and z[j] le 25.50 $
then printf,lun,myb.field001[j],i[j],z[j],i[j]-z[j],myz.field042[j],j
Endfor
free_lun, lun
End

——————————————–
.run
openw,lun,”nv_drop.txt”,/get_lun,width=105
For j=0u,39431u Do Begin
if (v[j]-i[j] gt 2.0 or v[j]-i[j]-0.890*(i[j]-z[j]) gt 1.470) $
and v[j]-i[j] gt 1.20 and i[j]-z[j] lt 1.30 $
and si[j]/ni[j] gt 5.0 and sb[j]/nb[j] lt 3.0 $
and (sb[j]/nb[j] lt 3.0 or b[j]-i[j]-v[j]+i[j] gt 1.0) $
and radi[j] lt 3.40 $
and flagi[j] le 8.0 $
and clsi[j] ge 0.90 $
and i[j] le 25.50 $
then printf,lun,myb.field001[j],v[j],i[j],z[j],v[j]-i[j],i[j]-z[j],j
Endfor
free_lun, lun
End