/[MITgcm]/MITgcm/model/src/external_forcing_surf.F
ViewVC logotype

Diff of /MITgcm/model/src/external_forcing_surf.F

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

revision 1.50 by jmc, Sun Sep 5 22:29:16 2010 UTC revision 1.51 by jmc, Sat Sep 11 21:27:13 2010 UTC
# Line 129  C     of surf.thickness (hFac): account Line 129  C     of surf.thickness (hFac): account
129            ENDDO            ENDDO
130           ENDDO           ENDDO
131  # endif /* DISABLE_RSTAR_CODE */  # endif /* DISABLE_RSTAR_CODE */
132            ELSEIF ( selectSigmaCoord.NE.0 ) THEN
133    # ifndef DISABLE_SIGMA_CODE
134             DO j=jMin,jMax
135              DO i=iMin,iMax
136                surfaceForcingT(i,j,bi,bj) = surfaceForcingT(i,j,bi,bj)
137         &        *( 1. _d 0 + dEtaHdt(i,j,bi,bj)*deltaTfreesurf
138         &                    *dBHybSigF(ks)*recip_drF(ks)
139         &                    *recip_hFacC(i,j,ks,bi,bj)
140         &         )
141                surfaceForcingS(i,j,bi,bj) = surfaceForcingS(i,j,bi,bj)
142         &        *( 1. _d 0 + dEtaHdt(i,j,bi,bj)*deltaTfreesurf
143         &                    *dBHybSigF(ks)*recip_drF(ks)
144         &                    *recip_hFacC(i,j,ks,bi,bj)
145         &         )
146              ENDDO
147             ENDDO
148    # endif /* DISABLE_SIGMA_CODE */
149          ELSE          ELSE
150           DO j=jMin,jMax           DO j=jMin,jMax
151            DO i=iMin,iMax            DO i=iMin,iMax
152             IF (ks.EQ.ksurfC(i,j,bi,bj)) THEN             IF (ks.EQ.kSurfC(i,j,bi,bj)) THEN
153              surfaceForcingT(i,j,bi,bj) = surfaceForcingT(i,j,bi,bj)              surfaceForcingT(i,j,bi,bj) = surfaceForcingT(i,j,bi,bj)
154       &             *_recip_hFacC(i,j,ks,bi,bj)*hFac_surfC(i,j,bi,bj)       &             *_recip_hFacC(i,j,ks,bi,bj)*hFac_surfC(i,j,bi,bj)
155              surfaceForcingS(i,j,bi,bj) = surfaceForcingS(i,j,bi,bj)              surfaceForcingS(i,j,bi,bj) = surfaceForcingS(i,j,bi,bj)
Legend:
Removed from v.1.50  
changed lines
  Added in v.1.51
  ViewVC Help
Powered by ViewVC 1.1.22