question about msftbarot in OM2 CM2

[paolap]

According to the APP4 mappings both psiu and tx_trans are applied to calculate msftbarot (ocean_barotropic_mass_streamfunction)

using this function:

def msftbarot(psiu,tx_trans):
drake_trans=transAcrossLine(tx_trans,212,212,32,49)
#loop over times
for i,trans in enumerate(drake_trans):
#offset psiu by the drake passage transport at that time
psiu[i,:]=psiu[i,:]+trans
return psiu

In theory psiu and msfbarot are equivalent and just using psiu should be sufficient, provided psiu is calculated correctly by the model. Does anyone know if the issue with the way the model calculate psiu is still current for the OM2 and CM2 more recent model runs?

@taimoorsohail @aekiss @rmholmes tagging you just in case you might have an answer, thanks!