Home > atmlab > scattering > box_pnd_mono_size_1d.m

box_pnd_mono_size_1d

PURPOSE ^

-----------------------------------------------------------------------------

SYNOPSIS ^

function pnd_field = box_pnd_mono_size_1d(z_limits, imc, r)

DESCRIPTION ^

-----------------------------------------------------------------------------
  PND_MONO_SIZE_1D  Generates simple particle number density field
 
  This function generates the particle number density field for a 
  very simple 1D cirrus cloud layer. It assumes that all particles are of 
  the same size.  As input it requires the altitude of the cloud
  layer, the ice mass content and the radius of the
  particles.
 
  FORMAT pnd_field = box_pnd_mono_size_1d(z_limits, imc, r)
 
  OUT    pnd        particle number density (ARTS data format GField3) 
  IN     z_limits   cloud altitude     [z1 z2]  [m]
         imc        ice mass content   [kg/m^3]
         r          radius             [m]
 
  Created by Claudia Emde

-----------------------------------------------------------------------------

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

box_pnd_mono_size_1d.m

SOURCE CODE ^

0001 %-----------------------------------------------------------------------------
0002 %  PND_MONO_SIZE_1D  Generates simple particle number density field
0003 %
0004 %  This function generates the particle number density field for a
0005 %  very simple 1D cirrus cloud layer. It assumes that all particles are of
0006 %  the same size.  As input it requires the altitude of the cloud
0007 %  layer, the ice mass content and the radius of the
0008 %  particles.
0009 %
0010 %  FORMAT pnd_field = box_pnd_mono_size_1d(z_limits, imc, r)
0011 %
0012 %  OUT    pnd        particle number density (ARTS data format GField3)
0013 %  IN     z_limits   cloud altitude     [z1 z2]  [m]
0014 %         imc        ice mass content   [kg/m^3]
0015 %         r          radius             [m]
0016 %
0017 %  Created by Claudia Emde
0018 %
0019 %-----------------------------------------------------------------------------
0020 
0021 function pnd_field = box_pnd_mono_size_1d(z_limits, imc, r)
0022 
0023 % Convert altitude to pressure
0024 p_limits = z2p_simple(z_limits);
0025 
0026 % Density of ice 917 kg/m^3
0027 rho = 917;
0028 
0029 % Calculation of particle number density for given IMC
0030 % and radius
0031 pnd = 3/(4*pi) * imc/(rho*r^3);
0032 
0033 % Define pressure grid. Two points are sufficient for box cloud
0034 % We need to pressures nearby the limits, to assure that the
0035 % particle number density decreases steaply to, take 1Pa to both
0036 % sides.
0037 p_grid = [1000e2 p_limits(1)+10 p_limits(1) p_limits(2) p_limits(2)-10 1];
0038 % Assign pnd values
0039 data(:,1,1) = [0 0 pnd pnd 0 0];
0040 
0041 % Create GField3
0042 %
0043 pnd_field.grids     = { p_grid, [0], [0] };
0044 pnd_field.gridnames = { 'Pressure', 'Latitude', 'Longitude' };
0045 %
0046 pnd_field.data  = data;
0047

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