Salvus Mesher

[Saikiran Tharimena]

Hi,

I am new to Salvus and have been trying out the tutorials to familiarise with the code. I managed to successfully run mesher tutorials but got stuck in tutorial 5 (https://salvus.io/tutorials/0_Salvus_Mesher/5_Topo_Mesh/). I tried to input the topography file but it doesn't seem to work. In the YAML file that I generated, it says "Filename for topography, generated with salvus_topo: https://gitlab.com/Salvus/salvus_topo"

This page seems to be unavailable. I did however try running the YAML file with a path to the ETOPO1_Ice topography file from NOAA. But I get an error:

"KeyError: Unable to open object (object 'longitude_ETOPO1_Ice__lmax_10800' doesn't exist)"

I used the ETOPO1 cell registered grid with longitude (-180 180) and latitude (-90 90) dimensions of 21600 x 10800. The grid comes in the GMT4 ".grd" format. When I used the ".grd" file I got an error that it didn't recognise the file (OSError: Unable to open file (file signature not found).

I then used GMT GRDCONVERT to convert from ".grd" to 32-bit integer, COARDS, CD-1.5 NETCDF format. Since this didn't work either, I modified the NETCDF variable and attribute names so that it has latitude and longitude as variable names rather than lat/lon and also changed file name to "ETOPO1.nc". This doesn't seem to work either - I get the unable to open object error.

Any idea as to what I might be doing wrong? Should the topography file be in a specific NETCDF format?

Many thanks for your help.

Cheers,

Sai.

[Lion Krischer]

Hi Sai,

`salvus_topo` is a package which we use to preprocess any internal or external topographies. It is currently non-public because of dependency reasons - we'll try to fix this.

In any case - here are two topography files for you to download:

* Earth's topography at an about 4km resolution (Earth2014) - 3 GB: https://www.dropbox.com/s/tswwzkpaugdmmvl/topography_wgs84_filtered.nc?dl=1

* Moho topography (Crust 1.0) - 7 MB: https://www.dropbox.com/s/l72y5255u196pdy/moho_topography_wgs84_filtered.nc?dl=1

These are files with the netcdf CF conventions and they can be visualized using paraview. Please remember to adapt the `lmax_` settings for each in the mesher's YAML input to your domain size and mesh period.

Cheers!

Martin + Lion

--
You received this message because you are subscribed to the Google Groups "Salvus" group.
To unsubscribe from this group and stop receiving emails from it, send an email to salvus+un...@googlegroups.com.
To post to this group, send email to sal...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/salvus/ee0c5164-4a20-44d5-85a8-15ca9798ba9b%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Saikiran Tharimena]

Hi Lion / Martin,

Thanks for the quick reply. Well it makes sense. Thanks for sharing the topography files, I'll work with these.

Cheers,

Sai.