divide Hs=0 or not and change qs*dt to qs_ave*dt