Changeset 4135 for trunk/LMDZ.COMMON/libf/evolution/backup.F90

Timestamp:

Mar 16, 2026, 6:09:43 PM (10 days ago)

Author:

jbclement

Message:

PEM:
Relocate Mars-specific parameters out of "planet" module and into the modules that own their physics domain.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/backup.F90 (modified) (4 diffs)

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

@@ -114 +114 @@
 ! LOCAL VARIABLES
 ! ---------------
-real(dp),    dimension(:,:),   allocatable :: h2o_ice4PCM, co2_ice4PCM, tsurf4PCM, flux_geo4PCM, albedo4PCM, emissivity4PCM
-real(dp),    dimension(:,:),   allocatable :: teta4PCM, air_mass4PCM
-real(dp),    dimension(:,:,:), allocatable :: tsoil4PCM, inertiesoil4PCM, q4PCM
-real(dp),    dimension(:),     allocatable :: ps4PCM
-real(dp)                                   :: pa4PCM, preff4PCM
-logical(k4), dimension(:),     allocatable :: is_h2o_perice
+real(dp),    dimension(ngrid,nslope)           :: h2o_ice4PCM, co2_ice4PCM, tsurf4PCM, flux_geo4PCM, albedo4PCM, emissivity4PCM
+real(dp),    dimension(ngrid,nlayer)           :: teta4PCM, air_mass4PCM
+real(dp),    dimension(ngrid,nsoil_PCM,nslope) :: tsoil4PCM, inertiesoil4PCM
+real(dp),    dimension(ngrid,nlayer,nq)        :: q4PCM
+real(dp),    dimension(ngrid)                  :: ps4PCM
+logical(k4), dimension(ngrid)                  :: is_h2o_perice
+real(dp)                                       :: pa4PCM, preff4PCM
 
 ! CODE
 ! ----
 ! Build ice for the PCM
-allocate(h2o_ice4PCM(ngrid,nslope),co2_ice4PCM(ngrid,nslope),is_h2o_perice(ngrid))
 call build4PCM_perice(h2o_ice,co2_ice,is_h2o_perice,h2o_ice4PCM,co2_ice4PCM)
 
 ! Build surface temperature for the PCM
-allocate(tsurf4PCM(ngrid,nslope))
 call build4PCM_tsurf(tsurf_avg,tsurf_dev,tsurf4PCM)
 
 ! Build soil for the PCM
-if (do_soil) then
-    allocate(tsoil4PCM(ngrid,nsoil_PCM,nslope),inertiesoil4PCM(ngrid,nsoil_PCM,nslope),flux_geo4PCM(ngrid,nslope))
-    call build4PCM_soil(tsoil_avg,tsoil_dev,inertiesoil4PCM,tsoil4PCM,flux_geo4PCM)
-end if
+if (do_soil) call build4PCM_soil(tsoil_avg,tsoil_dev,inertiesoil4PCM,tsoil4PCM,flux_geo4PCM)
 
 ! Build atmosphere for the PCM
-allocate(ps4PCM(ngrid),teta4PCM(ngrid,nlayer),air_mass4PCM(ngrid,nlayer))
 call build4PCM_atmosphere(ps_avg,ps_dev,ps_avg_glob,ps_avg_glob_ini,ps4PCM,pa4PCM,preff4PCM,teta4PCM,air_mass4PCM)
 
 ! Build tracers for the PCM
-allocate(q4PCM(ngrid,nlayer,nq))
 call build4PCM_tracers(ps4PCM,q4PCM)
 
 ! Build surface radiative properties state for the PCM
-allocate(albedo4PCM(ngrid,nslope),emissivity4PCM(ngrid,nslope))
 call build4PCM_surf_rad_prop(h2o_ice,co2_ice,albedo4PCM,emissivity4PCM)
 
@@ -158 +151 @@
 end if
 
-! Deallocation
-if (allocated(emissivity4PCM)) deallocate(emissivity4PCM)
-if (allocated(albedo4PCM)) deallocate(albedo4PCM)
-if (allocated(q4PCM)) deallocate(q4PCM)
-if (allocated(ps4PCM)) deallocate(ps4PCM)
-if (allocated(teta4PCM)) deallocate(teta4PCM)
-if (allocated(air_mass4PCM)) deallocate(air_mass4PCM)
-if (allocated(tsoil4PCM)) deallocate(tsoil4PCM)
-if (allocated(inertiesoil4PCM)) deallocate(inertiesoil4PCM)
-if (allocated(flux_geo4PCM)) deallocate(flux_geo4PCM)
-if (allocated(tsurf4PCM)) deallocate(tsurf4PCM)
-if (allocated(co2_ice4PCM)) deallocate(co2_ice4PCM)
-if (allocated(h2o_ice4PCM)) deallocate(h2o_ice4PCM)
-if (allocated(is_h2o_perice)) deallocate(is_h2o_perice)
-
 END SUBROUTINE save_clim_state
 !=======================================================================
 
 !=======================================================================
-SUBROUTINE backup_restarts(idt)
+SUBROUTINE backup_restarts(backup_idt)
 !-----------------------------------------------------------------------
 ! NAME
@@ -205 +183 @@
 ! ARGUMENTS
 ! ---------
-integer(di), intent(in) :: idt
+integer(di), intent(in) :: backup_idt
 
 ! LOCAL VARIABLES
@@ -213 +191 @@
 ! CODE
 ! ----
-suffix = '_ts'//int2str(idt)
-
-call print_msg('> Backup at dt = '//int2str(idt),LVL_NFO)
+suffix = '_ts'//int2str(backup_idt)
+
+call print_msg('> Backup of "restart" files at dt = '//int2str(backup_idt),LVL_NFO)
 call copy_restart_if_present('re'//startevo_name,suffix2filename(startevo_name,suffix))
 call copy_restart_if_present('re'//startfi_name,suffix2filename(startfi_name,suffix))