Changeset 2916

Timestamp:

Mar 16, 2023, 11:50:06 AM (22 months ago)

Author:

emillour

Message:

Mars PCM
Add a "diagsoil" flag to trigger outputing a diagsoil.nc file
(default is diagsoil=.false.)
EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/conf_phys.F (modified) (2 diffs)
• libf/phymars/physiq_mod.F (modified) (1 diff)
• libf/phymars/writediagsoil.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/README

@@ -3949 +3949 @@
 == 14/03/2023 == RV
 Move the initialisation of inertiedat in the subslope dimension for the case where we start without startfi.
+
+== 16/03/2023 == EM
+Add a "diagsoil" flag to trigger outputing a diagsoil.nc file
+(default is diagsoil=.false.)
+

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

@@ -45 +45 @@
       use datafile_mod, only: datadir
       use wstats_mod, only: callstats
+      use writediagsoil_mod, only: diagsoil
       use calchim_mod, only: ichemistry
       use co2condens_mod, only: scavco2cond
@@ -133 +134 @@
          write(*,*) " callstats = ",callstats
 
+         write(*,*) "Write sub-surface fields in file diagsoil.nc ?"
+         diagsoil=.false. ! default value
+         call getin_p("diagsoil",diagsoil)
+         write(*,*) " diagsoil = ",diagsoil
+         
          write(*,*) "Save EOF profiles in file 'profiles' for ",
      &              "Climate Database?"

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

@@ -90 +90 @@
       use phyetat0_mod, only: phyetat0, tab_cntrl_mod
       use wstats_mod, only: callstats, wstats, mkstats
+      use writediagsoil_mod, only: writediagsoil
       use eofdump_mod, only: eofdump
       USE vertical_layers_mod, ONLY: ap,bp,aps,bps,presnivs,pseudoalt

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

@@ +1 @@
+module writediagsoil_mod
+
+implicit none
+
+logical,save :: diagsoil ! global flag to trigger generating a diagsoil.nc
+                         ! file or not. Initialized in conf_phys()
+
+contains
+
 subroutine writediagsoil(ngrid,name,title,units,dimpx,px)
 
@@ -16 +25 @@
 use time_phylmdz_mod, only: ecritphy, day_step, iphysiq
 use mod_phys_lmdz_para, only : is_mpi_root, is_master, gather
-use mod_grid_phy_lmdz, only : klon_glo, Grid1Dto2D_glo, &
-                              nbp_lon, nbp_lat
+use mod_grid_phy_lmdz, only : klon_glo, Grid1Dto2D_glo, nbp_lon, nbp_lat
+use mod_grid_phy_lmdz, only : grid_type, unstructured
 
 implicit none
@@ -72 +81 @@
 real px2(ngrid)
 #endif
+
+! 0. Do we ouput a diagsoil.nc file? If not just bail out now.
+if(.not.diagsoil) return
+
+! additional check: one can only output diagsoil.nc files
+! in lon-lat case (or 1D)
+if (grid_type==unstructured) then
+  write(*,*) "writediagsoil: Error !!!"
+  write(*,*) "diagsoil.nc outputs not possible on unstructured grids!!"
+  call abort_physic("writediagsoil","impossible on unstructured grid",1)
+endif
 
 ! 1. Initialization step
@@ -398 +418 @@
 
 end subroutine writediagsoil
+
+end module writediagsoil_mod