[aspect-devel] Trying to customize dynamic topography

[Austermann, Jacqueline]

Hi Felipe,

Glad it works now! If you add dynamic topography to your list of visualization output variables it will not produce an additional file in the output folder but instead be added as an output variable to your visualization. That means if you open your solution with paraview (or visit) you should now be able to select dynamic topography (in addition to the default p and T) and plot it for your model. Not sure if that was your problem, let me know if I misunderstood.

Cheers,
Jacky



On May 8, 2017, at 4:52 PM, FELIPE ORELLANA ROVIROSA <f_orellana at berkeley.edu<mailto:f_orellana at berkeley.edu>> wrote:


 Hi Jacqueline,

     I followed your tips somehow. It works !  Thank you !
     I get the dynamic topography output files (at each time step) hopefully subjected to a non-zero density above. I think I had misplaced 'end' instructions, somehow the subsection dyn topo wasn't being ended properly.

       Regarding visualization: I am using prescribed velocities at the surface, therefore it is not free, and that's why I am using dynamic topography. Beyond this, I don't know if visualization produces the equivalent topography graphic output for a case like this. As of now, running this with and without visualization of dynamic topography, I do not see any difference in the output, nothing was added (or appended) in the output directory. I don't know if this is sort of a bug or is because dynamic topography does not have a proper graphic output.

cheers,
Felipe

On Sun, May 7, 2017 at 5:57 AM, Austermann, Jacqueline <jaustermann at fas.harvard.edu<mailto:jaustermann at fas.harvard.edu>> wrote:
Hi Felipe,

   i) I want to use a custom value for 'Density above', but I cannot do it. I have tried writing the instruction ''      set Density above = 1000   '' at different places in the section Postprocess.

   I have also tried doing it creating '' subsection Dynamic Topography
                                                            set Density above = 1000      ''

 or writing it as Visualization/Dynamic Topography..

 I also browsed github at no avail.

     where and how should this variable and instruction be properly written?

The dynamic topography is calculated twice, once for the text output and once for the visualization output. So if you want the density above to be 1000 in both I believe you’ll have to do something like this in your parameter file:

subsection Postprocess
set List of postprocessors = dynamic topography, visualization
subsection Dynamic Topography
set Density above = 1000
end

   subsection Visualization
    set List of output variables      = dynamic topography
     subsection Dynamic Topography
       set Density above = 1000
     end
   end
end

You can of course combine that with any other postprocessor. Give that a try.


ii) Is there any way to prevent ASPECT from producing topography files at every step, and instead to force it to do it with a stride, every certain number of steps?

Unfortunately, I don’t think there is an option to do that yet. You’d have to implement it.

Hope that helps!

Best,
Jacky


_______________________________________________
Aspect-devel mailing list
Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel<https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.geodynamics.org_cgi-2Dbin_mailman_listinfo_aspect-2Ddevel&d=CwMFaQ&c=WO-RGvefibhHBZq3fL85hQ&r=j5AqZvMsoErn2L-vpbdTErPRtyT4BhzQUwKzsenbbTc&m=62S9HMjWTs1QpoBFgye2AqPRYEixtIpR4ay-jJ8r9-M&s=OohLqsmkxvizSJZuXK29s07YYDRtwBI3MngWQ_2f1fc&e=>

_______________________________________________
Aspect-devel mailing list
Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20170509/5f084151/attachment-0001.html>