_bc(state, bc_type)      type(state_type), intent(in)     :: state       character(len=*), intent(in)     :: bc_type      type(vector_field), pointer      :: positions      real                             :: gravity_magnitude     integer                          :: i     real, allocatable, dimension(:)  :: z0s, z0b, u_taus_squared, u_taub_squared     