/[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.4 by dgoldberg, Tue Sep 25 03:03:21 2012 UTC
# Line 33  C     !INPUT/OUTPUT ARGUMENTS Line 33  C     !INPUT/OUTPUT ARGUMENTS
33    
34    
35  C     LOCAL VARIABLES  C     LOCAL VARIABLES
36        INTEGER i, j, bi, bj, k, l, umid, vmid        INTEGER i, j, bi, bj, k, l
37        _RL ux, uy, vx, vy, exx, eyy, exy, second_inv, unorm        _RL ux, uy, vx, vy, exx, eyy, exy, second_inv, unorm
38          _RL umid, vmid
39    
40        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
41         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 46  cph( Line 47  cph(
47  cph)  cph)
48            IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN            IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN
49        
50    !            streamice_taubx (i,j,bi,bj) =
51    !      &      tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)
52    !            streamice_tauby (i,j,bi,bj) =
53    !      &      tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)
54    
55               umid = 0.0
56               vmid = 0.0
57    
58               DO k=0,1
59                DO l=0,1
60                 umid = umid + 0.25 *
61         &        dxG(i,j+l,bi,bj)*dyG(i+k,j,bi,bj) *
62         &        recip_rA(i,j,bi,bj) *
63         &        U_streamice(i+k,j+l,bi,bj)
64                 vmid = vmid + 0.25 *
65         &        dxG(i,j+l,bi,bj)*dyG(i+k,j,bi,bj) *
66         &        recip_rA(i,j,bi,bj) *
67         &        V_streamice(i+k,j+l,bi,bj)
68                ENDDO
69               ENDDO
70      
71             streamice_taubx (i,j,bi,bj) =             streamice_taubx (i,j,bi,bj) =
72       &      tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)       &      tau_beta_eff_streamice(i,j,bi,bj) * umid
73             streamice_tauby (i,j,bi,bj) =             streamice_tauby (i,j,bi,bj) =
74       &      tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)       &      tau_beta_eff_streamice(i,j,bi,bj) * vmid
75                        
76            ENDIF            ENDIF
77           ENDDO           ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22