/[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.10 by dimitri, Sat Feb 4 02:11:07 2012 UTC
# Line 23  C     == Global variables == Line 23  C     == Global variables ==
23  #include "PARAMS.h"  #include "PARAMS.h"
24  #include "DYNVARS.h"  #include "DYNVARS.h"
25  #include "GRID.h"  #include "GRID.h"
26    #include "FFIELDS.h"
27  #ifdef ALLOW_EXF  #ifdef ALLOW_EXF
28  # include "EXF_OPTIONS.h"  # include "EXF_OPTIONS.h"
29  # include "EXF_FIELDS.h"  # include "EXF_FIELDS.h"
# Line 58  CEOP Line 59  CEOP
59        Real*8  xfer_bc_veloc(2*(Nx+Ny)-6)        Real*8  xfer_bc_veloc(2*(Nx+Ny)-6)
60        _RL     local(1:sNx,1:sNy,nSx,nSy)        _RL     local(1:sNx,1:sNy,nSx,nSy)
61    
62        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  
63        _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)
64        _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)
65        _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 113  C     Send ice area
113          DO bi=1,nSx          DO bi=1,nSx
114           DO j=1,sNy           DO j=1,sNy
115            DO i=1,sNx            DO i=1,sNx
116             local(i,j,bi,bj) = AREA(i,j,1,bi,bj)             local(i,j,bi,bj) = AREA(i,j,bi,bj)
117            ENDDO            ENDDO
118           ENDDO           ENDDO
119          ENDDO          ENDDO
# Line 135  C     Send ice thickness Line 135  C     Send ice thickness
135          DO bi=1,nSx          DO bi=1,nSx
136           DO j=1,sNy           DO j=1,sNy
137            DO i=1,sNx            DO i=1,sNx
138             local(i,j,bi,bj) = HEFF(i,j,1,bi,bj)             local(i,j,bi,bj) = HEFF(i,j,bi,bj)
139            ENDDO            ENDDO
140           ENDDO           ENDDO
141          ENDDO          ENDDO
# Line 216  C     Send boundary ice area Line 216  C     Send boundary ice area
216         DO bi=1,nSx         DO bi=1,nSx
217          DO j=1,sNy          DO j=1,sNy
218           DO i=1,sNx           DO i=1,sNx
219            local(i,j,bi,bj) = AREA(i,j,1,bi,bj)            local(i,j,bi,bj) = AREA(i,j,bi,bj)
220           ENDDO           ENDDO
221          ENDDO          ENDDO
222         ENDDO         ENDDO
# Line 254  C     Send boundary ice thickness Line 254  C     Send boundary ice thickness
254         DO bi=1,nSx         DO bi=1,nSx
255          DO j=1,sNy          DO j=1,sNy
256           DO i=1,sNx           DO i=1,sNx
257            local(i,j,bi,bj) = HEFF(i,j,1,bi,bj)            local(i,j,bi,bj) = HEFF(i,j,bi,bj)
258           ENDDO           ENDDO
259          ENDDO          ENDDO
260         ENDDO         ENDDO
# Line 362  C     Send boundary u ice Line 362  C     Send boundary u ice
362         DO bi=1,nSx         DO bi=1,nSx
363          DO j=1,sNy          DO j=1,sNy
364           DO i=1,sNx           DO i=1,sNx
365            local(i,j,bi,bj) = UICE(i,j,1,bi,bj)            local(i,j,bi,bj) = UICE(i,j,bi,bj)
366           ENDDO           ENDDO
367          ENDDO          ENDDO
368         ENDDO         ENDDO
# Line 398  C     Send boundary v ice Line 398  C     Send boundary v ice
398         DO bi=1,nSx         DO bi=1,nSx
399          DO j=1,sNy          DO j=1,sNy
400           DO i=1,sNx           DO i=1,sNx
401            local(i,j,bi,bj) = VICE(i,j,1,bi,bj)            local(i,j,bi,bj) = VICE(i,j,bi,bj)
402           ENDDO           ENDDO
403          ENDDO          ENDDO
404         ENDDO         ENDDO
# Line 801  C     Receive u surface stress Line 801  C     Receive u surface stress
801         DO bi=1,nSx         DO bi=1,nSx
802          DO j=1,sNy          DO j=1,sNy
803           DO i=1,sNx           DO i=1,sNx
804            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) +
805       &                (1.-AREA(i,j,bi,bj) * fu_tmp(i,j,bi,bj)       &                (1.-AREA(i,j,bi,bj)) * fu_tmp(i,j,bi,bj)
806           ENDDO           ENDDO
807          ENDDO          ENDDO
808         ENDDO         ENDDO
# Line 833  C     Receive v surface stress Line 833  C     Receive v surface stress
833         DO bi=1,nSx         DO bi=1,nSx
834          DO j=1,sNy          DO j=1,sNy
835           DO i=1,sNx           DO i=1,sNx
836            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) +
837       &                (1.-AREA(i,j,bi,bj) * fv_tmp(i,j,bi,bj)       &                (1.-AREA(i,j,bi,bj)) * fv_tmp(i,j,bi,bj)
838           ENDDO           ENDDO
839          ENDDO          ENDDO
840         ENDDO         ENDDO
# Line 866  C     Receive residual shortwave Line 866  C     Receive residual shortwave
866          DO j=1,sNy          DO j=1,sNy
867           DO i=1,sNx           DO i=1,sNx
868            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) +
869       &                  (1.-AREA(i,j,bi,bj) * Qsw_tmp(i,j,bi,bj)       &                  (1.-AREA(i,j,bi,bj)) * Qsw_tmp(i,j,bi,bj)
870           ENDDO           ENDDO
871          ENDDO          ENDDO
872         ENDDO         ENDDO
# Line 899  C     Receive heat flux Line 899  C     Receive heat flux
899           DO i=1,sNx           DO i=1,sNx
900            fv(i,j,bi,bj) = Qsw(i,j,bi,bj) -            fv(i,j,bi,bj) = Qsw(i,j,bi,bj) -
901       &                   AREA(i,j,bi,bj) * local(i,j,bi,bj) +       &                   AREA(i,j,bi,bj) * local(i,j,bi,bj) +
902       &               (1.-AREA(i,j,bi,bj) * Qnet_tmp(i,j,bi,bj)       &               (1.-AREA(i,j,bi,bj)) * Qnet_tmp(i,j,bi,bj)
903           ENDDO           ENDDO
904          ENDDO          ENDDO
905         ENDDO         ENDDO
# Line 930  C     Receive freshwater flux Line 930  C     Receive freshwater flux
930         DO bi=1,nSx         DO bi=1,nSx
931          DO j=1,sNy          DO j=1,sNy
932           DO i=1,sNx           DO i=1,sNx
933            EmPmR(i,j,bi,bj) = - rhoConstFresh *            EmPmR(i,j,bi,bj) = - AREA(i,j,bi,bj)  * local    (i,j,bi,bj) +
934       &                    AREA(i,j,bi,bj) * local(i,j,bi,bj) +       &                  ( 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.10
  ViewVC Help
Powered by ViewVC 1.1.22