Changeset 474

Timestamp:

09/12/16 10:05:57 (8 years ago)

Author:

ymipsl

Message:

Bug fix : on physic_column. The "own" variable is correctly defined only inside local domain and on firts halo.

YM

File:

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

codes/icosagcm/trunk/src/physics_interface.f90

@@ -125 +125 @@
        END DO
     CASE(3) ! Copy non-halo points only, one at a time (works, slow ?)
-       n=COUNT(own,1)
+       n=0
+       n(jj_begin:jj_end)=COUNT(own(ii_begin:ii_end,jj_begin:jj_end),1)
        ngrid=SUM(n)
        info%ngrid=ngrid
@@ -146 +147 @@
     CASE DEFAULT ! Copy non-halo points only, as contiguous segments (works)
        n=0
-       n=COUNT(own,1)
+       n(jj_begin:jj_end)=COUNT(own(ii_begin:ii_end,jj_begin:jj_end),1)
        ngrid=SUM(n)
        info%ngrid=ngrid
@@ -154 +155 @@
        ALLOCATE(info%ij(nseg))
        ALLOCATE(info%k(nseg))
-
+       info%n(:)=0
+       info%ij(:)=0
+       info%k(:)=0
+       
        jj=1
        k=1
-       DO j=1,jjm
+       DO j=jj_begin,jj_end
           IF(n(j)>0) THEN
              ! find first .TRUE. value in own(:,j)
-             DO i=1,iim
+             DO i=ii_begin,ii_end
                 IF(own(i,j)) THEN
                    info%n(jj)=n(j)