# Path: pyJHTDB/dbinfo.py
# def generate_temp_dbinfo(
#         field):


import os
import sys
import numpy as np
import ctypes
import inspect
import h5py
import io
import pyJHTDB
from wurlitzer import pipes, STDOUT
from pyJHTDB.dbinfo import interpolation_code
        dset = fh.create_dataset("xcoor", (shape[2],), maxshape=(shape[2],))
        dset[...]=xcoor
        dset = fh.create_dataset("ycoor", (shape[1],), maxshape=(shape[1],))
        dset[...]=ycoor
        dset = fh.create_dataset("zcoor", (shape[0],), maxshape=(shape[0],))
        dset[...]=zcoor

        nl = '\r\n'

        tf=open(filename+".xmf", "w")
        print(f"<?xml version=\"1.0\" ?>{nl}", file=tf)
        print(f"<!DOCTYPE Xdmf SYSTEM \"Xdmf.dtd\" []>{nl}", file=tf)
        print(f"<Xdmf Version=\"2.0\">{nl}", file=tf)
        print(f"  <Domain>{nl}", file=tf)

        return fh, tf, shape

    def hdf5_writing(
            self,
            filename,
            result,
            data_set,
            VarName,
            dim,
            time_step,
            fh,
            tf,
            shape):

        H5_ds_name='{0}_{1:04d}'.format(VarName,time_step)

# Path: pyJHTDB/dbinfo.py
# def generate_temp_dbinfo(
#         field):


# Path: pyJHTDB/libJHTDB.py
import os
import sys
import numpy as np
import ctypes
import inspect
import h5py
import io
import pyJHTDB
from wurlitzer import pipes, STDOUT
from pyJHTDB.dbinfo import interpolation_code

        dset = fh.create_dataset("xcoor", (shape[2],), maxshape=(shape[2],))
        dset[...]=xcoor
        dset = fh.create_dataset("ycoor", (shape[1],), maxshape=(shape[1],))
        dset[...]=ycoor
        dset = fh.create_dataset("zcoor", (shape[0],), maxshape=(shape[0],))
        dset[...]=zcoor

        nl = '\r\n'

        tf=open(filename+".xmf", "w")
        print(f"<?xml version=\"1.0\" ?>{nl}", file=tf)
        print(f"<!DOCTYPE Xdmf SYSTEM \"Xdmf.dtd\" []>{nl}", file=tf)
        print(f"<Xdmf Version=\"2.0\">{nl}", file=tf)
        print(f"  <Domain>{nl}", file=tf)

        return fh, tf, shape

    def hdf5_writing(
            self,
            filename,
            result,
            data_set,
            VarName,
            dim,
            time_step,
            fh,
            tf,
            shape):

        H5_ds_name='{0}_{1:04d}'.format(VarName,time_step)

        dset = fh.create_dataset(H5_ds_name,