DISPATCH
|
IC module for reading in legacy Stagger snapshots into existing patches. More...
Data Types | |
type | rs_t |
Functions/Subroutines | |
subroutine | read_stagger_snap (task_list) |
Read a legacy Stagger raw data file, with auxiliary grid file. | |
Variables | |
character(len=64), save | meshfile ='mhd63.msh' |
type(rs_t), public | read_snapshot |
IC module for reading in legacy Stagger snapshots into existing patches.
The purpose with this separate module is to avoid having to load the entire snapshot on each MPI rank – there isn't enough memory to do that on the largest experiments from the archive. In this module one can loop over either variable index iv, or even iv and depth index, reading only one xy-slice at a time.
To save startup timme, separate OMP tasks are spawned for each patch, to work in parallel on the interpolation of the ICs to the new patch mesh.
The separate reader also saves memory because it operates before the RT solver has started to run, allocating additional scratch memory.