Changeset 3983 for trunk/LMDZ.COMMON/libf/evolution/read_XIOS_data.F90

Timestamp:

Dec 8, 2025, 11:27:43 AM (7 days ago)

Author:

jbclement

Message:

PEM:

• Removing completely the ice metamorphism computed by a threshold at the end of the PCM (which was commented).
• Addition of a module "metamorphism" to compute the PCM frost at the PEM beginning and give it back to the PCM at the PEM end. The frost is considered as the ice given by the PCM "startfi.nc" which is above the yearly minimum. Thereby, metamorphism is performed through this operation.
• Ice reservoirs representation in the PEM is modernized.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/read_XIOS_data.F90 (modified) (7 diffs)

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

@@ -19 +19 @@
 contains
 !=======================================================================
-SUBROUTINE read_PCM_data(ngrid,nslope,nsoil_PCM,nsol,ps_avg,tsurf_avg,tsurf_avg_y1,tsoil_avg,tsoil_ts,watersurf_density_avg,d_h2oice,d_co2ice, &
-                         ps_ts,q_h2o_ts,q_co2_ts,watersoil_density_ts)
-
-use compute_tend_mod, only: compute_tend
+SUBROUTINE read_PCM_data(ngrid,nslope,nsoil_PCM,nsol,h2ofrost_PCM,co2frost_PCM,ps_avg,tsurf_avg,tsurf_avg_y1,tsoil_avg,tsoil_ts,watersurf_density_avg,d_h2oice,d_co2ice, &
+                         ps_ts,q_h2o_ts,q_co2_ts,watersoil_density_ts,min_h2oice,min_co2ice)
+
 use grid_conversion,  only: lonlat2vect
 use comsoil_h_PEM,    only: soil_pem
+use compute_tend_mod, only: compute_tend
+use metamorphism,     only: compute_frost
 
 implicit none
@@ -30 +31 @@
 ! Arguments
 !----------
-integer, intent(in) :: ngrid, nslope, nsoil_PCM, nsol
+integer,                       intent(in) :: ngrid, nslope, nsoil_PCM, nsol
+real, dimension(ngrid,nslope), intent(in) :: h2ofrost_PCM, co2frost_PCM
 real, dimension(ngrid),                       intent(out) :: ps_avg
-real, dimension(ngrid,nslope),                intent(out) :: tsurf_avg, tsurf_avg_y1, watersurf_density_avg, d_h2oice, d_co2ice
+real, dimension(ngrid,nslope),                intent(out) :: tsurf_avg, tsurf_avg_y1, watersurf_density_avg, d_h2oice, d_co2ice, min_h2oice, min_co2ice
 real, dimension(ngrid,nsoil_PCM,nslope),      intent(out) :: tsoil_avg
 real, dimension(ngrid,nsol),                  intent(out) :: ps_ts, q_h2o_ts, q_co2_ts
@@ -44 +46 @@
 real, dimension(:,:,:,:), allocatable :: var_read_4d
 character(:),             allocatable :: num ! For reading slope variables
-real, dimension(ngrid,nslope,2)       :: min_h2oice, min_co2ice, min_h2ofrost, min_co2frost
+real, dimension(ngrid,nslope,2)       :: min_h2operice, min_co2perice, min_h2ofrost, min_co2frost
 
 ! Code
 !-----
 ! Initialization
-min_h2oice = 0.
-min_co2ice = 0.
+min_h2operice = 0.
+min_co2perice = 0.
 min_h2ofrost = 0.
 min_co2frost = 0.
@@ -80 +82 @@
     call get_var('co2ice'//num,var_read_2d)          ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2frost(:,islope,1))
     call get_var('h2o_ice_s'//num,var_read_2d)       ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2ofrost(:,islope,1))
-    call get_var('watercap'//num,var_read_2d)        ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2oice(:,islope,1))
-    call get_var('perennial_co2ice'//num,var_read_2d); call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2ice(:,islope,1))
+    call get_var('watercap'//num,var_read_2d)        ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2operice(:,islope,1))
+    call get_var('perennial_co2ice'//num,var_read_2d); call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2perice(:,islope,1))
     call get_var('tsurf'//num,var_read_2d)           ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,tsurf_avg_y1(:,islope))
 enddo
@@ -95 +97 @@
 
 ! Allocate and read the variables
-call get_var('ps',var_read_2d) ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,ps_avg)
+call get_var('ps',var_read_2d); call lonlat2vect(nlon,nlat,ngrid,var_read_2d,ps_avg)
 do islope = 1,nslope
     if (nslope /= 1) then
@@ -105 +107 @@
     call get_var('co2ice'//num,var_read_2d)          ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2frost(:,islope,2))
     call get_var('h2o_ice_s'//num,var_read_2d)       ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2ofrost(:,islope,2))
-    call get_var('watercap'//num,var_read_2d)        ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2oice(:,islope,2))
-    call get_var('perennial_co2ice'//num,var_read_2d); call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2ice(:,islope,2))
+    call get_var('watercap'//num,var_read_2d)        ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2operice(:,islope,2))
+    call get_var('perennial_co2ice'//num,var_read_2d); call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2perice(:,islope,2))
     if (soil_pem) then
         call get_var('soiltemp'//num,var_read_3d)
@@ -167 +169 @@
 
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+! Compute frost from yearly minima
+call compute_frost(ngrid,nslope,h2ofrost_PCM,min_h2ofrost,co2frost_PCM,min_co2frost)
+
+!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 ! Compute ice tendencies from yearly minima
 write(*,*) '> Computing surface ice tendencies'
-call compute_tend(ngrid,nslope,min_h2ofrost + min_h2oice,d_h2oice)
+call compute_tend(ngrid,nslope,min_h2operice + min_h2ofrost,d_h2oice)
 write(*,*) 'H2O ice tendencies (min/max):', minval(d_h2oice), maxval(d_h2oice)
-call compute_tend(ngrid,nslope,min_co2frost + min_co2ice,d_co2ice)
+call compute_tend(ngrid,nslope,min_co2perice + min_co2frost,d_co2ice)
 write(*,*) 'CO2 ice tendencies (min/max):', minval(d_co2ice), maxval(d_co2ice)