double precision g(4,4) data g/ . -1d0,0d0 ,0d0 , 0d0, . 0d0 ,-1d0,0d0 , 0d0, . 0d0 , 0d0,-1d0, 0d0, . 0d0, 0d0, 0d0,+1d0/ save g