Skip to content

Allowing for a selection of FEEC coefficients saving

Byung Kyu Na requested to merge improve_post_processing into devel

Solve issue #194 (closed)

We can save FEEC data selectively by inserting these lines at parameters.yml

em_fields :
    init :
        type : [ModesSin, ModesCos] # type of initial condition (possible types seen below)
        ModesSin : 
            coords : 'logical' # in which coordinates (logical or physical)
            comps :
                e1 : [True, False, False]  # components to be initialized (for scalar fields: no list)
                b2 : [False, False, False] # components to be initialized (for scalar fields: no list)
            ls : [0] # Integer mode numbers in x or eta_1 (depending on coords)
            ms : [1] # Integer mode numbers in y or eta_2 (depending on coords)
            ns : [0] # Integer mode numbers in z or eta_3 (depending on coords)
            amps : [0.001] # amplitudes of each mode
        ModesCos :
            coords : 'logical' # in which coordinates (logical or physical)
            comps :
                e1 : [False, True, False]  # components to be initialized (for scalar fields: no list)
                b2 : [False, False, False] # components to be initialized (for scalar fields: no list)
            ls : [1] # Integer mode numbers in x or eta_1 (depending on coords)
            ms : [0] # Integer mode numbers in y or eta_2 (depending on coords)
            ns : [0] # Integer mode numbers in z or eta_3 (depending on coords)
            amps : [0.002] # amplitudes of each mode
    save_data :
        comps :
            e1 : False
            b2 : True
fluid :
    mhd :
        phys_params:
            A : 1 # mass number in units of proton mass
            Z : 1 # signed charge number in units of elementary charge
        mhd_u_space : H1vec # Hdiv | H1vec
        init :
            type : [ModesSin, ModesCos] # type of initial condition (possible types seen below)
            ModesSin :
                coords : 'logical' # in which coordinates (logical or physical)
                comps :
                    n3 : False              # components to be initialized (for scalar fields: no list)
                    uv : [True, True, True] # components to be initialized (for scalar fields: no list)
                    p3 : False              # components to be initialized (for scalar fields: no list)
                ls : [1] # Integer mode numbers in x or eta_1 (depending on coords)
                ms : [0] # Integer mode numbers in y or eta_2 (depending on coords)
                ns : [0] # Integer mode numbers in z or eta_3 (depending on coords)
                amps : [0.001] # amplitudes of each mode
            ModesCos :
                coords : 'logical' # in which coordinates (logical or physical)
                comps :
                    n3 : False               # components to be initialized (for scalar fields: no list)
                    uv : [False, False, False]  # components to be initialized (for scalar fields: no list)
                    p3 : True               # components to be initialized (for scalar fields: no list)
                ls : [0] # Integer mode numbers in x or eta_1 (depending on coords)
                ms : [1] # Integer mode numbers in y or eta_2 (depending on coords)
                ns : [0] # Integer mode numbers in z or eta_3 (depending on coords)
                amps : [0.002] # amplitudes of each mode  
        save_data :
            comps :
                n3 : False
                uv : True
                p3 : False

If these lines are not included, save all coefficients by default.

Merge request reports