Changeset 3532 for trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

Timestamp:

Dec 4, 2024, 4:04:54 PM (3 weeks ago)

Author:

jbclement

Message:

PEM:
Removing unecessary module/subroutine "interpol_TI_PEM2PCM.F90" + Few small corrections/cleanings.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90 (modified) (9 diffs)

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

@@ -45 +45 @@
 real,                           intent(in) :: global_avg_pressure ! mean average pressure on the planet [Pa]
 real, dimension(ngrid,nslope),  intent(in) :: tsurf_avg_yr1       ! surface temperature at the first year of PCM call [K]
-real, dimension(ngrid,nslope),  intent(in) :: tsurf_avg_yr2       ! surface temperature at the second  year of PCM call [K]
+real, dimension(ngrid,nslope),  intent(in) :: tsurf_avg_yr2       ! surface temperature at the second year of PCM call [K]
 real, dimension(ngrid,timelen), intent(in) :: q_co2               ! MMR tracer co2 [kg/kg]
 real, dimension(ngrid,timelen), intent(in) :: q_h2o               ! MMR tracer h2o [kg/kg]
@@ -281 +281 @@
                 call compute_tsoil_pem(ngrid,nsoil_PEM,.true.,TI_PEM(:,:,islope),timestep,tsurf_avg_yr2(:,islope),tsoil_PEM(:,:,islope))
             else
-! predictor corrector: restart from year 1 of the PCM and build the evolution of
-! tsoil at depth
+! predictor corrector: restart from year 1 of the PCM and build the evolution of tsoil at depth
                 tsoil_tmp_yr1(:,:,islope) = tsoil_startPEM(:,:,islope)
                 call compute_tsoil_pem(ngrid,nsoil_PEM,.true.,TI_PEM(:,:,islope),timestep,tsurf_avg_yr1(:,islope),tsoil_tmp_yr1(:,:,islope))
@@ -289 +288 @@
                 call compute_tsoil_pem(ngrid,nsoil_PEM,.false.,TI_PEM(:,:,islope),timestep,tsurf_avg_yr2(:,islope),tsoil_tmp_yr2(:,:,islope))
 
-                do iloop = nsoil_PCM+1,nsoil_PEM
+                do iloop = nsoil_PCM + 1,nsoil_PEM
                     tsoil_PEM(:,iloop,islope) = tsoil_tmp_yr2(:,iloop,islope)
                 enddo
@@ -295 +294 @@
 
             do it = 1,timelen
-                do isoil = nsoil_PCM+1,nsoil_PEM
-                    tsoil_inst(:,isoil,islope,it) = tsoil_PEM(:,isoil,islope)
-                enddo
-            enddo
-            do isoil = nsoil_PCM+1,nsoil_PEM
-                do ig = 1,ngrid
-                    watersoil_avg(ig,isoil,islope) = exp(beta_clap_h2o/tsoil_PEM(ig,isoil,islope) + alpha_clap_h2o)/tsoil_PEM(ig,isoil,islope)*mmol(igcm_h2o_vap)/(mugaz*r)
-                enddo
+                tsoil_inst(:,nsoil_PCM + 1:nsoil_PEM,islope,it) = tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope)
+            enddo
+            do isoil = nsoil_PCM + 1,nsoil_PEM
+                watersoil_avg(:,isoil,islope) = exp(beta_clap_h2o/tsoil_PEM(:,isoil,islope) + alpha_clap_h2o)/tsoil_PEM(:,isoil,islope)*mmol(igcm_h2o_vap)/(mugaz*r)
             enddo
         enddo ! islope
@@ -314 +309 @@
                 write(*,*)'PEM settings: failed loading <ice_table_depth>'
                 write(*,*)'will reconstruct the values of the ice table given the current state'
-                call computeice_table_equilibrium(ngrid,nslope,nsoil_PEM,watercaptag,watersurf_avg,watersoil_avg, TI_PEM(:,1,:),ice_table_depth,ice_table_thickness)
+                call computeice_table_equilibrium(ngrid,nslope,nsoil_PEM,watercaptag,watersurf_avg,watersoil_avg,TI_PEM(:,1,:),ice_table_depth,ice_table_thickness)
                 call update_soil_thermalproperties(ngrid,nslope,nsoil_PEM,d_h2oice,h2o_ice,global_avg_pressure,ice_table_depth,ice_table_thickness,ice_porefilling,icetable_equilibrium,icetable_dynamic,TI_PEM)
                 do islope = 1,nslope
@@ -322 +317 @@
             write(*,*) 'PEMETAT0: ICE TABLE done'
         else if (icetable_dynamic) then
+            call get_field("ice_porefilling",ice_porefilling,found)
+            if (.not. found) then
+                write(*,*)'PEM settings: failed loading <ice_porefilling>'
+                ice_porefilling = 0.
+            endif
             call get_field("ice_table_depth",ice_table_depth,found)
             if (.not. found) then
@@ -331 +331 @@
                     call ini_tsoil_pem(ngrid,nsoil_PEM,TI_PEM(:,:,islope),tsurf_avg_yr2(:,islope),tsoil_PEM(:,:,islope))
                 enddo
-            endif
-            call get_field("ice_porefilling",ice_porefilling,found)
-            if (.not. found) then
-                write(*,*)'PEM settings: failed loading <ice_porefilling>'
-                ice_porefilling = 1.
             endif
             write(*,*) 'PEMETAT0: ICE TABLE done'
@@ -510 +505 @@
             write(*,*) 'PEMETAT0: Ice table done'
         else if (icetable_dynamic) then
-            ice_porefilling = 1.
+            ice_porefilling = 0.
             ice_table_depth = -9999.
             call update_soil_thermalproperties(ngrid,nslope,nsoil_PEM,d_h2oice,h2o_ice,global_avg_pressure,ice_table_depth,ice_table_thickness,ice_porefilling,icetable_equilibrium,icetable_dynamic,TI_PEM)
@@ -518 +513 @@
             write(*,*) 'PEMETAT0: Ice table done'
         endif
-        
+
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !d) Regolith adsorbed