<gh_stars>1-10
# -*- coding: utf-8 -*-
"""
Specific Frozen Wave parameters

@author: <NAME>
"""

import numpy as np

THETA = {-15: 20.2816,
         -14: 20.149,
         -13: 20.0155,
         -12: 19.8812,
         -11: 19.7459,
         -10: 19.6098,
         -9: 19.4728,
         -8: 19.3348,
         -7: 19.1959,
         -6: 19.056,
         -5: 18.9151,
         -4: 18.7732,
         -3: 18.6302,
         -2: 18.4862,
         -1: 18.3411,
         0: 18.1949,
         1: 18.0475,
         2: 17.899,
         3: 17.7492,
         4: 17.5982,
         5: 17.446,
         6: 17.2924,
         7: 17.1376,
         8: 16.9813,
         9: 16.8237,
         10: 16.6646,
         11: 16.504,
         12: 16.3419,
         13: 16.1782,
         14: 16.0128,
         15: 15.8458}

COEFF = {-15: -2.37169E-17,
          -14: 0.0133641,
          -13: 0.023287,
          -12: 0.0252276,
          -11: 0.0170089,
          -10: 1.05964E-14,
          -9: -0.0207887,
          -8: -0.0378413,
          -7: -0.0432472,
          -6: -0.031183,
          -5: -1.05783E-14,
          -4: 0.0467745,
          -3: 0.10091,
          -2: 0.151365,
          -1: 0.187098,
          0: 0.2,
          1: 0.187098,
          2: 0.151365,
          3: 0.10091,
          4: 0.0467745,
          5: -1.05783E-14,
          6: -0.031183,
          7: -0.0432472,
          8: -0.0378413,
          9: -0.0207887,
          10: 1.05964E-14,
          11: 0.0170089,
          12: 0.0252276,
          13: 0.023287,
          14: 0.0133641,
          15: -2.37169E-17}

def axicon_omega(degree, axicon):
    return np.sin(np.deg2rad(axicon)) * (degree + 1/2)