Using "usingPCoords" rather than "buoyancyRelation" to find surface k + modify comments
Using "usingPCoords" rather than "buoyancyRelation" to find surface k
+ modify comments

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