|
Last change
on this file since 5441 was
5095,
checked in by abarral, 5 months ago
|
|
Revert cosp*/ from the trunk, as it's external code
Add missing bits from FCM2 source
|
|
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.
