BEGIN_PROVIDER[real*8, val,(natomax)]
    implicit none
    val=0d0
END_PROVIDER