Changeset 22653

Timestamp:

04/02/18 11:12:46 (7 years ago)

Author:

dmenemen

Message:

sending melt rate

File:

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

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

@@ -132 +132 @@
             _END_MASTER( myThid )
          ENDIF
+
+C     Send melt rate
+         print*,'this is yC',yC
+         DO bj=1,nSy
+            DO bi=1,nSx
+               DO j=1,sNy
+                  DO i=1,sNx
+                     local(i,j,bi,bj) = shelficeHeatFlux(i,j,bi,bj)
+                  ENDDO
+               ENDDO
+            ENDDO
+         ENDDO
+         print*,'this is local shelficeHeatFlux',local
+         CALL BAR2( myThid ) 
+         CALL GATHER_2D_R8( xfer_array, local, Nx, Ny,
+     &                      .FALSE., .FALSE., myThid )
+         print*,'this is global shelficeHeatFlux',xfer_array
+         IF( myProcId .EQ. 0 ) THEN
+            _BEGIN_MASTER( myThid )
+            buffsize = Nx*Ny
+            CALL MPI_SEND(xfer_array,buffsize,MPI_DOUBLE_PRECISION,
+     &           0,10001007,toissmcomm,mpistatus)
+            _END_MASTER( myThid )
+         ENDIF
+         CALL BAR2( myThid )
+         print*,'Done Sending shelficeHeatFlux array.'
+
       ENDIF