Commit 0e7f3570 authored by Pierre Cazenave's avatar Pierre Cazenave

Fix the identification of a uniform sponge layer width

parent e36b849b
......@@ -84,7 +84,7 @@ Mobj.sponge_nodes(Mobj.nSponge,1:npts) = SpongeList;
Mobj.sponge_name{Mobj.nSponge} = SpongeName;
Mobj.sponge_fac(Mobj.nSponge) = SpongeCoeff;
if max(size(SpongeRadius))==1 % if you have a constant sponge radius
if numel(unique(SpongeRadius)) == 1 % if you have a constant sponge radius
Mobj.sponge_rad(Mobj.nSponge) = SpongeRadius;
else % if you have a variable sponge radius
Mobj.sponge_rad(Mobj.nSponge,1:npts) = SpongeRadius;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment