Changeset 5095 for LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/predict_mom07.F90

Timestamp:

Jul 22, 2024, 9:46:57 AM (2 months ago)

Author:

abarral

Message:

Revert cosp*/ from the trunk, as it's external code
Add missing bits from FCM2 source

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/predict_mom07.F90 (modified) (2 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/predict_mom07.F90

@@ -7 +7 @@
             implicit none
         
-            real(kind=8) :: a1,a2,a3,b1,b2,b3,c1,c2,c3
-            real(kind=8) :: m2,tc,n,m,a_,b_,c_,A,B,C,n2
+            real*8 :: a1,a2,a3,b1,b2,b3,c1,c2,c3
+            real*8 :: m2,tc,n,m,a_,b_,c_,A,B,C,n2
         
             a1=      13.6078
@@ -30 +30 @@
             
         ! predict m from m2 and tc
-                if(m2/=-9999) then
+                if(m2.ne.-9999) then 
                 m=A*exp(B*tc)*m2**C
                 endif
         ! get m2 if mass-dimension relationship not proportional to D**2
-                if(m2==-9999) then
+                if(m2.eq.-9999) then 
                 m2=(m/(A*exp(B*tc)))**(1.0/C)    
                 endif