Home > atmlab > mie > miecoated_rain3.m

miecoated_rain3

PURPOSE ^

Melting rain

SYNOPSIS ^

function result = miecoated_rain3(fGHz, TK, nrain, pam, coat)

DESCRIPTION ^

 Melting rain
 Extinction, scattering, absorption, backscattering and 
 asymmetric scattering coefficients in 1/km versus rain rate, 
 for Marshall-Palmer (MP) drop-size distribution 
 see Sauvageot et al. (1992),
 using Mie Theory, and Liebe '91 dielectric model. Input:
 fGHz: frequency in GHz, TK: Temp. in K, 
 nrain: Number of rain rates between Rmin=0.1 and Rmax=100mm/h
 pam: 0 if no costeta data to be given, 1 if they are needed
 coat: thickness of water coating of ice sphere in mm
 C. M�zler, July 2002.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

miecoated_rain3.m

SOURCE CODE ^

0001 function result = miecoated_rain3(fGHz, TK, nrain, pam, coat)
0002 
0003 % Melting rain
0004 % Extinction, scattering, absorption, backscattering and
0005 % asymmetric scattering coefficients in 1/km versus rain rate,
0006 % for Marshall-Palmer (MP) drop-size distribution
0007 % see Sauvageot et al. (1992),
0008 % using Mie Theory, and Liebe '91 dielectric model. Input:
0009 % fGHz: frequency in GHz, TK: Temp. in K,
0010 % nrain: Number of rain rates between Rmin=0.1 and Rmax=100mm/h
0011 % pam: 0 if no costeta data to be given, 1 if they are needed
0012 % coat: thickness of water coating of ice sphere in mm
0013 % C. M�zler, July 2002.
0014 
0015 opt=1;     
0016 Rmin=0.1;   
0017 nsteps=201;
0018 m1=sqrt(epsice(fGHz, TK));
0019 m2=sqrt(epswater(fGHz, TK));
0020 N0=0.08/10000;              % original MP N0 in 1/mm^4
0021 fact=1000^(1/(nrain-0.99999));
0022 R=Rmin/fact;
0023 nx=(1:nsteps)';
0024 c0=299.793;
0025 for jr = 1:nrain
0026     R=R*fact;
0027     dD=0.025*R^(1/6)/fGHz^0.05;
0028     D=(nx-1)*dD;
0029     y=pi*D*fGHz/c0;
0030     x=max(0,pi*(D-coat)*fGHz/c0);
0031     sigmag=pi*D.*D/4;
0032     LA=4.1/R^0.21;
0033     NMP=N0*exp(-LA*D);
0034     sn=sigmag.*NMP*1000000;
0035     for j = 1:nsteps    
0036         a(j,:)=miecoated(m1,m2,x(j),y(j),opt);
0037     end;
0038     b(:,1)=D;
0039     b(:,2)=a(:,1).*sn;
0040     b(:,3)=a(:,2).*sn;
0041     b(:,4)=a(:,3).*sn;
0042     b(:,5)=a(:,4).*sn;
0043     b(:,6)=a(:,2).*a(:,5).*sn; 
0044     gext= sum(b(:,2))*dD;
0045     gsca= sum(b(:,3))*dD;
0046     gabs= sum(b(:,4))*dD;
0047     gb=   sum(b(:,5))*dD;
0048     gteta=sum(b(:,6))*dD;
0049     res(jr,:)=[R gext gsca gabs gb gteta];
0050 end;
0051     if pam==0
0052         output_parameters='Gext, Gsca, Gabs, Gb';
0053         loglog(res(:,1),res(:,2:5))
0054         legend('Gext','Gsca','Gabs','Gb')
0055         title(sprintf('Propagation Coefficients Versus Melting-Rain Rate at f=%gGHz, T=%gK, coat=%gmm',fGHz,TK,coat))
0056         xlabel('R (mm/h)');    ylabel('Gi(1/km)')
0057     elseif pam==1
0058         output_parameters='Gext, Gsca, Gabs, Gb, Gsca*<costeta>';
0059         loglog(res(:,1),res(:,2:6))
0060         legend('Gext','Gsca','Gabs','Gb','Gsca*<costeta>')
0061         title(sprintf('Propagation Coefficients Versus Melting-Rain Rate at f=%gGHz, T=%gK, coat=%gmm',fGHz,TK,coat))
0062         xlabel('R (mm/h)');    ylabel('Gi(1/km)')
0063     end;
0064 result=res;

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