Changeset 3532 for trunk/LMDZ.COMMON/libf/evolution/pem.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/pem.F90 (modified) (6 diffs)

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

@@ -62 +62 @@
 use compute_tend_mod,           only: compute_tend
 use info_PEM_mod,               only: info_PEM
-use interpol_TI_PEM2PCM_mod,    only: interpol_TI_PEM2PCM
 use nb_time_step_PCM_mod,       only: nb_time_step_PCM
 use pemetat0_mod,               only: pemetat0
@@ -298 +297 @@
     call execute_command_line('squeue -j '//trim(job_id)//' -h --Format TimeLimit > tmp_cmdout.txt',cmdstat = cstat)
     if (cstat /= 0) then
-        call execute_command_line('qstat -f '//trim(job_id)//' | grep "Walltime" | awk ''{print $3}'' > tmp_cmdout.txt', cmdstat = cstat)
+        call execute_command_line('qstat -f '//trim(job_id)//' | grep "Walltime" | awk ''{print $3}'' > tmp_cmdout.txt',cmdstat = cstat)
         if (cstat > 0) then
             error stop 'pem: command execution failed!'
@@ -939 +938 @@
             do ig = 1,ngrid
                 do islope = 1,nslope
-                    call dyn_ss_ice_m(icetable_depth(ig,islope),tsurf_avg(ig,islope),tsoil_PEM(ig,:,islope),nsoilmx_PEM,TI_PEM(ig,1,nslope),ps(ig),sum(q_h2o_PEM_phys(ig,:))/size(q_h2o_PEM_phys,2),ice_porefilling(ig,:,islope),porefill,ssi_depth)
+                    call dyn_ss_ice_m(icetable_depth(ig,islope),tsurf_avg(ig,islope),tsoil_PEM(ig,:,islope),nsoilmx_PEM,TI_PEM(ig,1,nslope),ps(ig),(/sum(q_h2o_PEM_phys(ig,:))/size(q_h2o_PEM_phys,2)/),ice_porefilling(ig,:,islope),porefill,ssi_depth)
                     icetable_depth(ig,islope) = ssi_depth
                     ice_porefilling(ig,:,islope) = porefill
@@ -961 +960 @@
             do ig = 1,ngrid
                 do islope = 1,nslope
-                    do l = 1,nsoilmx_PEM 
+                    do l = 1,nsoilmx_PEM
                         if (l == 1) then
                             totmassco2_adsorbded = totmassco2_adsorbded + co2_adsorbded_phys(ig,l,islope)*(layer_PEM(l))* &
@@ -1008 +1007 @@
                 call writediagsoilpem(ngrid,'co2_ads_slope'//str2,'co2_ads','K',3,co2_adsorbded_phys(:,:,islope))
                 call writediagsoilpem(ngrid,'h2o_ads_slope'//str2,'h2o_ads','K',3,h2o_adsorbded_phys(:,:,islope))
-            endif                       
+            endif
         endif
     enddo
@@ -1101 +1100 @@
 ! III_a.2 Tsoil update (for startfi)
 if (soil_pem) then
-    call interpol_TI_PEM2PCM(ngrid,nslope,nsoilmx_PEM,nsoilmx,TI_PEM,inertiesoil)
+    inertiesoil = TI_PEM(:,:nsoilmx,:)
     tsoil = tsoil_PEM(:,1:nsoilmx,:) + tsoil_anom(:,:,:,timelen)
 #ifndef CPP_STD