BUGFIX for SO=1 (it worked as SO=2). ctrlgnM1.py (Ge 3d is added)
    
 BUG was in hambls.F. I forgot to comment out a line h(1:ndimh,1:ndimh,3) = 0 in hambls.F