From 754967ed61730520be0b449d8e8323155ccb1f55 Mon Sep 17 00:00:00 2001
From: root <>
Date: Sun, 8 Feb 2009 18:35:53 +0000
Subject: *** empty log message ***

---
 test/ex2.f | 123 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 123 insertions(+)
 create mode 100644 test/ex2.f

(limited to 'test/ex2.f')

diff --git a/test/ex2.f b/test/ex2.f
new file mode 100644
index 0000000..fa6584d
--- /dev/null
+++ b/test/ex2.f
@@ -0,0 +1,123 @@
+ 	
+	program moocow
+	dimension a(20000)
+	
+	integer*4 jwg_create_handle
+
+ 	parameter (ix=10,iy=10,iw=150,icolor=1,width=0.0)
+
+	jx=ix/2
+	jy=iy/2
+	sx=5.0/(1.0*jx)
+	sy=5.0/(1.0*jy)
+
+	do j=1,iy
+	 y=(j-jy)*sy
+
+	 k=iw*(j-1)
+	 
+	do i=1,ix
+	 x=(i-jx)*sx
+ 	a(k+i)=x*y
+C		4.0*r(x,y)
+
+	end do
+	end do 
+
+
+	ihandle=jwg_create_handle(1,'ex1.cgm')
+
+
+	call jwg_def_xform(ihandle)
+
+	sx=300.0/(1.0*(ix-1))
+	sy=300.0/(1.0*(iy-1))
+
+	call jwg_push_state(ihandle)
+
+	call jwg_add_xform(ihandle,100.0,100.0,sx,0.0,0.0,sy)
+	call jwg_add_xform(ihandle,-1.0,-1.0,1.0,0.0,0.0,1.0)
+
+ 	call jwg_pass_data(ihandle,a,ix,iy,iw)
+
+C	do i=1,5
+	
+	i=0
+	
+         clevel= -45+(i*15)
+
+	clevel=15.0
+
+	 red=(i*1.0)/5.0
+	 green=red
+	 blue=red
+
+	 call jwg_set_pen(ihandle,red,green,blue,0.0)
+	 call jwg_contour(ihandle,clevel,-1,1)
+C	end do
+
+	call jwg_destroy_data(ihandle)
+	
+
+	sx=5.0/(1.0*jx)
+	sy=5.0/(1.0*jy)
+
+	do j=1,iy
+	 y=(j-jy)*sy
+
+	 k=iw*(j-1)
+
+	do i=1,ix
+
+	 x=(i-jx)*sx
+
+ 	 a(k+i)=x*y
+
+	end do
+	end do 
+
+
+
+C	call jwg_add_xform(ihandle,0.0,iy*1.0,1.0,0.0,0.0,1.0)
+C	call jwg_add_xform(ihandle,0.0,iy*0.1,1.0,0.0,0.0,1.0)
+
+C	call jwg_pass_data(ihandle,a,ix,iy,iw)
+
+
+	do i=1,5
+	
+         clevel= -45+(i*15)
+	 red=(i*1.0)/5.0
+	 green=red
+	 blue=red
+
+C	 call jwg_set_pen(ihandle,red,green,blue,0.0)
+C	 call jwg_contour(ihandle,clevel,-1,1)
+	end do
+
+C	call jwg_destroy_data(ihandle)
+
+
+
+	call jwg_pop_state(ihandle)
+
+
+	call jwg_destroy_handle(ihandle)
+
+	end
+
+	function r(x,y)
+	r=sqrt((x*x)+(y*y))
+	return
+ 	end
+
+	function sinc(a)
+
+	if(a.ne.0.)then
+		sinc=sin(a)/a
+	else
+		sinc=1.0
+	endif
+ 	return	
+	end
+	
-- 
cgit v1.2.3