Optional vertical interpolation

[jdha]

Need to have option to switch on/off vertical interpolation i.e. keep boundary conditions on the native grid or interpolate onto the destination grid

[jdha]

@anwiseNOCL looks like I opened this 4 years ago!

[jdha]

# TODO: put conditional here as we may want to keep data on paren
#              vertical grid

    DstCoord.depths = {"t": {}, "u": {}, "v": {}}

    for grd in ["t", "u", "v"]:
        DstCoord.depths[grd]["bdy_H"] = np.nanmax(z.zpoints["w" + grd], axis=0)
        DstCoord.depths[grd]["bdy_dz"] = np.diff(z.zpoints["w" + grd], axis=0)
        DstCoord.depths[grd]["bdy_dz"] = np.vstack(
            [DstCoord.depths[grd]["bdy_dz"], np.zeros((1, nbdy[grd]))]
        )
        DstCoord.depths[grd]["bdy_z"] = z.zpoints[grd]
    logger.info("Depths defined")

see profiler.py

A conditional can be placed here to pass the parent vertical gird to the child if interpolation isn't required.

NB in NEMO 4.2 deptht -> gdept likewise for u and v

e3[t/u/v] also need to be output into the BDY files