sukye

corners of CDFS GOODS

03 32 39.1, -27 38 51
03 31 56.4, -27 42 16
03 33 4.2, -27 54 15
03 32 21.5, -27 57 42

typo

myb=SB_ACS_D10Z_ASCII
myi=SI_ACS_D10Z_ASCII
myv=SV_ACS_D10Z_ASCII
myz=SZ_ACS_D10Z_ASCII

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

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

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

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

.run
openw,lun,’s_v_drop.txt’,/get_lun,width=105
For j=0u,33954u 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 26.5 $
then printf,lun,myi.field001[j],myi.field004[j],myv.field034[j],myi.field034[j],myz.field042[j],j
Endfor
free_lun, lun
End

synphot

5 parameters:
• Observation mode (obsmode) - Passband to be calculated
• Spectrum (spectrum) - Spectrum to be calculated
• Form (form) - Units of the output data; for example, counts or magnitudes
• Wavelength table (wavetab) - Wavelength grid on which passband
and spectrum will be calculated
• Reference data parameters (refdata) - Pointer to reference data
for instrument graph and component lookup tables

plband acs,hrc,f555w acs_f555w.fits
plspec “” crcalspec$hz44_stis_001.fits flam
calcphot wfpc2,1,f439w,a2d7 crcalspec$hz44_stis_001.fits counts
calcspec “rn(bb(5000),band(wfpc2,f555w),100,counts)” rnbb.fits
calcspec “rn(pl(5500,-2,fnu),band(v),1.e-14,flam)*ebmv(0.1)” pl.fits form=fnu
calcphot wfpc2,3,f439w “rn(crgridkc96$sb_template.fits,band(v),18,vegamag)” counts
??”So this galaxy produces a total countrate of 42.4 e–/sec, which means
we’ll need an exposure time of at least 235 seconds in order to approach an
accuracy of 1% in measuring its total light.”
calcphot “band(wfpc2,4,f439w)-band(wfpc2,4,f555w)” “bb(5000)*ebmv($0)” obmag vzero=”0,0.25,0.5″
calcphot “acs,wfc1,f850lp,aper#0.1″ “rn(icat(k93models,5250,0.0,5.0),band(v),24,vegamag)” counts func=effstim

sy> plband johnson,v
sy> plband wfpc2,3,f555w norm+ app+ ltype=dashed
sy> plband wfpc2,3,f569w norm+ app+ ltype=dotted

–> plspec “” crgridkc96$sb_template.fits flam
–> plspec “” crgridkc96$elliptical_template.fits flam app+ ltype=dotted

countrate
sgraph

calcphot
http://www.stsci.edu/hst/wfpc2/software/wfpc2_extended.html