Changeset 5282 for LMDZ6/trunk/libf/phylmd/ocean_cpl_mod.f90

Timestamp:

Oct 28, 2024, 1:11:48 PM (6 hours ago)

Author:

abarral

Message:

Turn iniprint.h clesphys.h into modules
Remove unused description.h

File:

• LMDZ6/trunk/libf/phylmd/ocean_cpl_mod.f90 (modified) (14 diffs)

LMDZ6/trunk/libf/phylmd/ocean_cpl_mod.f90

@@ -8 +8 @@
 !
 
-  IMPLICIT NONE
+  USE clesphys_mod_h
+    IMPLICIT NONE
   PRIVATE
 
@@ -40 +41 @@
 ! Initialize module cpl_init
     CALL cpl_init(dtime, rlon, rlat)
-    
+
   END SUBROUTINE ocean_cpl_init
 !
@@ -60 +61 @@
 !
 ! This subroutine treats the "open ocean", all grid points that are not entierly covered
-! by ice. The subroutine first receives fields from coupler, then some calculations at 
+! by ice. The subroutine first receives fields from coupler, then some calculations at
 ! surface is done and finally it sends some fields to the coupler.
 !
@@ -82 +83 @@
          cpl_send_ocean_fields
     use config_ocean_skin_m, only: activate_ocean_skin
-
-
-    INCLUDE "clesphys.h"
-!    
-! Input arguments  
+    USE clesphys_mod_h
+! Input arguments
 !****************************************************************************************
     INTEGER, INTENT(IN)                      :: itime, knon
@@ -132 +130 @@
     REAL, DIMENSION(klon), INTENT(INOUT)     :: snow
     REAL, DIMENSION(klon), INTENT(INOUT)     :: agesno
-  
+
 ! Output arguments
 !****************************************************************************************
@@ -139 +137 @@
     REAL, DIMENSION(klon), INTENT(OUT)       :: flux_u1, flux_v1
     REAL, DIMENSION(klon), INTENT(OUT)       :: tsurf_new
-    REAL, DIMENSION(klon), INTENT(OUT)       :: dflux_s, dflux_l      
+    REAL, DIMENSION(klon), INTENT(OUT)       :: dflux_s, dflux_l
     REAL, intent(out):: sens_prec_liq(:) ! (knon)
 
     REAL, INTENT(OUT):: sss(:) ! (klon)
     ! bulk salinity of the surface layer of the ocean, in ppt
-  
+
 
 ! Local variables
@@ -156 +154 @@
     REAL, DIMENSION(klon) :: u1_lay, v1_lay
     LOGICAL               :: check=.FALSE.
-    REAL sens_prec_sol(knon)  
-    REAL, DIMENSION(klon) :: lat_prec_liq, lat_prec_sol    
+    REAL sens_prec_sol(knon)
+    REAL, DIMENSION(klon) :: lat_prec_liq, lat_prec_sol
 
 ! End definitions
@@ -180 +178 @@
     agesno(:) = 0.
     lat_prec_liq = 0.; lat_prec_sol = 0.
-    
+
 
     DO i = 1, knon
@@ -203 +201 @@
 
     ! assertion: tsurf_new == tsurf_cpl
-    
+
     do j = 1, knon
       i = knindex(j)
@@ -213 +211 @@
 
 
-    
+
 ! - Flux calculation at first modele level for U and V
     CALL calcul_flux_wind(knon, dtime, &
@@ -219 +217 @@
          AcoefU, AcoefV, BcoefU, BcoefV, &
          p1lay, temp_air, &
-         flux_u1, flux_v1)  
+         flux_u1, flux_v1)
 
 !****************************************************************************************
@@ -226 +224 @@
 !****************************************************************************************
     fder_new(:) = fder_old(:) + dflux_s(:) + dflux_l(:)
-    
+
     iloc = MAXLOC(fder_new(1:klon))
     IF (check .AND. fder_new(iloc(1))> 0.) THEN
@@ -262 +260 @@
        tsurf_new, dflux_s, dflux_l, rhoa)
 !
-! This subroutine treats the ocean where there is ice. The subroutine first receives 
-! fields from coupler, then some calculations at surface is done and finally sends 
+! This subroutine treats the ocean where there is ice. The subroutine first receives
+! fields from coupler, then some calculations at surface is done and finally sends
 ! some fields to the coupler.
-!    
+!
     USE yomcst_mod_h, ONLY: RPI, RCLUM, RHPLA, RKBOL, RNAVO                   &
           , RDAY, REA, REPSM, RSIYEA, RSIDAY, ROMEGA                  &
@@ -285 +283 @@
 
 
-    INCLUDE "clesphys.h"
 
 ! Input arguments