c      main program
c      parameter(nx=101,ny=101,nx2=nx*2,nxy=nx*ny,ico=7)
       parameter(nx=41,ny=41,nx2=nx*2,nxy=nx*ny,ico=7)
       parameter(nx0=nx-1,ny0=ny-1,nxy0=nx0*ny0)
       parameter(nxa=nx0-1,nya=ny0-1)
       parameter(x0=0.0,y0=0.0,z0=0.0)
c      parameter(xll=8.0,yll=8.0,zll=8.0,ipx0=3)
       parameter(xll=8.0,yll=8.0,zll=8.0,ipx0=4)
       dimension v(nxy0),u(nxy),p(nx2)
       character ch1(10)*40
c      icc=1:x-y plane, 2:y-z plane, 3:x-z plane
c
       data ch1(1)/"Nagoya University                       "/
       data ch1(2)/"Solar-Terrestrial Environment Laboratory"/
       pi=3.1415926
       amp=1.0
c
       do 10 j=1,ny
       do 10 i=1,nx
       x=2.0*pi*float(i-1)/float(nx-1)
       y=2.0*pi*float(j-1)/float(ny-1)
       i1=i+nx*(j-1)
       u(i1)=amp*sin(x)*sin(y)
       u(i1)=0.5*(u(i1)+amp)
c      u(i1)=0.5*(u(i1)+1.0)
   10  continue
c
       do 12 j=1,ny0
       do 12 i=1,nx0
       x=2.0*pi*(float(i-1)+0.5)/float(nx-1)
       y=2.0*pi*(float(j-1)+0.5)/float(ny-1)
       i1=i+nx0*(j-1)
       v(i1)=amp*sin(x)*sin(y)
       v(i1)=0.5*(v(i1)+amp)
   12  continue
c
       call trans1(nx,ny,v,u)
c 
       xb=x0
       yb=y0
       xl=xll
       yl=yll
       icc=3
       zcc=0.0
c      
c
       call initvrml
c
       x=x0
       y=y0+yll*1.3
       z=z0
       h=0.06*xll
       n=40
       call symblv(x,y,z,h,ch1(1),n)
       y=y0+yll*1.2
       z=z0
       h=0.06*xll
       n=40
       call symblv(x,y,z,h,ch1(2),n)
c
c
       dx=0.00*xll
       dy=0.00*yll
       dz=0.00*zll
c
       r=0.0
       g=1.0
       b=1.0
       icc=1
       zcc=z0
       x1=x0-dx
       y1=y0-dy
       x2=x0+dx+xll
       y2=y0+dy+yll
       call rect1(icc,zcc,r,g,b,x1,y1,x2,y2)
       zcc=z0+zll
       call rect1(icc,zcc,r,g,b,x1,y1,x2,y2)
c
cc
       r=0.0
       g=1.0
       b=1.0
       icc=2
       zcc=x0+xll
       x1=y0-dy
       y1=z0-dz
       x2=y0+dy+yll
       y2=z0+dz+zll
       call rect1(icc,zcc,r,g,b,x1,y1,x2,y2)
       zcc=x0
       call rect1(icc,zcc,r,g,b,x1,y1,x2,y2)
c
cc
       r=0.0
       g=1.0
       b=1.0
       icc=3
       zcc=y0
       x1=x0-dx
       y1=z0-dz
       x2=x0+dx+xll
       y2=z0+dz+zll
       call rect1(icc,zcc,r,g,b,x1,y1,x2,y2)
       zcc=y0+yll
       call rect1(icc,zcc,r,g,b,x1,y1,x2,y2)
c
c
       xb=x0
       yb=y0
       zb=z0
       xl=xll
       yl=yll
       zl=zll
       icc=1
       zcc=0.0
       zcc=yb
c      call image1(nx,ny,xb,yb,xl,yl,ico,icc,zcc,u,p)
       call pixel2(nx,ny,xb,yb,xl,yl,ipx0,ico,icc,zcc,u)
c
       icc=2
       zcc=0.0
       zcc=xl
c      call image1(nx,ny,xb,yb,xl,yl,ico,icc,zcc,u,p)
       call pixel2(nx,ny,xb,yb,xl,yl,ipx0,ico,icc,zcc,u)
c
       icc=3
       zcc=0.0
       zcc=zb
c      call image1(nx,ny,xb,yb,xl,yl,ico,icc,zcc,u,p)
       call pixel2(nx,ny,xb,yb,xl,yl,ipx0,ico,icc,zcc,u)
c
       stop
       end