module constants character(50),parameter::readname="idd1n540.data" character(50),parameter::writename="out.data" real,parameter::delta=0.3 !mesh size[Re] real,parameter::eta=0.002 !resistivity integer,parameter::mx=268,my=134,mz=134 integer,parameter::mx2=mx+2,my2=my+2,mz2=mz+2 real::d=delta*2 real,dimension(mx2,my2,mz2,21)::f end module