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

Timestamp:

Dec 4, 2024, 4:04:54 PM (3 weeks ago)

Author:

jbclement

Message:

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

File:

• trunk/LMDZ.COMMON/libf/evolution/soil_settings_PEM_mod.F90 (modified) (3 diffs)

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

@@ -15 +15 @@
 
 !=======================================================================
-!  Author: LL, based on work by  Ehouarn Millour (07/2006)
+!  Author: LL, based on work by Ehouarn Millour (07/2006)
 !
 !  Purpose: Read and/or initialise soil depths and properties
@@ -53 +53 @@
 alpha = 2
 do iloop = 0,nsoil_PCM - 1
-    mlayer_PEM(iloop) = lay1*(1+iloop**2.9*(1-exp(-real(iloop)/20.)))
+    mlayer_PEM(iloop) = lay1*(1 + iloop**2.9*(1 - exp(-real(iloop)/20.)))
 enddo
 
-do iloop = 0, nsoil_PEM-1
-    if(lay1*(alpha**(iloop-0.5)) > mlayer_PEM(nsoil_PCM-1)) then
+do iloop = 0,nsoil_PEM - 1
+    if(lay1*(alpha**(iloop-0.5)) > mlayer_PEM(nsoil_PCM - 1)) then
         index_powerlaw = iloop
         exit
@@ -63 +63 @@
 enddo
 
-do iloop = nsoil_PCM, nsoil_PEM-1 
-    mlayer_PEM(iloop) = lay1*(alpha**(index_powerlaw + (iloop-nsoil_PCM)-0.5))
+do iloop = nsoil_PCM,nsoil_PEM - 1
+    mlayer_PEM(iloop) = lay1*(alpha**(index_powerlaw + (iloop - nsoil_PCM) - 0.5))
 enddo
 
 ! Build layer_PEM()
-do iloop=1,nsoil_PEM-1
-layer_PEM(iloop)=(mlayer_PEM(iloop)+mlayer_PEM(iloop-1))/2.
+do iloop = 1,nsoil_PEM - 1
+layer_PEM(iloop) = (mlayer_PEM(iloop) + mlayer_PEM(iloop - 1))/2.
 enddo
-layer_PEM(nsoil_PEM)=2*mlayer_PEM(nsoil_PEM-1) - mlayer_PEM(nsoil_PEM-2)
+layer_PEM(nsoil_PEM) = 2*mlayer_PEM(nsoil_PEM - 1) - mlayer_PEM(nsoil_PEM - 2)