Source code for pycif.plugins.models.wrfchem.io.outputs2native

from .outputs.read_sim import read_sim
from .outputs.endconcs import fetch_end

"""
VERSION HISTORY
2021-09-26   freum   Original code copied from chimere
- edited out what WRF doesn't use
"""




def outputs2native(
        self, data2dump, input_type, di, df, runsubdir, mode="fwd", dump=True,
        onlyinit=False, check_transforms=False, **kwargs
):
    """Reads outputs to pyCIF objects.

    If the mode is 'fwd' or 'tl', only observation-like outputs are extracted.
    For the 'adj' mode, all outputs relative to model sensitivity are extracted.

    Dumps to a NetCDF file with output concentrations if needed"""
    
    ddi = min(di, df)
    ddf = max(di, df)
    
    # If no data to extract, pass
    if data2dump == {}:
        return data2dump
    
    # Read simulations if input_type is "concs"
    # freum: only have concs, pressure and dpressure for now
    # outcomp = ["concs", "pressure", "dpressure", "airm", "hlay"]
    outcomp = ["concs", "pressure", "dpressure"]
    if input_type in outcomp and not onlyinit:
        return read_sim(self, data2dump, runsubdir, mode, ddi, ddf)

    # freum: comment, related to subsimu
    # Fetch end concentration if input_type is endconcs
    elif input_type == "endconcs":
        return fetch_end(self, data2dump, runsubdir, mode, ddi, ddf)
    
    else:
        return data2dump