---+ Example Atmosphere and related code provided by Friends

Code for Subtracting the atmosphere using PCA, from Tom Crawford (SPT) 2006.

This code was given to Tom from Sunil (Caltech) who uses it with Bolocam. It DOES NOT run on IDL, but instead uses the open source version called Octave. See tom's comments about porting it to IDL.
-------- Original Message --------
Subject: 	more atmo stuff
Date: 	Fri, 21 Apr 2006 19:23:44 -0500 (CDT)
From: 	Tom Crawford 
Reply-To: 	Tom Crawford 
To: 	Matt Dobbs 
Hi Matt -

Here's the PCA routine and some stuff you'll need to run the array-mode
projector I sent you in the other e-mail.  The PCA routine makes an
external call to an octave routine called "exteigen" (which I've attached)
to diagonalize the covariance matrix and find its eigenvectors and
eigenvalues.  This is only because octave is much faster than IDL at
linear algebra -- it may be easier to start out doing this the slow way in 
IDL, and you can do this by uncommenting some lines in pca_clean_tc.pro 
and commenting out some others (it should be clear which ones; if not, 
just let me know).

Have fun!
Tom


Code for simulating, measureing, and subtracting the atmosphere using the Triangular Fourier method, from Tom Crawford (SPT) 2004.

My array-mode-projection method is buried inside that 
"spt_scan_sim_paredown_temp.pro" routine.  Unfortunately, it doesn't exist 
in any much cleaner state.

-------- Original Message --------
Subject: 	atmo code
Date: 	Fri, 08 Oct 2004 18:24:33 -0500 (CDT)
From: 	Tom Crawford 
Reply-To: 	Tom Crawford 
To: 	Matt Dobbs 
Hi Matt -

I've attached the code I use to make the screen of atmospheric
fluctuations that get blown past / scanned by the detector array.  I've
also attached the huge, awful routine that does the observation and the
atmospheric subtraction, because buried somewhere in there is how I
incorporate the atmosphere screen into sky observations.  The atmospheric
screen creator is pretty simple and should make sense if you've read the
Lay & Halverson paper, but the other piece of kludged-up spaghetti code 
may be impenetrable, in which case feel free to ask questions.


Enjoy.

-- 
Tom Crawford
University of Chicago
Department of Astronomy & Astrophysics
773-834-2103

Topic attachments
I Attachment Action Size Date Who Comment
Unknown file formatEXT exteigen2 manage 0.4 K 2006-05-11 - 23:01 MattDobbs  
Unknown file formatpro make_atmo_strip.pro manage 4.5 K 2006-05-11 - 22:53 MattDobbs  
Unknown file formatpro make_projmodes.pro manage 3.3 K 2006-05-11 - 23:00 MattDobbs  
Unknown file formatpro one_on_f_noise.pro manage 5.3 K 2006-05-11 - 23:00 MattDobbs  
Unknown file formatpro pca_clean_tc.pro manage 8.8 K 2006-05-11 - 23:01 MattDobbs  
Unknown file formatpro spt_scan_sim_paredown_temp.pro manage 9.4 K 2006-05-11 - 22:53 MattDobbs  

This topic: APEX_SZ > WebHome > 2006_AtmosphereSubtraction > 2006_AtmosphereCode Topic revision: r1 - 2006-05-11 - MattDobbs
© 2020 Winterland Cosmology Lab, McGill University, Montréal, Québec, Canada