/[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.3 by dgoldberg, Thu Sep 20 15:45:13 2012 UTC
# Line 40  C     LOCAL VARIABLES Line 40  C     LOCAL VARIABLES
40         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
41          DO j=1,sNy          DO j=1,sNy
42           DO i=1,sNx           DO i=1,sNx
43    cph(
44              streamice_taubx (i,j,bi,bj) = 0.
45              streamice_tauby (i,j,bi,bj) = 0.
46    cph)
47            IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN            IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN
48        
49    !            streamice_taubx (i,j,bi,bj) =
50    !      &      tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)
51    !            streamice_tauby (i,j,bi,bj) =
52    !      &      tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)
53    
54               umid = 0.0
55               vmid = 0.0
56    
57               DO k=0,1
58                DO l=0,1
59                 umid = umid + 0.25 *
60         &        dxG(i,j+l,bi,bj)*dyG(i+k,j,bi,bj) *
61         &        recip_rA(i,j,bi,bj) *
62         &        U_streamice(i+k,j+l,bi,bj)
63                 vmid = vmid + 0.25 *
64         &        dxG(i,j+l,bi,bj)*dyG(i+k,j,bi,bj) *
65         &        recip_rA(i,j,bi,bj) *
66         &        V_streamice(i+k,j+l,bi,bj)
67                ENDDO
68               ENDDO
69      
70             streamice_taubx (i,j,bi,bj) =             streamice_taubx (i,j,bi,bj) =
71       &      tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)       &      tau_beta_eff_streamice(i,j,bi,bj) * umid
72             streamice_tauby (i,j,bi,bj) =             streamice_tauby (i,j,bi,bj) =
73       &      tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)       &      tau_beta_eff_streamice(i,j,bi,bj) * vmid
74                        
75            ENDIF            ENDIF
76           ENDDO           ENDDO
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22