Coarsening horizontal particle resolution currently unfeasible for high resolution

pwolfram / MPAS-Model

Repository for MPAS models and shared framework releases.

Other

1 stars 1 forks source link

Coarsening horizontal particle resolution currently unfeasible for high resolution #19

Closed bradyrx closed 5 years ago

bradyrx commented 5 years ago

The current use of pyAMG in make_particle_file.py to coarsen particle seeding in the horizontal is unfeasible for high resolution meshes. From memory, coarsening by 1 level in the 60to30km case takes ~45min.

I tried to coarsen 1 level in the 30to10km case (i.e. ./make_particle_file.py -i input.nc -g graph.9600 -p 9600 -o particles.nc --spatialfilter SouthernOceanXYZ -t 'passive','surface' -d 1) and the script is still running ~20 hours later. Is this an inherent limitation of pyAMG?

bradyrx commented 5 years ago

Closed to avoid duplicate work from https://github.com/pwolfram/MPAS-Model/issues/20. It seems like https://github.com/pwolfram/MPAS-Model/issues/20 relates directly to this issue.