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

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_visc_beta_hybrid.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.7 by dgoldberg, Wed Aug 27 19:29:15 2014 UTC
# Line 36  C     !INPUT/OUTPUT ARGUMENTS Line 36  C     !INPUT/OUTPUT ARGUMENTS
36  #ifdef STREAMICE_HYBRID_STRESS  #ifdef STREAMICE_HYBRID_STRESS
37    
38  C     LOCAL VARIABLES  C     LOCAL VARIABLES
39        INTEGER i, j, bi, bj, k, l, umid, vmid, m        INTEGER i, j, bi, bj, k, l, m
40        INTEGER ikey_1        INTEGER ikey_1
41        _RL ux, uy, vx, vy, exx, eyy, exy, unorm, second_inv        _RL ux, uy, vx, vy, exx, eyy, exy, unorm, second_inv
42        _RL ub, vb, fb, mean_u_shear, mean_v_shear        _RL ub, vb, fb, mean_u_shear, mean_v_shear, umid, vmid
43        _RL omega_temp (Nr+1), u_shear(Nr+1), v_shear(Nr+1)        _RL omega_temp (Nr+1), u_shear(Nr+1), v_shear(Nr+1)
44    
45        _RL STREAMICE_BSTRESS_EXPONENT        _RL STREAMICE_BSTRESS_EXP
46  !       _RL total_vol_out  !       _RL total_vol_out
47        external STREAMICE_BSTRESS_EXPONENT        external STREAMICE_BSTRESS_EXP
48    
49        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
50         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 133  CADJ &     = comlev1_stream_hybrid, key= Line 133  CADJ &     = comlev1_stream_hybrid, key=
133    
134              second_inv =              second_inv =
135       &       sqrt(exx**2+eyy**2+exx*eyy+exy**2+eps_glen_min**2+       &       sqrt(exx**2+eyy**2+exx*eyy+exy**2+eps_glen_min**2+
136       &            streamice_vert_shear_uz(m)**2 +       &            0.25 * streamice_vert_shear_uz(m)**2 +
137       &            streamice_vert_shear_vz(m)**2)       &            0.25 * streamice_vert_shear_vz(m)**2)
138    
139    #ifdef STREAMICE_3D_GLEN_CONST
140              visc_streamice_full(i,j,m,bi,bj) =              visc_streamice_full(i,j,m,bi,bj) =
141       &       .5 * A_glen(i,j,bi,bj)**(-1./n_glen) *       &       .5 * B_glen(i,j,m,bi,bj)**2 *
142       &        second_inv**((1-n_glen)/n_glen)       &        second_inv**((1-n_glen)/n_glen)
143    #else
144                visc_streamice_full(i,j,m,bi,bj) =
145         &       .5 * B_glen(i,j,bi,bj)**2 *
146         &        second_inv**((1-n_glen)/n_glen)
147    #endif
148    
149              visc_streamice (i,j,bi,bj) = visc_streamice (i,j,bi,bj) +              visc_streamice (i,j,bi,bj) = visc_streamice (i,j,bi,bj) +
150       &       H_streamice(i,j,bi,bj) * streamice_delsigma (m) *       &       H_streamice(i,j,bi,bj) * streamice_delsigma (m) *
# Line 183  CADJ &     = comlev1_stream_hybrid, key= Line 189  CADJ &     = comlev1_stream_hybrid, key=
189    
190             ub = umid - streamice_taubx(i,j,bi,bj) *             ub = umid - streamice_taubx(i,j,bi,bj) *
191       &      streamice_omega(i,j,bi,bj) / H_streamice(i,j,bi,bj)       &      streamice_omega(i,j,bi,bj) / H_streamice(i,j,bi,bj)
192    
193               streamice_u_bed (i,j,bi,bj) = ub
194    
195             vb = vmid - streamice_tauby(i,j,bi,bj) *             vb = vmid - streamice_tauby(i,j,bi,bj) *
196       &      streamice_omega(i,j,bi,bj) / H_streamice(i,j,bi,bj)       &      streamice_omega(i,j,bi,bj) / H_streamice(i,j,bi,bj)
197    
198             unorm = sqrt(ub**2+vb**2+eps_glen_min**2)             streamice_v_bed (i,j,bi,bj) = vb
199    
200               unorm = sqrt(ub**2+vb**2+eps_u_min**2)
201    
202             fb = C_basal_friction(i,j,bi,bj)**2 *             fb = C_basal_friction(i,j,bi,bj)**2 *
203       &      STREAMICE_BSTRESS_EXPONENT (unorm,n_basal_friction) *       &      STREAMICE_BSTRESS_EXP (unorm,n_basal_friction) *
204       &      streamice_basal_geom(i,j,bi,bj)       &      streamice_basal_geom(i,j,bi,bj) *
205         &      float_frac_streamice(i,j,bi,bj)
206    
207             tau_beta_eff_streamice(i,j,bi,bj) =             tau_beta_eff_streamice(i,j,bi,bj) =
208       &       fb /       &       fb /
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22