<div dir="ltr"><div><div><div><div>I am not sure about the logic with some of these Mathcad files, have a number to deal with in various forms.<br><br></div>The result for calc_grv is invalid because it should compute all values on the profiel to npts<br>
<br></div>The PDF I loaded before shows what it is meant to generate<br><br></div><div>Checked the loading and plotting of the inputs:<br><br><pre style="font-family:Monospaced;font-style:normal"><span style="color:rgb(100,174,100);font-style:italic">//Read in topography and free-air gravity files</span>
<span style="color:rgb(0,0,0)">npts</span><span style="color:rgb(92,92,92)">=</span><span style="color:rgb(188,143,143)">128</span>
<span style="color:rgb(0,0,0)">Topo</span> <span style="color:rgb(92,92,92)">=</span> <span style="color:rgb(50,185,185)">read</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">C:\Users\Lester\Documents\Scilab_functions\Geoscience\</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(92,92,92)">+</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">haw_topo.prn</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(92,92,92)">-</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">2</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(0,0,0)">Grav</span> <span style="color:rgb(92,92,92)">=</span> <span style="color:rgb(50,185,185)">read</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">C:\Users\Lester\Documents\Scilab_functions\Geoscience\</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(92,92,92)">+</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">haw_grv.prn</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(92,92,92)">-</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">2</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(0,0,0)">Topo</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">=</span><span style="color:rgb(0,0,0)">Topo</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">*</span><span style="color:rgb(188,143,143)">1000</span><span style="color:rgb(0,0,0)">;</span> 
<span style="color:rgb(0,0,0)">Grav</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">=</span><span style="color:rgb(0,0,0)">Grav</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">*</span><span style="color:rgb(188,143,143)">1000</span><span style="color:rgb(0,0,0)">;</span>

<span style="color:rgb(0,0,0)">dx</span> <span style="color:rgb(92,92,92)">=</span> <span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">Grav</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">$</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">-</span><span style="color:rgb(0,0,0)">Grav</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">/</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">npts</span><span style="color:rgb(92,92,92)">-</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(0,0,0)">XKint</span> <span style="color:rgb(92,92,92)">=</span> <span style="color:rgb(188,143,143)">2</span><span style="color:rgb(92,92,92)">*</span><span style="color:rgb(218,112,214)">%pi</span><span style="color:rgb(92,92,92)">/</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">npts</span><span style="color:rgb(92,92,92)">-</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">*</span><span style="color:rgb(0,0,0)">dx</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(0,0,0)">Xint</span> <span style="color:rgb(92,92,92)">=</span> <span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">Grav</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">$</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">-</span><span style="color:rgb(0,0,0)">Grav</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">/</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">npts</span><span style="color:rgb(92,92,92)">-</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>

<span style="color:rgb(160,32,240)">for</span> <span style="color:rgb(0,0,0)">i</span><span style="color:rgb(92,92,92)">=</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">npts</span>
    <span style="color:rgb(0,0,0)">x1</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(74,85,219)">)</span> <span style="color:rgb(92,92,92)">=</span> <span style="color:rgb(0,0,0)">Topo</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">+</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">i</span><span style="color:rgb(92,92,92)">-</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(92,92,92)">*</span><span style="color:rgb(0,0,0)">Xint</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(160,32,240)">end</span>

<span style="color:rgb(0,0,0)">obs_topo</span> <span style="color:rgb(92,92,92)">=</span> <span style="color:rgb(174,92,176);text-decoration:underline">interp1</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">Topo</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(0,0,0)">Topo</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">2</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(0,0,0)">x1</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">linear</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(0,0,0)">obs_grav</span> <span style="color:rgb(92,92,92)">=</span> <span style="color:rgb(174,92,176);text-decoration:underline">interp1</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">Grav</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(0,0,0)">Grav</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">2</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(0,0,0)">x1</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">linear</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>

<span style="color:rgb(0,0,0)">k</span> <span style="color:rgb(92,92,92)">=</span> <span style="color:rgb(188,143,143)">0</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(255,170,0)">:</span><span style="color:rgb(0,0,0)">npts</span><span style="color:rgb(92,92,92)">-</span><span style="color:rgb(188,143,143)">1</span><span style="color:rgb(0,0,0)">;</span>

<span style="color:rgb(100,174,100);font-style:italic">//Plot input files - ignore labels/scale</span>
<span style="color:rgb(174,92,176);text-decoration:underline">scf</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(188,143,143)">0</span><span style="color:rgb(74,85,219)">)</span>
<span style="color:rgb(174,92,176);text-decoration:underline">clf</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(188,143,143)">0</span><span style="color:rgb(74,85,219)">)</span>
<span style="color:rgb(174,92,176);text-decoration:underline">xlabel</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(188,143,143)">Distance (km)</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline">ylabel</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(188,143,143)">gravity</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span> 
<span style="color:rgb(50,185,185)">set</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(174,92,176);text-decoration:underline">gca</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(188,143,143)">auto_clear</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(188,143,143)">off</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(74,85,219)">)</span>
<span style="color:rgb(0,0,0)">a</span><span style="color:rgb(92,92,92)">=</span><span style="color:rgb(174,92,176);text-decoration:underline">gca</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(0,0,0)">a</span><span style="color:rgb(92,92,92)">.</span><span style="color:rgb(170,170,170)">x_location</span><span style="color:rgb(92,92,92)">=</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(188,143,143)">bottom</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(0,0,0)">a</span><span style="color:rgb(92,92,92)">.</span><span style="color:rgb(170,170,170)">axes_reverse</span><span style="color:rgb(92,92,92)">=</span><span style="color:rgb(74,85,219)">[</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(188,143,143)">off</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(188,143,143)">off</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(188,143,143)">off</span><span style="color:rgb(188,143,143)">"</span><span style="color:rgb(74,85,219)">]</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline">plot</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">k</span><span style="color:rgb(92,92,92)">.*</span><span style="color:rgb(0,0,0)">Xint</span><span style="color:rgb(92,92,92)">/</span><span style="color:rgb(188,143,143)">1000</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(0,0,0)">obs_grav</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">m</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline">plot</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(0,0,0)">k</span><span style="color:rgb(92,92,92)">.*</span><span style="color:rgb(0,0,0)">Xint</span><span style="color:rgb(92,92,92)">/</span><span style="color:rgb(188,143,143)">1000</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(0,0,0)">obs_topo</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">b</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline">legend</span><span style="color:rgb(74,85,219)">(</span><span style="color:rgb(74,85,219)">[</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">bathymetry</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(0,0,0)">;</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(188,143,143)">gravity</span><span style="color:rgb(188,143,143)">'</span><span style="color:rgb(74,85,219)">]</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(188,143,143)">4</span><span style="color:rgb(74,85,219)">)</span><span style="color:rgb(0,0,0)">;</span><span style="color:rgb(0,0,0)"></span></pre>
<br></div><div>The plots of the input data (attached) show the data is read in fine. The issue is somewhere in the computation, which at least narrows the problem.<br><br></div>Cheers<br></div>Lester<br></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On 26 August 2013 12:40, Antoine Monmayrant <span dir="ltr"><<a href="mailto:antoine.monmayrant@laas.fr" target="_blank">antoine.monmayrant@laas.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">On 26/08/2013 11:03, Lester Anderson wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Antoine,<br>
<br>
The prn files are just plain text files, no idea why they were given prn as<br>
an extension, that's how they came with the Mathcad code. They can be<br>
loaded into Notepad with no problem)<br>
<br>
The complex issue is coming from the forward and inverse fft as you say. I<br>
am kind of suspicious if it is actually loading the data files properly?<br>
<br>
c = fft(obs_topo1); =  COMPLEX<br>
k = 0:1:npts-1;<br>
phik = (D*(k*XKint)^4/(g*(rho_m-rho_<u></u>c))+1)^-1;Zflexure =<br>
(2*%pi*G*1e5)*(rho_c-rho_w).*<u></u>exp(-k*XKint*mean_depth).*(1-<u></u>phik.*exp(-XKint*ocean_thick))<u></u>;<br>
factork = Zflexure * c; = COMPLEX (single value eg<br>
-4.490857598202565E-10-8.<u></u>427771813321494E-11*%i)<br>
</blockquote></div>
Well, I don't know exactly what this means :<div class="im"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
calc_grv = ifft(factork) = COMPLEX - this should be a real set of values.<br>
</blockquote></div>
factork is a constant ie a 1x1 matrix and I don't see what it means to take the fft or the ifft of a constant!<br>
<br>
Antoine<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
<br>
<br>
Cheers<br>
Lester<br>
<br>
<br>
<br>
<br>
On 26 August 2013 09:24, Antoine Monmayrant <<a href="mailto:antoine.monmayrant@laas.fr" target="_blank">antoine.monmayrant@laas.fr</a>><u></u>wrote:<br>
<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
  On 25/08/2013 12:45, Lester Anderson wrote:<br>
<br>
Hi all,<br>
<br>
I am attempting to code up a Mathcad file into Scilab and it should be<br>
straightforward, but for some reason it is not working as it should. And<br>
checking the variable values, some are still complex as output from the<br>
fft/ifft runs - the inputs look right as far as I can see.<br>
<br>
  I cannot open the prn files and it's thus difficult for me to check your<br>
data.<br>
What do you mean by still complex?<br>
When you run fft or ifft on a real dataset, the output is complex, that's<br>
to be expected.<br>
Maybe (it's just a wild guess) you ended up with complex values that have<br>
an imaginary part that is zero.<br>
Scilab still consider these values to be complex.<br>
<br>
Antoine<br>
<br>
<br>
If someone can see what the issue is and a solution, that would be great.<br>
<br>
Attached my quick attempt at converting the Mathcad code (see PDF) and the<br>
input profile data (*.prn files)<br>
<br>
Cheers<br>
<br>
Lester<br>
<br>
<br>
<br>
______________________________<u></u>_________________<br></div>
users mailing listusers@lists.scilab.<u></u>orghttp://<a href="http://lists.scilab.org/mailman/listinfo/users" target="_blank">lists.scilab.org/<u></u>mailman/listinfo/users</a><div class="im"><br>
<br>
<br>
<br>
<br>
______________________________<u></u>_________________<br>
users mailing list<br>
<a href="mailto:users@lists.scilab.org" target="_blank">users@lists.scilab.org</a><br>
<a href="http://lists.scilab.org/mailman/listinfo/users" target="_blank">http://lists.scilab.org/<u></u>mailman/listinfo/users</a><br>
<br>
<br>
</div></blockquote></blockquote>
<br>
<br>
</blockquote></div><br></div>