1
 #define SAF_3DF(nx,ny,nz)       _saf_indexspec(3,nx,ny,nz,1,1,1,SAF_F_ORDER)