Hello. I am trying to use the OBCS package to implement a sponge layer. Instead of relaxing temperature to tRef, which is the default parameter, I want to default to the temperature set up by HydrogThetaFile. Does anyone know what the command should be in obcs_calc.F instead of, for example, OBWt(j,k,bi,bj) =tRef(k)? Thanks Ashley Stroman