& objf_vwindm, objf_vwindsmoo,
& objf_vwind, objf_vwindm, objf_vwindsmoo,
_RL objf_vwindsmoo(nsx,nsy)
& +mult_smooth_bc*objf_vwindsmoo(bi,bj) )
f_vwindsmoo = f_vwindsmoo + objf_vwindsmoo(bi,bj)