Home > atmlab > scattering > Asymmetry_parameter.m

Asymmetry_parameter

PURPOSE ^

Asymmetry_parameter extracts the asymmetry parameter from an array of

SYNOPSIS ^

function g = Asymmetry_parameter(D)

DESCRIPTION ^

 Asymmetry_parameter   extracts the asymmetry parameter from an array of
                       ARTS single scattering data

        Returns an array of asymmetry parameters depending on the length 
        of T-matrix's output. 
        The output of theT-matrix (D), is the input of the Asymmetry_parameter function.
     


 FORMAT   g = Asymmetry_parameter(D)
  

 OUT      g        Array of asymmetry parameters         [-] 
                   Its dimension is: g{length of D}(length of f_grid , length of t_grid)

 IN       D        Array of ARTS single scattering data(Tmatrix's output) 

 History: 2014-02-10      Created by Maryam Jamali

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

Asymmetry_parameter.m

SOURCE CODE ^

0001 % Asymmetry_parameter   extracts the asymmetry parameter from an array of
0002 %                       ARTS single scattering data
0003 %
0004 %        Returns an array of asymmetry parameters depending on the length
0005 %        of T-matrix's output.
0006 %        The output of theT-matrix (D), is the input of the Asymmetry_parameter function.
0007 %
0008 %
0009 %
0010 % FORMAT   g = Asymmetry_parameter(D)
0011 %
0012 %
0013 % OUT      g        Array of asymmetry parameters         [-]
0014 %                   Its dimension is: g{length of D}(length of f_grid , length of t_grid)
0015 %
0016 % IN       D        Array of ARTS single scattering data(Tmatrix's output)
0017 %
0018 % History: 2014-02-10      Created by Maryam Jamali
0019 
0020 
0021 
0022 function g = Asymmetry_parameter(D)
0023 
0024 for i = 1 : length(D)
0025     for j = 1 : length (D{1,i}.f_grid) 
0026         for k = 1 : length (D{1,i}.T_grid) 
0027  
0028           b  =  D{1,i}.pha_mat_data(j,k,1:end,1,1,1,1) ;
0029           norm_factor = 1/2 * trapz ( D{1,i}.za_grid , squeeze(b) .* sind(D{1,i}.za_grid) );
0030  
0031           g{i}(j,k) = 1/2 * trapz ( D{1,i}.za_grid , cosd(D{1,i}.za_grid) .* (squeeze(b) ./ norm_factor) .* sind(D{1,i}.za_grid) );
0032  
0033         end
0034     end
0035 end

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