/[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.7 by edhill, Thu Oct 9 04:19:19 2003 UTC revision 1.8 by heimbach, Thu Nov 6 22:11:03 2003 UTC
# Line 2  c $Header$ Line 2  c $Header$
2    
3  #include "EXF_OPTIONS.h"  #include "EXF_OPTIONS.h"
4    
   
5        subroutine exf_mapfields( mythid )        subroutine exf_mapfields( mythid )
6    
7  c     ==================================================================  c     ==================================================================
# Line 49  c     == global variables == Line 48  c     == global variables ==
48  #include "EEPARAMS.h"  #include "EEPARAMS.h"
49  #include "SIZE.h"  #include "SIZE.h"
50  #include "FFIELDS.h"  #include "FFIELDS.h"
51    #include "GRID.h"
52    
53  #include "exf_param.h"  #include "exf_param.h"
54  #include "exf_constants.h"  #include "exf_constants.h"
55  #include "exf_fields.h"  #include "exf_fields.h"
# Line 140  CADJ STORE ustress(:,:,bi,bj) = comlev1_ Line 141  CADJ STORE ustress(:,:,bi,bj) = comlev1_
141              enddo              enddo
142            enddo            enddo
143            do j = jmin,jmax            do j = jmin,jmax
144              do i = imin,imax              do i = imin+1,imax
145               fu(i,j,bi,bj)   = exf_outscal_ustress*ustress(i,j,bi,bj)  #if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND))
146    c     Shift wind stresses calculated at C-points to W/S points
147                  fu(i,j,bi,bj) = exf_outscal_ustress*
148         &              (ustress(i,j,bi,bj)+ustress(i-1,j,bi,bj))/2.*
149         &              maskW(i,j,1,bi,bj)
150    #else
151                  fu(i,j,bi,bj) = exf_outscal_ustress*ustress(i,j,bi,bj)
152    #endif
153              enddo              enddo
154            enddo            enddo
155    
# Line 166  CADJ STORE vstress(:,:,bi,bj) = comlev1_ Line 174  CADJ STORE vstress(:,:,bi,bj) = comlev1_
174                endif                endif
175              enddo              enddo
176            enddo            enddo
177            do j = jmin,jmax            do j = jmin+1,jmax
178              do i = imin,imax              do i = imin,imax
179               fv(i,j,bi,bj)   = exf_outscal_vstress*vstress(i,j,bi,bj)  #if (defined (ALLOW_ATM_TEMP) || defined (ALLOW_ATM_WIND))
180    c     Shift wind stresses calculated at C-points to W/S points
181                  fv(i,j,bi,bj) = exf_outscal_vstress*
182         &              (vstress(i,j,bi,bj)+vstress(i,j-1,bi,bj))/2.*
183         &              maskS(i,j,1,bi,bj)
184    #else
185                  fv(i,j,bi,bj) = exf_outscal_vstress*vstress(i,j,bi,bj)
186    #endif
187              enddo              enddo
188            enddo            enddo
189    

Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8

  ViewVC Help
Powered by ViewVC 1.1.22