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

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

@@ -32 +32 @@
   real(8) equildepth  ! =zdepthE
   !real(8), external :: zint  ! defined in allinterfaces.mod
-  
+
   typeE = -9; equildepth = -9999.
-  do i=1,nz 
+  do i=1,nz
      if (rhosatav(i) <= avrho1) then
         typeE=i
@@ -67 +67 @@
   integer, intent(OUT) :: typeF  ! index of depth below which filling occurs
   real(8), intent(INOUT) :: zdepthF
-  real(8), intent(IN) :: B 
+  real(8), intent(IN) :: B
   real(8), intent(OUT) :: ypp(nz), zdepthG
   integer, intent(INOUT) :: typeG  ! positive on input when Fgeotherm>0
@@ -112 +112 @@
 
 !-depth to shallowest perennial ice
-  typeP = -9 
+  typeP = -9
   do i=1,nz
      if (porefill(i)>0.) then
@@ -157 +157 @@
   endif
   if (typeG>0 .and. typeT<0) then
-     call colint(porefill(:)/eta(:),z,nz,typeG-1,nz,cumfillabove) 
+     call colint(porefill(:)/eta(:),z,nz,typeG-1,nz,cumfillabove)
      newtypeG = -9
      do i=typeG,nz
@@ -178 +178 @@
      if (newtypeG>0) typeG=newtypeG
   end if
-  ! if typeG>0, then all ice at and below typeG should be erased 
+  ! if typeG>0, then all ice at and below typeG should be erased
 end subroutine depths_avmeth
 
@@ -204 +204 @@
 
   ! advance ice table, avdrho>0 is retreat
-  if (zdepthT>=0. .and. avdrho>0.) then 
+  if (zdepthT>=0. .and. avdrho>0.) then
      typeP=-9999; typeT=-9999
      do j=1,nz
@@ -226 +226 @@
   endif
   if (zdepthT>z(nz)) zdepthT=-9999.
-  
+
   ! advance interface, avdrhoP>0 is loss from zdepthP
   if (avdrhoP>0.) then
@@ -232 +232 @@
      do j=1,nz
         if (typeF>0 .and. j>=typeF) exit ! don't retreat beyond typeF
-        if (zdepthT>=0. .and. z(j)>zdepthT) exit 
+        if (zdepthT>=0. .and. z(j)>zdepthT) exit
         call colint(porefill(1:nz)*z(1:nz),z(1:nz),nz,1,j,integ)
         erase = j
@@ -241 +241 @@
 
   ! new depth
-  newtypeP = -9 
+  newtypeP = -9
   do j=1,nz
      if (zdepthT>=0. .and. z(j)>zdepthT) exit
@@ -253 +253 @@
   ub = typeF
   if (newtypeP>0 .and. typeF>0 .and. newtypeP<ub) ub=newtypeP
-  if (ub>0) then  
+  if (ub>0) then
      do j=ub,nz
         porefill(j) = porefill(j) + B*ypp(j)