mpas_tools.mesh.creation.mesh_definition_tools.EC_CellWidthVsLat

mpas_tools.mesh.creation.mesh_definition_tools.EC_CellWidthVsLat(lat, cellWidthEq=30.0, cellWidthMidLat=60.0, cellWidthPole=35.0, latPosEq=15.0, latPosPole=73.0, latTransition=40.0, latWidthEq=6.0, latWidthPole=9.0)[source]

Create Eddy Closure spacing as a function of lat. This is intended as part of the workflow to make an MPAS global mesh.

Parameters:

  • lat (ndarray) – vector of length n, with entries between -90 and 90, degrees

  • cellWidthEq (float, optional) – Cell width in km at the equator

  • cellWidthMidLat (float, optional) – Cell width in km at mid latitudes

  • cellWidthPole (float, optional) – Cell width in km at the poles

  • latPosEq (float, optional) – Latitude in degrees of center of the equatorial transition region

  • latPosPole (float, optional) – Latitude in degrees of center of the polar transition region

  • latTransition (float, optional) – Latitude in degrees of the change from equatorial to polar function

  • latWidthEq (float, optional) – Width in degrees latitude of the equatorial transition region

  • latWidthPole (float, optional) – Width in degrees latitude of the polar transition region

Returns:

cellWidthOut (ndarray) – 1D array of same length as lat with entries that are cell width as a function of lat

Examples

Default

>>> EC60to30 = EC_CellWidthVsLat(lat)

Half the default resolution:

>>> EC120to60 = EC_CellWidthVsLat(lat, cellWidthEq=60., cellWidthMidLat=120., cellWidthPole=70.)