Adjust sediment temperature model

models/modules/nivafjord/basin.txt

@@ -240,7 +240,9 @@ Authors: Magnus D. Norling
 		(in_flux(sw_vert, layer.water.heat)->>)*(1 - attn) * A[vert.below]/A
 	}
 
+	sed_alb : constant("Sediment albedo", [], 0.3)
+
 	flux(layer.water.heat, sw_sed, [J, day-1], "Shortwave to sediments") {
-		(in_flux(sw_vert, layer.water.heat)->>)*(1 - attn) * (A - A[vert.below])/A
+		(in_flux(sw_vert, layer.water.heat)->>)*(1 - attn)*(1 - sed_alb) * (A - A[vert.below])/A
 	}
 }

models/modules/nivafjord/sediment.txt

@@ -41,8 +41,8 @@ module("NIVAFjord sediments", version(0, 0, 3),
 
 	# Set these as parameters: ?
 	sed_rho : constant("Sediment density", [k g, m-3], 2600)
-	sed_C   : constant("Sediment specific heat capacity", [J, k g-1, K-1], 5000) #TODO: Find a good value
-	sed_k   : constant("Sediment heat conductivity coefficient", [W, m-1, K-1], 20) #TODO: This is 10 times too high, but it probably must black-body radiate??
+	sed_C   : constant("Sediment specific heat capacity", [J, k g-1, K-1], 4000) #TODO: Find a good value
+	sed_k   : constant("Sediment heat conductivity coefficient", [W, m-1, K-1], 2) #TODO: This is 10 times too high, but otherwise heat just builds up...
 
 	par_group("Sediment temperature") {
 		dzh : par_real("Thickness of thermally active sediment layer", [m], 0.5)
@@ -59,10 +59,11 @@ module("NIVAFjord sediments", version(0, 0, 3),
 		V := area*dzh,
 		T_avg := heat / (sed_C*sed_rho*area*dzh) -> [deg_c],
 		2*T_avg - T_bot
+		#(T_avg - T_bot/4)*(4/3)
 	} @initial { layer.water.temp } # Or T_bot?
 
 	flux(layer.water.heat, sediment.heat, [J, day-1], "Water-sediment heat conduction") {
-		dz_ := 0.5*(layer.dz + dzh),
+		dz_ := 0.5[m],#layer.dz,#0.5*(layer.dz + dzh),
 		sediment.area*(layer.water.temp->[K] - sediment.temp->[K])*sed_k/dz_ ->>
 	}