anchor.diff:                  Kn = lambda/rwet
                   GG = 1.+ Kn*(A+B*exp(-C/Kn))
                   vc = terminal_vel(rwet,pdn,adn(k,i,j),avis,GG,flag)
  -                
  +
  +
  +                ! This algorithm breaks down if the fall velocity is large enough to make the fall 
  +                ! distance greater than the grid level thickness (mainly an issue with very high 
  +                ! vertical resolutions). As a simple solution, limit the fall velocity based on grid 
  +                ! box thickness and timestep (should not cause any further issues).
  +                vc = MIN( vc, MIN(0.5*(1./dzt(k))/dt, 2.0) )
   
                   IF (k==2) THEN ! The level just above surface
                       ! Particle diffusitivity  (15.29) in jacobson book
  
positive.code: 
        IF (sed_ice) THEN
  
         CALL DepositionSlow(n1,n2,n3,n4,nice,tk,a_dn,spec%rhoic,ustar,nicep,micep,dzt,nlim,indiv,imdiv,indep,remice,4)
  
         nicet = nicet - indiv 
         micet = micet - imdiv 
  
negative.code: 
  
  
    SUBROUTINE DepositionSlow(n1,n2,n3,n4,nn,tk,adn,pdn,ustar,numc,mass,dzt,dt,clim,flxdivn,flxdivm,depflxn,depflxm,flag)
      IMPLICIT NONE
  
      INTEGER, INTENT(in) :: n1,n2,n3,n4       ! Grid numbers, number of chemical species
  