code: def cas2tas(cas, h):
      
      p, rho, T = atmos(h)
      qdyn = p0*((1.+rho0*cas*cas/(7.*p0))**3.5-1.)
      tas = np.sqrt(7.*p/rho*((1.+qdyn/p)**(2./7.)-1.))
      tas = -1 * tas if cas < 0 else tas
      return tas
docs: cas2tas conversion both m/s h in m