Changeset 4071 for trunk/LMDZ.COMMON/libf/evolution/xios_data.F90

Timestamp:

Feb 16, 2026, 10:28:56 AM (13 days ago)

Author:

jbclement

Message:

PEM:

• Making the computation of ice tendencies more reliable by doing it after 'read_startpem' to know exactly where perennial ice is (no matter if there is a "startpem.nc" or not). Moreover, when there is no "startpem.nc", location of perennial ice depends now on 'watercaptag' and on huge amount of frost. This prevents negative ice tendency while there is no ice which can happen with weird PCM inputs (i.e. 'watercaptag' = F & 'watercap' /= 0 & no "stratpem.nc").
• Few small safeguards throughout the code.

JBC

File:

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

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

@@ -27 +27 @@
 !=======================================================================
 SUBROUTINE load_xios_data(ps_avg,ps_ts,tsurf_avg,tsurf_avg_yr1,tsoil_avg,tsoil_ts,h2o_surfdensity_avg,h2o_soildensity_avg, &
-                          q_h2o_ts,q_co2_ts,d_h2oice,d_co2ice)
+                          q_h2o_ts,q_co2_ts,minPCM_h2operice,minPCM_co2perice,minPCM_h2ofrost,minPCM_co2frost)
 !-----------------------------------------------------------------------
 ! NAME
@@ -47 +47 @@
 use io_netcdf,  only: xios_day_name2, xios_yr_name1, xios_yr_name2, open_nc, close_nc, get_var_nc, get_dim_nc
 use soil,       only: do_soil
-use tendencies, only: compute_tend
 use frost,      only: compute_frost4PCM
 use stoppage,   only: stop_clean
@@ -60 +59 @@
 ! ---------
 real(dp), dimension(:),       intent(out) :: ps_avg
-real(dp), dimension(:,:),     intent(out) :: tsurf_avg, tsurf_avg_yr1, h2o_surfdensity_avg, d_h2oice, d_co2ice, ps_ts, q_h2o_ts, q_co2_ts
+real(dp), dimension(:,:),     intent(out) :: tsurf_avg, tsurf_avg_yr1, h2o_surfdensity_avg, ps_ts, q_h2o_ts, q_co2_ts
 real(dp), dimension(:,:,:),   intent(out) :: tsoil_avg, h2o_soildensity_avg
 real(dp), dimension(:,:,:,:), intent(out) :: tsoil_ts
+real(dp), dimension(:,:,:),   intent(out) :: minPCM_h2operice, minPCM_co2perice, minPCM_h2ofrost, minPCM_co2frost
 
 ! LOCAL VARIABLES
@@ -72 +72 @@
 real(dp), dimension(:,:,:,:), allocatable    :: var_read_4d
 character(:),                 allocatable    :: num ! To read slope variables
-real(dp), dimension(ngrid,nslope,2)          :: min_h2operice, min_co2perice, min_h2ofrost, min_co2frost
 real(dp), dimension(ngrid,nsoil,nslope,nday) :: h2o_soildensity_ts
 ! CODE
 ! ----
 ! Initialization
-min_h2operice(:,:,:) = 0._dp
-min_co2perice(:,:,:) = 0._dp
-min_h2ofrost(:,:,:) = 0._dp
-min_co2frost(:,:,:) = 0._dp
+minPCM_h2operice(:,:,:) = 0._dp
+minPCM_co2perice(:,:,:) = 0._dp
+minPCM_h2ofrost(:,:,:) = 0._dp
+minPCM_co2frost(:,:,:) = 0._dp
 if (nslope == 1) then ! No slope
     allocate(character(0) :: num)
@@ -103 +102 @@
         num = '_slope'//num
     end if
-    call get_var_nc('co2ice'//num,var_read_2d)          ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2frost(:,islope,1))
-    call get_var_nc('h2o_ice_s'//num,var_read_2d)       ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2ofrost(:,islope,1))
-    call get_var_nc('watercap'//num,var_read_2d)        ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2operice(:,islope,1))
-    call get_var_nc('perennial_co2ice'//num,var_read_2d); call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2perice(:,islope,1))
+    call get_var_nc('co2ice'//num,var_read_2d)          ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,minPCM_co2frost(:,islope,1))
+    call get_var_nc('h2o_ice_s'//num,var_read_2d)       ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,minPCM_h2ofrost(:,islope,1))
+    call get_var_nc('watercap'//num,var_read_2d)        ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,minPCM_h2operice(:,islope,1))
+    call get_var_nc('perennial_co2ice'//num,var_read_2d); call lonlat2vect(nlon,nlat,ngrid,var_read_2d,minPCM_co2perice(:,islope,1))
     call get_var_nc('tsurf'//num,var_read_2d)           ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,tsurf_avg_yr1(:,islope))
 end do
@@ -129 +128 @@
     end if
     call get_var_nc('tsurf'//num,var_read_2d)           ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,tsurf_avg(:,islope))
-    call get_var_nc('co2ice'//num,var_read_2d)          ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2frost(:,islope,2))
-    call get_var_nc('h2o_ice_s'//num,var_read_2d)       ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2ofrost(:,islope,2))
-    call get_var_nc('watercap'//num,var_read_2d)        ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_h2operice(:,islope,2))
-    call get_var_nc('perennial_co2ice'//num,var_read_2d); call lonlat2vect(nlon,nlat,ngrid,var_read_2d,min_co2perice(:,islope,2))
+    call get_var_nc('co2ice'//num,var_read_2d)          ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,minPCM_co2frost(:,islope,2))
+    call get_var_nc('h2o_ice_s'//num,var_read_2d)       ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,minPCM_h2ofrost(:,islope,2))
+    call get_var_nc('watercap'//num,var_read_2d)        ; call lonlat2vect(nlon,nlat,ngrid,var_read_2d,minPCM_h2operice(:,islope,2))
+    call get_var_nc('perennial_co2ice'//num,var_read_2d); call lonlat2vect(nlon,nlat,ngrid,var_read_2d,minPCM_co2perice(:,islope,2))
     if (do_soil) then
         call get_var_nc('soiltemp'//num,var_read_3d)
@@ -204 +203 @@
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 ! Compute frost from yearly minima
-call compute_frost4PCM(min_h2ofrost(:,:,2),min_co2frost(:,:,2))
-
-!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-! Compute ice tendencies from yearly minima
-call print_msg('> Computing surface ice tendencies')
-call compute_tend(min_h2operice + min_h2ofrost,d_h2oice)
-call print_msg('H2O ice tendencies [kg/m2/yr] (min|max): '//real2str(minval(d_h2oice))//' | '//real2str(maxval(d_h2oice)))
-call compute_tend(min_co2perice + min_co2frost,d_co2ice)
-call print_msg('CO2 ice tendencies [kg/m2/yr] (min|max): '//real2str(minval(d_co2ice))//' | '//real2str(maxval(d_co2ice)))
+call compute_frost4PCM(minPCM_h2ofrost(:,:,2),minPCM_co2frost(:,:,2))
 
 END SUBROUTINE load_xios_data