getPoci model coefficient for Northern Hemisphere

[Wenbo-Duan]

When generating tracks, TCRM calculates modified pressure for the outermost closed isobar based on a model of daily long-term mean SLP values, central pressure, latitude and day of year. The default coefficients of the model are based on Southern Hemisphere data. Would it be possible to provide some references on the model and the ways to fit the model base on Northern Hemisphere data?

def getPoci(penv, pcentre, lat, jdays, eps, coeffs=[2324.1564738613392, -0.6539853183796136, -1.3984456535888878, 0.00074072928008818927, 0.0044469231429346088, -1.4337623534206905], missingValue=sys.maxsize):

    """
    Calculate a modified pressure for the outermost closed isobar, based
    on a model of daily long-term mean SLP values, central pressure,
    latitude and day of year.

    :param penv: environmental pressure estimate (from long term mean pressure
                 dataset, hPa).
    :param pcentre: Central pressure of storm (hPa).
    :param lat: Latitude of storm (degrees).
    :param jdays: Julian day (day of year).
    :param eps: random variate. Retained as a constant for a single storm.
    :param list coeffs: Coefficients of the model. Defaults based on 
                        **Southern Hemisphere data** 
                        (IBTrACS v03r06, 1981-2014). 

    :returns: Revised estimate for the pressure of outermost closed isobar.
    """

[wcarthur]

The coefficients in the Poci model were derived using this notebook: https://github.com/GeoscienceAustralia/tcrm/blob/c7447089be8dfe056a81a0231f5322c6af4f966b/Validating%20environmental%20pressure%20values.ipynb

You could re-run the notebook and use Northern Hemisphere data - note we used JTWC data, which uses b-deck format for observed tracks.