add subroutine drawband, independent file drawband.f90. However, none of above versions works fine