/[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.6 by dgoldberg, Thu Mar 7 15:23:19 2013 UTC revision 1.7 by dgoldberg, Wed Aug 27 19:29:15 2014 UTC
# Line 42  C     LOCAL VARIABLES Line 42  C     LOCAL VARIABLES
42        _RL ub, vb, fb, mean_u_shear, mean_v_shear, umid, vmid        _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 136  CADJ &     = comlev1_stream_hybrid, key= Line 136  CADJ &     = comlev1_stream_hybrid, key=
136       &            0.25 * streamice_vert_shear_uz(m)**2 +       &            0.25 * streamice_vert_shear_uz(m)**2 +
137       &            0.25 * 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) =
141         &       .5 * B_glen(i,j,m,bi,bj)**2 *
142         &        second_inv**((1-n_glen)/n_glen)
143    #else
144              visc_streamice_full(i,j,m,bi,bj) =              visc_streamice_full(i,j,m,bi,bj) =
145       &       .5 * B_glen(i,j,bi,bj)**2 *       &       .5 * B_glen(i,j,bi,bj)**2 *
146       &        second_inv**((1-n_glen)/n_glen)       &        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               streamice_v_bed (i,j,bi,bj) = vb
199    
200             unorm = sqrt(ub**2+vb**2+eps_u_min**2)             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)       &      float_frac_streamice(i,j,bi,bj)
206    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22