Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/neutral.f90

Timestamp:

Jul 24, 2024, 4:23:34 PM (2 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/neutral.f90 (modified) (3 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/neutral.f90

@@ -43 +43 @@
   do i = 1, klon
 
-    if (u10_mps(i) < 0.) u10_mps(i) = 0.0
+    IF (u10_mps(i) < 0.) u10_mps(i) = 0.0
 
-    if  (obklen_m(i) < 0.) THEN
+    IF  (obklen_m(i) < 0.) THEN
       phi = (1. - 160. / obklen_m(i))**(-0.25)
       phi_inv = 1. / phi
@@ -55 +55 @@
       f3 = atan(dum1)
       psi = 2. * log(f1) + log(f2) - 2. * f3 + pi / 2.
-    else if (obklen_m(i) > 0.) THEN
+    ELSE IF (obklen_m(i) > 0.) THEN
       psi = -50. / obklen_m(i)
     end if
@@ -61 +61 @@
     u10n_mps(i) = u10_mps(i) + (ustar_mps(i) * psi / von_karman)
     ! u10n set to 0. if -1 < obklen < 20
-    if ((obklen_m(i)>-1.).and.(obklen_m(i)<20.)) THEN
+    IF ((obklen_m(i)>-1.).AND.(obklen_m(i)<20.)) THEN
       u10n_mps(i) = 0.
     endif
-    if (u10n_mps(i) < 0.) u10n_mps(i) = 0.0
+    IF (u10n_mps(i) < 0.) u10n_mps(i) = 0.0
 
   enddo