Changeset 2091

Timestamp:

Jul 14, 2014, 10:56:11 PM (10 years ago)

Author:

fhourdin

Message:

Correction
Bug fixing

File:

• LMDZ5/trunk/libf/phymar/PHY________INI.f90 (modified) (1 diff)

LMDZ5/trunk/libf/phymar/PHY________INI.f90

@@ -153 +153 @@
 ! Adapted for MAR/LMDZ coupling:
 
-      ii__AP(1)=ixp1
-      jj__AP(1)=jyp1
-      PRINT*,'jyp1=',jyp1
-      PRINT*,'jj__AP(1)=',jj__AP(1)
-      ii__AP(kcolp)=ixp1
-      jj__AP(kcolp)=mypp
-
+!     ii__AP(1)=ixp1
+!     jj__AP(1)=jyp1
+!     PRINT*,'jyp1=',jyp1
+!     PRINT*,'jj__AP(1)=',jj__AP(1)
+!     ii__AP(kcolp)=ixp1
+!     jj__AP(kcolp)=mypp
+!
+!     DO i=ixp1,mxpp
+!     DO j=jyp1+1,mypp-1
+!
+!               ikl    = (j-(jyp1+1)) *mxpp +1 + (i-ixp1+1) ! Tout est décalé de 1 à cause du point isolé au pole dans la grille physique LMD
+!        !       ikl    = (j-jyp1) *mxpp + i-ixp1+1
+!        !       ikl    = (j-jyp1) *(mxp-1) + i-ixp1+1
+!        ii__AP(ikl)   =                 i
+!        jj__AP(ikl)   =  j
+!        ikl_AP(i,j)   =  ikl
+!        PRINT*,'ii__AP(',ikl,')=',ii__AP(ikl)
+!        PRINT*,'jj__AP(',ikl,')=',jj__AP(ikl)
+!     ENDDO
+!     ENDDO
+
+! Modification Gilles Delaygue 2014/07/14 !
+      ikl=1
+      ii__AP(ikl)=ixp1
+      jj__AP(ikl)=jyp1
+      ikl_AP(:,jyp1) = ikl
+
+      DO j=jyp1+1,mypp-1
       DO i=ixp1,mxpp
-      DO j=jyp1+1,mypp-1
-
-                ikl    = (j-(jyp1+1)) *mxpp +1 + (i-ixp1+1) ! Tout est décalé de 1 à cause du point isolé au pole dans la grille physique LMD
-         !       ikl    = (j-jyp1) *mxpp + i-ixp1+1
-         !       ikl    = (j-jyp1) *(mxp-1) + i-ixp1+1
-         ii__AP(ikl)   =                 i
+         ikl=ikl+1
+         ii__AP(ikl)   =  i
          jj__AP(ikl)   =  j
          ikl_AP(i,j)   =  ikl
-!        PRINT*,'ii__AP(',ikl,')=',ii__AP(ikl)
-!        PRINT*,'jj__AP(',ikl,')=',jj__AP(ikl)
       ENDDO
       ENDDO
+
+      ikl=ikl+1
+      ii__AP(ikl)=ixp1
+      jj__AP(ikl)=mypp
+      ikl_AP(:,mypp) = ikl
+
+
+