Changeset 2887 for trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90

Timestamp:

Feb 1, 2023, 9:34:37 AM (22 months ago)

Author:

llange

Message:

Mars PCM
Include geothermal flux as boundary conditions when solving the
diffusion equation for the soil temperature
By default, the geothermal flux is set to 0. To have geothermal flux >
0., it must be specified in the startfi. (I suggest to use a value of
30e-3 in this case)
Will be use with the PEM for studies of "paleosubsurfaces"
LL

File:

• trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90

@@ -25 +25 @@
   real,save,allocatable :: beta(:,:)        ! beta_k coefficients
   real,save :: mu
+  real,save,allocatable :: flux_geo(:)       ! Geothermal Flux (W/m^2)
 
-!$OMP THREADPRIVATE(tsoil,mthermdiff,thermdiff,coefq,coefd,alph,beta,mu )
+!$OMP THREADPRIVATE(tsoil,mthermdiff,thermdiff,coefq,coefd,alph,beta,mu,flux_geo)
 
 contains
@@ -47 +48 @@
     allocate(alph(ngrid,nsoilmx-1))
     allocate(beta(ngrid,nsoilmx-1))
-
+    allocate(flux_geo(ngrid))
  
   end subroutine ini_comsoil_h
@@ -66 +67 @@
     if (allocated(alph)) deallocate(alph)
     if (allocated(beta)) deallocate(beta)
-
+    if (allocated(flux_geo)) deallocate(flux_geo)
   end subroutine end_comsoil_h