      subroutine dummy(a,b)
      use nrtype
      real(SP) :: a
      real(DP) :: b
     
      write(6,*) 'a=', a, 'b=', b
     
      return 

      end subroutine dummy