/[MITgcm]/MITgcm_contrib/MPMice/beaufort/code/cpl_mpmice.F
ViewVC logotype

Diff of /MITgcm_contrib/MPMice/beaufort/code/cpl_mpmice.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by dimitri, Sun May 31 03:41:36 2009 UTC revision 1.2 by dimitri, Mon Jun 22 19:40:11 2009 UTC
# Line 57  CEOP Line 57  CEOP
57        Real*8  xfer_bc_tracer(2*(Nx+Ny)-4)        Real*8  xfer_bc_tracer(2*(Nx+Ny)-4)
58        Real*8  xfer_bc_veloc(2*(Nx+Ny)-6)        Real*8  xfer_bc_veloc(2*(Nx+Ny)-6)
59        _RL     local(1:sNx,1:sNy,nSx,nSy)        _RL     local(1:sNx,1:sNy,nSx,nSy)
60        _RL ScatArray(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
61    #ifdef CPL_DEBUG
62         _RL ScatArray(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
63          DO bj=1,nSy
64           DO bi=1,nSx
65            DO j=1-OLy,sNy+OLy
66             DO i=1-OLx,sNx+OLx
67              ScatArray(i,j,bi,bj) = 0.0 _d 0
68             ENDDO
69            ENDDO
70           ENDDO
71          ENDDO
72    #endif
73    
74        IF( myTime .EQ. startTime ) THEN        IF( myTime .EQ. startTime ) THEN
75    
# Line 633  C     Receive ice area Line 645  C     Receive ice area
645        ENDIF        ENDIF
646        _END_MASTER( myThid )        _END_MASTER( myThid )
647        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
648    #ifdef CPL_DEBUG
649        DO bj=1,nSy        DO bj=1,nSy
650         DO bi=1,nSx         DO bi=1,nSx
651          DO j=1,sNy          DO j=1,sNy
# Line 642  C     Receive ice area Line 655  C     Receive ice area
655          ENDDO          ENDDO
656         ENDDO         ENDDO
657        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
658        CALL PLOT_FIELD_XYRL( ScatArray, 'ice area', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'ice area', myIter, myThid )
659  #endif  #endif
660    
# Line 655  C     Receive ice thickness Line 667  C     Receive ice thickness
667        ENDIF        ENDIF
668        _END_MASTER( myThid )        _END_MASTER( myThid )
669        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
670    #ifdef CPL_DEBUG
671        DO bj=1,nSy        DO bj=1,nSy
672         DO bi=1,nSx         DO bi=1,nSx
673          DO j=1,sNy          DO j=1,sNy
# Line 664  C     Receive ice thickness Line 677  C     Receive ice thickness
677          ENDDO          ENDDO
678         ENDDO         ENDDO
679        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
680        CALL PLOT_FIELD_XYRL( ScatArray, 'ice thickness', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'ice thickness', myIter, myThid )
681  #endif  #endif
682    
# Line 677  C     Receive ice salinity Line 689  C     Receive ice salinity
689        ENDIF        ENDIF
690        _END_MASTER( myThid )        _END_MASTER( myThid )
691        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
692    #ifdef CPL_DEBUG
693        DO bj=1,nSy        DO bj=1,nSy
694         DO bi=1,nSx         DO bi=1,nSx
695          DO j=1,sNy          DO j=1,sNy
# Line 686  C     Receive ice salinity Line 699  C     Receive ice salinity
699          ENDDO          ENDDO
700         ENDDO         ENDDO
701        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
702        CALL PLOT_FIELD_XYRL( ScatArray, 'ice salinity', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'ice salinity', myIter, myThid )
703  #endif  #endif
704    
# Line 699  C     Receive snow thickness Line 711  C     Receive snow thickness
711        ENDIF        ENDIF
712        _END_MASTER( myThid )        _END_MASTER( myThid )
713        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
714    #ifdef CPL_DEBUG
715        DO bj=1,nSy        DO bj=1,nSy
716         DO bi=1,nSx         DO bi=1,nSx
717          DO j=1,sNy          DO j=1,sNy
# Line 708  C     Receive snow thickness Line 721  C     Receive snow thickness
721          ENDDO          ENDDO
722         ENDDO         ENDDO
723        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
724        CALL PLOT_FIELD_XYRL( ScatArray, 'ice thickness', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'ice thickness', myIter, myThid )
725  #endif  #endif
726    
# Line 721  C     Receive u surface stress Line 733  C     Receive u surface stress
733        ENDIF        ENDIF
734        _END_MASTER( myThid )        _END_MASTER( myThid )
735        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
736    #ifdef CPL_DEBUG
737        DO bj=1,nSy        DO bj=1,nSy
738         DO bi=1,nSx         DO bi=1,nSx
739          DO j=1,sNy          DO j=1,sNy
# Line 730  C     Receive u surface stress Line 743  C     Receive u surface stress
743          ENDDO          ENDDO
744         ENDDO         ENDDO
745        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
746        CALL PLOT_FIELD_XYRL( ScatArray, 'u stress', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'u stress', myIter, myThid )
747  #endif  #endif
748    
# Line 743  C     Receive v surface stress Line 755  C     Receive v surface stress
755        ENDIF        ENDIF
756        _END_MASTER( myThid )        _END_MASTER( myThid )
757        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
758    #ifdef CPL_DEBUG
759        DO bj=1,nSy        DO bj=1,nSy
760         DO bi=1,nSx         DO bi=1,nSx
761          DO j=1,sNy          DO j=1,sNy
# Line 752  C     Receive v surface stress Line 765  C     Receive v surface stress
765          ENDDO          ENDDO
766         ENDDO         ENDDO
767        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
768        CALL PLOT_FIELD_XYRL( ScatArray, 'v stress', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'v stress', myIter, myThid )
769  #endif  #endif
770    
# Line 765  C     Receive residual shortwave Line 777  C     Receive residual shortwave
777        ENDIF        ENDIF
778        _END_MASTER( myThid )        _END_MASTER( myThid )
779        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
780        DO bj=1,nSy   #ifdef CPL_DEBUG
781         DO bj=1,nSy
782         DO bi=1,nSx         DO bi=1,nSx
783          DO j=1,sNy          DO j=1,sNy
784           DO i=1,sNx           DO i=1,sNx
# Line 774  C     Receive residual shortwave Line 787  C     Receive residual shortwave
787          ENDDO          ENDDO
788         ENDDO         ENDDO
789        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
790        CALL PLOT_FIELD_XYRL( ScatArray, 'shortwave', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'shortwave', myIter, myThid )
791  #endif  #endif
792    
# Line 787  C     Receive heat flux Line 799  C     Receive heat flux
799        ENDIF        ENDIF
800        _END_MASTER( myThid )        _END_MASTER( myThid )
801        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
802    #ifdef CPL_DEBUG
803        DO bj=1,nSy        DO bj=1,nSy
804         DO bi=1,nSx         DO bi=1,nSx
805          DO j=1,sNy          DO j=1,sNy
# Line 796  C     Receive heat flux Line 809  C     Receive heat flux
809          ENDDO          ENDDO
810         ENDDO         ENDDO
811        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
812        CALL PLOT_FIELD_XYRL( ScatArray, 'heat flux', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'heat flux', myIter, myThid )
813  #endif  #endif
814    
# Line 809  C     Receive freshwater flux Line 821  C     Receive freshwater flux
821        ENDIF        ENDIF
822        _END_MASTER( myThid )        _END_MASTER( myThid )
823        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
824    #ifdef CPL_DEBUG
825        DO bj=1,nSy        DO bj=1,nSy
826         DO bi=1,nSx         DO bi=1,nSx
827          DO j=1,sNy          DO j=1,sNy
# Line 818  C     Receive freshwater flux Line 831  C     Receive freshwater flux
831          ENDDO          ENDDO
832         ENDDO         ENDDO
833        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
834        CALL PLOT_FIELD_XYRL( ScatArray, 'freshwater', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'freshwater', myIter, myThid )
835  #endif  #endif
836    
# Line 831  C     Receive salt flux Line 843  C     Receive salt flux
843        ENDIF        ENDIF
844        _END_MASTER( myThid )        _END_MASTER( myThid )
845        CALL SCATTER_2D( xfer_array, local, myThid )        CALL SCATTER_2D( xfer_array, local, myThid )
846    #ifdef CPL_DEBUG
847        DO bj=1,nSy        DO bj=1,nSy
848         DO bi=1,nSx         DO bi=1,nSx
849          DO j=1,sNy          DO j=1,sNy
# Line 840  C     Receive salt flux Line 853  C     Receive salt flux
853          ENDDO          ENDDO
854         ENDDO         ENDDO
855        ENDDO        ENDDO
 #ifdef CPL_DEBUG  
856        CALL PLOT_FIELD_XYRL( ScatArray, 'salt flux', myIter, myThid )        CALL PLOT_FIELD_XYRL( ScatArray, 'salt flux', myIter, myThid )
857  #endif  #endif
858    

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2

  ViewVC Help
Powered by ViewVC 1.1.22