Integrating across 2 grids

3 次查看（过去 30 天）

Tim Holgate 2021-2-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/750214-integrating-across-2-grids

评论： darova 2021-2-20

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?