/* * $Id: c_tdex01.c,v 1.4 2004/07/23 01:00:33 kennison Exp $ */ /* * Include standard I/O and math libraries. */ #include #include /* * Include NCAR Graphics function prototypes. */ #include /* * Define the workstation type and identifier. */ #define IWTYPE 1 #define WKID 1 /* * Define the maximum dimensions of the data arrays. */ #define IMAX 41 #define JMAX 41 #define KMAX 41 /* * Define the size of the triangle list. */ #define MTRI 200000 /* * Define a constant to use in converting angles from degrees to radians. */ #define DTOR .017453292519943 /* * Declare local dimensioned variables to hold data defining a simple * surface and an isosurface. */ float u[IMAX],v[JMAX],w[KMAX],s[JMAX][IMAX],f[KMAX][JMAX][IMAX]; /* * Declare the triangle-list array and a couple of temporary arrays to * be used in sorting the list. */ float rtri[MTRI][10],rtwk[2][MTRI]; int itwk[MTRI]; /* * Start of program. */ main() { /* * Set the desired values of the dimensions of the data arrays. Note that * "idim" must not exceed IMAX, that "jdim" must not exceed JMAX, and that * "kdim" must not exceed KMAX. */ int idim=31; int jdim=31; int kdim=31; /* * Declare the desired minimum and maximum values of U, V, and W. */ float umin=-1.; float vmin=-1.; float wmin=-1.; float umax= 1.; float vmax= 1.; float wmax= 1.; /* * Set the desired values of parameters determining the eye position. * "ang1" is a bearing angle, "ang2" is an elevation angle, and "rmul" * is a multiplier of the length of the diagonal of the data box, * specifying the distance from the center of the box to the eye. */ float ang1=-31.; float ang2= 35.; float rmul= 2.9; /* * Set the desired value of the flag that says whether the basic color * scheme will be white on black ("ibow"=0) or black on white ("ibow"=1). */ int ibow=1; /* * Define the range of color indices to be used for a gray scale ranging * all the way from pure white to pure black. (Following blocks of color * indices of the same length will be used for shades of other colors.) */ int iofc= 8; int iolc=27; /* * Set the desired values of the shading parameters. Values of "shde" near * 0 give brighter colors and values near 1 give pastel shades. Values of * "SHDR" near 0 give a narrow range of shades and values near 1 give a wide * range of shades. */ float shde=.1; float shdr=.8; /* * The flag "iste" says whether to do a simple image ("iste"=0), a * one-frame stereo image ("iste"=-1), or a two-frame stereo image * ("iste"=+1). */ int iste=-1; /* * The value of "aste" is the desired angle (in degrees) between the lines * of sight for a pair of stereo views. */ float aste=4.; /* * The value of "wosw" is the width of the stereo windows to be used in * one-frame stereo images; the width is stated as a fraction of the width * of the plotter frame. (The windows are centered vertically; horizontally, * they are placed as far apart as possible in the plotter frame.) The value * used must be positive and non-zero; it may be slightly greater than .5, if * it is desired that the stereo windows should overlap slightly. */ float wosw=.5; /* * Define labels for the edges of the box. */ char* unlb=" -1 -.8 -.6 -.4 -.2 0 .2 .4 .6 .8 1 "; char* vnlb=" -1 -.8 -.6 -.4 -.2 0 .2 .4 .6 .8 1 "; char* wnlb=" -1 -.8 -.6 -.4 -.2 0 .2 .4 .6 .8 1 "; char* uilb="Values of U Coordinate"; char* vilb="Values of V Coordinate"; char* wilb="Values of W Coordinate"; /* * Declare other variables to be used below. */ int i,irst,j,k,nshd,ntri; float otep,r,ueye,veye,weye,umid,vmid,wmid,uthi,vthi,wthi; float xvpl,xvpr,yvpb,yvpt,xwdl,xwdr,ywdb,ywdt; int lnlg; float xzro,yzro,xone,yone,xzbk,yzbk,xobk,yobk; float csm1,csm2,save_csm2; /* * Fill the data arrays. */ for (i=0;i0.) c_frame(); } /* * Order the triangles in the triangle list. */ c_tdotri(&rtri[0][0],MTRI,&ntri,&rtwk[0][0],&itwk[0],1); /* * Draw the sides of the box that could be hidden. */ c_tdgrds(umin,vmin,wmin,umax,vmax,wmax, .1*(umax-umin),.1*(vmax-vmin),.1*(wmax-wmin),03,1); /* * Add a couple of "equatorial plane" lines on the far sides of the box. */ gset_line_colr_ind((Gint)(iofc+7*nshd+1*nshd/8)); gset_linewidth((Gdouble)2.); c_tdline( 1., 1.,0.,-1., 1.,0.); c_tdline(-1., 1.,0.,-1.,-1.,0.); gset_linewidth((Gdouble)1.); gset_line_colr_ind((Gint)1); /* * Draw the triangles in the triangle list. */ c_tddtri(&rtri[0][0],MTRI,&ntri,&itwk[0]); /* * Draw the sides of the box that cannot be hidden. */ c_tdgrds(umin,vmin,wmin,umax,vmax,wmax, .1*(umax-umin),.1*(vmax-vmin),.1*(wmax-wmin),03,0); /* * Add a couple of "equatorial plane" lines on the near sides of the box. */ gset_line_colr_ind((Gint)(iofc+7*nshd+1*nshd/8)); gset_linewidth((Gdouble)2.); c_tdline(-1.,-1.,0., 1.,-1.,0.); c_tdline( 1.,-1.,0., 1., 1.,0.); gset_linewidth((Gdouble)1.); gset_line_colr_ind((Gint)1); /* * Draw labels for the axes. */ c_tdlbls(umin,vmin,wmin,umax,vmax,wmax,unlb,vnlb,wnlb, uilb,vilb,wilb,1); /* * Put a label above the data box, with a double-width arrow from the * label to the feature that it labels. */ gset_line_colr_ind((Gint)(iofc+6*nshd+1*nshd/8)); c_tdpara(0.,.25,1.25,0.,1.,0.,0.,0.,1.); c_pcseti("CC - CHARACTER COLOR",iofc+6*nshd+1*nshd/8); c_tdplch(0.,0.,"Maximum Z Value",.08,0.,-1.05); c_pcseti("CC - CHARACTER COLOR",-1); gset_linewidth((Gdouble)2.); c_tdline(0.,0.,1.,0.,.25,1.25); c_tdline(0.,0.,1.,0.,.09,1.05); c_tdline(0.,0.,1.,0.,.05,1.09); gset_linewidth((Gdouble)1.); gset_line_colr_ind((Gint)1); /* * Add a label for the "equatorial plane". */ gset_line_colr_ind((Gint)(iofc+7*nshd+1*nshd/8)); c_tdpara(1.1,-1.1,0.,0.,1.,0.,-1.,0.,0.); c_pcseti("CC - CHARACTER COLOR",iofc+7*nshd+1*nshd/8); c_tdplch(0.,0.,"\"Equatorial Plane\"",.08,0.,-1.15); c_pcseti("CC - CHARACTER COLOR",-1); gset_line_colr_ind((Gint)1); /* * If a left-eye view has just been done, loop back for a right-eye view. */ } while ((otep=-otep)>0.); /* * Advance the frame. */ c_frame (); /* * Close gks. */ gdeactivate_ws(WKID); gclose_ws(WKID); gclose_gks(); }