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

Timestamp:

Dec 11, 2025, 12:56:05 PM (5 weeks ago)

Author:

jbclement

Message:

PEM:
Massive structural refactor of the PEM codebase for improved readability, consistency and maintainability. The goal is to modernize, standardize and consolidate the code while removing legacy complexity. In detail, this change:

• Performs large-scale cleanup removing unused code, obsolete routines, duplicated functionality and deprecated initialization logic;
• Removes "*_mod" wrappers;
• Replaces mixed naming conventions with clearer variable names, domain-based file/module names and purpose-based routine names;
• Adds native reading/writing capabilities to the PEM removing the cumbersome dependency on Mars PCM subroutines;
• Centralizes soil/ice/orbital/PEM configuration logic into dedicated modules;
• Simplifies routines and updates calls/interfaces to match the new structure.

This refactor significantly clarifies the codebase and provides a cleaner foundation for forthcoming developments.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/xios_data.F90 (moved) (moved from trunk/LMDZ.COMMON/libf/evolution/read_XIOS_data.F90) (10 diffs)

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

@@ -1 @@
-MODULE read_XIOS_data
+MODULE xios_data
 
 use netcdf, only: nf90_open, nf90_close, nf90_inquire_dimension, nf90_inq_dimid, nf90_noerr, nf90_nowrite, nf90_get_var, nf90_inq_varid
@@ -19 +19 @@
 contains
 !=======================================================================
-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, &
+
+!=======================================================================
+SUBROUTINE load_xios_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
+use grid_conversion, only: lonlat2vect
+use soil,            only: do_soil
+use tendencies,      only: compute_tend
+use metamorphism,    only: compute_frost
 
 implicit none
@@ -109 +111 @@
     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
+    if (do_soil) then
         call get_var('soiltemp'//num,var_read_3d)
         do isoil = 1,nsoil_PCM
@@ -141 +143 @@
     call lonlat2vect(nlon,nlat,ngrid,var_read_3d(:,:,isol),q_co2_ts(:,isol))
 enddo
-if (soil_pem) then
+if (do_soil) then
     do islope = 1,nslope
         if (nslope /= 1) then
@@ -180 +182 @@
 write(*,*) 'CO2 ice tendencies (min/max):', minval(d_co2ice), maxval(d_co2ice)
 
-END SUBROUTINE read_PCM_data
+END SUBROUTINE load_xios_data
+!=======================================================================
+
+!=======================================================================
+SUBROUTINE get_timelen(filename,timelen)
+
+use netcdf
+
+implicit none
+
+! Arguments
+! ---------
+character(*), intent(in)  :: filename
+integer,      intent(out) :: timelen
+
+! Local variables
+! ---------------
+integer :: ncid  ! File ID
+integer :: dimid ! Dimension ID
+integer :: ierr  ! Return codes
+
+! Code
+! ----
+! Open the NetCDF file
+ierr = nf90_open(trim(filename),NF90_NOWRITE,ncid)
+if (ierr /= nf90_noerr) then
+    write(*,*) "Error opening file:", trim(nf90_strerror(ierr))
+    error stop
+endif
+
+! Get the dimension ID for 'time_counter'
+ierr = nf90_inq_dimid(ncid,"time_counter",dimid)
+if (ierr /= nf90_noerr) then
+    write(*,*) "Error getting dimid 'time_counter':", trim(nf90_strerror(ierr))
+    error stop
+endif
+
+! Get the size of the dimension 'time_counter'
+ierr = nf90_inquire_dimension(ncid,dimid,len = timelen)
+if (ierr /= nf90_noerr) then
+    write(*,*) "Error getting dimension length:", trim(nf90_strerror(ierr))
+    error stop
+endif
+
+! Close the file
+ierr = nf90_close(ncid)
+if (ierr /= nf90_noerr) then
+    write(*,*) "Error closing file:", trim(nf90_strerror(ierr))
+    error stop
+endif
+
+END SUBROUTINE get_timelen
+!=======================================================================
 
 !=======================================================================
@@ -205 +259 @@
 
 END SUBROUTINE error_msg
+!=======================================================================
 
 !=======================================================================
@@ -218 +273 @@
 
 END SUBROUTINE get_var_1d
+!=======================================================================
 
 !=======================================================================
@@ -231 +287 @@
 
 END SUBROUTINE get_var_2d
+!=======================================================================
 
 !=======================================================================
@@ -244 +301 @@
 
 END SUBROUTINE get_var_3d
+!=======================================================================
 
 !=======================================================================
@@ -257 +315 @@
 
 END SUBROUTINE get_var_4d
-
-END MODULE read_XIOS_data
+!=======================================================================
+
+END MODULE xios_data