/[MITgcm]/MITgcm/pkg/seaice/seaice_ocean_stress.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_ocean_stress.F

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

revision 1.16 by mlosch, Tue Apr 24 18:38:15 2007 UTC revision 1.17 by jmc, Mon May 14 14:29:49 2007 UTC
# Line 56  c     introduce turning angle (default i Line 56  c     introduce turning angle (default i
56        SINWIN=SIN(SEAICE_airTurnAngle*deg2rad)        SINWIN=SIN(SEAICE_airTurnAngle*deg2rad)
57        COSWIN=COS(SEAICE_airTurnAngle*deg2rad)        COSWIN=COS(SEAICE_airTurnAngle*deg2rad)
58    
 C--   Update overlap regions  
       CALL EXCH_UV_XY_RL(WINDX, WINDY, .TRUE., myThid)  
   
 #ifndef SEAICE_EXTERNAL_FLUXES  
 C--   Interpolate wind stress (N/m^2) from C-points of C-grid  
 C     to U and V points of C-grid for forcing the ocean model.  
       DO bj=myByLo(myThid),myByHi(myThid)  
        DO bi=myBxLo(myThid),myBxHi(myThid)  
         DO j=1,sNy  
          DO i=1,sNx  
           fu(I,J,bi,bj)=0.5*(WINDX(I,J,bi,bj) + WINDX(I-1,J,bi,bj))  
           fv(I,J,bi,bj)=0.5*(WINDY(I,J,bi,bj) + WINDY(I,J-1,bi,bj))  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
 #endif /* ifndef SEAICE_EXTERNAL_FLUXES */  
   
59        IF ( useHB87StressCoupling ) THEN        IF ( useHB87StressCoupling ) THEN
60  C  C
61  C     use an intergral over ice and ocean surface layer to define  C     use an intergral over ice and ocean surface layer to define
# Line 188  C     stress and internal ice stresses t Line 170  C     stress and internal ice stresses t
170           ENDIF           ENDIF
171          ENDDO          ENDDO
172         ENDDO         ENDDO
173    
174        ELSE        ELSE
175    C     else: useHB87StressCoupling=F
176    
177  C--   Compute ice-affected wind stress (interpolate to U/V-points)  C--   Compute ice-affected wind stress (interpolate to U/V-points)
178  C     by averaging wind stress and ice-ocean stress according to  C     by averaging wind stress and ice-ocean stress according to
# Line 232  C     ice cover Line 216  C     ice cover
216        ENDIF        ENDIF
217        CALL EXCH_UV_XY_RS(fu, fv, .TRUE., myThid)        CALL EXCH_UV_XY_RS(fu, fv, .TRUE., myThid)
218    
219  #endif /* not SEAICE_CGRID */  #endif /* SEAICE_CGRID */
220    
221        RETURN        RETURN
222        END        END
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22