← Back to branch summary

~fluidity-core/fluidity/sea-ice-branch

« back to all changes in this revision

Viewing changes to feseaice/ice_evp.F90

  • Committer: Simon Mouradian
  • Date: 2012-04-04 21:26:15 UTC
  • Revision ID: mouradian@gmail.com-20120404212615-2kc4roq1tvwc32tf
introduce viscosity limiting

Show diffs side-by-side

added added

removed removed

Lines of Context:
feseaice/ice_evp.F90
70
70
     pressure=0.5_8*pressure
71
71
     ! ===== viscosity zeta should exceed zeta_min
72
72
     ! (done via limiting delta from above)
73
 
     !if(delta>pressure/zeta_min) delta=pressure/zeta_min
 
73
     if(delta>pressure/zeta_min) delta=pressure/zeta_min
74
74
     ! ===== if viscosity is too big, it is limited too
75
75
     ! (done via correcting delta_inv)
76
76
     delta_inv=1.0_8/max(delta,delta_min)