Home > atmlab > h2o > parametrisations > n_ice_warren84.m

n_ice_warren84

PURPOSE ^

N_ICE_WARREN84 Complex refractive index for ice

SYNOPSIS ^

function n = n_ice_warren84( f, t )

DESCRIPTION ^

 N_ICE_WARREN84   Complex refractive index for ice

   Refractive index of ice following the Warren84 parameterization.

   Calculates complex refractive index of Ice 1H for wavelengths
   between 45 nm and 8.6 m.
   For wavelengths above 167 microns, temperature dependence is
   included for temperatures between 213 and 272K.
   
   This function requires that ARTS is at hand. For more information about
   the ARTS workspace method that is used by this Atmlab function:
      arts -d complex_refr_indexIceWarren84

 FORMAT   n = n_ice_warren84( f, t )
        
 OUT   n   Complex refractive index, for each combination of f and t.
 IN    f   Frequency
       t   Temperature

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

n_ice_warren84.m

SOURCE CODE ^

0001 % N_ICE_WARREN84   Complex refractive index for ice
0002 %
0003 %   Refractive index of ice following the Warren84 parameterization.
0004 %
0005 %   Calculates complex refractive index of Ice 1H for wavelengths
0006 %   between 45 nm and 8.6 m.
0007 %   For wavelengths above 167 microns, temperature dependence is
0008 %   included for temperatures between 213 and 272K.
0009 %
0010 %   This function requires that ARTS is at hand. For more information about
0011 %   the ARTS workspace method that is used by this Atmlab function:
0012 %      arts -d complex_refr_indexIceWarren84
0013 %
0014 % FORMAT   n = n_ice_warren84( f, t )
0015 %
0016 % OUT   n   Complex refractive index, for each combination of f and t.
0017 % IN    f   Frequency
0018 %       t   Temperature
0019 
0020 % 2014-05-27   Created by Patrick Eriksson.
0021 
0022 
0023 function n = n_ice_warren84( f, t )
0024 %
0025 rqre_datatype( f, @istensor1 );
0026 rqre_datatype( t, @istensor1 );
0027 
0028 
0029 %- Create workfolder?
0030 %
0031 workfolder = create_tmpfolder;
0032 cu = onCleanup( @()delete_tmpfolder( workfolder ) );
0033 
0034 
0035 %- Store input arguments to files
0036 %
0037 xmlStore( fullfile(workfolder,'f.xml'), f, 'Vector' );
0038 xmlStore( fullfile(workfolder,'t.xml'), t, 'Vector' );
0039 
0040 
0041 %- Define control file
0042 %
0043 S{1}     = 'Arts2{';
0044 S{end+1} = 'VectorCreate( f )';
0045 S{end+1} = 'VectorCreate( t )';
0046 S{end+1} = 'ReadXML( f, "f.xml" )';
0047 S{end+1} = 'ReadXML( t, "t.xml" )';
0048 S{end+1} = 'complex_refr_indexIceWarren84( complex_refr_index, f, t )';
0049 S{end+1} = 'WriteXML( "ascii", complex_refr_index, "n.xml" )';
0050 S{end+1} = '}';
0051 
0052 
0053 %- Run ARTS
0054 %
0055 cfile = fullfile( workfolder, 'cfile.arts' );
0056 strs2file( cfile, S );
0057 arts( cfile );
0058 
0059 
0060 %- Load data
0061 %
0062 N = xmlLoad( fullfile( workfolder, 'n.xml' ) );
0063 n = N.data(:,:,1) + i*N.data(:,:,2);

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