bug fixed in roughness_re

srcF/module_z0t_lake.f90

@@ -138,10 +138,12 @@ module module_z0t_lake
        
       
        
-        ! --- define roughness [thermal]
-        z0_t = z0_m (-0.56*(4.0*sqtr(Re)-3.4))
+        
+        B=alog(-0.56*(4.0*(Re)**(0.5)-3.4))     !Repina, 2023
+
 
-        B=log(z0_m/z0_t)      !Repina, 2023
+       ! --- define roughness [thermal]
+        z0_t = z0_m / exp(B)   
        
     end subroutine