Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/lmdz_1dutils.f90

Timestamp:

Aug 2, 2024, 2:12:03 PM (3 months ago)

Author:

abarral

Message:

Add missing klon on strataer_emiss_mod.F90
Correct various missing explicit declarations
Replace tabs by spaces (tabs are not part of the fortran charset)
Continue cleaning modules
Removed unused arguments and variables

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/lmdz_1dutils.f90 (modified) (12 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/lmdz_1dutils.f90

@@ -4 +4 @@
           disvert0, advect_vert, advect_va, lstendh, nudge_rht_init, nudge_uv_init, &
           nudge_rht, nudge_uv, interp2_case_vertical
+
 CONTAINS
   REAL FUNCTION fq_sat(kelvin, millibar)
@@ -67 +68 @@
     USE lmdz_fcs_gcssold, ONLY: imp_fcg_gcssold, ts_fcg_gcssold, Tp_fcg_gcssold, Tp_ini_gcssold, xTurb_fcg_gcssold
     USE lmdz_tsoilnudge, ONLY: nudge_tsoil, isoil_nudge, Tsoil_nudge, tau_soil_nudge
+    USE lmdz_compar1d
+
+    IMPLICIT NONE
     !-----------------------------------------------------------------------
     !     Auteurs :   A. Lahellec  .
-
-    !   Declarations :
-    !   --------------
-
-    include "compar1d.h"
-    include "fcg_racmo.h"
-
-
-    !   local:
-    !   ------
-
-    !      CHARACTER ch1*72,ch2*72,ch3*72,ch4*12
 
     !  -------------------------------------------------------------------
@@ -1148 +1140 @@
     !---------------------------------------------------------------
 
-    do l = 1, llm
+    DO l = 1, llm
       zw(l) = dt * w(l)
       zm(l) = plev(l) - plev(l + 1)
@@ -1156 +1148 @@
     zw(llm + 1) = 0.
 
-    do l = 1, llm
+    DO l = 1, llm
       qold = q(l)
       q(l) = (q(l) * zm(l) + zwq(l + 1) - zwq(l)) / (zm(l) + zw(l + 1) - zw(l))
@@ -1186 +1178 @@
     REAL alpha, omgdown, omgup
 
-    do l = 1, llm
+    DO l = 1, llm
       IF(l==1) THEN
         !si omgup pour la couche 1, alors tendance nulle
@@ -1271 +1263 @@
     cor(:) = rkappa * temp * (1. + q(:, 1) * rv / rd) / (play * (1. + q(:, 1)))
 
-    do k = 2, llm - 1
+    DO k = 2, llm - 1
 
       dph  (k - 1) = (play(k) - play(k - 1))
@@ -1286 +1278 @@
     dtdp (llm) = dtdp (llm - 1)
 
-    do k = 2, llm - 1
+    DO k = 2, llm - 1
       omdn = max(0.0, omega(k + 1))
       omup = min(0.0, omega(k))
@@ -1325 +1317 @@
 
     IMPLICIT NONE
- INCLUDE "FCTTRE.h"
+    INCLUDE "FCTTRE.h"
 
     ! ========================================================
@@ -1399 +1391 @@
 
     IMPLICIT NONE
- INCLUDE "FCTTRE.h"
+    INCLUDE "FCTTRE.h"
 
     ! ========================================================
@@ -1593 +1585 @@
     !      enddo
 
-    do l = 1, llm
+    DO l = 1, llm
 
       IF (play(l)>=plev_prof_cas(nlev_cas)) THEN
@@ -1602 +1594 @@
 
         IF (play(l)<=plev_prof_cas(1)) THEN
-          do k = 1, nlev_cas - 1
+          DO k = 1, nlev_cas - 1
             IF (play(l)<=plev_prof_cas(k).AND. play(l)>plev_prof_cas(k + 1)) THEN
               k1 = k
@@ -1612 +1604 @@
             WRITE(*, *) 'PB! k1, k2 = ', k1, k2
             WRITE(*, *) 'l,play(l) = ', l, play(l) / 100
-            do k = 1, nlev_cas - 1
+            DO k = 1, nlev_cas - 1
               WRITE(*, *) 'k,plev_prof_cas(k) = ', k, plev_prof_cas(k) / 100
             enddo