     if (have_option("/material_phase[0]/subgridscale_parameterisations/GLS/relax_diffusivity")) then         call get_option("/material_phase[0]/subgridscale_parameterisations/GLS/relax_diffusivity", relax