c c c ===================================================== subroutine ic(maxmx,maxmy,maxmz,meqn,mbc,mx,my,mz, & x,y,z,dx,dy,dz,q) c ===================================================== c implicit double precision (a-h,o-z) c include "cuser.i" c dimension q(meqn, 1-mbc:maxmx+mbc, 1-mbc:maxmy+mbc, & 1-mbc:maxmz+mbc) dimension x(1-mbc:maxmx+mbc),y(1-mbc:maxmy+mbc), & z(1-mbc:maxmz+mbc) c idisc = 1 do 60 k = 1, mz do 60 j = 1, my x0 = zmax0 y0 = 0.d0 alf = 1.d0 beta = 0.d0 call cellave(z(k)-dz/2.d0,y(j)-dy/2.d0,dz,dy,wlz1) x0 = zmin0 y0 = 0. alf = -1.d0 beta = 0.d0 call cellave(z(k)-dz/2.d0,y(j)-dy/2.d0,dz,dy,wlz2) wlz=wlz1*wlz2 do 60 i = 1, mx q(1,i,j,k) = wlz*rhos + (1.d0-wlz)*rhoamb q(2,i,j,k) = 0.d0 q(3,i,j,k) = 0.d0 q(4,i,j,k) = 0.d0 q(5,i,j,k) = wlz*(ps/gamma1) + (1.d0-wlz)*(pamb/gamma1) 60 continue return end