CosmoCoffee

Posted: **October 10 2004**

I'm unfamiliar with dealing with CMB maps.
I would like to be able to:

a) Take say one of the publicly available WMAP map files, and
extract the sky temperature as a function of pixel (and I'd really like
the pixel coordinates too..)

b) Take the FITS format output of say the Healpix routine anafast
and read the C_l's and alm's. I realize there are f90 routines
provided with Healpix which do this, but i'm using OS X and was
never able to get the relevant f90 routines to compile with all the
dependencies in the makefile.

What is the best way to go about doing this (short of being a FITS guru master)?

Posted: **October 11 2004**

Try something like this:

Code: Select all

SUBROUTINE LOAD_HEAL_FILE&#40;filename, DATAR, nside, npixtot&#41;
use healpix_types
use fitstools
use pix_tools
use alm_tools
use utilities 

implicit none
CHARACTER&#40;LEN=*&#41;, PARAMETER &#58;&#58; code = 'tlike'
CHARACTER &#40;LEN=*&#41; filename
INTEGER, INTENT &#40;OUT&#41;, OPTIONAL &#58;&#58; nside, npixtot
INTEGER&#40;I4B&#41; &#58;&#58; npixtoto, nmapso, info, status
INTEGER&#40;I4B&#41; &#58;&#58; ordering, obs_npix, nsideo, mlpol, type, polarisation

REAL&#40;SP&#41;, DIMENSION&#40;&#58;,&#58;&#41;, POINTER &#58;&#58; DATAR


     IF &#40;debug.gt.0&#41; PRINT *,'Reading file ', filename
     npixtoto = getsize_fits&#40;filename, nmapso, ordering, obs_npix, nsideo, mlpol, type, polarisation &#41;
     IF &#40;debug.gt.1&#41;   PRINT *,'nsideo=', nsideo
     ALLOCATE&#40;DATAR&#40;0&#58;12*nsideo**2-1,1&#58;2&#41;, stat = status&#41;
     IF &#40;status /= 0&#41; CALL die_alloc&#40;code,'alloc '//filename&#41;
     CALL input_map&#40;filename, DATAR, npixtoto, nmapso, fmissval=-1.6375e30_sp&#41;

     IF &#40;ordering == 2 &#41; then 
        IF &#40;debug.gt.0&#41; PRINT *,'Converting NEST to RING...'
        CALL convert_nest2ring&#40;nsideo,DATAR&#40;&#58;,1&#41;&#41;
        CALL convert_nest2ring&#40;nsideo,DATAR&#40;&#58;,2&#41;&#41;
        ordering =1 
     ENDIF
     if &#40;present&#40;npixtot&#41;&#41; npixtot=npixtoto
     if &#40;present&#40;nside&#41;&#41; nside = nsideo
end

with

Code: Select all

   REAL&#40;SP&#41;, DIMENSION&#40;&#58;,&#58;&#41;, POINTER &#58;&#58; map_t

and to get positions you should do something like:

Code: Select all

   DO II=0,NPIXTOT-1
         CALL PIX2VEC_RING&#40;nside, ii, VEC1&#41;
           !! Now VEC1 contains point on a sphere oorresponding to the
           !!   current pixel

Posted: **October 12 2004**

If you're desperate for a FORTRAN 95 compiler for MacOS X, you can get a binary for one here:

http://www.g95.org/

I haven't used it much but CAMB at least compiles and seems to run OK. If you want to wait for a less experimental compiler, gcc 4.0.0 is scheduled to include FORTRAN 95; it's tentatively scheduled for release early next year.

Or am I misinterpreting your compiler difficulties?

CosmoCoffee

Reading Healpix and FITS files

Reading Healpix and FITS files