DISPATCH
read_snapshot_mod Module Reference

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
 

Detailed Description

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.