Last change
on this file since 3595 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
|
Rev | Line | |
---|
[2432] | 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.