0.5e0    kappa    # thermal conductivity of fluid [W/mK]
1e0      kappa_1  # thermal conductivity of solid 1 [W/mK]
1e0      kappa_2  # thermal conductivity of solid 2 [W/mK]
1e0      kappa_m1 # thermal conductivity of thin film 1 [W/mK]
1e0      kappa_m2 # thermal conductivity of thin film 2 [W/mK]
1e3      rho      # density of a fluid [kg/m^3]
4.2e3    C_p      # ** not used ** 
2.4e3    Aa       # absorption coefficient of a fluid [1/m]
0e0      Ab1      # absorption coefficient of solid 1 [1/m]
0e0      Ab2      # absorption coefficient of solid 2 [1/m]
50e-3    P        # laser power
10e-6    w_0      # beam waist at the focal plane [m]
0e-6     z_0      # z-position of the focal plane [m]
1480e-9  lambea   # wavelength of the laser [m] (affects z_R)
9.8e0    g        # gravity [m/s^2]
298e0    T_0      # reference temperature [K]
2.4e-4   betac    # thermal expansion coefficient [1/K]
40e-6    H        # height of a fluid domain [m]
80e-6    H_1      # thickness of solid 1 [m]
80e-6    H_2      # thickness of solid 2 [m]
0e-6     H_m1     # thickness of thin film 1 [m]
0e-6     H_m2     # thickness of thin film 2 [m]
1e-6     nu       # kinematic viscosity [m^2/s]
20       Nr       # the number of grid points in r direction
1e0      Aspect   # Aspect ratio of the domain
20       Nz       # the number of grid points in z direction in the fluid
10       Nz_1     # the number of grid points in z direction in the solid 1
10       Nz_2     # the number of grid points in z direction in the solid 2
28       Nm       # the number of grid points in the numerical integration
1        case     # case 1 -> use simpson's rule; case 2 -> use trapezoidal rule
2e0      s_max0   # the parameter used to determine the end point of the numerical integration
100      save1    # ** not used ** 
1        threads  # ** not used **  
1e-12    D_T      # thermophoretic mobility [m^2/sK]  (not used in the python code)
0.5e-6   raeius   # the radius of dispersed particles (not used in the python code)
-2e-12   Kd1      # thermal slip coefficient at the surface 1 [m^2/sK]
-2e-12   Kd2      # thermal slip coefficient at the surface 2 [m^2/sK]
ref      memo     # ** not used **