Turning on field-dependent Dirichlet BC when Dirichlet data are passed from
CISM.  The Dirichlet boundary is defined according to a mask passed from CISM and
the Dirichlet values are set to the initial u and v velocities passed from the CISM *.nc file.

Thanks to Mauro for adding "dirichlet_field" to FELIX to enable this!