Changeset 4121 for trunk/LMDZ.GENERIC

Timestamp:

Mar 12, 2026, 3:17:58 PM (4 weeks ago)

Author:

cmetz

Message:

GENERIC PCM - CM
Added subsurface outputs for XIOS, and outputting tsoil for XIOS.

Location:

trunk/LMDZ.GENERIC

Files:

• changelog.txt (modified) (1 diff)
• deftank/xios/context_pcm_physics.xml (modified) (5 diffs)
• deftank/xios/field_def_physics.xml (modified) (2 diffs)
• libf/phygeneric/physiq_mod.F90 (modified) (1 diff)
• libf/phygeneric/xios_output_mod.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/changelog.txt

@@ -2205 +2205 @@
 must be a call to getin("nsoilmx",nsoilmx) in rcm1d before the call to routine
 phys_state_var_init.
+
+== 12/03/2026 == CM
+Added subsurface outputs for XIOS, and outputting tsoil for XIOS.

trunk/LMDZ.GENERIC/deftank/xios/context_pcm_physics.xml

@@ -6 +6 @@
     <domain_definition>
         <!-- NB: dom_glo data_dim is 1 for DYNAMICO and 2 for LMDZ.COMMON -->
-        <domain_group id="dom_glo" data_dim="1">
+        <domain_group id="dom_glo" data_dim="2">
           <domain id="dom_glo" />
         </domain_group>
@@ -49 +49 @@
 
         <!-- NB: dom_out data_ref should be "dom_glo" with LMDZ.COMMON -->
-        <domain id="dom_out" domain_ref="dom_128_96"/>
+        <domain id="dom_out" domain_ref="dom_glo"/>
     </domain_definition>
 
@@ -64 +64 @@
               standard_name="Pseudo-altitude of model vertical levels"
               unit="km">
-        </axis>
+      </axis>
+        <axis id="soildepth"
+              standard_name="Soil mid-layer depth"
+              unit="m">
+      </axis>
+
         <axis id="controle_axe"
               standard_name="Controle"
@@ -93 +98 @@
             <domain domain_ref="dom_glo" />
             <axis axis_ref="altitude" />
+        </grid>
+        <grid id="grid_soil">
+            <domain domain_ref="dom_glo" />
+            <axis axis_ref="soildepth" />
         </grid>
         <grid id="grid_2D">
@@ -119 +128 @@
             <domain domain_ref="dom_out" />
             <axis axis_ref="altitude" />
+        </grid>
+        <grid id="grid_soil_out">
+            <domain domain_ref="dom_out" />
+            <axis axis_ref="soildepth" />
         </grid>
         <grid id="grid_2D_out">

trunk/LMDZ.GENERIC/deftank/xios/field_def_physics.xml

@@ -337 +337 @@
             <field id="temperature"
                    long_name="Atmospheric temperature"
-                   unit="K" />
+                   unit="K" />
             <field id="teta"
                    long_name="Potential temperature"
@@ -555 +555 @@
                    long_name="zdqconv_prec"
                    unit="" />
-
-        </field_group>
+   
+   </field_group>
+
+
+   <!--Subsurface 3D variables-->
+        <field_group id="soil_fields_3D" grid_ref="grid_soil">
+            <field id="tsoil"
+                   long_name="Soil temperature"
+                   unit="K" />
+        </field_group>
 
         <!--Spectral IR 3D variables-->

trunk/LMDZ.GENERIC/libf/phygeneric/physiq_mod.F90

@@ -2778 +2778 @@
       endif
 
+      if(callsoil)then
+         call write_output("tsoil","Soil temperature","K",tsoil)
+      endif
+
 #else
    comm_HR_SW(1:ngrid,1:nlayer) = zdtsw(1:ngrid,1:nlayer)

trunk/LMDZ.GENERIC/libf/phygeneric/xios_output_mod.F90

@@ -31 +31 @@
   USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
   USE nrtype, ONLY: pi
+  USE comsoil_h, ONLY: mlayer
+
 #ifdef CPP_XIOS
   USE xios
@@ -62 +64 @@
     CALL xios_set_axis_attr("altitude", n_glo=size(pseudoalt), value=pseudoalt,&
                             unit="km",positive="up")
+    IF (prt_level>=10) WRITE(lunout,*) "initialize_xios_output: call xios_set_axis_attr for soildepth"
+    CALL xios_set_axis_attr("soildepth", n_glo=size(mlayer), value=mlayer, &
+                            unit="m", positive="down")
     if (prt_level >=10) write(lunout,*) "initialize_xios_output: call xios_set_axis_attr for IR_Wavenumber"
     write(lunout,*) "writing IR_Wavenumber now in initialize_xios_output"