Ieri ho importato le mappe da seaborg sul mio laptop. Ho cominciato a scrivere il programma che le doveva portare a bassa risoluzione (30') come il canale a 33GHz.
Le mappe fatte da springtide per ora sono in coordinate Ecliptiche. Potremmo trovarci a dover guardare mappe in coordinate diverse dalle galattiche per cui un possibile improuvement per CCA potrebbe essere quello di lavorare con mappe in qualsiasi tipo di coordinate. Siccome CCA estae patches a partire da coordinate galattiche sipotrebbe mettere un if iniziale e nel caso coord ne 'G' fa la conversione.
Penso che sia piu' facile e corretto passare per gli alm. Dopo averli calcolati si riproietta la mappa nelle nuove coordinate. Non e' una procedura perfetta perche' bisognerebbe rifare la mappa da capo col mapmaking ma questo non e' sempre possibile.
Questi sarebbero gli steps:
calcolare gli alm con synfast o con map2alm
call coordsys2euler_zyz(2000.0_dp, 2000.0_dp, 'E', 'G', psi, theta, phi)
call rotate_alm(nside, alm_TGC, psi, theta, phi)
Oggi sto facendo la rotazione scrivendo un programmino IDL che mi cambia le coordinate in coda al programma che abbassa l;a risoluzione. Federico mi ha detto come lo fa di solito lui e sto tentando di riprodurre la sua strategia.
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment