/[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.2 by heimbach, Thu Sep 20 02:04:45 2012 UTC revision 1.3 by dgoldberg, Thu Sep 20 15:45:13 2012 UTC
# Line 46  cph( Line 46  cph(
46  cph)  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.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22