/[MITgcm]/MITgcm/pkg/exf/exf_mapfields.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_mapfields.F

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

revision 1.23 by jmc, Tue Apr 28 18:15:33 2009 UTC revision 1.24 by jmc, Thu Dec 22 19:03:41 2011 UTC
# Line 69  c     mythid - thread number for this in Line 69  c     mythid - thread number for this in
69    
70  c     == local variables ==  c     == local variables ==
71    
72        integer bi,bj        INTEGER bi,bj
73        integer i,j,k        INTEGER i,j,ks
74        INTEGER imin, imax        INTEGER imin, imax
75        INTEGER jmin, jmax        INTEGER jmin, jmax
76        PARAMETER ( imin = 1-OLx , imax = sNx+OLx )        PARAMETER ( imin = 1-OLx , imax = sNx+OLx )
# Line 78  c     == local variables == Line 78  c     == local variables ==
78    
79  c     == end of interface ==  c     == end of interface ==
80    
81    C--   set surface level index:
82          ks = 1
83    
84        DO bj = myByLo(myThid),myByHi(myThid)        DO bj = myByLo(myThid),myByHi(myThid)
85          DO bi = myBxLo(myThid),myBxHi(myThid)          DO bi = myBxLo(myThid),myBxHi(myThid)
86    
# Line 160  CADJ STORE ustress(:,:,bi,bj) = comlev1_ Line 163  CADJ STORE ustress(:,:,bi,bj) = comlev1_
163  c     Shift wind stresses calculated at Grid-center to W/S points  c     Shift wind stresses calculated at Grid-center to W/S points
164                fu(i,j,bi,bj) = exf_outscal_ustress*                fu(i,j,bi,bj) = exf_outscal_ustress*
165       &              (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj))       &              (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj))
166       &              *exf_half*maskW(i,j,1,bi,bj)       &              *exf_half*maskW(i,j,ks,bi,bj)
167              enddo              enddo
168             enddo             enddo
169            ENDIF            ENDIF
# Line 198  CADJ STORE vstress(:,:,bi,bj) = comlev1_ Line 201  CADJ STORE vstress(:,:,bi,bj) = comlev1_
201  c     Shift wind stresses calculated at C-points to W/S points  c     Shift wind stresses calculated at C-points to W/S points
202                fv(i,j,bi,bj) = exf_outscal_vstress*                fv(i,j,bi,bj) = exf_outscal_vstress*
203       &              (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj))       &              (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj))
204       &              *exf_half*maskS(i,j,1,bi,bj)       &              *exf_half*maskS(i,j,ks,bi,bj)
205              enddo              enddo
206             enddo             enddo
207            ENDIF            ENDIF
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.22