Home > atmlab > h2o > parametrisations > ice_particleshape_Hogan2012.m

ice_particleshape_Hogan2012

PURPOSE ^

ice_particleshape_Hogan2012 Dimension parameters(long and short diameters)

SYNOPSIS ^

function [diameter_max diameter_short aspect_ratio mixfrac rho]= ice_particleshape_Hogan2012(d)

DESCRIPTION ^

 ice_particleshape_Hogan2012   Dimension parameters(long and short diameters)
                               and density of a non-spherical and non-solid 
                               ice particle that composes ice matrix with air inclusion.
                         

             The equivalent mass of a spherical solid ice is calculated upon the 
             distribution of mass equivalent spheres m=(power(d,3)*pi*rhoice)/6;  
             then according to Brown and Francis(1995)relationship between
             particle mass and size, the mean and max(long) diameters, 
             and inclusion media fraction of a non-spherical(spheroidal) 
             particle which has the *same mass*, are computed.
             
             Note that all of the outputs are as a function of mass-equivalent 
             diameter (d).
             
             The parameterization is taken from R. Hogan et al(2012).
             "Radar Scattering from Ice Aggregates Using the Horizontally 
             Aligned Oblate Spheroid Approximation".
             
                        
 
 FORMAT         [diameter_max diameter_short aspect_ratio mixfrac rho]= ice_particleshape_Hogan2012(d)
        
 OUT  diameter_max   longest diameter of an aligned oblate spheroid particle  [m]      
      diameter_short shortest diameter of an aligned oblate spheroid particle [m]  
      aspect_ratio   d_short / d_long
      mixfrac        Fraction of inclusion media (air) in ice matrix.
      rho            Density of a sheroid of non-solid ice particle      [kg/m^3]
                     (mixture of ice and air)               

 IN   d              mass equivalent diameter       [m]

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

ice_particleshape_Hogan2012.m

SOURCE CODE ^

0001 % ice_particleshape_Hogan2012   Dimension parameters(long and short diameters)
0002 %                               and density of a non-spherical and non-solid
0003 %                               ice particle that composes ice matrix with air inclusion.
0004 %
0005 %
0006 %             The equivalent mass of a spherical solid ice is calculated upon the
0007 %             distribution of mass equivalent spheres m=(power(d,3)*pi*rhoice)/6;
0008 %             then according to Brown and Francis(1995)relationship between
0009 %             particle mass and size, the mean and max(long) diameters,
0010 %             and inclusion media fraction of a non-spherical(spheroidal)
0011 %             particle which has the *same mass*, are computed.
0012 %
0013 %             Note that all of the outputs are as a function of mass-equivalent
0014 %             diameter (d).
0015 %
0016 %             The parameterization is taken from R. Hogan et al(2012).
0017 %             "Radar Scattering from Ice Aggregates Using the Horizontally
0018 %             Aligned Oblate Spheroid Approximation".
0019 %
0020 %
0021 %
0022 % FORMAT         [diameter_max diameter_short aspect_ratio mixfrac rho]= ice_particleshape_Hogan2012(d)
0023 %
0024 % OUT  diameter_max   longest diameter of an aligned oblate spheroid particle  [m]
0025 %      diameter_short shortest diameter of an aligned oblate spheroid particle [m]
0026 %      aspect_ratio   d_short / d_long
0027 %      mixfrac        Fraction of inclusion media (air) in ice matrix.
0028 %      rho            Density of a sheroid of non-solid ice particle      [kg/m^3]
0029 %                     (mixture of ice and air)
0030 %
0031 % IN   d              mass equivalent diameter       [m]
0032        
0033 % 2013-08-09    Created by Maryam Jamali
0034 
0035 
0036 function  [diameter_max diameter_short aspect_ratio mixfrac rho]= ice_particleshape_Hogan2012(d)
0037 
0038 
0039 rhoice=0.917*1e3;     % kg/m^3
0040 
0041 %To calculate diameter_max and diameter_short :
0042 for i=1:length(d)
0043     m=(power(d(i),3)*pi*rhoice)/6;   %kg
0044     D_0=power((m./(pi*rhoice/6)),(1/3));
0045     
0046     if D_0 < 97e-6
0047     diameter_mean(i)=D_0;
0048     else
0049     diameter_mean(i)=power((m./0.0185), (1/1.9));
0050     end
0051 
0052     if  D_0 < 66e-6
0053     diameter_max(i)=D_0;
0054     else
0055     diameter_max(i)=power((m./0.0121), (1/1.9));
0056     end
0057     
0058     diameter_long(i) =diameter_max(i);
0059     diameter_short(i)=(2*diameter_mean(i))-diameter_max(i);
0060     aspect_ratio(i)=diameter_short(i)./diameter_long(i);
0061     
0062     volume_ice(i)= (pi/6).* power(d(i),3);
0063     
0064     if diameter_short(i)==diameter_long(i)
0065     rho(i)=rhoice;
0066     volume_spheroid(i)=volume_ice(i);
0067     else
0068     volume_spheroid(i)=(pi/6).*power(diameter_long(i),2).*diameter_short(i);
0069     rho(i)=m./volume_spheroid(i); % kg/m^3
0070     end
0071 end
0072 
0073 mixfrac=(volume_spheroid - volume_ice)./ volume_spheroid;

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