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

Timestamp:

Jan 29, 2024, 11:14:07 PM (10 months ago)

Author:

llange

Message:

PEM
Fixing bug when recomputing Tsoil for the startfi. It is now done with: Tsoil averaged + Delta T where delta T is tthe difference between the instantaneous soil temperature and the yearly averaged soil temperature in the original startfi.
LL

File:

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

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

@@ -202 +202 @@
 real, dimension(:,:),     allocatable :: tsurf_ave                          ! Physic x SLOPE field: Averaged Surface Temperature [K]
 real, dimension(:,:,:),   allocatable :: tsoil_ave                          ! Physic x SOIL x SLOPE field: Averaged Soil Temperature [K]
+real, dimension(:,:,:),   allocatable :: tsoil_anom                         ! Amplitude between instataneous and yearly average soil temperature [K]
 real, dimension(:,:,:),   allocatable :: tsurf_PCM_timeseries               ! ngrid x SLOPE XTULES field: Surface Temperature in timeseries [K]
 real, dimension(:,:,:,:), allocatable :: tsoil_phys_PEM_timeseries          ! IG x SLOPE XTULES field: Non averaged Soil Temperature [K]
@@ -476 +477 @@
 
 allocate(tsoil_ave(ngrid,nsoilmx,nslope))
+allocate(tsoil_anom(ngrid,nsoilmx,nslope))
 allocate(watersoil_density_PEM_ave(ngrid,nsoilmx_PEM,nslope))
 allocate(vmr_co2_PCM(ngrid,timelen))
@@ -523 +525 @@
 call ini_ice_table_porefilling(ngrid,nslope)
 
+tsoil_anom(:,:,:) = tsoil(:,:,:) - tsoil_ave(:,:,:) ! compute anomaly between Tsoil(t) in the startfi - <Tsoil> to recompute properly tsoil in the restart
 if (soil_pem) then
     call soil_settings_PEM(ngrid,nslope,nsoilmx_PEM,nsoilmx,inertiesoil,TI_PEM)
@@ -985 +988 @@
 if (soil_pem) then
     call interpol_TI_PEM2PCM(ngrid,nslope,nsoilmx_PEM,nsoilmx,TI_PEM,inertiesoil)
-    tsoil = tsoil_phys_PEM_timeseries(:,1:nsoilmx,:,timelen)
+    tsoil = tsoil_PEM(:,1:nsoilmx,:) + tsoil_anom
 #ifndef CPP_STD
     flux_geo = fluxgeo