zonca

nside = 128
import healpy as hp
import numpy as np
tmap = np.random.normal(scale=1, size=12*nside**2)
qmap = np.random.normal(scale=1.2, size=12*nside**2)
umap = np.random.normal(scale=.8, size=12*nside**2)

cltt, clee, clbb, clte, cleb, cltb = hp.anafast( [tmap, qmap, umap], lmax=2*nside )

exp = 4.*np.pi/12./nside**2

zonca

zonca

zonca

zonca

import numpy as np;import cPickle;from math import sin, cos, sqrt

dpcconv = np.load("dipole_27M_od100.npy")
dpcpnt = cPickle.load(open("pointing_DPC_LFI27M_100_DX9.pkl", 'rb'))

# sample 20000 of OD 100 LFI28S
p = 20000
scet = dpcpnt['sampleSCET'][p] # 1629538470557322

# relativistic add of:

zonca

#!/usr/bin/env python

#Specify this:
username = ''


"""
Get an IPython notebook running on Carver. 

Usage:

zonca

import numpy as np
import cPickle
import pycfitsio
from collections import defaultdict
import healpy

from planck import Planck

from testenv import remix, todtools
from planck.pointingtools import AHF_btw_OBT

zonca

See the references below directly

zonca

from mayavi import mlab
import numpy as np
import healpy as hp

# load a WMAP map
m = hp.read_map("data/wmap_band_iqumap_r9_7yr_W_v4.fits", 0) * 1e3 # muK
nside = hp.npix2nside(len(m))

vmin = -1e3; vmax = 1e3

zonca

#! /bin/bash
# Generated by configure.
# Run this file to recreate the current configuration.
# Compiler output produced by configure, useful for debugging
# configure, is in config.log if it exists.

debug=false
ac_cs_recheck=false
ac_cs_silent=false