Integrating across 2 grids

I'm having issues with the integral function requiring floating point scalars instead of gridded data. In the code below, I wish to integrate between sea_surface and january_mld_int (everything is gridded in 38X1388 arrays)

%% Mixed Layer Freshwater Function (prior to integrating across mixed layer)
sea_surface=zeros(38,1388); %makes a grid of zeroes to integrate from/to
jan_fun=@(x)((january_ref_sal-january_mld_int)/january_ref_sal)*x;  %makes the formulae as defined in parampil 2010 jan_fwf=integral(jan_fun,sea_surface,january_mld_int,'ArrayValued',true);%attempts to integrate

where january_mld_int is the interpolated MLD data, and jan_ref_sal is the reference salinity (both 38x1388 grids)

However, Integral requires floating point vectors, not arrays!

Any help for this issue?

[darova]

Can't you add zeros after dot?

2.0

3.0

...

[darova]

Can't you add zeros after dot?

2.0

3.0

...