Home > atmlab > forwardmodel > isotopologue_data.m

isotopologue_data

PURPOSE ^

ISOTOPOLUGUE_DATA Returns data for different isotopologues

SYNOPSIS ^

function data = isotopologue_data(isos,fields)

DESCRIPTION ^

 ISOTOPOLUGUE_DATA   Returns data for different isotopologues

    Isotopologues are named following the ARTS scheme, such as O3-666
    for standard ozone.

    The following data fields are defined:
       isofrac   Standard values for relative fraction of isotopologue.
       weight    Molecular weight, such as 47.98... for standard ozone.
       hitran    HITRAN number, with molecular isotopologue numbers joined
                 (such as 11 for H2O-161). Present numbers are valid for
                 HITRAN 2000.
       jpl       JPL tag number(s). This could be a vector (when JPL 
                 includes vibrational levels).

    NaN is returned for missing data.

    An example:  isotopologue_data( {'02-66','O3-666'}, 'isofrac' )

    If no input argumnets are given, all data are returned. A field with
    name *name* is then added with name of all defined isotopolugues.
    A list of all isotopolugues is obtained from output as
       data(:).name 

 FORMAT   data = isotopologue_data(isos,fields)
        
 OUT   data   Extracted data. A struct array, with fields following *fields*.
 IN    isos   Name of isotopolugues for which data shall be extracted.
              A cell string array, or a string for a single isotopologue.
       fields Data fields to extract.
              A cell string array, or a string for a single field.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

isotopologue_data.m

SOURCE CODE ^

0001 % ISOTOPOLUGUE_DATA   Returns data for different isotopologues
0002 %
0003 %    Isotopologues are named following the ARTS scheme, such as O3-666
0004 %    for standard ozone.
0005 %
0006 %    The following data fields are defined:
0007 %       isofrac   Standard values for relative fraction of isotopologue.
0008 %       weight    Molecular weight, such as 47.98... for standard ozone.
0009 %       hitran    HITRAN number, with molecular isotopologue numbers joined
0010 %                 (such as 11 for H2O-161). Present numbers are valid for
0011 %                 HITRAN 2000.
0012 %       jpl       JPL tag number(s). This could be a vector (when JPL
0013 %                 includes vibrational levels).
0014 %
0015 %    NaN is returned for missing data.
0016 %
0017 %    An example:  isotopologue_data( {'02-66','O3-666'}, 'isofrac' )
0018 %
0019 %    If no input argumnets are given, all data are returned. A field with
0020 %    name *name* is then added with name of all defined isotopolugues.
0021 %    A list of all isotopolugues is obtained from output as
0022 %       data(:).name
0023 %
0024 % FORMAT   data = isotopologue_data(isos,fields)
0025 %
0026 % OUT   data   Extracted data. A struct array, with fields following *fields*.
0027 % IN    isos   Name of isotopolugues for which data shall be extracted.
0028 %              A cell string array, or a string for a single isotopologue.
0029 %       fields Data fields to extract.
0030 %              A cell string array, or a string for a single field.
0031 
0032 % 2006-11-26   Created by Patrick Eriksson.
0033 
0034 
0035 function data = isotopologue_data(isos,fields)
0036 
0037 if nargin == 1  |  nargin > 2
0038   error( 'The function requires 0 or 2 input arguments.' );
0039 end
0040 
0041 
0042 %- Modification of input data
0043 %
0044 if nargin
0045   if ischar( isos )
0046     isos = { isos };
0047   end
0048   %
0049   if ischar( fields )
0050     fields = { fields };
0051   end
0052 end
0053 
0054 
0055 
0056 %-----------------------------------------------------------------------------
0057 %
0058 % Isotopologue data
0059 %
0060 % The data were copied from the ARTS-1.0 file species_data.cc 2006-11-26
0061 %
0062 % Iso. names must here use '_' instead of '-' as Matlab will take '-' as a
0063 % minus sign. This for seperator between molecule name and isotopologue
0064 % number.
0065 % The + in ions are here replaced by 'x'.
0066 % All fields must be specified, for all isotopologues, to get nargin=0
0067 % version to work.
0068 %
0069 %-----------------------------------------------------------------------------
0070 
0071 %D..isofrac = ;
0072 %D..weight  = ;
0073 %D..hitran  = ;
0074 %D..jpl     = [   ];
0075 
0076 
0077 
0078 %---
0079 %--- BrO
0080 %---
0081 %
0082 D.BrO_96.isofrac   = 0.50582466;
0083 D.BrO_96.weight    = 95;
0084 D.BrO_96.hitran    = NaN;
0085 D.BrO_96.jpl       = [ 95001 ];
0086 %
0087 D.BrO_16.isofrac   = 0.49431069;
0088 D.BrO_16.weight    = 97;
0089 D.BrO_16.hitran    = NaN;
0090 D.BrO_16.jpl       = [ 97001 ];
0091 
0092 
0093 
0094 %---
0095 %--- C2H2
0096 %---
0097 %
0098 D.C2H2_1221.isofrac   = .977599E+00;     
0099 D.C2H2_1221.weight    = 26.015650;
0100 D.C2H2_1221.hitran    = 261;
0101 D.C2H2_1221.jpl       = [ NaN ];
0102 %
0103 D.C2H2_1231.isofrac   = 2.19663E-02;     
0104 D.C2H2_1231.weight    = 27.019005;
0105 D.C2H2_1231.hitran    = 261;
0106 D.C2H2_1231.jpl       = [ NaN ];
0107 
0108 
0109 
0110 %---
0111 %--- C2H4
0112 %---
0113 %
0114 D.C2H4_221.isofrac   = .977294E+00;     
0115 D.C2H4_221.weight    = 28.031300;
0116 D.C2H4_221.hitran    = 381;
0117 D.C2H4_221.jpl       = [ NaN ];
0118 %
0119 D.C2H4_231.isofrac   = .219595E-01;     
0120 D.C2H4_231.weight    = 29.034655;
0121 D.C2H4_231.hitran    = 382;
0122 D.C2H4_231.jpl       = [ NaN ];
0123 
0124 
0125 
0126 %---
0127 %--- C2H6
0128 %---
0129 %
0130 D.C2H6_1221.isofrac   = .976990E+00;     
0131 D.C2H6_1221.weight    = 30.046950;
0132 D.C2H6_1221.hitran    = 271;
0133 D.C2H6_1221.jpl       = [ NaN ];
0134 
0135 
0136 
0137 %---
0138 %--- CH3Cl
0139 %---
0140 %
0141 D.CH3Cl_215.isofrac   = .748937E+00;     
0142 D.CH3Cl_215.weight    = 49.992328;
0143 D.CH3Cl_215.hitran    = 241;
0144 D.CH3Cl_215.jpl       = [ 50007 ];
0145 %
0146 D.CH3Cl_217.isofrac   = .239491E+00;     
0147 D.CH3Cl_217.weight    = 51.989379;
0148 D.CH3Cl_217.hitran    = 242;
0149 D.CH3Cl_217.jpl       = [ 52009 ];
0150 
0151 
0152 
0153 %---
0154 %--- CH3CN
0155 %---
0156 %
0157 D.CH3CN_211124.isofrac   = 0.97366840;
0158 D.CH3CN_211124.weight    = 41;
0159 D.CH3CN_211124.hitran    = NaN;
0160 D.CH3CN_211124.jpl       = [ 41001 ];
0161 %
0162 D.CH3CN_311124.isofrac   = 0.011091748;
0163 D.CH3CN_311124.weight    = 42;
0164 D.CH3CN_311124.hitran    = NaN;
0165 D.CH3CN_311124.jpl       = [ 42006 ];
0166 %
0167 D.CH3CN_211134.isofrac   = 0.011091748;
0168 D.CH3CN_211134.weight    = 42;
0169 D.CH3CN_211134.hitran    = NaN;
0170 D.CH3CN_211134.jpl       = [ 42007 ];
0171 %
0172 D.CH3CN_211125.isofrac   = 0.0036982817;
0173 D.CH3CN_211125.weight    = 42;
0174 D.CH3CN_211125.hitran    = NaN;
0175 D.CH3CN_211125.jpl       = [ 42001 ];
0176 %
0177 D.CH3CN_211224.isofrac   = 0.00044977985;
0178 D.CH3CN_211224.weight    = 42;
0179 D.CH3CN_211224.hitran    = NaN;
0180 D.CH3CN_211224.jpl       = [ 42008 ];
0181 
0182 
0183 
0184 %---
0185 %--- CH4
0186 %---
0187 %
0188 D.CH4_211.isofrac   = .988274E+00;     
0189 D.CH4_211.weight    = 16.031300;
0190 D.CH4_211.hitran    = 61;
0191 D.CH4_211.jpl       = [ NaN ];
0192 %
0193 D.CH4_311.isofrac   = 1.11031E-02;     
0194 D.CH4_311.weight    = 17.034655;
0195 D.CH4_311.hitran    = 62;
0196 D.CH4_311.jpl       = [ NaN ];
0197 %
0198 D.CH4_212.isofrac   = 6.15751E-04;     
0199 D.CH4_212.weight    = 17.037475;
0200 D.CH4_212.hitran    = 63;
0201 D.CH4_212.jpl       = [ 17003 ];
0202 
0203 
0204 %---
0205 %--- Cl2O2
0206 %---
0207 %
0208 D.Cl2O2_565.isofrac   = 0.57016427;     
0209 D.Cl2O2_565.weight    = 102;
0210 D.Cl2O2_565.hitran    = NaN;
0211 D.Cl2O2_565.jpl       = [ 102001 ];
0212 %
0213 D.Cl2O2_765.isofrac   = 0.36982818;     
0214 D.Cl2O2_765.weight    = 104;
0215 D.Cl2O2_765.hitran    = NaN;
0216 D.Cl2O2_765.jpl       = [ 104001 ];
0217 
0218 
0219 
0220 %---
0221 %--- ClO
0222 %---
0223 %
0224 D.ClO_56.isofrac   = .755908E+00;     
0225 D.ClO_56.weight    = 50.963768;
0226 D.ClO_56.hitran    = 181;
0227 D.ClO_56.jpl       = [ 51002, 51003 ];
0228 %
0229 D.ClO_76.isofrac   = .241720E+00;     
0230 D.ClO_76.weight    = 52.960819;
0231 D.ClO_76.hitran    = 182;
0232 D.ClO_76.jpl       = [ 53002, 53006 ];
0233 
0234 
0235 
0236 %---
0237 %--- ClONO2
0238 %---
0239 %
0240 D.ClONO2_5646.isofrac   = .749570E+00; 
0241 D.ClONO2_5646.weight    = 96.956672;
0242 D.ClONO2_5646.hitran    = 351;
0243 D.ClONO2_5646.jpl       = [ 97002 ];
0244 %
0245 D.ClONO2_7646.isofrac   = .239694E+00;     
0246 D.ClONO2_7646.weight    = 98.953723;
0247 D.ClONO2_7646.hitran    = 352;
0248 D.ClONO2_7646.jpl       = [ 99001 ];
0249 
0250 
0251          
0252 %---
0253 %--- CO
0254 %---
0255 %
0256 D.CO_26.isofrac   = .986544E+00;     
0257 D.CO_26.weight    = 27.994915;
0258 D.CO_26.hitran    = 51;
0259 D.CO_26.jpl       = [ 28001 ];
0260 %
0261 D.CO_36.isofrac   = 1.10836E-02;     
0262 D.CO_36.weight    = 28.998270;
0263 D.CO_36.hitran    = 52;
0264 D.CO_36.jpl       = [ 29001 ];
0265 %
0266 D.CO_28.isofrac   = 1.97822E-03;     
0267 D.CO_28.weight    = 29.999161;
0268 D.CO_28.hitran    = 53;
0269 D.CO_28.jpl       = [ 30001 ];
0270 %
0271 D.CO_27.isofrac   = 3.67867E-04;     
0272 D.CO_27.weight    = 28.999130;
0273 D.CO_27.hitran    = 54;
0274 D.CO_27.jpl       = [ 29006 ];
0275 %
0276 D.CO_38.isofrac   = 2.22250E-05;     
0277 D.CO_38.weight    = 31.002516;
0278 D.CO_38.hitran    = 55;
0279 D.CO_38.jpl       = [ NaN ];
0280 %
0281 D.CO_37.isofrac   = 4.13292E-06;     
0282 D.CO_37.weight    = 30.002485;
0283 D.CO_37.hitran    = 56;
0284 D.CO_37.jpl       = [ NaN ];
0285 
0286 
0287 
0288 %---
0289 %--- CO2
0290 %---
0291 %
0292 D.CO2_626.isofrac = 0.984204;     
0293 D.CO2_626.weight  = 43.989830;
0294 D.CO2_626.hitran  = 21;
0295 D.CO2_626.jpl     = [ NaN ];
0296 %
0297 D.CO2_636.isofrac = 1.10574E-02;     
0298 D.CO2_636.weight  = 44.993185;
0299 D.CO2_636.hitran  = 22;
0300 D.CO2_636.jpl     = [ NaN ];
0301 %
0302 D.CO2_628.isofrac = 3.94707E-03;     
0303 D.CO2_628.weight  = 45.994076;
0304 D.CO2_628.hitran  = 23;
0305 D.CO2_628.jpl     = [ 46013 ];
0306 %
0307 D.CO2_627.isofrac = 7.33989E-04;     
0308 D.CO2_627.weight  = 44.994045;
0309 D.CO2_627.hitran  = 24;
0310 D.CO2_627.jpl     = [ 45012 ];
0311 %
0312 D.CO2_638.isofrac = 4.43446E-05;     
0313 D.CO2_638.weight  = 46.997431;
0314 D.CO2_638.hitran  = 25;
0315 D.CO2_638.jpl     = [ NaN ];
0316 %
0317 D.CO2_637.isofrac = 8.24623E-06;     
0318 D.CO2_637.weight  = 45.997400;
0319 D.CO2_637.hitran  = 26;
0320 D.CO2_637.jpl     = [ NaN ];
0321 %
0322 D.CO2_828.isofrac = 3.95734E-06;     
0323 D.CO2_828.weight  = 47.998322;
0324 D.CO2_828.hitran  = 27;
0325 D.CO2_828.jpl     = [ NaN ];
0326 %
0327 D.CO2_728.isofrac = 1.47180E-06;     
0328 D.CO2_728.weight  = 46.998291;
0329 D.CO2_728.hitran  = 28;
0330 D.CO2_728.jpl     = [ NaN ];
0331 
0332 
0333 
0334 %---
0335 %--- COF2
0336 %---
0337 %
0338 D.COF2_269.isofrac = .986544E+00;     
0339 D.COF2_269.weight  = 65.991722;
0340 D.COF2_269.hitran  = 291;
0341 D.COF2_269.jpl     = [ 66001 ];
0342 
0343 
0344 
0345 %---
0346 %--- H2CO
0347 %---
0348 %
0349 D.H2CO_1126.isofrac = .986237E+00;     
0350 D.H2CO_1126.weight  = 30.010565;
0351 D.H2CO_1126.hitran  = 201;
0352 D.H2CO_1126.jpl     = [ 30004 ];
0353 %
0354 D.H2CO_1136.isofrac = 1.10802E-02;     
0355 D.H2CO_1136.weight  = 31.013920;
0356 D.H2CO_1136.hitran  = 202;
0357 D.H2CO_1136.jpl     = [ 31002 ];
0358 %
0359 D.H2CO_1128.isofrac = 1.97761E-03;     
0360 D.H2CO_1128.weight  = 32.014811;
0361 D.H2CO_1128.hitran  = 203;
0362 D.H2CO_1128.jpl     = [ 32004 ];
0363 %
0364 D.H2CO_1226.isofrac = 0.00029578940;     
0365 D.H2CO_1226.weight  = 31;
0366 D.H2CO_1226.hitran  = NaN;
0367 D.H2CO_1226.jpl     = [ 31003 ];
0368 %
0369 D.H2CO_2226.isofrac = 2.2181076E-08;     
0370 D.H2CO_2226.weight  = 32;
0371 D.H2CO_2226.hitran  = NaN;
0372 D.H2CO_2226.jpl     = [ 32006 ];
0373 
0374 
0375 
0376 %---
0377 %--- H2O
0378 %---
0379 %
0380 D.H2O_161.isofrac = 0.99731702;     
0381 D.H2O_161.weight  = 18.010565;
0382 D.H2O_161.hitran  = 11;
0383 D.H2O_161.jpl     = [ 18003, 18005 ];
0384 %
0385 D.H2O_181.isofrac = 0.00199983;     
0386 D.H2O_181.weight  = 20.014811;
0387 D.H2O_181.hitran  = 12;
0388 D.H2O_181.jpl     = [ 20003 ];
0389 %
0390 D.H2O_171.isofrac = 0.00037200;     
0391 D.H2O_171.weight  = 19.014780;
0392 D.H2O_171.hitran  = 13;
0393 D.H2O_171.jpl     = [ 19003 ];
0394 %
0395 D.H2O_162.isofrac = 0.000310693; 
0396 D.H2O_162.weight  = 19.016740;     
0397 D.H2O_162.hitran  = 14;
0398 D.H2O_162.jpl     = [ 19002 ];
0399 %
0400 D.H2O_182.isofrac = 6.23003E-07; 
0401 D.H2O_182.weight  = 21.020985;     
0402 D.H2O_182.hitran  = 15;
0403 D.H2O_182.jpl     = [ 21001 ];
0404 %
0405 D.H2O_172.isofrac = 1.15853E-07;     
0406 D.H2O_172.weight  = 20.020956;
0407 D.H2O_172.hitran  = 16;
0408 D.H2O_172.jpl     = [ NaN ];
0409 %
0410 D.H2O_262.isofrac = 2.2430204E-08;     
0411 D.H2O_262.weight  = 20;
0412 D.H2O_262.hitran  = NaN;
0413 D.H2O_262.jpl     = [ 20001 ];
0414 
0415 
0416 
0417 %---
0418 %--- H2O2
0419 %---
0420 %
0421 D.H2O2_1661.isofrac = .994952E+00;
0422 D.H2O2_1661.weight  = 34.005480;
0423 D.H2O2_1661.hitran  = 251;
0424 D.H2O2_1661.jpl     = [ 34004 ];
0425 
0426 
0427 
0428 %---
0429 %--- H2S
0430 %---
0431 %
0432 D.H2S_121.isofrac = .949884E+00;
0433 D.H2S_121.weight  = 33.987721;
0434 D.H2S_121.hitran  = 311;
0435 D.H2S_121.jpl     = [ 34002 ];
0436 %
0437 D.H2S_141.isofrac = 4.21369E-02;
0438 D.H2S_141.weight  = 35.983515;
0439 D.H2S_141.hitran  = 312;
0440 D.H2S_141.jpl     = [ NaN ];
0441 %
0442 D.H2S_131.isofrac = 7.49766E-03;
0443 D.H2S_131.weight  = 34.987105;
0444 D.H2S_131.hitran  = 313;
0445 D.H2S_131.jpl     = [ NaN ];
0446 %
0447 D.H2S_122.isofrac = 0.00029991625;
0448 D.H2S_122.weight  = 35;
0449 D.H2S_122.hitran  = NaN;
0450 D.H2S_122.jpl     = [ 35001 ];
0451 
0452 
0453 
0454 %---
0455 %--- H2SO4
0456 %---
0457 %
0458 D.H2SO4_126.isofrac = 0.95060479;
0459 D.H2SO4_126.weight  = 98;
0460 D.H2SO4_126.hitran  = NaN;
0461 D.H2SO4_126.jpl     = [ 98001 ];
0462 
0463 
0464 
0465 %---
0466 %--- HBr
0467 %---
0468 %
0469 D.HBr_19.isofrac = 0.506781;
0470 D.HBr_19.weight  = 79.926160;
0471 D.HBr_19.hitran  = 161;
0472 D.HBr_19.jpl     = [ 80001 ];
0473 %
0474 D.HBr_11.isofrac = 0.493063;
0475 D.HBr_11.weight  = 81.924115;
0476 D.HBr_11.hitran  = 162;
0477 D.HBr_11.jpl     = [ 82001 ];
0478 
0479 
0480 
0481 %---
0482 %--- HCl
0483 %---
0484 %
0485 D.HCl_15.isofrac = 0.757587;
0486 D.HCl_15.weight  = 35.976678;
0487 D.HCl_15.hitran  = 151;
0488 D.HCl_15.jpl     = [ 36001 ];
0489 %
0490 D.HCl_17.isofrac = 0.242257;
0491 D.HCl_17.weight  = 37.973729;
0492 D.HCl_17.hitran  = 152;
0493 D.HCl_17.jpl     = [ 38001 ];
0494 %
0495 D.HCl_25.isofrac = 0.00011324004;
0496 D.HCl_25.weight  = 37;
0497 D.HCl_25.hitran  = NaN;
0498 D.HCl_25.jpl     = [ 37001 ];
0499 %
0500 D.HCl_27.isofrac = 3.6728230E-05;
0501 D.HCl_27.weight  = 39;
0502 D.HCl_27.hitran  = NaN;
0503 D.HCl_27.jpl     = [ 39004 ];
0504 
0505 
0506 
0507 %---
0508 %--- HCN
0509 %---
0510 %
0511 D.HCN_124.isofrac = .985114E+00;
0512 D.HCN_124.weight  = 27.010899;
0513 D.HCN_124.hitran  = 231;
0514 D.HCN_124.jpl     = [ 27001, 27003 ];
0515 %
0516 D.HCN_134.isofrac = 1.10676E-02;
0517 D.HCN_134.weight  = 28.014254;
0518 D.HCN_134.hitran  = 232;
0519 D.HCN_134.jpl     = [ 28002 ];
0520 %
0521 D.HCN_125.isofrac = 3.62174E-03;
0522 D.HCN_125.weight  = 28.007933;
0523 D.HCN_125.hitran  = 233;
0524 D.HCN_125.jpl     = [ 28003 ];
0525 %
0526 D.HCN_224.isofrac = 0.00014773545;
0527 D.HCN_224.weight  = 28;
0528 D.HCN_224.hitran  = NaN;
0529 D.HCN_224.jpl     = [ 28004 ];
0530 
0531 
0532 
0533 %---
0534 %--- HCOOH
0535 %---
0536 %
0537 D.HCOOH_1261.isofrac = 0.983898E+00;
0538 D.HCOOH_1261.weight  = 46.005480;
0539 D.HCOOH_1261.hitran  = 321;
0540 D.HCOOH_1261.jpl     = [ 46005 ];
0541 %
0542 D.HCOOH_1361.isofrac = 0.010913149;
0543 D.HCOOH_1361.weight  = 47;
0544 D.HCOOH_1361.hitran  = NaN;
0545 D.HCOOH_1361.jpl     = [ 47002 ];
0546 %
0547 D.HCOOH_2261.isofrac = 0.00014755369;
0548 D.HCOOH_2261.weight  = 47;
0549 D.HCOOH_2261.hitran  = NaN;
0550 D.HCOOH_2261.jpl     = [ 47003 ];
0551 %
0552 D.HCOOH_1262.isofrac = 0.00014755369;
0553 D.HCOOH_1262.weight  = 47;
0554 D.HCOOH_1262.hitran  = NaN;
0555 D.HCOOH_1262.jpl     = [ 47004 ];
0556 
0557 
0558 
0559 %---
0560 %--- HF
0561 %---
0562 %
0563 D.HF_19.isofrac = 0.99984425;
0564 D.HF_19.weight  = 20.006229;
0565 D.HF_19.hitran  = 141;
0566 D.HF_19.jpl     = [ 20001 ];
0567 %
0568 D.HF_29.isofrac = 0.00014994513;
0569 D.HF_29.weight  = 21;
0570 D.HF_29.hitran  = NaN;
0571 D.HF_29.jpl     = [ 21001 ];
0572 
0573 
0574 
0575 %---
0576 %--- HI
0577 %---
0578 %
0579 D.HI_17.isofrac = 0.99984425;
0580 D.HI_17.weight  = 127.912297;
0581 D.HI_17.hitran  = 171;
0582 D.HI_17.jpl     = [ NaN ];
0583 
0584 
0585 
0586 %---
0587 %--- HNC
0588 %---
0589 %
0590 D.HNC_142.isofrac = 0.98505998;
0591 D.HNC_142.weight  = 27;
0592 D.HNC_142.hitran  = NaN;
0593 D.HNC_142.jpl     = [ 27002 ];
0594 %
0595 D.HNC_143.isofrac = 0.011091748;
0596 D.HNC_143.weight  = 28;
0597 D.HNC_143.hitran  = NaN;
0598 D.HNC_143.jpl     = [ 28005 ];
0599 %
0600 D.HNC_152.isofrac = 0.0036982817;
0601 D.HNC_152.weight  = 28;
0602 D.HNC_152.hitran  = NaN;
0603 D.HNC_152.jpl     = [ 28006 ];
0604 %
0605 D.HNC_242.isofrac = 0.00014996849;
0606 D.HNC_242.weight  = 28;
0607 D.HNC_242.hitran  = NaN;
0608 D.HNC_242.jpl     = [ 28007 ];
0609 
0610 
0611 
0612 %---
0613 %--- HNO3
0614 %---
0615 %
0616 D.HNO3_146.isofrac = 0.989110;
0617 D.HNO3_146.weight  = 62.995644;
0618 D.HNO3_146.hitran  = 121;
0619 D.HNO3_146.jpl     = [ 63001, 63002, 63003, 63004, 63005, 63006 ];
0620 
0621 
0622 
0623 %---
0624 %--- HO2
0625 %---
0626 %
0627 D.HO2_166.isofrac = 0.995107;
0628 D.HO2_166.weight  = 32.997655;
0629 D.HO2_166.hitran  = 331;
0630 D.HO2_166.jpl     = [ 33001 ];
0631 
0632 
0633 
0634 %---
0635 %--- HOBr
0636 %---
0637 %
0638 D.HOBr_169.isofrac = .505579E+00;
0639 D.HOBr_169.weight  = 95.921076;
0640 D.HOBr_169.hitran  = 371;
0641 D.HOBr_169.jpl     = [ 96001 ];
0642 %
0643 D.HOBr_169.isofrac = .491894E+00;
0644 D.HOBr_169.weight  = 97.919027;
0645 D.HOBr_169.hitran  = 372;
0646 D.HOBr_169.jpl     = [ 98001 ];
0647 
0648 
0649 
0650 %---
0651 %--- HOCl
0652 %---
0653 %
0654 D.HOCl_165.isofrac = .755790E+00;
0655 D.HOCl_165.weight  = 51.971593;
0656 D.HOCl_165.hitran  = 211;
0657 D.HOCl_165.jpl     = [ 52006 ];
0658 %
0659 D.HOCl_167.isofrac = .241683E+00;
0660 D.HOCl_167.weight  = 53.968644;
0661 D.HOCl_167.hitran  = 212;
0662 D.HOCl_167.jpl     = [ 54005 ];
0663 
0664 
0665 
0666 %---
0667 %--- N2
0668 %---
0669 %
0670 D.N2_44.isofrac = 0.9926874;
0671 D.N2_44.weight  = 28.006147;
0672 D.N2_44.hitran  = 221;
0673 D.N2_44.jpl     = [ NaN ];
0674 
0675 
0676 
0677 %---
0678 %--- N2O
0679 %---
0680 %
0681 D.N2O_446.isofrac = .990333E+00;
0682 D.N2O_446.weight  = 44.001062;
0683 D.N2O_446.hitran  = 41;
0684 D.N2O_446.jpl     = [ 44004, 44009, 44012 ];
0685 %
0686 D.N2O_456.isofrac = 3.64093E-03;
0687 D.N2O_456.weight  = 44.998096;
0688 D.N2O_456.hitran  = 42;
0689 D.N2O_456.jpl     = [ 45007 ];
0690 %
0691 D.N2O_546.isofrac = 3.64093E-03;
0692 D.N2O_546.weight  = 44.998096;
0693 D.N2O_546.hitran  = 43;
0694 D.N2O_546.jpl     = [ 45008 ];
0695 %
0696 D.N2O_448.isofrac = 1.98582E-03;
0697 D.N2O_448.weight  = 46.005308;
0698 D.N2O_448.hitran  = 44;
0699 D.N2O_448.jpl     = [ 46007 ];
0700 %
0701 D.N2O_447.isofrac = 3.69280E-04;
0702 D.N2O_447.weight  = 45.005278;
0703 D.N2O_447.hitran  = 45;
0704 D.N2O_447.jpl     = [ NaN ];
0705 
0706 
0707 
0708 %---
0709 %--- NH3
0710 %---
0711 %
0712 D.NH3_4111.isofrac   = .9958715E+00;
0713 D.NH3_4111.weight    = 17.026549;
0714 D.NH3_4111.hitran    = 111;
0715 D.NH3_4111.jpl       = [ 17002, 17004 ];
0716 %
0717 D.NH3_5111.isofrac   = 3.66129E-03;
0718 D.NH3_5111.weight    = 18.023583;
0719 D.NH3_5111.hitran    = 112;
0720 D.NH3_5111.jpl       = [ 18002 ];
0721 %
0722 D.NH3_4112.isofrac   = 0.00044792294;
0723 D.NH3_4112.weight    = 18;
0724 D.NH3_4112.hitran    = NaN;
0725 D.NH3_4112.jpl       = [ 18004 ];
0726 
0727 
0728 
0729 %---
0730 %--- NO
0731 %---
0732 %
0733 D.NO_46.isofrac   = .993974E+00;
0734 D.NO_46.weight    = 29.997989;
0735 D.NO_46.hitran    = 81;
0736 D.NO_46.jpl       = [ 30008 ];
0737 %
0738 D.NO_56.isofrac   = 3.65431E-03;
0739 D.NO_56.weight    = 30.995023;
0740 D.NO_56.hitran    = 82;
0741 D.NO_56.jpl       = [ NaN ];
0742 %
0743 D.NO_48.isofrac   = 1.99312E-03;
0744 D.NO_48.weight    = 32.002234;
0745 D.NO_48.hitran    = 83;
0746 D.NO_48.jpl       = [ NaN ];
0747 
0748 
0749 
0750 %---
0751 %--- NO+
0752 %---
0753 %
0754 D.NOx_46.isofrac   = .993974E+00;
0755 D.NOx_46.weight    = 29.997989;
0756 D.NOx_46.hitran    = 361;
0757 D.NOx_46.jpl       = [ 30011 ];
0758 
0759          
0760 %---
0761 %--- NO2
0762 %---
0763 %
0764 D.NO2_646.isofrac   = .991616E+00;
0765 D.NO2_646.weight    = 45.992904;
0766 D.NO2_646.hitran    = 101;
0767 D.NO2_646.jpl       = [ 46006 ];
0768 
0769 
0770 
0771 %---
0772 %--- O
0773 %---
0774 %
0775 D.O_6.isofrac   = 0.997628;
0776 D.O_6.weight    = 15.994915;
0777 D.O_6.hitran    = 341;
0778 D.O_6.jpl       = [ 16001 ];
0779 
0780 
0781 
0782 %---
0783 %--- O2
0784 %---
0785 %
0786 D.O2_66.isofrac   = .995262E+00;
0787 D.O2_66.weight    = 31.989830;
0788 D.O2_66.hitran    = 71;
0789 D.O2_66.jpl       = [ 32001, 32002 ];
0790 %
0791 D.O2_68.isofrac   = 3.99141E-03;
0792 D.O2_68.weight    = 33.994076;
0793 D.O2_68.hitran    = 72;
0794 D.O2_68.jpl       = [ 34001 ];
0795 %
0796 D.O2_67.isofrac   = 7.42235E-04;
0797 D.O2_67.weight    = 32.994045;
0798 D.O2_67.hitran    = 73;
0799 D.O2_67.jpl       = [ 330002 ];
0800 
0801 
0802 
0803 %---
0804 %--- O3
0805 %---
0806 %
0807 D.O3_666.isofrac = .992901E+00;     
0808 D.O3_666.weight  = 47.984745;
0809 D.O3_666.hitran  = 31;
0810 D.O3_666.jpl     = [ 48004, 48005, 48006, 48007, 48008 ];
0811 %
0812 D.O3_668.isofrac = 3.98194E-03;     
0813 D.O3_668.weight  = 49.988991;
0814 D.O3_668.hitran  = 32;
0815 D.O3_668.jpl     = [ 50004, 50006 ];
0816 %
0817 D.O3_686.isofrac = 1.99097E-03;     
0818 D.O3_686.weight  = 49.988991;
0819 D.O3_686.hitran  = 33;
0820 D.O3_686.jpl     = [ 50003, 50005 ];
0821 %
0822 D.O3_667.isofrac = 7.40475E-04;     
0823 D.O3_667.weight  = 48.988960;
0824 D.O3_667.hitran  = 34;
0825 D.O3_667.jpl     = [ 49002 ];
0826 %
0827 D.O3_676.isofrac = 3.70237E-04;     
0828 D.O3_676.weight  = 48.988960;
0829 D.O3_676.hitran  = 35;
0830 D.O3_676.jpl     = [ 49001 ];
0831 
0832 
0833 
0834 %---
0835 %--- OClO
0836 %---
0837 %
0838 D.OClO_656.isofrac   = 0.75509223;
0839 D.OClO_656.weight    = 67;
0840 D.OClO_656.hitran    = NaN;
0841 D.OClO_656.jpl       = [ 67001 ];
0842 %
0843 D.OClO_656.isofrac   = 0.24490632;
0844 D.OClO_656.weight    = 69;
0845 D.OClO_656.hitran    = NaN;
0846 D.OClO_656.jpl       = [ 69001 ];
0847 
0848 
0849 
0850 %---
0851 %--- OCS
0852 %---
0853 %
0854 % OCS-623 introduced in Hitran 2000, with an iso. number already used for 822
0855 %
0856 D.OCS_622.isofrac   = .937395E+00;
0857 D.OCS_622.weight    = 59.966986;
0858 D.OCS_622.hitran    = 191;
0859 D.OCS_622.jpl       = [ 60001 ];
0860 %
0861 D.OCS_624.isofrac   = 4.15828E-02;
0862 D.OCS_624.weight    = 61.962780;
0863 D.OCS_624.hitran    = 192;
0864 D.OCS_624.jpl       = [ 62001 ];
0865 %
0866 D.OCS_632.isofrac   = 1.05315E-02;
0867 D.OCS_632.weight    = 60.970341;
0868 D.OCS_632.hitran    = 193;
0869 D.OCS_632.jpl       = [ 61001 ];
0870 %
0871 D.OCS_623.isofrac   = 7.39908E-03;
0872 D.OCS_623.weight    = 60.966371;
0873 D.OCS_623.hitran    = 194;
0874 D.OCS_623.jpl       = [ NaN ];
0875 %
0876 D.OCS_822.isofrac   = 1.87967E-03;
0877 D.OCS_822.weight    = 61.971231;
0878 D.OCS_822.hitran    = 195;
0879 D.OCS_822.jpl       = [ 62002 ];
0880 
0881 
0882 
0883 %---
0884 %--- OH
0885 %---
0886 %
0887 D.OH_61.isofrac   = .997473E+00;
0888 D.OH_61.weight    = 17.002740;
0889 D.OH_61.hitran    = 131;
0890 D.OH_61.jpl       = [ 17001 ];
0891 %
0892 D.OH_81.isofrac   = 2.00014E-03;
0893 D.OH_81.weight    = 19.006986;
0894 D.OH_81.hitran    = 132;
0895 D.OH_81.jpl       = [ 19001 ];
0896 %
0897 D.OH_62.isofrac   = 1.55371E-04;
0898 D.OH_62.weight    = 18.008915;
0899 D.OH_62.hitran    = 133;
0900 D.OH_62.jpl       = [ 18001 ];
0901 
0902 
0903 
0904 %---
0905 %--- PH3
0906 %---
0907 %
0908 D.PH3_1111.isofrac   = 0.99953283;
0909 D.PH3_1111.weight    = 33.997238;
0910 D.PH3_1111.hitran    = 281;
0911 D.PH3_1111.jpl       = [ 34003 ];
0912 
0913 
0914 
0915 %---
0916 %--- SF6
0917 %---
0918 %
0919 D.SF6_29.isofrac   = 0.95018;
0920 D.SF6_29.weight    = 145.962492;
0921 D.SF6_29.hitran    = 301;
0922 D.SF6_29.jpl       = [ NaN ];
0923 
0924 
0925 
0926 %---
0927 %--- SO2
0928 %---
0929 %
0930 D.SO2_626.isofrac = .945678E+00;     
0931 D.SO2_626.weight  = 63.961901;
0932 D.SO2_626.hitran  = 91;
0933 D.SO2_626.jpl     = [ 64002, 64005 ];
0934 %
0935 D.SO2_646.isofrac = 4.19503E-02;     
0936 D.SO2_646.weight  = 65.957695;
0937 D.SO2_646.hitran  = 92;
0938 D.SO2_646.jpl     = [ 66002 ];
0939 %
0940 D.SO2_636.isofrac = 0.0074989421;     
0941 D.SO2_636.weight  = 65;
0942 D.SO2_636.hitran  = 93;
0943 D.SO2_636.jpl     = [ 65001 ];
0944 %
0945 D.SO2_628.isofrac = 0.0020417379;     
0946 D.SO2_628.weight  = 66;
0947 D.SO2_628.hitran  = 94;
0948 D.SO2_628.jpl     = [ 66004 ];
0949 
0950 
0951 
0952 
0953 %-----------------------------------------------------------------------------
0954 %
0955 % Extarct selected data
0956 %
0957 %-----------------------------------------------------------------------------
0958 
0959 
0960 data = [];
0961 
0962 
0963 if nargin == 0
0964   %
0965   f = fieldnames( D );
0966   %
0967   data(1).isofrac = NaN;
0968   data(1).weight  = NaN;
0969   data(1).hitran  = NaN;
0970   data(1).jpl     = NaN;
0971   data(1).name    = NaN;
0972   %
0973   for i = 1:length(f)
0974     s             = f{i};
0975     ind           = find( s == '_' );
0976     s(ind)        = '-';
0977     ind           = find( s == 'x' );
0978     s(ind)        = '+';
0979     D.(f{i}).name = s;
0980     data(i)       = D.(f{i});
0981   end
0982   return
0983 end
0984 
0985 
0986 for i = 1 : length( isos )
0987 
0988   s      = isos{i};
0989   ind    = find( s == '-' );
0990   s(ind) = '_';
0991   ind    = find( s == '+' );
0992   s(ind) = 'x';
0993   
0994   if ~isfield( D, s )
0995     error( sprintf('No data defined for %s.',isos{i}) );
0996   end
0997     
0998   d = getfield( D, s );
0999   
1000   for f = 1 : length( fields )
1001 
1002     if ~isfield( d, fields{f} )
1003       error( sprintf( 'Field *%s* not defined for %s.', fields{f}, isos{i} ) );
1004     end
1005     
1006     data(i).(fields{f}) = getfield( d, fields{f} );
1007     
1008   end
1009 end
1010 
1011 
1012 
1013 
1014 
1015 
1016

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