Last change
on this file since 3312 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.