Changeset 3189 for trunk/LMDZ.COMMON/libf/evolution
- Timestamp:
- Jan 29, 2024, 11:14:07 PM (11 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/evolution/pem.F90
r3188 r3189 202 202 real, dimension(:,:), allocatable :: tsurf_ave ! Physic x SLOPE field: Averaged Surface Temperature [K] 203 203 real, dimension(:,:,:), allocatable :: tsoil_ave ! Physic x SOIL x SLOPE field: Averaged Soil Temperature [K] 204 real, dimension(:,:,:), allocatable :: tsoil_anom ! Amplitude between instataneous and yearly average soil temperature [K] 204 205 real, dimension(:,:,:), allocatable :: tsurf_PCM_timeseries ! ngrid x SLOPE XTULES field: Surface Temperature in timeseries [K] 205 206 real, dimension(:,:,:,:), allocatable :: tsoil_phys_PEM_timeseries ! IG x SLOPE XTULES field: Non averaged Soil Temperature [K] … … 476 477 477 478 allocate(tsoil_ave(ngrid,nsoilmx,nslope)) 479 allocate(tsoil_anom(ngrid,nsoilmx,nslope)) 478 480 allocate(watersoil_density_PEM_ave(ngrid,nsoilmx_PEM,nslope)) 479 481 allocate(vmr_co2_PCM(ngrid,timelen)) … … 523 525 call ini_ice_table_porefilling(ngrid,nslope) 524 526 527 tsoil_anom(:,:,:) = tsoil(:,:,:) - tsoil_ave(:,:,:) ! compute anomaly between Tsoil(t) in the startfi - <Tsoil> to recompute properly tsoil in the restart 525 528 if (soil_pem) then 526 529 call soil_settings_PEM(ngrid,nslope,nsoilmx_PEM,nsoilmx,inertiesoil,TI_PEM) … … 985 988 if (soil_pem) then 986 989 call interpol_TI_PEM2PCM(ngrid,nslope,nsoilmx_PEM,nsoilmx,TI_PEM,inertiesoil) 987 tsoil = tsoil_ phys_PEM_timeseries(:,1:nsoilmx,:,timelen)990 tsoil = tsoil_PEM(:,1:nsoilmx,:) + tsoil_anom 988 991 #ifndef CPP_STD 989 992 flux_geo = fluxgeo
Note: See TracChangeset
for help on using the changeset viewer.