source: LMDZ5/trunk/libf/phylmd/cosp/predict_mom07.F90 @ 3982

Last change on this file since 3982 was 2432, checked in by idelkadi, 9 years ago

Correction : rajout de fichiers manquants pour les simulateurs Modis, Cloudsat, Misr et Isccp

File size: 1.3 KB
Line 
1! $Revision: 88 $, $Date: 2013-11-13 15:08:38 +0100 (mer. 13 nov. 2013) $
2! $URL: http://cfmip-obs-sim.googlecode.com/svn/stable/v1.4.0/predict_mom07.F90 $
3        subroutine predict_mom07(m2,tc,n,m)
4       
5        ! subroutine to predict nth moment m given m2, tc, n
6        ! if m2=-9999 then the routine will predict m2 given m,n,tc
7            implicit none
8       
9            real*8 :: a1,a2,a3,b1,b2,b3,c1,c2,c3
10            real*8 :: m2,tc,n,m,a_,b_,c_,A,B,C,n2
11       
12            a1=      13.6078
13            a2=     -7.76062
14            a3=     0.478694
15            b1=   -0.0360722
16            b2=    0.0150830
17            b3=   0.00149453
18            c1=     0.806856
19            c2=   0.00581022
20            c3=    0.0456723
21       
22            n2=n*n
23            a_=a1+a2*n+a3*n2
24            b_=b1+b2*n+b3*n2
25            c_=c1+c2*n+c3*n2
26       
27            A=exp(a_)
28            B=b_
29            C=c_
30           
31        ! predict m from m2 and tc
32                if(m2.ne.-9999) then
33                m=A*exp(B*tc)*m2**C
34                endif
35        ! get m2 if mass-dimension relationship not proportional to D**2
36                if(m2.eq.-9999) then
37                m2=(m/(A*exp(B*tc)))**(1.0/C)   
38                endif
39       
40                return
41        end subroutine predict_mom07
Note: See TracBrowser for help on using the repository browser.