Changeset 2897 for trunk/LMDZ.COMMON/libf/evolution/compute_tendencies_mod_slope.F90

Timestamp:

Feb 16, 2023, 5:29:48 PM (22 months ago)

Author:

romain.vande

Message:

Mars PEM:
Deep cleaning of variables name and allocate.
All the "dyn to phys" grid change is done in subroutines and not in the main program.

File:

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

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

@@ -2 +2 @@
 ! $Id $
 !
-SUBROUTINE compute_tendencies_slope(tendencies_h2o_ice,min_h2o_ice_Y1,&
-     min_h2o_ice_Y2,iim_input,jjm_input,ngrid,tendencies_h2o_ice_phys,nslope)
+SUBROUTINE compute_tendencies_slope(ngrid,nslope,min_ice_Y1,&
+     min_ice_Y2,tendencies_ice)
 
       IMPLICIT NONE
@@ -18 +18 @@
 !   INPUT
 
-     INTEGER, intent(in) :: iim_input,jjm_input,ngrid  ,nslope                           ! # of grid points along longitude/latitude/ total
-     REAL, intent(in) , dimension(iim_input+1,jjm_input+1,nslope):: min_h2o_ice_Y1       ! LON x LAT field : minimum of water ice at each point for the first year
-     REAL, intent(in) , dimension(iim_input+1,jjm_input+1,nslope):: min_h2o_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, 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
 
 !   OUTPUT
-     REAL, intent(out) , dimension(iim_input+1,jjm_input+1,nslope) :: tendencies_h2o_ice ! LON x LAT field : difference between the minima = evolution of perenial ice
-     REAL, intent(out) , dimension(ngrid,nslope)   :: tendencies_h2o_ice_phys            ! physical point field : difference between the minima = evolution of perenial ice
+     REAL, intent(out) , dimension(ngrid,nslope)   :: tendencies_ice            ! physical point field : difference between the minima = evolution of perenial ice
 
 !   local:
 !   ------
-
-     INTEGER :: i,j,ig0,islope                                                           ! loop variable
+     INTEGER :: ig,islope                                                           ! loop variable
 
 !=======================================================================
@@ -35 +33 @@
 !  We compute the difference
 
-  DO j=1,jjm_input+1
-    DO i = 1, iim_input
-       DO islope = 1, nslope
-         tendencies_h2o_ice(i,j,islope)=min_h2o_ice_Y2(i,j,islope)-min_h2o_ice_Y1(i,j,islope)
-       enddo
-    ENDDO
+  DO ig=1,ngrid
+    DO islope = 1, nslope
+      tendencies_ice(ig,islope)=min_ice_Y2(ig,islope)-min_ice_Y1(ig,islope)
+    enddo
   ENDDO
 
 !  If the difference is too small; there is no evolution
-  DO j=1,jjm_input+1
-    DO i = 1, iim_input
-       DO islope = 1, nslope
-         if(abs(tendencies_h2o_ice(i,j,islope)).LT.1.0E-10) then
-            tendencies_h2o_ice(i,j,islope)=0.
-         endif
-       enddo
-    ENDDO
-  ENDDO
-
-  DO islope = 1,nslope
-    CALL gr_dyn_fi(1,iim_input+1,jjm_input+1,ngrid,tendencies_h2o_ice(:,:,islope),tendencies_h2o_ice_phys(:,islope))
+  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