/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_taub.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_taub.F

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

revision 1.1 by dgoldberg, Tue Sep 18 17:04:40 2012 UTC revision 1.5 by dgoldberg, Mon Dec 10 02:34:45 2012 UTC
# Line 30  C     !INPUT/OUTPUT ARGUMENTS Line 30  C     !INPUT/OUTPUT ARGUMENTS
30  !       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  !       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
31    
32  #ifdef ALLOW_STREAMICE  #ifdef ALLOW_STREAMICE
33    #ifdef STREAMICE_HYBRID_STRESS
34    
35    
36  C     LOCAL VARIABLES  C     LOCAL VARIABLES
37        INTEGER i, j, bi, bj, k, l, umid, vmid        INTEGER i, j, bi, bj, k, l
38        _RL ux, uy, vx, vy, exx, eyy, exy, second_inv, unorm        _RL ux, uy, vx, vy, exx, eyy, exy, second_inv, unorm
39          _RL umid, vmid
40    
41        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
42         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
43          DO j=1,sNy          DO j=1,sNy
44           DO i=1,sNx           DO i=1,sNx
45    cph(
46              streamice_taubx (i,j,bi,bj) = 0.
47              streamice_tauby (i,j,bi,bj) = 0.
48    cph)
49            IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN            IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN
50        
51    !            streamice_taubx (i,j,bi,bj) =
52    !      &      tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)
53    !            streamice_tauby (i,j,bi,bj) =
54    !      &      tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)
55    
56               umid = 0.0
57               vmid = 0.0
58    
59               DO k=0,1
60                DO l=0,1
61                 umid = umid + 0.25 *
62         &        dxG(i,j+l,bi,bj)*dyG(i+k,j,bi,bj) *
63         &        recip_rA(i,j,bi,bj) *
64         &        U_streamice(i+k,j+l,bi,bj)
65                 vmid = vmid + 0.25 *
66         &        dxG(i,j+l,bi,bj)*dyG(i+k,j,bi,bj) *
67         &        recip_rA(i,j,bi,bj) *
68         &        V_streamice(i+k,j+l,bi,bj)
69                ENDDO
70               ENDDO
71      
72             streamice_taubx (i,j,bi,bj) =             streamice_taubx (i,j,bi,bj) =
73       &      tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)       &      tau_beta_eff_streamice(i,j,bi,bj) * umid
74             streamice_tauby (i,j,bi,bj) =             streamice_tauby (i,j,bi,bj) =
75       &      tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)       &      tau_beta_eff_streamice(i,j,bi,bj) * vmid
76                        
77            ENDIF            ENDIF
78           ENDDO           ENDDO
# Line 55  C     LOCAL VARIABLES Line 82  C     LOCAL VARIABLES
82    
83    
84  #endif  #endif
85    #endif
86        RETURN        RETURN
87        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22