Hi Michael
You should first interpolate your runoff data onto your mesh. Assuming that the data you have is on a regular grid on and uses the same projection as your ISSM model, you can do this:
%Allocate matrix
runoff = zeros(md.mesh.numberofvertices,365);`
%Interpolate data on model mesh
for i=1:365
runoff(:,i) = InterpFromGridToMesh(X(1,: ), Y(:,1), data,md.mesh.x, md.mesh.y, 0);
end
%Construct time vector (daily data in yr)
time = linspace(0, 1, 365);
%Insert in model as melting rate
md.basalforcings.groundedice_melting_rate = [runoff; time];
%If you use GlaDS, you can do
md.hydrology = hydrologyglads();
md.hydrology.melt_flag = 1; %Tell GlaDS to use md.basalforcings.groundedice_melting_rate
%Force time stepping to be daily? (you probably need shorter time steps though...)
md.timestepping.start_time=0.;
md.timestepping.final_time=1.; %1 year
md.timestepping.time_step=1./365; %Daily time step
I hope this helps!
Cheers
Mathieu