Changeset 4136 for trunk/LMDZ.COMMON/libf/evolution/geometry.F90

Timestamp:

Mar 17, 2026, 11:12:37 AM (12 days ago)

Author:

jbclement

Message:

PEM:

• Temporary memory load reduced in "xios_data".
• Rename local variables in procedures to avoid masking variables in parent scope.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/geometry.F90 (modified) (16 diffs)

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

@@ -153 +153 @@
 
 !=======================================================================
-SUBROUTINE lonlat2vect(nlon,nlat,ngrid,v_ll,v_vect)
+SUBROUTINE lonlat2vect(nlon_in,nlat_in,ngrid_in,v_ll,v_vect)
 !-----------------------------------------------------------------------
 ! NAME
@@ -180 +180 @@
 ! ARGUMENTS
 !----------
-integer(di),                    intent(in)  :: nlon, nlat, ngrid
-real(dp), dimension(nlon,nlat), intent(in)  :: v_ll
-real(dp), dimension(ngrid),     intent(out) :: v_vect
+integer(di),                          intent(in)  :: nlon_in, nlat_in, ngrid_in
+real(dp), dimension(nlon_in,nlat_in), intent(in)  :: v_ll
+real(dp), dimension(ngrid_in),        intent(out) :: v_vect
 
 ! LOCAL VARIABLES
@@ -190 +190 @@
 ! CODE
 !-----
-if (ngrid == 1) then ! 1D case
+if (ngrid_in == 1) then ! 1D case
     v_vect(1) = v_ll(1,1)
     return
 else ! 3D
     ! Check
-    if (ngrid /= nlon*(nlat - 2) + 2) call stop_clean(__FILE__,__LINE__,'problem of dimensions!',1)
+    if (ngrid_in /= nlon_in*(nlat_in - 2) + 2) call stop_clean(__FILE__,__LINE__,'problem of dimensions!',1)
 
     ! Initialization
@@ -202 +202 @@
     ! Treatment of the poles
     v_vect(1) = v_ll(1,1)
-    v_vect(ngrid) = v_ll(1,nlat)
+    v_vect(ngrid_in) = v_ll(1,nlat_in)
 
     ! Treatment of regular points
-    do j = 2,nlat - 1
-        do i = 1,nlon
-            v_vect(1 + i + (j - 2)*nlon) = v_ll(i,j)
+    do j = 2,nlat_in - 1
+        do i = 1,nlon_in
+            v_vect(1 + i + (j - 2)*nlon_in) = v_ll(i,j)
         end do
     end do
@@ -216 +216 @@
 
 !=======================================================================
-SUBROUTINE vect2lonlat(nlon,nlat,ngrid,v_vect,v_ll)
+SUBROUTINE vect2lonlat(nlon_in,nlat_in,ngrid_in,v_vect,v_ll)
 !-----------------------------------------------------------------------
 ! NAME
@@ -243 +243 @@
 ! ARGUMENTS
 !----------
-integer(di),                    intent(in)  :: nlon, nlat, ngrid
-real(dp), dimension(ngrid),     intent(in)  :: v_vect
-real(dp), dimension(nlon,nlat), intent(out) :: v_ll
+integer(di),                          intent(in)  :: nlon_in, nlat_in, ngrid_in
+real(dp), dimension(ngrid_in),        intent(in)  :: v_vect
+real(dp), dimension(nlon_in,nlat_in), intent(out) :: v_ll
 
 ! LOCAL VARIABLES
@@ -253 +253 @@
 ! CODE
 !-----
-if (ngrid == 1) then ! 1D case
+if (ngrid_in == 1) then ! 1D case
     v_ll(1,1) = v_vect(1)
     return
 else ! 3D
     ! Check
-    if (ngrid /= nlon*(nlat - 2) + 2) call stop_clean(__FILE__,__LINE__,'problem of dimensions!',1)
+    if (ngrid_in /= nlon_in*(nlat_in - 2) + 2) call stop_clean(__FILE__,__LINE__,'problem of dimensions!',1)
 
     ! Initialization
@@ -265 +265 @@
     ! Treatment of the poles
     v_ll(:,1) = v_vect(1)
-    v_ll(:,nlat) = v_vect(ngrid)
+    v_ll(:,nlat_in) = v_vect(ngrid_in)
 
     ! Treatment of regular points
-    do j = 2,nlat - 1
-        do i = 1,nlon
-            v_ll(i,j) = v_vect(1 + i + (j - 2)*nlon)
+    do j = 2,nlat_in - 1
+        do i = 1,nlon_in
+            v_ll(i,j) = v_vect(1 + i + (j - 2)*nlon_in)
         end do
     end do
@@ -279 +279 @@
 
 !=======================================================================
-SUBROUTINE dyngrd2vect(ni,nj,ngrid,v_dyn,v_vect)
+SUBROUTINE dyngrd2vect(ni,nj,ngrid_in,v_dyn,v_vect)
 !-----------------------------------------------------------------------
 ! NAME
@@ -306 +306 @@
 ! ARGUMENTS
 !----------
-integer(di),                intent(in)  :: ni, nj, ngrid
-real(dp), dimension(ni,nj), intent(in)  :: v_dyn
-real(dp), dimension(ngrid), intent(out) :: v_vect
+integer(di),                   intent(in)  :: ni, nj, ngrid_in
+real(dp), dimension(ni,nj),    intent(in)  :: v_dyn
+real(dp), dimension(ngrid_in), intent(out) :: v_vect
 
 ! LOCAL VARIABLES
@@ -316 +316 @@
 ! CODE
 !-----
-if (ngrid == 1) then ! 1D case
+if (ngrid_in == 1) then ! 1D case
     v_vect(1) = v_dyn(1,1)
     return
 else ! 3D
     ! Check
-    if (ngrid /= (ni - 1)*(nj - 2) + 2) call stop_clean(__FILE__,__LINE__,'problem of dimensions!',1)
+    if (ngrid_in /= (ni - 1)*(nj - 2) + 2) call stop_clean(__FILE__,__LINE__,'problem of dimensions!',1)
 
     ! Initialization
@@ -328 +328 @@
     ! Treatment of the poles
     v_vect(1) = v_dyn(1,1)
-    v_vect(ngrid) = v_dyn(1,nj)
+    v_vect(ngrid_in) = v_dyn(1,nj)
 
     ! Treatment of regular points
@@ -342 +342 @@
 
 !=======================================================================
-SUBROUTINE vect2dyngrd(ni,nj,ngrid,v_vect,v_dyn)
+SUBROUTINE vect2dyngrd(ni,nj,ngrid_in,v_vect,v_dyn)
 !-----------------------------------------------------------------------
 ! NAME
@@ -369 +369 @@
 ! ARGUMENTS
 !----------
-integer(di),                intent(in)  :: ni, nj, ngrid
-real(dp), dimension(ngrid), intent(in)  :: v_vect
-real(dp), dimension(ni,nj), intent(out) :: v_dyn
+integer(di),                   intent(in)  :: ni, nj, ngrid_in
+real(dp), dimension(ngrid_in), intent(in)  :: v_vect
+real(dp), dimension(ni,nj),    intent(out) :: v_dyn
 
 ! LOCAL VARIABLES
@@ -379 +379 @@
 ! CODE
 !-----
-if (ngrid == 1) then ! 1D case
+if (ngrid_in == 1) then ! 1D case
     v_dyn(1,1) = v_vect(1)
     return
 else ! 3D
     ! Check
-    if (ngrid /= (ni - 1)*(nj - 2) + 2) call stop_clean(__FILE__,__LINE__,'problem of dimensions!',1)
+    if (ngrid_in /= (ni - 1)*(nj - 2) + 2) call stop_clean(__FILE__,__LINE__,'problem of dimensions!',1)
 
     ! Initialization
@@ -391 +391 @@
     ! Treatment of the poles
     v_dyn(:,1) = v_vect(1)
-    v_dyn(:,nj) = v_vect(ngrid)
+    v_dyn(:,nj) = v_vect(ngrid_in)
 
     ! Treatment of regular points