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

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

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

revision 1.3 by dgoldberg, Tue Sep 4 21:11:44 2012 UTC revision 1.7 by dgoldberg, Sat Jun 8 22:15:34 2013 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 92  C     LOCAL VARIABLES Line 93  C     LOCAL VARIABLES
93    
94    
95             visc_streamice(i,j,bi,bj) =             visc_streamice(i,j,bi,bj) =
96       &      .5 * A_glen(i,j,bi,bj)**(-1./n_glen) *       &      .5 * (B_glen(i,j,bi,bj))**2 *
97       &      second_inv**((1-n_glen)/n_glen) * H_streamice(i,j,bi,bj)       &      second_inv**((1-n_glen)/n_glen) * H_streamice(i,j,bi,bj)
98    
99             unorm = sqrt(umid**2+vmid**2+eps_glen_min**2)             unorm = sqrt(umid**2+vmid**2+eps_u_min**2)
100             tau_beta_eff_streamice(i,j,bi,bj) =             tau_beta_eff_streamice(i,j,bi,bj) =
101       &       C_basal_friction(i,j,bi,bj)**2 *       &       C_basal_friction(i,j,bi,bj)**2 *
102       &       unorm ** (n_basal_friction-1) *       &       unorm ** (n_basal_friction-1.0) *
103       &       float_frac_streamice (i,j,bi,bj)       &       float_frac_streamice (i,j,bi,bj)
104                        
105            ENDIF            ENDIF
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22