_RL CT(irun), USTAR(irun)
CT(I) = CTSAVE(I)
SSDEV(I,NLEV)=CT(I)*(AA(I,NLEV)*DSH(I,NLEV)
9 CT,RIB,ZETA,WS,
DTHV(I,NLEV) = CT(I) * YY(I) *
CTIN(I) = CT(I)
CTSAVE(I) = CT(I)
ct_diff = gsw_ct_from_pt(sa,pt) - CT
_RL FUNCTION GSW_PT_FROM_CT(SA,CT)
_RL SA,CT
a5ct = a5*CT
b3ct = b3*CT
pt_num = a0 + s1*(a1 + a2*s1) + CT*ct_factor
pt_den = b0 + b1*s1 + CT*(b2 + b3ct)