Home > atmlab > demos > qarts_abstable_demo.m

qarts_abstable_demo

PURPOSE ^

QARTS_ABSTABLE_DEMO A simple demonstration calculation of abs. table

SYNOPSIS ^

function A = qarts_abstable_demo

DESCRIPTION ^

 QARTS_ABSTABLE_DEMO   A simple demonstration calculation of abs. table

    An absorption table is created. This example case gives the absorption
    for the centre channel of the SSM/I channels. Absorption for N2, O2 and
    H2O is provided, where some well known complete absorption models are
    used. The table covers rougly 0 to 25 km, in steps of 250 m. The
    temperature covers +-35K (in steps of 10K) around the Fascod sub-arctic
    winter scenario. 

 FORMAT   A = qarts_abstable_demo
        
 OUT   A  Absorption table.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

qarts_abstable_demo.m

SOURCE CODE ^

0001 % QARTS_ABSTABLE_DEMO   A simple demonstration calculation of abs. table
0002 %
0003 %    An absorption table is created. This example case gives the absorption
0004 %    for the centre channel of the SSM/I channels. Absorption for N2, O2 and
0005 %    H2O is provided, where some well known complete absorption models are
0006 %    used. The table covers rougly 0 to 25 km, in steps of 250 m. The
0007 %    temperature covers +-35K (in steps of 10K) around the Fascod sub-arctic
0008 %    winter scenario.
0009 %
0010 % FORMAT   A = qarts_abstable_demo
0011 %
0012 % OUT   A  Absorption table.
0013 
0014 % 2010-04-15   Created by Patrick Eriksson.
0015 
0016 
0017 function A = qarts_abstable_demo
0018 
0019 arts_xmldata_path = atmlab( 'ARTS_XMLDATA_PATH' );
0020 %
0021 if isnan( arts_xmldata_path )
0022   error('You need to ARTS_XMLDATA_PATH to run this example.');
0023 end
0024 %
0025 fascod = fullfile( arts_xmldata_path, 'planets', 'Earth', 'Fascod' );
0026 
0027 
0028 %= Init Q structures
0029 %
0030 Q  = qarts;
0031 %
0032 Q.INCLUDES              = { fullfile( 'ARTS_INCLUDES', 'general.arts' ), ...
0033                             fullfile( 'ARTS_INCLUDES', 'agendas.arts' ), ...
0034                             fullfile( 'ARTS_INCLUDES', 'continua.arts' ), ...
0035                             fullfile( 'ARTS_INCLUDES', 'planet_earth.arts' ) };
0036 
0037 
0038 %= Atmospheric species and absorption models
0039 %
0040 Q.ABS_SPECIES(1).TAG{1}  = 'N2-SelfContStandardType';
0041 Q.ABS_SPECIES(2).TAG{1}  = 'O2-PWR93';
0042 Q.ABS_SPECIES(3).TAG{1}  = 'H2O-PWR98';
0043 
0044 
0045 %= Set-up atmosphere and surface
0046 %
0047 Q.ATMOSPHERE_DIM         = 1;
0048 %
0049 Q.P_GRID                 = z2p_simple( [0:250:40e3]' );
0050 
0051 
0052 %= Reference atmosphere
0053 %
0054 Q.RAW_ATMOSPHERE         = fullfile( fascod, 'subarctic-winter', ...
0055                                              'subarctic-winter' );
0056 %
0057 % Needed, but dummy variables
0058 Q.REFELLIPSOID           = ellipsoidmodels( 'SphericalEarth' );
0059 Q.Z_SURFACE              = 500;
0060 
0061 
0062 %= Absorption settings
0063 %
0064 % For 22 GHz channel, the centre of the H2O transition is avoided by taking
0065 % a frequency half way to band edge.
0066 Q.F_GRID                 = [19.35 22.235-0.06 37.0 85.5]'*1e9;
0067 Q.STOKES_DIM             = 1;
0068 %
0069 Q.ABS_LINES_FORMAT       = 'None';
0070 Q.ABS_LINES              = [];
0071 Q.ABS_NLS                = [];
0072 
0073 
0074 %= Calculate absorption
0075 %
0076 Q                        = qarts_abstable( Q, 8, 35 ); 
0077 A                        = arts_abstable( Q );

Generated on Mon 15-Sep-2014 13:31:28 by m2html © 2005