Changeset 3482

Timestamp:

Oct 24, 2024, 1:10:51 PM (4 weeks ago)

Author:

afalco

Message:

Pluto PCM: renamed soil.F to soil.F90.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/soil.F90 (moved) (moved from trunk/LMDZ.PLUTO/libf/phypluto/soil.F) (8 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/soil.F90

@@ -16 +16 @@
 !
 !  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.h
 !-----------------------------------------------------------------------
@@ -25 +25 @@
 !  ---------
 !  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 
+      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
       logical,intent(in) :: lastcall
       real,intent(in) :: therm_i(ngrid,nsoil) ! thermal inertia
@@ -46 +46 @@
       real,save :: mu
 !$OMP THREADPRIVATE(mthermdiff,thermdiff,coefq,coefd,alph,beta,mu)
-            
+
 ! local variables:
       integer ig,ik
@@ -56 +56 @@
       if (firstcall) then
       ! note: firstcall is set to .true. or .false. by the caller
-      !       and not changed by soil.F 
+      !       and not changed by soil.F
 
       ALLOCATE(mthermdiff(ngrid,0:nsoil-1)) ! mid-layer thermal diffusivity
@@ -100 +100 @@
           coefd(ig,ik)=thermdiff(ig,ik)/(mlayer(ik)-mlayer(ik-1))
         enddo
-        
+
         ! alph_{N-1}
         alph(ig,nsoil-1)=coefd(ig,nsoil-1)/
@@ -120 +120 @@
       !write(*,*)'soil: ig=',ig,' capcal(ig)=',capcal(ig)
       enddo ! of do ig=1,ngrid
-      
+
       ! Additional checks: is the vertical discretization sufficient
       ! to resolve diurnal and annual waves?
@@ -153 +153 @@
         endif
       enddo ! of do ig=1,ngrid
-      
+
       else ! of if (firstcall)
 
@@ -171 +171 @@
         enddo
       enddo
-      
+
       endif! of if (firstcall)