subroutine gr_gazrare(nelem,jz,rarflag)
! identifie une configuration de gaz rare (G.Alecian, 26/06/2001)
	implicit none
	integer, dimension (6)      :: zgrar
	integer                     :: nelem,j,jz,rarflag
	data zgrar/2,10,18,36,54,86/
	do j=1,6
		if(zgrar(j).eq.(nelem-jz)) then
			rarflag=1
			exit
		end if
	end do
	end subroutine gr_gazrare