#!/bin/env python
import smallfield
import healmap
import healpy as hp
import numpy as np
if __name__ == "__main__":
use_split = 'D'
use_nside = 512
out_nside = 32
outfname = 'smallfield_0016.fits'
dl_lcdm,lc = smallfield.calc_lcdm_dl()
print dl_lcdm
h1,h2 = healmap.get_input_maps(use_split,use_nside=use_nside,dir='/n/bicepfs3/general/input_maps/planck_pr2/')
lat,lon,ipix = smallfield.gen_field_centers(nside=out_nside,minlat=-1.0)
i0 = 0
if i0==0:
rb = np.zeros([12*out_nside*out_nside])
sb = 0.0 * rb
re = 0.0 * rb
se = 0.0 * rb
else:
rb = hp.read_map(outfname,field=0)
sb = hp.read_map(outfname,field=1)
re = hp.read_map(outfname,field=2)
se = hp.read_map(outfname,field=3)
for i in xrange(i0,lat.size):
if (rb[ipix[i]]!=0.) or (sb[ipix[i]]!=0.):
#!/opt/local/bin/python2.7 import smallfield import healmap import healpy as hp import numpy as np if __name__ == "__main__": use_split = 'D' use_nside = 512 h1,h2 = healmap.get_input_maps(use_split,use_nside=use_nside,dir='..') sfac = 0.0395 squ = 10 st = 200 rot=[0,90,0] # hp.orthview(h1.map[0,:]*sfac,rot=rot,nest=False,title="Split 1 T",min=0,max=st) hp.orthview(h1.map[1,:]*sfac,rot=rot,nest=False,title="Split 1 Q",min=-squ,max=squ) hp.orthview(h1.map[2,:]*sfac,rot=rot,nest=False,title="Split 1 U",min=-squ,max=squ) # hp.orthview(h2.map[0,:]*sfac,rot=rot,nest=False,title="Split 2 T",min=0,max=st) hp.orthview(h2.map[1,:]*sfac,rot=rot,nest=False,title="Split 2 Q",min=-squ,max=squ) hp.orthview(h2.map[2,:]*sfac,rot=rot,nest=False,title="Split 2 U",min=-squ,max=squ) # plt.show()
