/[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.8 by heimbach, Thu Nov 6 22:11:03 2003 UTC revision 1.12 by heimbach, Wed Jun 29 07:11:18 2005 UTC
# Line 125  CADJ STORE ustress(:,:,bi,bj) = comlev1_ Line 125  CADJ STORE ustress(:,:,bi,bj) = comlev1_
125            do j = jmin,jmax            do j = jmin,jmax
126              do i = imin,imax              do i = imin,imax
127  c             Zonal wind stress.  c             Zonal wind stress.
128                if (ustress(i,j,bi,bj).gt.2.0D0) then                if (ustress(i,j,bi,bj).gt.windstressmax) then
129                  ustress(i,j,bi,bj)=2.0D0                  ustress(i,j,bi,bj)=windstressmax
130                endif                endif
131              enddo              enddo
132            enddo            enddo
# Line 135  CADJ STORE ustress(:,:,bi,bj) = comlev1_ Line 135  CADJ STORE ustress(:,:,bi,bj) = comlev1_
135  #endif  #endif
136            do j = jmin,jmax            do j = jmin,jmax
137              do i = imin,imax              do i = imin,imax
138                if (ustress(i,j,bi,bj).lt.-2.0D0) then                if (ustress(i,j,bi,bj).lt.-windstressmax) then
139                  ustress(i,j,bi,bj)=-2.0D0                  ustress(i,j,bi,bj)=-windstressmax
140                endif                endif
141              enddo              enddo
142            enddo            enddo
143            do j = jmin,jmax            do j = jmin,jmax
144              do i = imin+1,imax              do i = imin+1,imax
145  #if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND))  #if (defined (ALLOW_BULKFORMULAE) || defined (USE_EXF_INTERPOLATION))
146  c     Shift wind stresses calculated at C-points to W/S points  c     Shift wind stresses calculated at C-points to W/S points
147                fu(i,j,bi,bj) = exf_outscal_ustress*                fu(i,j,bi,bj) = exf_outscal_ustress*
148       &              (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj))/2.*       &              (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj))/2.*
# Line 159  CADJ STORE vstress(:,:,bi,bj) = comlev1_ Line 159  CADJ STORE vstress(:,:,bi,bj) = comlev1_
159            do j = jmin,jmax            do j = jmin,jmax
160              do i = imin,imax              do i = imin,imax
161  c             Meridional wind stress.  c             Meridional wind stress.
162                if (vstress(i,j,bi,bj).gt.2.0D0) then                if (vstress(i,j,bi,bj).gt.windstressmax) then
163                  vstress(i,j,bi,bj)=2.0D0                  vstress(i,j,bi,bj)=windstressmax
164                endif                endif
165              enddo              enddo
166            enddo            enddo
# Line 169  CADJ STORE vstress(:,:,bi,bj) = comlev1_ Line 169  CADJ STORE vstress(:,:,bi,bj) = comlev1_
169  #endif  #endif
170            do j = jmin,jmax            do j = jmin,jmax
171              do i = imin,imax              do i = imin,imax
172                if (vstress(i,j,bi,bj).lt.-2.0D0) then                if (vstress(i,j,bi,bj).lt.-windstressmax) then
173                  vstress(i,j,bi,bj)=-2.0D0                  vstress(i,j,bi,bj)=-windstressmax
174                endif                endif
175              enddo              enddo
176            enddo            enddo
177            do j = jmin+1,jmax            do j = jmin+1,jmax
178              do i = imin,imax              do i = imin,imax
179  #if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND))  #if (defined (ALLOW_BULKFORMULAE) || defined (USE_EXF_INTERPOLATION))
180  c     Shift wind stresses calculated at C-points to W/S points  c     Shift wind stresses calculated at C-points to W/S points
181                fv(i,j,bi,bj) = exf_outscal_vstress*                fv(i,j,bi,bj) = exf_outscal_vstress*
182       &              (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj))/2.*       &              (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj))/2.*
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22