_RL bMfresh, bMsalt, bMpres, BulkMod
bMsalt = 0. _d 0
bMsalt =
bulkMod = bMfresh + bMsalt + bMpres
_RL bMfresh, bMsalt, bMpres
bulkMod(i,j) = bMfresh + bMsalt + bMpres