Changeset 2887 for trunk/LMDZ.MARS/libf/phymars/soil_settings.F

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/soil_settings.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/soil_settings.F

@@ -2 +2 @@
 
 !      use netcdf
-      use comsoil_h, only: layer, mlayer, inertiedat, volcapa
+      use comsoil_h, only: layer, mlayer, inertiedat, volcapa,flux_geo
       use iostart, only: inquire_field_ndims, get_var, get_field,
      &                   inquire_field, inquire_dimension_length
@@ -435 +435 @@
         deallocate(oldtsoil)
       endif ! of if (interpol)
-      
-! 6. Report min and max values of soil temperatures and thermal inertias
+
+
+      
+! 6. Load Geothermal Flux
+! ----------------------------------------------------------------------
+
+
+      call get_field("flux_geo",flux_geo,found,timeindex=indextime)
+      if (.not.found) then
+           write(*,*) "soil_settings: Failed loading <flux_geo>"
+           write(*,*) "soil_settings: Will put  <flux_geo> to 0."
+           flux_geo(:) = 0.
+      endif
+
+
+      
+! 7. Report min and max values of soil temperatures and thermal inertias
 ! ----------------------------------------------------------------------