110 likes | 214 Views
3rd ACS Workshop and advanced course ESO Garching Headquarter, January 15-19, 2006. A tmospheric T ransmission at M icrowaves ( ATM ) C++ implementation within ALMA TelCal Subsystem. Juan R. Pardo 1 (1) Consejo superior de Investigaciones Científicas (Spain).
E N D
3rd ACS Workshop and advanced course ESO Garching Headquarter, January 15-19, 2006 Atmospheric Transmission at Microwaves (ATM) C++ implementation within ALMA TelCal Subsystem Juan R. Pardo1 (1) Consejo superior de Investigaciones Científicas (Spain) Physical model to implement (described last Monday) Current C++ implementation Architecture and design review
1. Physical model to implement: Atmospheric Refractivity at ALMA frequencies 1a. Imaginary Part (absorption) Chajnantor zenith transmission for 0.5 mm H2O / Water lines / Oxygen lines / ozone lines / H2O-foreign / N2-N2 + N2-O2 + O2-O2 1b. Real Part (phase delay)
1a. ALMA needs related to this component • Able to Guess atmospheric T, P, gas profiles from available site data(Tground, Pground, local humidity, complementary soundings) • Given an atmospheric profile is known, then able to obtain: • Atmospheric opacities sorted out by component (dry, wet, individual molecules, continuum-like terms, etc..). • Phase delays sorted out by component, but specially those related to H2O. • Simulate atmospheric brightness temperatures along a given propagation path (the easy-to-get observable) • Inversion capabilities. Given atmospheric brightness temperatures can be measured, then it should be able to: • Retrieve atmospheric parameters (mainly the H2O column). • Retrieve hardware implementation parameters such as the coupling to the sky of water vapor radiometers. • Use the retrieved information to provide correction parameters at the frequencies of the current astronomical observation. At all ALMA 's
2. C++ implimentation history • Starting point: Fortran code developed during 20 years. • Contract with ESO to develop ATM for alma. • Fortran library created to encapsulate "fundamental" physics of the problem. • C++ interface developed specifically for ALMA needs. • Updated via CVS (TelCal subsystem) • Doxygen documented. • Test examples provided using real FTS and WVR data. • Work done within the TelCal working group.
3. Architecture and design review: Old ATM fortran code Ground Temperature Tropospheric Lapse Rate Ground Pressure Rel. humidity (ground) Water vapor scale height Primary pressure step Pressure step factor Altitude of site Top of atmospheric profile ATM_telluric 1st guess of water vapor column Layer thickness (NPP) Number of atmospheric layers, NPP Layer pressure (NPP) Layer temperature (NPP) Layer water vapor (NPP) Layer_O3 (NPP) Layer CO (NPP) Layer N2O (NPP) User Other sources to obtain atmospheric parameters ATM_st76 ATM_rwat DATA_xx_lines DATA_xx_index xx all species ABSORPTION COEFS. kh2o_lines (NPP) kh2o_cont (NPP) ko2_lines (NPP) kdry_cont (NPP) kminor_gas (NPP) PHASE FACTORS total_dispersive_phase (NPP) total_nondisp_phase (NPP) INI_singlefreq ABS_h2oABS_o2 ABS_co ABS_o3_161618 ABS_h2o_v2ABS_o2_vib ABS_n2o ABS_o3_161816 ABS_hdoABS_16o17o ABS_o3 ABS_o3_nu1 ABS_hh17oABS_16o18o ABS_o3_161617 ABS_o3_nu2 ABS_hh18oABS_cont ABS_o3_161716 ABS_o3_nu3 kabs (NPP) Frequency Air mass Bgr. Temp. Applications: INV_telluric for WVR and FTS data (available in current release) RT_telluric Atmospheric radiance INTERFACE LEVEL I/O Parameters DEEP LEVEL(LIBRARY)
3. Architecture and design review: Basic C++ structure: Collaboration diagram between the most important classes Class AtmProfile: Profiles of physical conditions & chemical abundances Class AbsorptionPhaseProfile: Profiles of refractive index for an array of frequencies Class WaterVaporRadiometer: Water Vapor Radiometer system in place for phase correction Class SkyStatus: Relevant atmospheric information for antenna operations
3. Architecture and design review: Comments Class SpectralGrid: replicates essentially what would a Class of the receiver or autocorrelator components. Class Temperautre Class Pressure Class NumberDensity, Class Length... Classes for the physical parameters relevant to ATM have been created within this component with a namespace atm to avoid confussion. Interfaces with other components need to take this into account. Standard classes for this physical parameters within the whole software? Let's have a look to the documentation
3. Architecture and design review: Tests Using real atmospheric transmission curves measured with a Fourier Transform Spectrometer
3. Architecture and design review: Tests 183 GHz water vapor radiometry (preferred method for phase correction in ALMA) Taking an average Precipitable Water Vapor amount of 0.5 mm, we have: 23 mk/μm60 mk/μm 173 mk/μm
GOES-10 Water vapor 350 m opacity meter
3. Architecture and design review: Tests 183 GHz water vapor radiometry with actual phase correction Observing Frequency: 230.5 GHz, PWV to phase conversion factor: 1.944 deg/m a Time since the beginning of the observation (min) on Nov. 25, 2001