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

import copy
import pandas as pd

from .....utils.datastores.empty import init_empty

""" Per advise from Antoine, Copied from flexpart


VERSION HISTORY
2021-10-14   freum   Replaced whole function with current one from flexpart on wrf branch, commented self.reset_obs
2021-09-26   freum   Original code copied from https://git.nilu.no/VERIFY/CIF/-/blob/TM5/pycif/plugins/models/flexpart/io/outputs2native_adj.py
"""




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

       Does nothing for now as we instead read FLEXPART output 
       inside loop over observations in obsoper.py 

    """

    ddi = min(di, df)

    # Otherwise, crop the datastore to active species
    if not hasattr(self, "dataobs"):
        self.dataobs = {ddi: {spec: init_empty()
                        for spec in self.chemistry.acspecies.attributes}}

    if self.reset_obs[ddi]:
        self.dataobs[ddi] = {
            spec: init_empty()
            for spec in self.chemistry.acspecies.attributes}
        self.reset_obs[ddi] = False

    for trid in data2dump:
        mod_input = trid[0]
        trcr = trid[1]

        if mod_input != "endconcs":
            if onlyinit:
                self.dataobs[ddi][trcr] = \
                    pd.concat([self.dataobs[ddi][trcr],
                               data2dump[trid][ddi]],
                              axis=0, sort=False)