multiply sigmaR by maskC(k)*maskC(k-1) (instead of just maskC(k))
multiply sigmaR by maskC(k)*maskC(k-1) (instead of just maskC(k))

altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h