Changeset 3532 for trunk/LMDZ.COMMON/libf/evolution/compute_soiltemp_mod.F90

Timestamp:

Dec 4, 2024, 4:04:54 PM (7 months ago)

Author:

jbclement

Message:

PEM:
Removing unecessary module/subroutine "interpol_TI_PEM2PCM.F90" + Few small corrections/cleanings.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/compute_soiltemp_mod.F90 (modified) (11 diffs)

trunk/LMDZ.COMMON/libf/evolution/compute_soiltemp_mod.F90

@@ -4 +4 @@
 !-----------------------------------------------------------------------
 !  Author: LL
-!  Purpose: This module gathers the different routines used in the PEM to compute the soil temperature evolution and initialisation. 
-!  
-!  Note: depths of layers and mid-layers, soil thermal inertia and 
+!  Purpose: This module gathers the different routines used in the PEM to compute the soil temperature evolution and initialisation.
+!
+!  Note: depths of layers and mid-layers, soil thermal inertia and
 !        heat capacity are commons in comsoil_PEM.h
 !-----------------------------------------------------------------------
@@ -12 +12 @@
 !=======================================================================
 
-
-
 SUBROUTINE compute_tsoil_pem(ngrid,nsoil,firstcall,therm_i,timestep,tsurf,tsoil)
 
@@ -24 +22 @@
 !  Author: LL
 !  Purpose: Compute soil temperature using an implict 1st order scheme
-!  
-!  Note: depths of layers and mid-layers, soil thermal inertia and 
+!
+!  Note: depths of layers and mid-layers, soil thermal inertia and
 !        heat capacity are commons in comsoil_PEM.h
 !-----------------------------------------------------------------------
@@ -31 +29 @@
 #include "dimensions.h"
 
-!-----------------------------------------------------------------------
-!  arguments
-!  ---------
-!  inputs:
-integer,                      intent(in) :: ngrid     ! number of (horizontal) grid-points 
-integer,                      intent(in) :: nsoil     ! number of soil layers  
-logical,                      intent(in) :: firstcall ! identifier for initialization call 
+! Inputs:
+! -------
+integer,                      intent(in) :: ngrid     ! number of (horizontal) grid-points
+integer,                      intent(in) :: nsoil     ! number of soil layers
+logical,                      intent(in) :: firstcall ! identifier for initialization call
 real, dimension(ngrid,nsoil), intent(in) :: therm_i   ! thermal inertia [SI]
 real,                         intent(in) :: timestep  ! time step [s]
 real, dimension(ngrid),       intent(in) :: tsurf     ! surface temperature [K]
- 
-! outputs:
+! Outputs:
+!---------
 real, dimension(ngrid,nsoil), intent(inout) :: tsoil ! soil (mid-layer) temperature [K]
-! local variables:
-integer :: ig, ik    
+! Local:
+!-------
+integer :: ig, ik
 
 ! 0. Initialisations and preprocessing step
@@ -52 +49 @@
     do ig = 1,ngrid
         do ik = 0,nsoil - 1
-            mthermdiff_PEM(ig,ik) = therm_i(ig,ik + 1)*therm_i(ig,ik + 1)/volcapa    
+            mthermdiff_PEM(ig,ik) = therm_i(ig,ik + 1)*therm_i(ig,ik + 1)/volcapa
         enddo
     enddo
@@ -64 +61 @@
     enddo
 
-! 0.3 Build coefficients mu_PEM, q_{k+1/2}, d_k, alph_PEMa_k and capcal
+! 0.3 Build coefficients mu_PEM, q_{k+1/2}, d_k, alph_PEM
     ! mu_PEM
     mu_PEM = mlayer_PEM(0)/(mlayer_PEM(1) - mlayer_PEM(0))
@@ -99 +96 @@
 ! Other layers:
         do ik = 1,nsoil - 1
-                tsoil(ig,ik + 1) = alph_PEM(ig,ik)*tsoil(ig,ik) + beta_PEM(ig,ik) 
+                tsoil(ig,ik + 1) = alph_PEM(ig,ik)*tsoil(ig,ik) + beta_PEM(ig,ik)
         enddo
     enddo
@@ -139 +136 @@
 #include "dimensions.h"
 
-!-----------------------------------------------------------------------
-!  arguments
-!  ---------
-!  inputs:
+! Inputs:
+!--------
 integer,                      intent(in) :: ngrid   ! number of (horizontal) grid-points
 integer,                      intent(in) :: nsoil   ! number of soil layers
 real, dimension(ngrid,nsoil), intent(in) :: therm_i ! thermal inertia [SI]
 real, dimension(ngrid),       intent(in) :: tsurf   ! surface temperature [K]
-
-! outputs:
+! Outputs:
+!---------
 real, dimension(ngrid,nsoil), intent(inout) :: tsoil ! soil (mid-layer) temperature [K]
-! local variables:
+! Local:
+!-------
 integer :: ig, ik, iloop
 
@@ -169 +165 @@
 enddo
 
-! 0.3 Build coefficients mu_PEM, q_{k+1/2}, d_k, alph_PEMa_k and capcal
+! 0.3 Build coefficients mu_PEM, q_{k+1/2}, d_k, alph_PEM
 ! mu_PEM
 mu_PEM = mlayer_PEM(0)/(mlayer_PEM(1) - mlayer_PEM(0))
@@ -179 +175 @@
 do ig = 1,ngrid
     ! d_k
-    do ik = 1,nsoil-1
+    do ik = 1,nsoil - 1
         coefd_PEM(ig,ik) = thermdiff_PEM(ig,ik)/(mlayer_PEM(ik) - mlayer_PEM(ik - 1))
     enddo
@@ -207 +203 @@
 !  2. Compute soil temperatures
 do iloop = 1,10 !just convergence
-    ! First layer:
-    do ig = 1,ngrid
-        tsoil(ig,1)=(tsurf(ig) + mu_PEM*beta_PEM(ig,1)*thermdiff_PEM(ig,1)/mthermdiff_PEM(ig,0))/ &
-                   (1. + mu_PEM*(1. - alph_PEM(ig,1))*thermdiff_PEM(ig,1)/mthermdiff_PEM(ig,0))
-    ! Other layers:
-    do ik = 1,nsoil - 1
-        tsoil(ig,ik + 1) = alph_PEM(ig,ik)*tsoil(ig,ik) + beta_PEM(ig,ik)
-    enddo
-enddo
-
+    do ig = 1,ngrid
+        ! First layer:
+        tsoil(ig,1) = (tsurf(ig) + mu_PEM*beta_PEM(ig,1)*thermdiff_PEM(ig,1)/mthermdiff_PEM(ig,0))/ &
+                      (1. + mu_PEM*(1. - alph_PEM(ig,1))*thermdiff_PEM(ig,1)/mthermdiff_PEM(ig,0))
+        ! Other layers:
+        do ik = 1,nsoil - 1
+            tsoil(ig,ik + 1) = alph_PEM(ig,ik)*tsoil(ig,ik) + beta_PEM(ig,ik)
+        enddo
+    enddo
 enddo ! iloop