/[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.5 by dimitri, Fri Feb 3 14:08:28 2012 UTC revision 1.8 by dimitri, Fri Feb 3 18:12:59 2012 UTC
# Line 58  CEOP Line 58  CEOP
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    
61        COMMON /FFIELDS_tmp/        COMMON /FFIELDS_tmp/ fu_tmp, fv_tmp, Qnet_tmp, Qsw_tmp, EmPmR_tmp
      & fu_tmp, fv_tmp, Qnet_tmp, Qsw_tmp, EmPmR_tmp, saltFlux_tmp  
62        _RS  fu_tmp       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  fu_tmp       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
63        _RS  fv_tmp       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  fv_tmp       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
64        _RS  Qnet_tmp     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  Qnet_tmp     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 113  C     Send ice area Line 112  C     Send ice area
112          DO bi=1,nSx          DO bi=1,nSx
113           DO j=1,sNy           DO j=1,sNy
114            DO i=1,sNx            DO i=1,sNx
115             local(i,j,bi,bj) = AREA(i,j,1,bi,bj)             local(i,j,bi,bj) = AREA(i,j,bi,bj)
116            ENDDO            ENDDO
117           ENDDO           ENDDO
118          ENDDO          ENDDO
# Line 135  C     Send ice thickness Line 134  C     Send ice thickness
134          DO bi=1,nSx          DO bi=1,nSx
135           DO j=1,sNy           DO j=1,sNy
136            DO i=1,sNx            DO i=1,sNx
137             local(i,j,bi,bj) = HEFF(i,j,1,bi,bj)             local(i,j,bi,bj) = HEFF(i,j,bi,bj)
138            ENDDO            ENDDO
139           ENDDO           ENDDO
140          ENDDO          ENDDO
# Line 216  C     Send boundary ice area Line 215  C     Send boundary ice area
215         DO bi=1,nSx         DO bi=1,nSx
216          DO j=1,sNy          DO j=1,sNy
217           DO i=1,sNx           DO i=1,sNx
218            local(i,j,bi,bj) = AREA(i,j,1,bi,bj)            local(i,j,bi,bj) = AREA(i,j,bi,bj)
219           ENDDO           ENDDO
220          ENDDO          ENDDO
221         ENDDO         ENDDO
# Line 254  C     Send boundary ice thickness Line 253  C     Send boundary ice thickness
253         DO bi=1,nSx         DO bi=1,nSx
254          DO j=1,sNy          DO j=1,sNy
255           DO i=1,sNx           DO i=1,sNx
256            local(i,j,bi,bj) = HEFF(i,j,1,bi,bj)            local(i,j,bi,bj) = HEFF(i,j,bi,bj)
257           ENDDO           ENDDO
258          ENDDO          ENDDO
259         ENDDO         ENDDO
# Line 362  C     Send boundary u ice Line 361  C     Send boundary u ice
361         DO bi=1,nSx         DO bi=1,nSx
362          DO j=1,sNy          DO j=1,sNy
363           DO i=1,sNx           DO i=1,sNx
364            local(i,j,bi,bj) = UICE(i,j,1,bi,bj)            local(i,j,bi,bj) = UICE(i,j,bi,bj)
365           ENDDO           ENDDO
366          ENDDO          ENDDO
367         ENDDO         ENDDO
# Line 398  C     Send boundary v ice Line 397  C     Send boundary v ice
397         DO bi=1,nSx         DO bi=1,nSx
398          DO j=1,sNy          DO j=1,sNy
399           DO i=1,sNx           DO i=1,sNx
400            local(i,j,bi,bj) = VICE(i,j,1,bi,bj)            local(i,j,bi,bj) = VICE(i,j,bi,bj)
401           ENDDO           ENDDO
402          ENDDO          ENDDO
403         ENDDO         ENDDO
# Line 801  C     Receive u surface stress Line 800  C     Receive u surface stress
800         DO bi=1,nSx         DO bi=1,nSx
801          DO j=1,sNy          DO j=1,sNy
802           DO i=1,sNx           DO i=1,sNx
803            fu(i,j,bi,bj) = AREA(i,j,bi,bj) * local(i,j,bi,bj)            fu(i,j,bi,bj) = AREA(i,j,bi,bj) * local(i,j,bi,bj) +
804       &                (1.-AREA(i,j,bi,bj) * fu_tmp(i,j,bi,bj)       &                (1.-AREA(i,j,bi,bj)) * fu_tmp(i,j,bi,bj)
805           ENDDO           ENDDO
806          ENDDO          ENDDO
807         ENDDO         ENDDO
# Line 833  C     Receive v surface stress Line 832  C     Receive v surface stress
832         DO bi=1,nSx         DO bi=1,nSx
833          DO j=1,sNy          DO j=1,sNy
834           DO i=1,sNx           DO i=1,sNx
835            fv(i,j,bi,bj) = AREA(i,j,bi,bj) * local(i,j,bi,bj)            fv(i,j,bi,bj) = AREA(i,j,bi,bj) * local(i,j,bi,bj) +
836       &                (1.-AREA(i,j,bi,bj) * fv_tmp(i,j,bi,bj)       &                (1.-AREA(i,j,bi,bj)) * fv_tmp(i,j,bi,bj)
837           ENDDO           ENDDO
838          ENDDO          ENDDO
839         ENDDO         ENDDO
# Line 866  C     Receive residual shortwave Line 865  C     Receive residual shortwave
865          DO j=1,sNy          DO j=1,sNy
866           DO i=1,sNx           DO i=1,sNx
867            Qsw(i,j,bi,bj) = -AREA(i,j,bi,bj) * local(i,j,bi,bj) +            Qsw(i,j,bi,bj) = -AREA(i,j,bi,bj) * local(i,j,bi,bj) +
868       &                  (1.-AREA(i,j,bi,bj) * Qsw_tmp(i,j,bi,bj)       &                  (1.-AREA(i,j,bi,bj)) * Qsw_tmp(i,j,bi,bj)
869           ENDDO           ENDDO
870          ENDDO          ENDDO
871         ENDDO         ENDDO
# Line 899  C     Receive heat flux Line 898  C     Receive heat flux
898           DO i=1,sNx           DO i=1,sNx
899            fv(i,j,bi,bj) = Qsw(i,j,bi,bj) -            fv(i,j,bi,bj) = Qsw(i,j,bi,bj) -
900       &                   AREA(i,j,bi,bj) * local(i,j,bi,bj) +       &                   AREA(i,j,bi,bj) * local(i,j,bi,bj) +
901       &               (1.-AREA(i,j,bi,bj) * Qnet_tmp(i,j,bi,bj)       &               (1.-AREA(i,j,bi,bj)) * Qnet_tmp(i,j,bi,bj)
902           ENDDO           ENDDO
903          ENDDO          ENDDO
904         ENDDO         ENDDO
# Line 932  C     Receive freshwater flux Line 931  C     Receive freshwater flux
931           DO i=1,sNx           DO i=1,sNx
932            EmPmR(i,j,bi,bj) = - rhoConstFresh *            EmPmR(i,j,bi,bj) = - rhoConstFresh *
933       &                    AREA(i,j,bi,bj) * local(i,j,bi,bj) +       &                    AREA(i,j,bi,bj) * local(i,j,bi,bj) +
934       &                (1.-AREA(i,j,bi,bj) * EmPmR_tmp(i,j,bi,bj)       &                (1.-AREA(i,j,bi,bj)) * EmPmR_tmp(i,j,bi,bj)
935           ENDDO           ENDDO
936          ENDDO          ENDDO
937         ENDDO         ENDDO
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22