#pragma once // template<typename T> // __device__ void get_charnock_roughness(T &z0_m, T &u_dyn0, // const T h, const T U, // const T kappa, // const T h_charnock, const T c1_charnock, const T c2_charnock, // const int maxiters);