integer smooth3Dnbt(smoothOpNbMax),
& smooth3Dnbt, smooth2Dnbt,
nbt_in=smooth3Dnbt(smoothOpNb)
nbt_in=smooth3Dnbt(smoothOpNb)/2
smooth3DtotTime=smooth3Dnbt(smoothOpNb)*smooth3DdelTime
& smoothOpNb,smooth3Dnbt(smoothOpNb),
& smooth3Dnbt,
smooth3Dnbt(smoothOpNb)=0