Ug = 8.0; Vg = 0.0; hbl_0 = 100.0; Tgrad = 0.01; cooling_rate = -0.25 * 3600.0; time{ begin = 0.0; end = 20.0; #9.0 * 3600.0; dt = 20.0; out_dt = 360.0; } fluid { pref= 1013.4; #hPa tref = 265.0; #K g = 9.81; # m/s2 beta = g /tref; kappa = 0.4; # von Karman constant } grid { type = "uniform"; nz = 8; h_bot = 0.0; h_top = 400.0; } surface_model { id = "esm"; # optional: "esm" (default), "sheba", "most", "log" type = "land"; z0_m = 0.01; # aerodynamic roughness [m] z0_h = -1; # no prescribed value # -> using scheme assigned by surface type } output{ dt = 0.2 *3600.0; } closure { type = "FOM"; name = "Louis"; # "Esau"; # "EFB"; }