ERF_WeatherDataInterpolation.cpp File Reference

#include <filesystem>
#include <stdexcept>
#include "ERF.H"
#include "ERF_ReadCustomBinaryIC.H"
#include "ERF_Interpolation_Bilinear.H"

Include dependency graph for ERF_WeatherDataInterpolation.cpp:

Macros
#define	ERF_WEATHERDATAINTERPOLATION_H_

Enumerations
enum class	MultiFabType { CC , NC }

Functions
void	PlotMultiFab (const MultiFab &mf, const Geometry &geom_mf, const std::string plotfilename, MultiFabType mftype)

Macro Definition Documentation

ERF_WEATHERDATAINTERPOLATION_H_

#define ERF_WEATHERDATAINTERPOLATION_H_

Enumeration Type Documentation

MultiFabType

enum MultiFabType

strong

Enumerator
CC
NC

{ CC, NC };

Function Documentation

PlotMultiFab()

void PlotMultiFab	(	const MultiFab &	mf,
		const Geometry &	geom_mf,
		const std::string	plotfilename,
		MultiFabType	mftype
	)

{
 
    Vector<std::string> varnames = {
    "rho", "uvel", "vvel", "wvel", "theta", "qv", "qc", "qr", "latitude", "longitude"
    }; // Customize variable names
 
    const Real time = 0.0;
 
 
    // Assume weather_mf is nodal in all directions
    if(mftype == MultiFabType::NC) {
        BoxArray cba = mf.boxArray();
        cba = amrex::convert(mf.boxArray(), IntVect::TheCellVector());
 
        MultiFab cc_mf(cba, mf.DistributionMap(),
               mf.nComp(), 0);
 
        amrex::average_node_to_cellcenter(cc_mf, 0, mf, 0, mf.nComp());
 
        WriteSingleLevelPlotfile(
            plotfilename,
            cc_mf,
            varnames,
            geom_mf,
            time,
            0 // level
        );
    } else {
        WriteSingleLevelPlotfile(
            plotfilename,
            mf,
            varnames,
            geom_mf,
            time,
            0 // level
        );
    }
}