Changeset 3498 for trunk/LMDZ.COMMON/libf/evolution/compute_tend_mod.F90

Timestamp:

Nov 7, 2024, 2:48:08 PM (3 weeks ago)

Author:

jbclement

Message:

PEM:

• Correction of the variable name for the ice table depth in "pemetat0.F90". So it is now got as intended from the "startpem.nc" file;
• Renaming of the tendencies in the PEM with the prefix 'd_' instead of 'tend_';
• Modification of the PEM time step type from integer to real. As a consequence, all time variables are now of real type. This change adds the possibility to consider fractions of year as time step.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/compute_tend_mod.F90 (modified) (2 diffs)

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

@@ -7 +7 @@
 !=======================================================================
 
-SUBROUTINE compute_tend(ngrid,nslope,min_ice,tendencies_ice)
+SUBROUTINE compute_tend(ngrid,nslope,min_ice,d_ice)
 
 implicit none
@@ -25 +25 @@
 
 !   OUTPUT
-real, dimension(ngrid,nslope), intent(out) :: tendencies_ice ! Difference between the minima = evolution of perennial ice
+real, dimension(ngrid,nslope), intent(out) :: d_ice ! Difference between the minima = evolution of perennial ice
 !=======================================================================
 ! We compute the difference
-tendencies_ice = min_ice(:,:,2) - min_ice(:,:,1)
+d_ice = min_ice(:,:,2) - min_ice(:,:,1)
 
 ! If the difference is too small, then there is no evolution
-where (abs(tendencies_ice) < 1.e-10) tendencies_ice = 0.
+where (abs(d_ice) < 1.e-10) d_ice = 0.
 
 ! If the minimum over the last year is 0, then we have no perennial ice
-where (abs(min_ice(:,:,2)) < 1.e-10) tendencies_ice = 0.
+where (abs(min_ice(:,:,2)) < 1.e-10) d_ice = 0.
 
 END SUBROUTINE compute_tend