Changeset 3070 for trunk/LMDZ.COMMON/libf/evolution/compute_tendencies_slope.F90

Timestamp:

Oct 5, 2023, 11:29:45 AM (14 months ago)

Author:

jbclement

Message:

PEM:
Correction of a bug: the variable 'g' was not correctly initialized in 1D. A little of code cleaning.
JBC

File:

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

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

@@ -1 @@
-!
-! $Id $
-!
-SUBROUTINE compute_tendencies_slope(ngrid,nslope,min_ice_Y1,&
-     min_ice_Y2,tendencies_ice)
+SUBROUTINE compute_tendencies_slope(ngrid,nslope,min_ice_Y1,min_ice_Y2,tendencies_ice)
 
-      IMPLICIT NONE
+implicit none
 
 !=======================================================================
@@ -15 +11 @@
 !   arguments:
 !   ----------
-
 !   INPUT
-
-     INTEGER, intent(in) :: ngrid, nslope                           ! # of grid points along longitude/latitude/ total
-     REAL, intent(in) , dimension(ngrid,nslope):: min_ice_Y1       ! LON x LAT field : minimum of water ice at each point for the first year
-     REAL, intent(in) , dimension(ngrid,nslope):: min_ice_Y2       ! LON x LAT field : minimum of water ice at each point for the second year
+integer,                       intent(in) :: ngrid, nslope ! # of grid points along longitude/latitude/ total
+real, dimension(ngrid,nslope), intent(in) :: min_ice_Y1    ! LON x LAT field : minimum of water ice at each point for the first year
+real, dimension(ngrid,nslope), intent(in) :: min_ice_Y2    ! LON x LAT field : minimum of water ice at each point for the second year
 
 !   OUTPUT
-     REAL, intent(out) , dimension(ngrid,nslope)   :: tendencies_ice            ! physical point field : difference between the minima = evolution of perenial ice
-
-!   local:
-!   ------
-     INTEGER :: ig,islope                                                           ! loop variable
+real, dimension(ngrid,nslope), intent(out) :: tendencies_ice ! physical point field : difference between the minima = evolution of perenial ice
 
 !=======================================================================
 
 !  We compute the difference
-
-  DO ig=1,ngrid
-    DO islope = 1, nslope
-      tendencies_ice(ig,islope)=min_ice_Y2(ig,islope)-min_ice_Y1(ig,islope)
-    enddo
-  ENDDO
+tendencies_ice = min_ice_Y2 - min_ice_Y1
 
 !  If the difference is too small; there is no evolution
-  DO ig=1,ngrid
-    DO islope = 1, nslope
-      if(abs(tendencies_ice(ig,islope)).LT.1.0E-10) then
-        tendencies_ice(ig,islope)=0.
-      endif
-    enddo
-  ENDDO
+where (abs(tendencies_ice) < 1.e-10) tendencies_ice = 0.
 
 END SUBROUTINE compute_tendencies_slope