Just a quick query. Now that I have the plot generated with the wavenumber (k) on the x-axis (bottom), is it possible to have a different scale on the top x-axis which is the wavelength? So we would have a linear wavenumber scale on the bottom and the corresponding wavelength on the top:

wavelength = 2*Pi / wavenumber(k)

I can alter the properties of the graph from the menu options, but the secondary axis only applies to Y and not X

Will be interested to know how this would work and if it is easy or tricky to do.

Cheers

FAA_isostatic1e.mw

The plots now all work! Also spotted that I needed a different flexp parameter for the Zsubload, so put that in and the plot scale was all sorted. Bit fiddly but got there in the end :) Thanks for the help on this, still pretty new with using Maple so it takes a little getting used to.

The equations relate to different isostatic models of the Earth's lithosphere (the Crust and uppermost mantle) to see how they vary with wavelength and the admittance function (Z(k) )=fft(gravity)/fft(topography). The actual equations simplify to more usable structures.

I am a geologist by training and currently completing a PhD in geodynamics.

The equations are found in the book by Tony Watts - Isostasy and flexure of the lithosphere (2001)

Have some more complex functions to work on in Maple so will give those a try later on. Early days with Maple and I am sure I will have more queries as I learn it more.

Thanks again for the help

I added the changes you spotted, but the last plot is incorrect compared to the SCILAB output, so need to investigate that. Otherwise so far so good it is all getting there. One thing I needed to remember was that restart clears the variables and so they need to be reinitialised.

Spotted a typo DC should be 113700m (113.7km) - still does not solve the problem