glen_coefficients

PURPOSE ^

GLEN_COEFFICIENTS - compute A and B, given the temperature and Ao

SYNOPSIS ^

function [A ,B]=glen_coefficients(Ao,T);

DESCRIPTION ^

GLEN_COEFFICIENTS - compute A and B, given the temperature and Ao

   Ao,T according to patterson,Ao in s-1kPa-3.
   values range between 3.9 and 5.3*10^-16;
   Return: A in s-1kPa-3. B in Pa.a1/3

   Usage:
      [A ,B]=glen_coefficients(Ao,T)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function [A ,B]=glen_coefficients(Ao,T);
0002 %GLEN_COEFFICIENTS - compute A and B, given the temperature and Ao
0003 %
0004 %   Ao,T according to patterson,Ao in s-1kPa-3.
0005 %   values range between 3.9 and 5.3*10^-16;
0006 %   Return: A in s-1kPa-3. B in Pa.a1/3
0007 %
0008 %   Usage:
0009 %      [A ,B]=glen_coefficients(Ao,T)
0010 
0011 To=263;
0012 R=8.314;
0013 E1=139000;
0014 E2=60000;
0015 
0016 %For T>263,
0017 Ac1=Ao*exp(E1/R/To);
0018 %For T<263,
0019 Ac2=Ao*exp(E2/R/To);
0020 
0021 A=Ac1*exp(-E1/R./T);
0022 pos=find(T<263);
0023 A(pos)=Ac2*exp(-E2/R./T(pos));
0024 
0025 B=(A*365*24*3600).^(-1/3);
0026

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003