Changeset 4140 for trunk/LMDZ.COMMON/libf/evolution/tendencies.F90

Timestamp:

Mar 18, 2026, 1:56:02 PM (2 weeks ago)

Author:

jbclement

Message:

PEM:
Fix to avoid building large temporary arrays in the arguments evaluation of "compute_tendice" which caused memory to stall.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/tendencies.F90 (modified) (3 diffs)

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

@@ -28 +28 @@
 
 !=======================================================================
-SUBROUTINE compute_tendice(min_ice,is_perice,d_ice)
+SUBROUTINE compute_tendice(min_perice,min_frost,perice,d_ice)
 !-----------------------------------------------------------------------
 ! NAME
@@ -51 +51 @@
 ! ARGUMENTS
 ! ---------
-real(dp),    dimension(:,:,:), intent(in)  :: min_ice
-logical(k4), dimension(:,:),   intent(in)  :: is_perice
-real(dp),    dimension(:,:),   intent(out) :: d_ice
+real(dp), dimension(:,:,:), intent(in)  :: min_perice
+real(dp), dimension(:,:,:), intent(in)  :: min_frost
+real(dp), dimension(:,:),   intent(in)  :: perice
+real(dp), dimension(:,:),   intent(out) :: d_ice
 
 ! CODE
 ! ----
 ! We compute the difference to get the tendency
-d_ice(:,:) = min_ice(:,:,2) - min_ice(:,:,1)
+d_ice(:,:) = min_perice(:,:,2) + min_frost(:,:,2) - (min_perice(:,:,1) + min_frost(:,:,1))
 
 ! If the difference is too small, then there is no evolution
@@ -64 +65 @@
 
 ! If the tendency is negative but there is no ice reservoir for the PEM
-where (abs(d_ice) < 0._dp .and. .not. is_perice) d_ice = 0._dp
+where (d_ice < 0._dp .and. abs(perice) < minieps) d_ice = 0._dp
 
 END SUBROUTINE compute_tendice