Changeset 22546

Timestamp:

03/16/18 12:09:04 (7 years ago)

Author:

dmenemen

Message:

MITgcm sends XC/YC to ISSM.

File:

• issm/trunk-jpl/test/MITgcm/code/cpl_issm.F (modified) (1 diff)

issm/trunk-jpl/test/MITgcm/code/cpl_issm.F

@@ -48 +48 @@
      &     mpiMyWid, toissmcomm
       INTEGER mpiMyWid, toissmcomm
-      integer mpistatus(MPI_STATUS_SIZE), mpiRC
+      integer mpistatus(MPI_STATUS_SIZE)
       real*8 CouplingTime, IceModelTime
+      integer i, j, bi, bj, buffsize
+      Real*8  xfer_array(Nx,Ny)
+      _RL     local(1:sNx,1:sNy,nSx,nSy)
 
-      print*,'choubi 1'
-      
       IF( myTime .EQ. startTime ) THEN
-      print*,'choubi 2'
 
 C     Send deltatimestep
-       _BEGIN_MASTER( myThid )
-
-       print*,'choubi 3'
-
+         _BEGIN_MASTER( myThid )
          call MPI_Recv(CouplingTime,1,MPI_DOUBLE,0,10001000,toissmcomm,
-     &        mpistatus,mpiRC)
+     &        mpistatus)
          print*, 'Ocean received CouplingTime: ', CouplingTime
          call MPI_Recv(IceModelTime,1,MPI_DOUBLE,0,10001001,toissmcomm,
-     &        mpistatus,mpiRC)
+     &        mpistatus)
          print*, 'Ocean received IceModelTime: ', IceModelTime
          call MPI_Send(myTime,1,MPI_DOUBLE,0,10001002,toissmcomm,
-     &        mpistatus,mpiRC)
+     &        mpistatus)
          call MPI_Send(Nx,1,MPI_INT,0,10001003,toissmcomm,
-     &        mpistatus,mpiRC)
+     &        mpistatus)
          call MPI_Send(Ny,1,MPI_INT,0,10001004,toissmcomm,
-     &        mpistatus,mpiRC)
-      
-cdb      call MPI_Bcast(dummy1,1,MPI_INT,0,MPI_COMM_MODEL,mpiRC)
-cdb      print*, 'Ocean : dummy received:',dummy1
+     &        mpistatus)
+         _END_MASTER( myThid )
 
-      _END_MASTER( myThid )
+C     Send longitude East of center of cell
+         DO bj=1,nSy
+            DO bi=1,nSx
+               DO j=1,sNy
+                  DO i=1,sNx
+                     local(i,j,bi,bj) = xC(i,j,bi,bj)
+                  ENDDO
+               ENDDO
+            ENDDO
+         ENDDO
+         CALL GATHER_2D_R8( xfer_array, local, myThid )
+         _BEGIN_MASTER( myThid )
+         buffsize = Nx*Ny
+         CALL MPI_SEND(xfer_array,buffsize,MPI_DOUBLE_PRECISION,
+     &        0,10001005,toissmcomm,mpistatus)
+         _END_MASTER( myThid )
+         
+C     Send latitude North of center of cell
+         DO bj=1,nSy
+            DO bi=1,nSx
+               DO j=1,sNy
+                  DO i=1,sNx
+                     local(i,j,bi,bj) = yC(i,j,bi,bj)
+                  ENDDO
+               ENDDO
+            ENDDO
+         ENDDO
+         CALL GATHER_2D_R8( xfer_array, local, myThid )
+         _BEGIN_MASTER( myThid )
+         buffsize = Nx*Ny
+         CALL MPI_SEND(xfer_array,buffsize,MPI_DOUBLE_PRECISION,
+     &        0,10001006,toissmcomm,mpistatus)
+         _END_MASTER( myThid )
+
       ENDIF