44 real(dp),
intent(in) :: temp_val
50 #if defined(FRAC_DUST)
51 c_ice = c(floor(temp_val)) &
52 +(c(floor(temp_val)+1)-c(floor(temp_val))) &
53 *(temp_val-
real(floor(temp_val),dp))
55 c_val = c(floor(temp_val)) &
56 +(c(floor(temp_val)+1)-c(floor(temp_val))) &
57 *(temp_val-
real(floor(temp_val),dp))
63 #if defined(FRAC_DUST)
64 c_val = rho_inv * ( (1.0_dp-frac_dust)*rho_i*c_ice + frac_dust*rho_c*c_c )