anchor.diff:          write(*,'(a,l1)') 'lhf     ', bnd%lhf
           write(*,'(a,l1)') 'shf     ', bnd%shf
           write(*,'(a,l1)') 'subl    ', bnd%subl
  -        write(*,'(a,l1)') 'massbal ', bnd%massbal
  +        write(*,'(a,l1)') 'smb     ', bnd%smb
       end subroutine
   
   
  
positive.code:                     bnd%melt    = .TRUE. 
                  case("refr")
                      bnd%refr    = .TRUE. 
                  case("massbal")
                      bnd%massbal = .TRUE. 
                  case("acc")
                      bnd%acc     = .TRUE. 
                  case("lhf")
  
negative.code:         allocate(now%wind(npts))
          allocate(now%rhoa(npts))
          allocate(now%qq(npts))
          allocate(now%mask(npts))
  
          return 
      end subroutine surface_alloc 
  