altitude.py 1.92 KiB

"""
ZS was extracted from a netcdf file
"""
import numpy as np

ZS = """[ 300.  600.  900. 1200. 1500. 1800. 2100. 2400. 2700. 3000. 3300.  300.
600.  900. 1200. 1500. 1800. 2100. 2400. 2700. 3000.  300.  600.  900.
1200. 1500. 1800. 2100. 2400. 2700. 3000. 3300. 3600. 3900. 4200. 4500.
4800.    0.  300.  600.  900. 1200. 1500. 1800. 2100. 2400.  300.  600.
900. 1200. 1500. 1800. 2100. 2400. 2700. 3000.  600.  900. 1200. 1500.
1800. 2100. 2400. 2700. 3000. 3300. 3600. 3900.    0.  300.  600.  900.
1200. 1500. 1800. 2100.    0.  300.  600.  900. 1200. 1500. 1800. 2100.
2400. 2700. 3000.    0.  300.  600.  900. 1200. 1500. 1800. 2100. 2400.
2700. 3000. 3300.  300.  600.  900. 1200. 1500. 1800. 2100. 2400. 2700.
3000. 3300. 3600. 3900.  900. 1200. 1500. 1800. 2100. 2400. 2700. 3000.
3300. 3600. 3900.  600.  900. 1200. 1500. 1800. 2100. 2400. 2700. 3000.
3300. 3600.  900. 1200. 1500. 1800. 2100. 2400. 2700. 3000. 3300.    0.
300.  600.  900. 1200. 1500. 1800. 2100. 2400.    0.  300.  600.  900.
1200. 1500. 1800. 2100. 2400. 2700. 3000. 3300. 3600. 3900. 4200.  900.
1200. 1500. 1800. 2100. 2400. 2700. 3000. 3300. 3600. 3900. 4200.  900.
1200. 1500. 1800. 2100. 2400. 2700. 3000. 3300.  300.  600.  900. 1200.
1500. 1800. 2100. 2400. 2700. 3000.  600.  900. 1200. 1500. 1800. 2100.
2400. 2700. 3000. 3300. 3600.  600.  900. 1200. 1500. 1800. 2100. 2400.
2700. 3000. 3300.  600.  900. 1200. 1500. 1800. 2100. 2400. 2700. 3000.
3300. 3600.    0.  300.  600.  900. 1200. 1500. 1800. 2100. 2400. 2700.
3000. 3300.  300.  600.  900. 1200. 1500. 1800. 2100. 2400. 2700. 3000.
3300.  300.  600.  900. 1200. 1500. 1800. 2100. 2400. 2700. 3000.]"""

ZS_INT = [int(float(e)) for e in ZS[1:-1].split()]
ALTITUDES = sorted(set(ZS_INT))


# Create a ZS_INT with only the 23 first massifs
ZS_INT_23 = ZS_INT[:-10].copy()

# Create a ZS_INT with np.nan all altitudes corresponding to the 24th massif
ZS_INT_MASK = np.array(ZS_INT)
ZS_INT_MASK[-10:] = np.nan