问题描述
我需要执行两个积分，以便基于密度分布和分布之间的宽度来获取区域内水的累积传输。
首先，速度计算为；
请注意，V在表面处可以为非零，并且在深度H处必须为零。
第二，运输量计算为；
现在，我已经使用numpy.ma.cumsum编码了上述方程式
# Compute veLocity
V = constant*numpy.ma.cumsum((density/width*dz)[::-1],axis=0)[::-1]
# Compute transport in sverdrup (10**6 m/s^2)
T = width * numpy.ma.cumsum((V[::-1]-1/H*numpy.ma.cumsum((V*dz)[::-1],axis=0))*dz[::-1],axis=0)[::-1] / 10**6
这些方程式给我以下结果；
如果您对如何强制运输工具T在地面和深度H处为零有任何想法，请告诉我。
这是我的数据集：
constant = -100
density = [0.08162372168703769,-0.04471999854399655,-0.3110550684352802,-0.4783609962985338,-0.2606431434960541,-0.1734227201536669,-0.20306011069737906,-0.22652204210817217,-0.24526932281150948,-0.2899378668189456,-0.2968024166812029,-0.30554869432489795,-0.3145306210531089,0.04810820373199931,0.08374671599813155,0.11754102251029508,0.14828751252093753,0.17433953329941687,0.19374034901125015,0.20420143703017857,0.2026197673765182,0.1848936664282519,0.131398194555004,0.3103135616245254,0.3938907065328294,0.4770966560975012,0.5276251000107095,0.48519935265176173,0.41467987590567645,0.34373320113900263,0.12321804603743658,0.05310501021676828,0.006689604194434651,-0.02221187988652673,-0.03928717776148005,-0.06361049772976912,-0.06657255480854474,-0.07648006375347904,-0.0671293538559894,-0.04341857793974668,-0.04089445488693855,-0.025967806515836855,0.010263609448429634,0.016120765154255423,0.016567011411552812,0.02229294316521191,0.02531710921425656,0.027188241842850402,0.027440846175068145,--]
width = [5370872.025103143,5370872.025103143,5286952.149710908,5203032.274318671,5035192.523534197,4867352.772749724,4783432.897357487,4699513.02196525,4615593.146573015,4531673.271180778,4447753.395788541,4363833.520396304,4279913.645004068,4195993.769611831,4028154.018827358,3776394.3926506485,--]
dz = [-5.03355,-10.0671,-10.118699999999999,-10.139100000000006,-10.217899999999986,-10.276900000000012,-10.40849999999999,-10.541100000000029,-10.774099999999947,-11.046300000000059,-11.472899999999953,-12.010500000000036,-12.801899999999975,-13.838099999999997,-15.307100000000048,-17.25869999999992,-19.95050000000009,-23.514299999999935,-28.290900000000022,-34.48270000000002,-42.456899999999905,-52.35730000000012,-64.36129999999991,-78.2561,-93.72070000000008,-109.98889999999983,-126.30150000000026,-141.72169999999983,-155.66089999999986,-167.62050000000045,-177.5652999999993,-185.49530000000095,-191.7308999999991,-196.46390000000042,-200.07110000000011,-202.71990000000005,-204.71409999999923,-206.14090000000124,-207.21909999999843,-207.96890000000167,-208.54709999999795,-208.93390000000272,-209.24569999999676,-209.43890000000374,-209.61189999999624,-209.70250000000306,-209.80249999999796,-209.84110000000146,-209.90109999999913,--]
H = 5185.11075
integral of dx = width
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）