Changeset 1011 for codes/icosagcm/devel/src/physics/physics_interface.f90

Timestamp:

01/21/20 00:16:21 (4 years ago)

Author:

dubos

Message:

devel : bugfix, missing halo exchange before call to physics_column

File:

• codes/icosagcm/devel/src/physics/physics_interface.f90 (modified) (3 diffs)

codes/icosagcm/devel/src/physics/physics_interface.f90

@@ -7 +7 @@
   TYPE t_physics_inout
      ! Input, time-independent
-     INTEGER :: ngrid
+     INTEGER :: ngrid, it
      REAL(rstd) :: dt_phys
-     REAL(rstd), DIMENSION(:), POINTER :: Ai, lon, lat, phis
+     REAL(rstd), DIMENSION(:), ALLOCATABLE :: Ai, lon, lat, phis
      ! Input, time-dependent
-     REAL(rstd), DIMENSION(:,:), POINTER :: geopot, p, pk, Temp, ulon, ulat
-     REAL(rstd), DIMENSION(:,:,:), POINTER :: q
+     REAL(rstd), DIMENSION(:,:), ALLOCATABLE :: geopot, p, pk, Temp, ulon, ulat
+     REAL(rstd), DIMENSION(:,:,:), ALLOCATABLE :: q
      ! Output arrays
-     REAL(rstd), DIMENSION(:,:), POINTER :: dTemp, dulon, dulat
-     REAL(rstd), DIMENSION(:,:,:), POINTER :: dq
+     REAL(rstd), DIMENSION(:,:), ALLOCATABLE :: dTemp, dulon, dulat
+     REAL(rstd), DIMENSION(:,:,:), ALLOCATABLE :: dq
   END TYPE t_physics_inout
 
@@ -165 +165 @@
        
     CASE DEFAULT ! Copy non-halo points only, as contiguous segments (works)
-       n=0
+       n(:)=0
        n(jj_begin:jj_end)=COUNT(own(ii_begin:ii_end,jj_begin:jj_end),1)
        ngrid=SUM(n)
@@ -185 +185 @@
              DO i=ii_begin,ii_end
                 IF(own(i,j)) THEN
-                   info%n(jj)=n(j)
-                   info%k(jj)=k
-                   info%ij(jj) = iim*(j-1)+i
+                   info%n(jj)=n(j)              ! size of segment
+                   info%k(jj)=k                 ! start index in packed array
+                   info%ij(jj) = iim*(j-1)+i    ! start index in unpacked array
                    IF(COUNT(own(i:i+n(j)-1,j)) /= n(j)) STOP
                    EXIT