_RL rfresh, rsalt, rhoP0
rsalt = 0. _d 0
rsalt =
rhoP0 = rfresh + rsalt
_RL rfresh, rsalt
rhoP0(i,j) = rfresh + rsalt