_RL smooth2DdelTime,smooth3DdelTime
PARAMETER ( smooth2DdelTime = 1. _d 0 )
smooth2DtotTime=nbt_in*smooth2DdelTime
& -gt_in(i,j,bi,bj)*recip_rA(i,j,bi,bj)*smooth2DdelTime
smooth2DtotTime=smooth2Dnbt(smoothOpNb)*smooth2DdelTime