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