Sky model.pro

ডিপথট থেকে
; *** Create a n*n grid and find the RA & DEC of the centres of every grid element ***
pro SkyModel,n,gap,ra,dec

x=dblarr(n,n)
y=dblarr(n,n)
c=rnd(n/2,/down)

radectodeg,16,14,0.044,54,29,59.90,radec

for i=0,c do begin
  x[c+i,*]=radec[0]+gap*i
  x[c-i,*]=radec[0]-gap*i
  y[*,c+i]=radec[1]+gap*i
  y[*,c-i]=radec[1]-gap*i
endfor

ra=dblarr(n*3,n)
dec=dblarr(n*3,n)

for i=0,n-1 do begin
  radec, x[i,0],y[0,i],ihr,imin,xsec,ideg,imn,xsc
  print,ihr,imin,xsec,ideg,imn,xsc
  ra[3*i,*]=ihr
  ra[3*i+1,*]=imin
  ra[3*i+2,*]=xsec
  for j=0,n-1 do begin
    dec[3*j,i]=ideg
    dec[3*j+1,i]=imn
    dec[3*j+2,i]=xsc
  endfor
endfor

end

Associated programs