- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/predict_mom07.F90
r2432 r5082 7 7 implicit none 8 8 9 real *8:: a1,a2,a3,b1,b2,b3,c1,c2,c310 real *8:: m2,tc,n,m,a_,b_,c_,A,B,C,n29 real(kind=8) :: a1,a2,a3,b1,b2,b3,c1,c2,c3 10 real(kind=8) :: m2,tc,n,m,a_,b_,c_,A,B,C,n2 11 11 12 12 a1= 13.6078 … … 30 30 31 31 ! predict m from m2 and tc 32 if(m2 .ne.-9999) then32 if(m2/=-9999) then 33 33 m=A*exp(B*tc)*m2**C 34 34 endif 35 35 ! get m2 if mass-dimension relationship not proportional to D**2 36 if(m2 .eq.-9999) then36 if(m2==-9999) then 37 37 m2=(m/(A*exp(B*tc)))**(1.0/C) 38 38 endif
Note: See TracChangeset
for help on using the changeset viewer.