Changeset 3543 for trunk/LMDZ.COMMON/libf/evolution/pem.F90

Timestamp:

Dec 10, 2024, 2:16:29 PM (13 days ago)

Author:

jbclement

Message:

PEM:
Using the correct subroutine to output 'ice_porefilling'.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pem.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/evolution/pem.F90

@@ -1003 +1003 @@
             call writediagsoilpem(ngrid,'tsoil_PEM_slope'//str2,'tsoil_PEM','K',3,tsoil_PEM(:,:,islope))
             call writediagsoilpem(ngrid,'inertiesoil_PEM_slope'//str2,'TI_PEM','K',3,TI_PEM(:,:,islope))
-            if (icetable_dynamic) call writediagpem(ngrid,'ice_porefilling'//str2,'ice pore filling','-',3,ice_porefilling(:,:,islope))
+            if (icetable_dynamic) call writediagsoilpem(ngrid,'ice_porefilling'//str2,'ice pore filling','-',3,ice_porefilling(:,:,islope))
             if (adsorption_pem) then
                 call writediagsoilpem(ngrid,'co2_ads_slope'//str2,'co2_ads','K',3,co2_adsorbded_phys(:,:,islope))