The description below in psource_ncfile part of Section 16.2 of the documentation is:
“qbardir: Orientation: South=0 or North=0, East=0 or West=1”
But the codes in the relevant part (get_vbc.F) in line 610 is :
“Qbar(is)=qbardir(is)*Qbar(is)”
Should the index be -1 instead of 0 in the description?