anchor.diff:              endif
           endif
       end subroutine
  -    
  +
       !>----------------------------------------------------------
       !!  Compute sedimentation in a column
  -    !! 
  +    !!
       !!  Takes a mixing ratio of some species along with the fall velocities,
  -    !!  air densities and layer thicknesses.  
  -    !! 
  +    !!  air densities and layer thicknesses.
  +    !!
       !!  @param  q   1D microphysical species to sediment (snow, rain, etc) [kg/kg]
       !!  @param  v   1D vertical velocity [m/s]
       !!  @param  rho 1D air density [kg/m^3]
       !!  @param  dz  1D layer thickness [m]
  -    !!  @param  n   0D number of layers 
  +    !!  @param  n   0D number of layers
       !!
       !!----------------------------------------------------------
       real function sediment(q,v,rho,dz,n)
  
positive.code:         real,intent(inout) :: qc,q
          real,intent(in) :: conversion, qcmin
          real::delta
          
          if (qc > qcmin) then
              delta = qc-(qc*conversion)
          else
              delta = 0
          endif
          
          if (delta<qc) then
              qc = qc-delta
              q = q+delta
  
negative.code:     !!  @param qc           cloud water mixing ratio [kg/kg]
      !!  @param qvsat        saturated mixing ratio   [kg/kg]
      !!  @param dt           time step                [sec.]
      !!
      !!----------------------------------------------------------
      subroutine cloud_conversion(pressure,temperature,qv,qc,qvsat,dt)
          implicit none
  