/* FAN (C) Copyright 2012
Källkoderna tillhör IDC (Industrial Development Center i USA) och Evelina Rockefeller
Dom är författad av Bengt-Olof Drugge och är sensat ändrad 2012-07-09 kl 18.00 
Det här programmet är AutCADs del i ett expertsystem för att bygga optimala fläktar
Jag kallar mitt expert system för FANnik och rätten till det expert systemet ägs av IDC
och jag som jobbar för IDC ger källkoderna till IDC och Evelina Rockefeller som chefs ingenjör på
IDC eller (Industrial Development Center i USA).

Bengt-Olof Drugge

*/


#include  <stdio.h>
#include  <math.h>
#include  <string.h>
#include <dbsymtb.h>
#include <dbidmap.h>
#include <dbents.h>
#include <adslib.h>
#include <adsdlg.h>
#include <ads.h>
#include <acmem.h>
#include <rxregsvc.h>

#define   fcount 6
#define   err 1e-5
#define   err2 1e-10

#define ELEMENTS(array) (sizeof(array)/sizeof((array)[0]))

static double pi = 3.141592653589793;
static int loadfuncs();
ads_point ptup1[15001],ptdw1[15001],ptup2[15001],ptdw2[15001];
int spiral();
int hus();
int inlp();
int hus1();
int prop1();
int funcload();
int funcunload();
int fan1();
int fnik(); 
int hol(ads_name ent1,ads_name ent2);
void cross();
void subb();
void add();
void trans();
void midp();
void trim(int n1, int n2);
    

/**************************************************************************/
/* ADS Function Table */
typedef struct {
    char    *name;
    int    (*fptr)();
} ftblent;


/* Table of ADS functions */
ftblent ftable[] = {
	{ "C:3DSPIRAL",	spiral},
	{ "C:FANHOUSE", hus1},
	{ "C:FAN",   	fan1},
    { "C:PROP",    prop1}, 
    { "C:FANNIK",    fnik},
};

int	funcload()
{
	int i;
    
    
	for (i = 0; i < ELEMENTS(ftable); i++) {
        if (!ads_defun(ftable[i].name, i))
            return RTERROR;
    /*ds_regfunc(ftable[i].fptr,i);*/

	}

    return RTNORM;
}

int	funcunload()
{
    int i;
    
	/* Undefine each function we defined */

	for (i = 0; i < ELEMENTS(ftable); i++) {
		ads_undef(ftable[i].name,i);
    }

    return RTNORM;
    
	
   

}

int	dofun()
{
    int    val;
	int    rc;
     
	

	if ((val = ads_getfuncode()) < 0 || val > ELEMENTS(ftable))
        return RTERROR;
 
    rc = (*ftable[val].fptr)();
 
    return ((rc == RTNORM) ? RSRSLT:RSERR);

}

extern "C" AcRx::AppRetCode
acrxEntryPoint(AcRx::AppMsgCode msg, void* ptr) {

   if (ptr != NULL) {
        // We have been handed some kind of object
		// but we aren't going to do anything with it.
	}

	switch(msg) {
	case AcRx::kInitAppMsg:
            acrxUnlockApplication(ptr);
            acrxRegisterAppMDIAware(ptr);
			break;
        case AcRx::kInvkSubrMsg:
            dofun();
            break;
        case AcRx::kLoadADSMsg:
	
			{                  
			 funcload();
			 ads_printf("\nFAN by B.O Drugge  (C) Copyright 2012\n");
           
   
			}
			break;
		case AcRx::kUnloadADSMsg:
			funcunload();
			ads_printf("Unloading.\n");
			break;
	case AcRx::kUnloadAppMsg:
        default:
	    break;
    }
    return AcRx::kRetOK;
}





ads_real dabs(ads_real a)
{return sqrt(a*a);}




void 
cross (ads_point v,ads_point w,ads_point n)
{
n[X]=v[Y]*w[Z]-v[Z]*w[Y];
n[Y]=v[Z]*w[X]-v[X]*w[Z];
n[Z]=v[X]*w[Y]-v[Y]*w[X];
}

void 
subb (ads_point v,ads_point w,ads_point n)
{
n[X]=v[X]-w[X];
n[Y]=v[Y]-w[Y];
n[Z]=v[Z]-w[Z];
}

void
add (ads_point v,ads_point w,ads_point n)
{
n[X]=v[X]+w[X];
n[Y]=v[Y]+w[Y];
n[Z]=v[Z]+w[Z];
}


ads_real dotprod(ads_point v,ads_point w)
{
ads_real back;
back = v[X]*w[X]+v[Y]*w[Y]+v[Z]*w[Z];
return back;
}

ads_real dabsprod(ads_point v,ads_point w)
{
ads_real back;
back=sqrt(v[X]*v[X]+v[Y]*v[Y]+v[Z]*v[Z])*sqrt(w[X]*w[X]+w[Y]*w[Y]+w[Z]*w[Z]);
return back;
}

ads_real angle(ads_point v,ads_point w)
{
ads_real back;
back = acos((dotprod(v,w)/dabsprod(v,w)));
return back;
}                                                                


void 
trans (ads_point p1,ads_point p2,ads_real a1,ads_real a2)
{
ads_real v,u,l1;
u=atan(p1[Y]/(p1[X]+err2));
v=atan(p1[Z]/(sqrt(p1[X]*p1[X]+p1[Y]*p1[Y])+err2));
if (p1[X]<0) 
 a1=a1+pi;
a1=a1+u;
a2=a2+v;
l1=sqrt(p1[X]*p1[X]+p1[Y]*p1[Y]+p1[Z]*p1[Z]);
p2[X]=l1*cos(a1)*cos(a2);
p2[Y]=l1*sin(a1)*cos(a2);
p2[Z]=l1*sin(a2);
}





int spiral()
{
ads_point pins,pt,result;
ads_real r1,r2,h,dr,dh,k1,k2;
struct resbuf *ebuf,rb,fromrb,torb;
int i1,n,test;
 
ads_getvar("SURFTAB1",&rb);
n=rb.resval.rint;
fromrb.restype = RTSHORT;
fromrb.resval.rint = 0;
torb.restype = RTSHORT;
torb.resval.rint = 1;
  
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nEnter radius1:",&r1);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nEnter radius2:",&r2);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nEnter hight of spiral:",&h);
if (test == RTCAN)
 return RTNORM;
test=ads_getpoint(NULL,"\nEnter insertion point:",pins);
if (test == RTCAN)
 return RTNORM;

k2=2*pi/n;
dr=(r2-r1)/n;
dh=h/n;
ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,8,NULL);
ads_entmake(ebuf);
k1=0;
n++;
for (i1=0;i1<n;i1++)
  {
  pt[X]=(r1+dr*i1)*cos(k1);
  pt[Y]=(r1+dr*i1)*sin(k1);
  pt[Z]=dh*i1;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,32,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ads_relrb(ebuf);
return RTNORM;
}

int hus()
{
ads_point pins,pt,result;
ads_real q,dq,vhus,r1,r,h,k1,k2,A;
struct resbuf *ebuf,rb,fromrb,torb;
int i1,n,test;
 
ads_getvar("SURFTAB1",&rb);
n=rb.resval.rint;
fromrb.restype = RTSHORT;
fromrb.resval.rint = 0;
torb.restype = RTSHORT;
torb.resval.rint = 1;
  
vhus=40.0;


ads_initget(5,NULL);
test=ads_getdist(NULL,"\nEnter the flow:",&q);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nAnge grundradie:",&r1);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nAnge utloppshöjd:",&h);
if (test == RTCAN)
 return RTNORM;
test=ads_getpoint(NULL,"\nEnter insertion point:",pins);
if (test == RTCAN)
 return RTNORM;


dq=q/n;
k2=2*pi/n;
n--;

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=r1*cos(k1);
  pt[Y]=r1*sin(k1);
  pt[Z]=0.0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+r+r*sin(pi/8))*cos(k1);
  pt[Y]=(r1+r+r*sin(pi/8))*sin(k1);
  pt[Z]=0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+r+r*sin(pi/8))*cos(k1);
  pt[Y]=(r1+r+r*sin(pi/8))*sin(k1);
  pt[Z]=0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+2*r)*cos(k1);
  pt[Y]=(r1+2*r)*sin(k1);
  pt[Z]=r-r*sin(pi/8);
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+2*r)*cos(k1);
  pt[Y]=(r1+2*r)*sin(k1);
  pt[Z]=r-r*sin(pi/8);
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+2*r)*cos(k1);
  pt[Y]=(r1+2*r)*sin(k1);
  pt[Z]=r+r*sin(pi/8);
  if (pt[Z]<h)
   pt[Z]=h;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+2*r)*cos(k1);
  pt[Y]=(r1+2*r)*sin(k1);
  pt[Z]=r+r*sin(pi/8);
  if (pt[Z]<h)
   pt[Z]=h;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+r+r*sin(pi/8))*cos(k1);
  pt[Y]=(r1+r+r*sin(pi/8))*sin(k1);
  pt[Z]=2*r;
  if (pt[Z]<h)
   pt[Z]=h;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+r+r*sin(pi/8))*cos(k1);
  pt[Y]=(r1+r+r*sin(pi/8))*sin(k1);
  pt[Z]=2*r;
  if (pt[Z]<h)
   pt[Z]=h;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+r-r*sin(pi/8))*cos(k1);
  pt[Y]=(r1+r-r*sin(pi/8))*sin(k1);
  pt[Z]=2*r;
  if (pt[Z]<h)
   pt[Z]=h;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=(r1+r-r*sin(pi/8))*cos(k1);
  pt[Y]=(r1+r-r*sin(pi/8))*sin(k1);
  pt[Z]=2*r;
  if (pt[Z]<h)
   pt[Z]=h;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=r1*cos(k1);
  pt[Y]=r1*sin(k1);
  pt[Z]=r+r*sin(pi/8);
  if (pt[Z]<h)
   pt[Z]=h;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=r1*cos(k1);
  pt[Y]=r1*sin(k1);
  pt[Z]=r+r*sin(pi/8);
  if (pt[Z]<h)
   pt[Z]=h;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A/pi)*1000;
  pt[X]=r1*cos(k1);
  pt[Y]=r1*sin(k1);
  pt[Z]=h;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ads_relrb(ebuf);
return RTNORM;
}

int inlp()
{
ads_point pins,pt,result;
ads_real r1,r2,k1,k2,rstep,h;
struct resbuf *ebuf,rb,fromrb,torb;
int i1,i2,n1,n2,test;
 
ads_getvar("SURFTAB1",&rb);
n1=rb.resval.rint;
ads_getvar("SURFTAB2",&rb);
n2=rb.resval.rint;

fromrb.restype = RTSHORT;
fromrb.resval.rint = 0;
torb.restype = RTSHORT;
torb.resval.rint = 1;
  


ads_initget(5,NULL);
test=ads_getdist(NULL,"\nInlopps diameter:",&r1);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nUtlopps diameter:",&r2);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nUtlopps höjd:",&h);
if (test == RTCAN)
 return RTNORM;
test=ads_getpoint(NULL,"\nEnter insertion point:",pins);
if (test == RTCAN)
 return RTNORM;

k2=2*pi/n1;
r1=r1/2;
r2=r2/2;
rstep=(r2-r1)/n2;

for (i2=0;i2<n2;i2++)
  {
  ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,48,71,2,72,n1,NULL);
  ads_entmake(ebuf);
  k1=0;
  for (i1=0;i1<n1;i1++)
    {
    pt[X]=r1*cos(k1);
    pt[Y]=r1*sin(k1);
    pt[Z]=h*r2/r1;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    k1=k1+k2;
    }
  k1=0;
  for (i1=0;i1<n1;i1++)
    {
    pt[X]=(r1+rstep)*cos(k1);
    pt[Y]=(r1+rstep)*sin(k1);
    pt[Z]=h*r2/(r1+rstep);
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    k1=k1+k2;
    }
  ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
  ads_entmake(ebuf);
  r1=r1+rstep;
  }


return RTNORM;
}

int hus1()
{
ads_point pins,pt,result,p1,p2,p3,p4;
ads_real q,dq,vhus,l,l1,r,r1,h,k1,k2,A,fdiam,kang, idiam,difdist,difbd, diabd;
struct resbuf *ebuf,rb,fromrb,torb;
int i1,n,test;
 
ads_getvar("SURFTAB1",&rb);
n=rb.resval.rint;
fromrb.restype = RTSHORT;
torb.restype = RTSHORT;


ads_initget(5,NULL);
test=ads_getreal("\nEnter the flow (m3/s):",&q);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getreal("\nScroll speed (m/s):",&vhus);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nInlet diameter (mm):",&idiam);
if (test == RTCAN)
 return RTNORM;
test=ads_getdist(NULL,"\nFan diameter (mm):",&fdiam);
if (test == RTCAN)
 return RTNORM;

test=ads_getdist(NULL,"\nDiameter on backshield (mm):",&diabd);
if (test == RTCAN)
 return RTNORM;

ads_initget(5,NULL);
test=ads_getdist(NULL,"\nDistance to diffusor (mm):",&difdist);
if (test == RTCAN)
 return RTNORM;
test=ads_getdist(NULL,"\nDistance to the backshield (mm):",&difbd);
if (test == RTCAN)
 return RTNORM;
 

ads_initget(5,NULL);
test=ads_getangle(NULL,"\nEnter the concentricity (deg):",&kang);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nLenght of diffusor (mm):",&l1);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getdist(NULL,"\nHight of outer fanwheel (mm):",&h);
if (test == RTCAN)
 return RTNORM;
test=ads_getpoint(NULL,"\nEnter insertion point:",pins);
if (test == RTCAN)
 return RTNORM;

dq=q/n;

k2=2*pi/n;
n++;


ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=-difbd;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=diabd/2.0*cos(k1);
  pt[Y]=diabd/2.0*sin(k1);
  pt[Z]=-difbd;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=0.0;
  pt[Y]=0.0;
  pt[Z]=-difbd;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=h/cos(kang);
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=h/cos(kang)+difdist;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);


ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=h+difdist;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=idiam/2.0*cos(k1);
  pt[Y]=idiam/2.0*sin(k1);
  pt[Z]=h+difdist+idiam/4.0-idiam*idiam/(fdiam*4);
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  pt[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  pt[Z]=-(l1+r)*sin(kang);
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);








ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  pt[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  pt[Z]=-(l1+r)*sin(kang);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;

  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;

  p1[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  p1[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  p1[Z]=-(l1+r)*sin(kang);
  subb(pt,p1,p2);

  p3[X]=0.0;
  p3[Y]=0.0;
  p3[Z]=1.0;
  cross(p3,p2,p4);
  cross(p2,p4,p3);
  
  if(r<h)
   r=h;
 
  l=sqrt(p3[X]*p3[X]+p3[Y]*p3[Y]+p3[Z]*p3[Z]);
  p3[X]=p3[X]/l*r;
  p3[Y]=p3[Y]/l*r;
  p3[Z]=p3[Z]/l*r;
  add(p3,p1,pt);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;

  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;

  p1[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  p1[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  p1[Z]=-(l1+r)*sin(kang);
  subb(pt,p1,p2);

  p3[X]=0.0;
  p3[Y]=0.0;
  p3[Z]=1.0;
  cross(p3,p2,p4);
  cross(p2,p4,p3);
  
  if(r<h)
   r=h;
 
  l=sqrt(p3[X]*p3[X]+p3[Y]*p3[Y]+p3[Z]*p3[Z]);
  p3[X]=p3[X]/l*r;
  p3[Y]=p3[Y]/l*r;
  p3[Z]=p3[Z]/l*r;
  add(p3,p1,pt);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;

  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;

  p1[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  p1[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  p1[Z]=-(l1+r)*sin(kang);
  subb(pt,p1,p2);
 
  p3[X]=0.0;
  p3[Y]=0.0;
  p3[Z]=1.0;
  cross(p3,p2,p4);
  cross(p2,p4,p3);
  r1=r;
  if(r<h)
   r1=h;
  l=sqrt(p3[X]*p3[X]+p3[Y]*p3[Y]+p3[Z]*p3[Z]);
  p3[X]=p3[X]/l*r1;
  p3[Y]=p3[Y]/l*r1;
  p3[Z]=p3[Z]/l*r1;

  l=sqrt(p2[X]*p2[X]+p2[Y]*p2[Y]+p2[Z]*p2[Z]);
  p2[X]=p2[X]/l*r;
  p2[Y]=p2[Y]/l*r;
  p2[Z]=p2[Z]/l*r;
  add(p3,p1,p4);
  add(p4,p2,pt);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;


  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;

  p1[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  p1[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  p1[Z]=-(l1+r)*sin(kang);
  subb(pt,p1,p2);
 
  p3[X]=0.0;
  p3[Y]=0.0;
  p3[Z]=1.0;
  cross(p3,p2,p4);
  cross(p2,p4,p3);
  r1=r;
  if(r<h)
   r1=h;
  l=sqrt(p3[X]*p3[X]+p3[Y]*p3[Y]+p3[Z]*p3[Z]);
  p3[X]=p3[X]/l*r1;
  p3[Y]=p3[Y]/l*r1;
  p3[Z]=p3[Z]/l*r1;

  l=sqrt(p2[X]*p2[X]+p2[Y]*p2[Y]+p2[Z]*p2[Z]);
  p2[X]=p2[X]/l*r;
  p2[Y]=p2[Y]/l*r;
  p2[Z]=p2[Z]/l*r;
  add(p3,p1,p4);
  add(p4,p2,pt);


  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0+l1*cos(kang)+h*sin(kang))*cos(k1);
  pt[Y]=(fdiam/2.0+l1*cos(kang)+h*sin(kang))*sin(k1);
  pt[Z]=-l1*sin(kang)+h*cos(kang);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0+l1*cos(kang)+h*sin(kang))*cos(k1);
  pt[Y]=(fdiam/2.0+l1*cos(kang)+h*sin(kang))*sin(k1);
  pt[Z]=-l1*sin(kang)+h*cos(kang);


  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0)*cos(k1);
  pt[Y]=(fdiam/2.0)*sin(k1);
  pt[Z]=h/cos(kang);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);
return RTNORM;
}



int prop1()
{
ads_point pins,pt,result;
ads_real alpha,beta,ndia,vrad,r1,r2,zstig;
ads_real k1,k2,l;
struct resbuf *ebuf,rb,fromrb,torb;
int sf1,sf2,i1,i2,i3,test,aving;
 
ads_getvar("SURFTAB1",&rb);
sf1=rb.resval.rint;
ads_getvar("SURFTAB2",&rb);
sf2=1;
fromrb.restype = RTSHORT;
fromrb.resval.rint = 0;
torb.restype = RTSHORT;
torb.resval.rint = 1;
  
ads_initget(7,NULL);
test=ads_getint("\nNumber of wings (st):",&aving);
if (test == RTCAN)
 return RTNORM;

ads_initget(7,NULL);
test=ads_getreal("\nVingvinkel (°):",&alpha);
if (test == RTCAN)
 return RTNORM;

ads_initget(7,NULL);
test=ads_getreal("\nNav diameter (mm):",&ndia);
if (test == RTCAN)
 return RTNORM;

ads_initget(7,NULL);
test=ads_getreal("\nVing radie (mm):",&vrad);
if (test == RTCAN)
 return RTNORM;

ads_initget(7,NULL);
test=ads_getreal("\nBladvikel (ø):",&beta);
if (test == RTCAN)
 return RTNORM;

test=ads_getpoint(NULL,"\nInsättnings punkt:",pins);
if (test == RTCAN)
 return RTNORM;

alpha=alpha*pi/180.0;
beta=beta*pi/180.0;
r1=ndia/2.0;
r2=vrad;

for (i1=0;i1<aving;i1++)
  {
  ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,sf2+1,72,sf1+1,NULL);
  ads_entmake(ebuf);
  for (i2=0;i2<sf2+1;i2++)
    {
    l=(r1+(r2-r1)/sf2*i2)*sin(alpha/sf1)/sin((pi-alpha/sf1)/2);
    k2=i1*2*pi/aving;
    k1=k2;
    zstig=tan(beta)*r1/(r1+(r2-r1)/sf2*i2);
    for (i3=0;i3<sf1+1;i3++)
      {
      pt[X]=(r1+(r2-r1)/sf2*i2)*cos(k1);
      pt[Y]=(r1+(r2-r1)/sf2*i2)*sin(k1);
      pt[Z]=zstig*l*i3;
      add(pt,pins,pt);
      ads_trans(pt,&torb,&fromrb,0,result);
      ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
      ads_entmake(ebuf);
      k1=k2+alpha/sf1*(i3+1);
      }
    }
  ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
  ads_entmake(ebuf);
  }

return RTNORM;
}

int fan1()
{
ads_point pins,pt,result;
ads_real u1,u2,kslip,beta1,beta2,press,flow,di,dy,n,dp,w;
ads_real dens,alpha1,alpha2,dalp,vi,h1,h2,r1,r2,k1,rstep,stodola;
struct resbuf *ebuf,rb,fromrb,torb;
int sf1,sf2,i1,i2,test,nv;
 
ads_getvar("SURFTAB1",&rb);
sf1=rb.resval.rint;
ads_getvar("SURFTAB2",&rb);
sf2=rb.resval.rint;
fromrb.restype = RTSHORT;
fromrb.resval.rint = 0;
torb.restype = RTSHORT;
torb.resval.rint = 1;
  
ads_initget(5,NULL);
test=ads_getreal("\nPressure (Pa):",&press);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getreal("\nFlow (m3/s):",&flow);
if (test == RTCAN)
 return RTNORM;

dens=1.2;
ads_initget(4,NULL);
test=ads_getreal("\nDensity (kg/m3) <1.2>:",&dens);
if (test == RTCAN)
 return RTNORM;
stodola=1.0;
ads_initget(4,NULL);
test=ads_getreal("\nStodola coef <1.0>:",&stodola);
if (test == RTCAN)
 return RTNORM;
ads_initget(5,NULL);
test=ads_getint("\nNumber of wings (st):",&nv);
if (test == RTCAN)
 return RTNORM;
n=0.0;
ads_initget(4,NULL);
test=ads_getreal("\nRotating speed (rpm):",&n);
if (test == RTCAN)
 return RTNORM;
if (n>0)
{
 ads_initget(4,NULL);
 test=ads_getreal("\nInlet speed (m/s):",&vi);
 if (test == RTCAN)
  return RTNORM;
}
if (n == 0.0)
 {
 ads_initget(5,NULL);
 test=ads_getdist(NULL,"\nInlet diameter (mm):",&di);
 if (test == RTCAN)
  return RTNORM;
 ads_initget(5,NULL);
 test=ads_getdist(NULL,"\nOutlet diameter (mm):",&dy);
 if (test == RTCAN)
  return RTNORM;
 }
test=ads_getpoint(NULL,"\nInsertion point:",pins);
if (test == RTCAN)
 return RTNORM;

if (n > 0.0)
 {
 r1=sqrt(flow/(vi*pi));
 w=2*pi*n/60;
 r2=r1;
 u1=w*r1;
 dp=0.0;
 while (dp < 2*press)
  {
  beta2=atan(2*vi/(w*r2));
  /*kslip=1/(1+1.5/nv*(1+beta2/(pi/3.0))/(1-r1*r1/(r2*r2)));*/
  u2=w*r2;
  dp=(u2*u2*(1-stodola*pi*sin(beta2)/nv)-u1*u1)*dens;
  r2=r2+.000001;
  }
 }



if (n == 0)
 {
 r1=di/2000;
 r2=dy/2000;
 vi=flow/(r1*r1*pi);
 dp=0.0;
 w=2*pi/60;
 while (dp < 2*press)
  {
  beta2=atan(2*vi/(w*r2));
  /*kslip=1/(1+1.5/nv*(1+beta2/(pi/3.0))/(1-di*di/(dy*dy)));*/
  dp=((w*r2)*(w*r2)*(1-stodola*pi*sin(beta2)/nv)-(w*r1)*(w*r1))*dens;
  w=w+2*pi/60000;
  }
 }
beta1=atan(2*vi/(w*r1));


h1=r1*r1/(2*r1);
h2=r1*r1/(2*r2);
alpha1=r1/(2*vi)*w;
alpha2=r2/(2*vi)*w;
dalp=(alpha2-alpha1)/sf1;
if ((2*vi/(w*r1))>1.0)
  ads_printf("\nWARNING !!! Inlet blade angle greater than 45 degrees");

ads_printf("\nInlett blade angle:%f\nOutlet blade angle:%f",beta1*180/pi,beta2*180/pi);

for (i2=0;i2<nv;i2++)
  {
  ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,sf1+1,NULL);
  ads_entmake(ebuf);
  k1=-2*pi/nv*i2;
  for (i1=0;i1<sf1+1;i1++)
    {
    pt[X]=(r1+(r2-r1)/sf1*i1)*cos(k1)*1000;
    pt[Y]=(r1+(r2-r1)/sf1*i1)*sin(k1)*1000;
    pt[Z]=0.0;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    k1=k1-dalp;
    }
  k1=-2*pi/nv*i2;
  for (i1=0;i1<sf1+1;i1++)
    {
    pt[X]=(r1+(r2-r1)/sf1*i1)*cos(k1)*1000;
    pt[Y]=(r1+(r2-r1)/sf1*i1)*sin(k1)*1000;
    pt[Z]=(h1+(h2-h1)*i1/sf1)*1000;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    k1=k1-dalp;
    }
  ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
  ads_entmake(ebuf);
  }
ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,48,71,2,72,sf1,NULL);
ads_entmake(ebuf);
for (i1=0;i1<sf1;i1++)
  {
  pt[X]=0.0;
  pt[Y]=0.0;
  pt[Z]=0.0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  }
for (i1=0;i1<sf1;i1++)
  {
  pt[X]=r2*cos(2*pi/sf1*i1)*1000;
  pt[Y]=r2*sin(2*pi/sf1*i1)*1000;
  pt[Z]=0.0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);
sf2=1;
rstep=(r2-r1)/sf2;
for (i2=0;i2<sf2;i2++)
  {
  ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,48,71,2,72,sf1,NULL);
  ads_entmake(ebuf);
  for (i1=0;i1<sf1;i1++)
    {
    pt[X]=r1*cos(2*pi/sf1*i1)*1000;
    pt[Y]=r1*sin(2*pi/sf1*i1)*1000;
    pt[Z]=h2*r2/r1*1000;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    }
  for (i1=0;i1<sf1;i1++)
    {
    pt[X]=(r1+rstep)*cos(2*pi/sf1*i1)*1000;
    pt[Y]=(r1+rstep)*sin(2*pi/sf1*i1)*1000;
    pt[Z]=h2*r2/(r1+rstep)*1000;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    }
  ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
  ads_entmake(ebuf);
  r1=r1+rstep;
  }

ads_printf("\nRotating speed:%f\n",w*30/pi);
return RTNORM;
}


int fnik()
{

FILE *fp;	
float lres[25];/*antalet siffror i *.FAN filen ökar du slisar så ökar med antalet slitsar*2+17
			för närvarande är det 4 st slitsar 4*2+17 =25*/
ads_real u1,u2,kslip,beta1,beta2,press,flow,di,dy,dp,w;
ads_real dens,alpha1,alpha2,dalp,vi,h1,h2,r2,k1,k2,rstep,stodola;
struct resbuf *ebuf,rb,fromrb,torb,*res;
int sf1,sf2,i1,i2,test,nv,i,ii,n,slitz,cont;
ads_point pins,pt,result,p1,p2,p3,p4;
ads_real q,dq,vhus,l,l1,r,r1,h,A,fdiam,kang, idiam,difdist,difbd, diabd,rr1[9]; 
short surf;
ads_name ent1,ent2,rent;





res=ads_newrb(RTSTR);
test=ads_getfiled("FANnik an expertsystem for buliding fans by B.O Drugge","*","FAN",8,res);
if (test == RTERROR)
 return RTNORM;



(fp=fopen(res->resval.rstring,"r"));

 for (i=0; i<26; i++) /*26 anger antal siffror i fan filen ökar du den siffran så ökar du
					  antalet slitsar, om du vill ha 20 slitsar så blir talet 58*/
	   fscanf(fp,"%f",&lres[i]); 

		 

surf=lres[13];
 

ads_command(RTSTR,"SURFTAB1",RTSHORT,surf,NULL);

ads_getvar("SURFTAB1",&rb);
sf1=rb.resval.rint;
ads_getvar("SURFTAB2",&rb);
sf2=rb.resval.rint;
fromrb.restype = RTSHORT;
fromrb.resval.rint = 0;
torb.restype = RTSHORT;
torb.resval.rint = 1;
  

flow=lres[2]/2.0;
dens=lres[4];
stodola=lres[6];
nv=lres[5];
di=lres[0];
dy=lres[1];



test=ads_getpoint(NULL,"\nInsertion point:",pins);
if (test == RTCAN)
 return RTNORM;





r1=di/2000;
r2=dy/2000; 
vi=flow/(r1*r1*pi);

w=2*pi*lres[3]/60;
 
beta1=atan(2*vi/(w*r1));


h1=r1*r1/(2*r1);
h2=r1*r1/(2*r2);
alpha1=r1/(2*vi)*w;
alpha2=r2/(2*vi)*w;
dalp=(alpha2-alpha1)/sf1;


for (i2=0;i2<nv;i2++)
  {
  ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,sf1+1,NULL);
  ads_entmake(ebuf);
  k1=-2*pi/nv*i2;
  for (i1=0;i1<sf1+1;i1++)
    {
    pt[X]=(r1+(r2-r1)/sf1*i1)*cos(k1)*1000;
    pt[Y]=(r1+(r2-r1)/sf1*i1)*sin(k1)*1000;
    pt[Z]=0.0;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    k1=k1-dalp;
    }
  k1=-2*pi/nv*i2;
  for (i1=0;i1<sf1+1;i1++)
    {
    pt[X]=(r1+(r2-r1)/sf1*i1)*cos(k1)*1000;
    pt[Y]=(r1+(r2-r1)/sf1*i1)*sin(k1)*1000;
    pt[Z]=(h1+(h2-h1)*i1/sf1)*1000;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    k1=k1-dalp;
    }
  ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
  ads_entmake(ebuf);
  }
ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,48,71,2,72,sf1,NULL);
ads_entmake(ebuf);
for (i1=0;i1<sf1;i1++)
  {
  pt[X]=0.0;
  pt[Y]=0.0;
  pt[Z]=0.0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  }
for (i1=0;i1<sf1;i1++)
  {
  pt[X]=r2*cos(2*pi/sf1*i1)*1000;
  pt[Y]=r2*sin(2*pi/sf1*i1)*1000;
  pt[Z]=0.0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);
sf2=1;
rstep=(r2-r1)/sf2;
for (i2=0;i2<sf2;i2++)
  {
  ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,48,71,2,72,sf1,NULL);
  ads_entmake(ebuf);
  for (i1=0;i1<sf1;i1++)
    {
    pt[X]=r1*cos(2*pi/sf1*i1)*1000;
    pt[Y]=r1*sin(2*pi/sf1*i1)*1000;
    pt[Z]=h2*r2/r1*1000;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    }
  for (i1=0;i1<sf1;i1++)
    {
    pt[X]=(r1+rstep)*cos(2*pi/sf1*i1)*1000;
    pt[Y]=(r1+rstep)*sin(2*pi/sf1*i1)*1000;
    pt[Z]=h2*r2/(r1+rstep)*1000;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    }
  ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
  ads_entmake(ebuf);
  r1=r1+rstep;
  }

ads_getvar("SURFTAB1",&rb);
n=rb.resval.rint;
fromrb.restype = RTSHORT;
torb.restype = RTSHORT;


q=lres[2];
vhus=lres[2]/(lres[0]*lres[0]*pi/4000000.0);
idiam=lres[0];
fdiam=lres[1]+lres[10]*2.0;
diabd=lres[7];
difdist=lres[9];
difbd=lres[8];
kang=lres[11]*pi/180.0;
l1=lres[12];
h=idiam/fdiam*idiam/4.0;


test=ads_getpoint(NULL,"\nEnter insertion point:",pins);
if (test == RTCAN)
 return RTNORM;

dq=q/n;

k2=2*pi/n;
n++;


ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=-difbd;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=diabd/2.0*cos(k1);
  pt[Y]=diabd/2.0*sin(k1);
  pt[Z]=-difbd;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=0.0;
  pt[Y]=0.0;
  pt[Z]=-difbd;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=h/cos(kang);
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=h/cos(kang)+difdist;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);


ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=h+difdist;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=idiam/2.0*cos(k1);
  pt[Y]=idiam/2.0*sin(k1);
  pt[Z]=h+difdist+idiam/4.0-idiam*idiam/(fdiam*4);
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  pt[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  pt[Z]=-(l1+r)*sin(kang);
  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);








ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  pt[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  pt[Z]=-(l1+r)*sin(kang);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;

  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;

  p1[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  p1[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  p1[Z]=-(l1+r)*sin(kang);
  subb(pt,p1,p2);

  p3[X]=0.0;
  p3[Y]=0.0;
  p3[Z]=1.0;
  cross(p3,p2,p4);
  cross(p2,p4,p3);
  
  if(r<h)
   r=h;
 
  l=sqrt(p3[X]*p3[X]+p3[Y]*p3[Y]+p3[Z]*p3[Z]);
  p3[X]=p3[X]/l*r;
  p3[Y]=p3[Y]/l*r;
  p3[Z]=p3[Z]/l*r;
  add(p3,p1,pt);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;

  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;

  p1[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  p1[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  p1[Z]=-(l1+r)*sin(kang);
  subb(pt,p1,p2);

  p3[X]=0.0;
  p3[Y]=0.0;
  p3[Z]=1.0;
  cross(p3,p2,p4);
  cross(p2,p4,p3);
  
  if(r<h)
   r=h;
 
  l=sqrt(p3[X]*p3[X]+p3[Y]*p3[Y]+p3[Z]*p3[Z]);
  p3[X]=p3[X]/l*r;
  p3[Y]=p3[Y]/l*r;
  p3[Z]=p3[Z]/l*r;
  add(p3,p1,pt);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);
  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;

  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;

  p1[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  p1[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  p1[Z]=-(l1+r)*sin(kang);
  subb(pt,p1,p2);
 
  p3[X]=0.0;
  p3[Y]=0.0;
  p3[Z]=1.0;
  cross(p3,p2,p4);
  cross(p2,p4,p3);
  r1=r;
  if(r<h)
   r1=h;
  l=sqrt(p3[X]*p3[X]+p3[Y]*p3[Y]+p3[Z]*p3[Z]);
  p3[X]=p3[X]/l*r1;
  p3[Y]=p3[Y]/l*r1;
  p3[Z]=p3[Z]/l*r1;

  l=sqrt(p2[X]*p2[X]+p2[Y]*p2[Y]+p2[Z]*p2[Z]);
  p2[X]=p2[X]/l*r;
  p2[Y]=p2[Y]/l*r;
  p2[Z]=p2[Z]/l*r;
  add(p3,p1,p4);
  add(p4,p2,pt);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;


  pt[X]=fdiam/2.0*cos(k1);
  pt[Y]=fdiam/2.0*sin(k1);
  pt[Z]=0.0;

  p1[X]=(fdiam/2.0+(l1+r)*cos(kang))*cos(k1);
  p1[Y]=(fdiam/2.0+(l1+r)*cos(kang))*sin(k1);
  p1[Z]=-(l1+r)*sin(kang);
  subb(pt,p1,p2);
 
  p3[X]=0.0;
  p3[Y]=0.0;
  p3[Z]=1.0;
  cross(p3,p2,p4);
  cross(p2,p4,p3);
  r1=r;
  if(r<h)
   r1=h;
  l=sqrt(p3[X]*p3[X]+p3[Y]*p3[Y]+p3[Z]*p3[Z]);
  p3[X]=p3[X]/l*r1;
  p3[Y]=p3[Y]/l*r1;
  p3[Z]=p3[Z]/l*r1;

  l=sqrt(p2[X]*p2[X]+p2[Y]*p2[Y]+p2[Z]*p2[Z]);
  p2[X]=p2[X]/l*r;
  p2[Y]=p2[Y]/l*r;
  p2[Z]=p2[Z]/l*r;
  add(p3,p1,p4);
  add(p4,p2,pt);


  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0+l1*cos(kang)+h*sin(kang))*cos(k1);
  pt[Y]=(fdiam/2.0+l1*cos(kang)+h*sin(kang))*sin(k1);
  pt[Z]=-l1*sin(kang)+h*cos(kang);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0+l1*cos(kang)+h*sin(kang))*cos(k1);
  pt[Y]=(fdiam/2.0+l1*cos(kang)+h*sin(kang))*sin(k1);
  pt[Z]=-l1*sin(kang)+h*cos(kang);


  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  A=dq*i1/vhus;
  r=sqrt(A)*1000;
  pt[X]=(fdiam/2.0)*cos(k1);
  pt[Y]=(fdiam/2.0)*sin(k1);
  pt[Z]=h/cos(kang);

  add(pt,pins,pt);
  ads_trans(pt,&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);























test=ads_getpoint(NULL,"\nInsertion point:",pins);
if (test == RTCAN)
 return RTNORM;
ads_command(RTSTR,"SURFTAB1",RTSHORT,surf*2,NULL);/*Här gångar jag surftab med 2 (surf*2)
												  och får då dubbelt surftab på beräkningskonan*/
ads_getvar("SURFTAB1",&rb);
n=rb.resval.rint;
k2=2*pi/n;

ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,48,71,2,72,n,NULL);
ads_entmake(ebuf);
k1=0;
beta1=atan((idiam*idiam/(4*idiam)-idiam*idiam/(4*fdiam))/(0.5*fdiam-0.5*idiam));
for (i1=0;i1<n;i1++)
  {
  
  
  ptup2[i1][X]=(lres[1]/2+cos(pi/2-beta1)*lres[14]/2)*cos(k1);
  ptup2[i1][Y]=(lres[1]/2.0+cos(pi/2.0-beta1)*lres[14]/2.0)*sin(k1);
  ptup2[i1][Z]=idiam*idiam/(4*lres[1])+lres[14]*cos(beta1)/2.0;


  add(pins,ptup2[i1],ptup2[i1]);
  ads_trans(ptup2[i1],&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);

  k1=k1+k2;
  }
k1=0;
for (i1=0;i1<n;i1++)
  {
  ptdw2[i1][X]=(idiam/2-cos(pi/2+beta1)*lres[14]/2)*cos(k1);
  ptdw2[i1][Y]=(idiam/2.0-cos(pi/2.0+beta1)*lres[14]/2.0)*sin(k1);
  ptdw2[i1][Z]=idiam/4.0+lres[14]*cos(beta1)/2.0;

  
  add(pins,ptdw2[i1],ptdw2[i1]);
  ads_trans(ptdw2[i1],&torb,&fromrb,0,result);

  ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
  ads_entmake(ebuf);
  k1=k1+k2;
  }
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
ads_entmake(ebuf);
ptup2[i1][X]=ptup2[0][X];
ptup2[i1][Y]=ptup2[0][Y];
ptup2[i1][Z]=ptup2[0][Z];
ptdw2[i1][X]=ptdw2[0][X];
ptdw2[i1][Y]=ptdw2[0][Y];
ptdw2[i1][Z]=ptdw2[0][Z];



sf1=sf1/4;/*Här dividerar jag surftab1 med 4 så att jag får fyra ggr mindre indelningar på slitsar
		  och vinge för surftab1=60 blir det 16 indelningar i vingbitarna*/

r1=di/2000;
r2=dy/2000; 
slitz=lres[17]*2+2;
rr1[0]=r1;
for (i=1; i<slitz-1; i++)
 rr1[i]=lres[17+i]/1000.0;
rr1[i]=r2;





k2=2*pi;
for (i=0;i<slitz-1;i++)
{
vi=flow/(r1*r1*pi);

w=2*pi*lres[3]/60;
 
beta1=atan(2*vi/(w*rr1[i]));


h1=rr1[0]*rr1[0]/(2*rr1[0]);
h2=rr1[0]*rr1[0]/(2*rr1[slitz-1]);
alpha1=rr1[i]/(2*vi)*w;
alpha2=rr1[i+1]/(2*vi)*w;
dalp=(alpha2-alpha1)/sf1;


  
  ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,16,71,2,72,sf1+1,NULL);
  ads_entmake(ebuf);
  
  
  k1=k2;
  for (i1=0;i1<sf1+1;i1++)
    {
    pt[X]=(rr1[i]+(rr1[i+1]-rr1[i])/sf1*i1)*cos(k1)*1000;
    pt[Y]=(rr1[i]+(rr1[i+1]-rr1[i])/sf1*i1)*sin(k1)*1000;
    pt[Z]=0.0;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    k1=k1-dalp;
    }
  k1=k2;
  for (i1=0;i1<sf1+1;i1++)
    {
    pt[X]=(rr1[i]+(rr1[i+1]-rr1[i])/sf1*i1)*cos(k1)*1000;
    pt[Y]=(rr1[i]+(rr1[i+1]-rr1[i])/sf1*i1)*sin(k1)*1000;
    kslip=(h2-h1)/(rr1[slitz-1]-rr1[0]); 
	pt[Z]=(2*h2+h1+rr1[0]*kslip+kslip*rr1[i]+kslip*(rr1[i+1]-rr1[i])*i1/sf1)*1000;
    add(pt,pins,pt);
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,70,64,NULL);
    ads_entmake(ebuf);
    k1=k1-dalp;
    }
  
	  k2=k2-dalp*sf1;
  
  
  ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
  ads_entmake(ebuf);
  
}



















for (ii=0;ii<nv;ii++)
{
vi=flow/(r1*r1*pi);

w=2*pi*lres[3]/60;
 



for (i=1;i<slitz-2;i=i+2)
{



alpha1=rr1[0]/(2*vi)*w;
alpha2=rr1[i]/(2*vi)*w;
dalp=(alpha2-alpha1)/sf1;

beta1=alpha2-alpha1;
k2=2*pi-2*pi/nv*ii-dalp*sf1;



	
	
	
	
	
	
h1=rr1[0]*rr1[0]/(2*rr1[0]);
h2=rr1[0]*rr1[0]/(2*rr1[slitz-1]);
alpha1=rr1[i]/(2*vi)*w;
alpha2=rr1[i+1]/(2*vi)*w;
dalp=(alpha2-alpha1)/sf1;




  
  ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,1,NULL);
  ads_entmake(ebuf);
  
  test=0;
  k1=k2;
  for (i1=0;i1<sf1+1;i1++)
    {
    
	   kslip=(rr1[i]+(rr1[i+1]-rr1[i])/sf1*i1);
    
	pt[X]=kslip*cos(k1)*1000.0;
    pt[Y]=kslip*sin(k1)*1000.0;
    pt[Z]=0.0;
 	kslip=(rr1[i]+(rr1[i+1]-rr1[i])/sf1*(i1+1));
	p1[X]=kslip*cos(k1-dalp)*1000.0;
    p1[Y]=kslip*sin(k1-dalp)*1000.0;
    p1[Z]=0.0;

	p2[X]=p1[X];
	p2[Y]=p1[Y];
    p2[Z]=100.0;
subb(p2,pt,p3);
subb (p1,pt,p4);
cross (p3,p4,p1);
p2[X]=0.0;
p2[Y]=0.0;
p2[Z]=0.0;

l1=ads_distance(p1,p2);     


p1[X]=p1[X]/l1*lres[15]/2.0;
p1[Y]=p1[Y]/l1*lres[15]/2.0;
p1[Z]=0.0;


add(p1,pt,p2);
add(p2,pins,pt);
ptdw1[i1][X]=pt[X];
ptdw1[i1][Y]=pt[Y];
ptdw1[i1][Z]=lres[0];
      

  ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,NULL);
    ads_entmake(ebuf);
	k1=k1-dalp;  
  
  }
 test=i1; 
k1=k1+dalp;
    for (i1=sf1 ;i1>-1 ;i1--)
    {
		
    kslip=(rr1[i]+(rr1[i+1]-rr1[i])/sf1*i1);
    
	pt[X]=kslip*cos(k1)*1000.0;
    pt[Y]=kslip*sin(k1)*1000.0;
    pt[Z]=0.0;
     kslip=(rr1[i]+(rr1[i+1]-rr1[i])/sf1*(i1-1)); 
	p1[X]=kslip*cos(k1+dalp)*1000.0;
    p1[Y]=kslip*sin(k1+dalp)*1000.0;
    p1[Z]=0.0;

	p2[X]=p1[X];
	p2[Y]=p1[Y];
    p2[Z]=100.0;
subb(p2,pt,p3);
subb (p1,pt,p4);
cross (p3,p4,p1);
p2[X]=0.0;
p2[Y]=0.0;
p2[Z]=0.0;

l1=ads_distance(p1,p2);     
p1[X]=p1[X]/l1*lres[15]/2.0;
p1[Y]=p1[Y]/l1*lres[15]/2.0;
p1[Z]=0.0;


add(p1,pt,p2);
add(p2,pins,pt);
ptdw1[test+sf1-i1][X]=pt[X];
ptdw1[test+sf1-i1][Y]=pt[Y];
ptdw1[test+sf1-i1][Z]=lres[0];
    
	
	
	
	
	
    ads_trans(pt,&torb,&fromrb,0,result);
    ebuf=ads_buildlist(RTDXF0,"VERTEX",10,result,NULL);
    ads_entmake(ebuf);
    k1=k1+dalp;
    

	
	}


  
	  
  
  
  ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
  ads_entmake(ebuf);
  






for (i1=0; i1<test+sf1+1;i1++)
{ 

ptup1[i1][X]=ptdw1[i1][X];
ptup1[i1][Y]=ptdw1[i1][Y];
ptup1[i1][Z]=0.0;

}

trim(test+sf1+1,n+1);


ebuf=ads_buildlist(RTDXF0,"POLYLINE",66,1,70,9,NULL);
  ads_entmake(ebuf);
for (i1=0; i1<test+sf1+1; i1++)
{
	ebuf=ads_buildlist(RTDXF0,"VERTEX",10,ptdw1[i1],70,32,NULL);
    ads_entmake(ebuf);
  
}  
  
ebuf=ads_buildlist(RTDXF0,"SEQEND",NULL);
  ads_entmake(ebuf);  
  
}





 
 

}  
  
ads_command(RTSTR,"CIRCLE",RTPOINT,pins,RTREAL,lres[1]/2.0,NULL);
ads_command(RTSTR,"SURFTAB1",RTSHORT,surf,NULL); 

ads_relrb(ebuf);
  return RTNORM;

}



void
trim(int n1,int n2)
{
int i1,i2;
ads_point pt,pt1,al,bl,ap,ap1,bp,bp1,norm;
ads_real sumt,t,alpha,alpha1,alpha2,beta,beta1,beta2;

for (i1=0;i1<n1;i1++)                                       /*Step lines of the trimming object */
  {
  for (i2=0;i2<n2-1;i2++)   /* Do a loop of the all surfaces the line ar intesected*/
    {
    subb(ptup1[i1],ptdw2[i2],al);
    subb(ptdw1[i1],ptdw2[i2],bl);
    subb(ptup2[i2],ptdw2[i2],ap);
    subb(ptup2[i2+1],ptdw2[i2],bp);
    subb(al,bl,pt);
    cross(ap,bp,norm);
    sumt=pt[X]*norm[X]+pt[Y]*norm[Y]+pt[Z]*norm[Z];
    if (sumt != 0)
     {
     t=-(bl[X]*norm[X]+bl[Y]*norm[Y]+bl[Z]*norm[Z])/sumt;
     if ((t>0) && (t<1))
      {
      pt[X]=pt[X]*t+bl[X];
      pt[Y]=pt[Y]*t+bl[Y];
      pt[Z]=pt[Z]*t+bl[Z];
      alpha = angle(ap,bp);
      alpha1= angle(ap,pt);
      alpha2= angle(bp,pt);
      if ( (dabs(pt[X])<err2 && dabs(pt[Y])<err2 && dabs(pt[Z])<err2) ||
           dabs(alpha1+alpha2-alpha) < err)
       {
       subb(pt,bp,pt1);
       subb(ap,bp,ap1);
       bp1[X]=-bp[X];
       bp1[Y]=-bp[Y];
       bp1[Z]=-bp[Z];
       beta = angle(ap1,bp1);
       beta1= angle(ap1,pt1);
       beta2= angle(bp1,pt1);
       if ( (dabs(pt1[X])<err2 && dabs(pt1[Y])<err2 && dabs(pt1[Z])<err2) ||
           dabs(beta1+beta2-beta) < err)
        {
        add(pt,ptdw2[i2],ptdw1[i1]);
        }
       }
      }
     }
    }
  for (i2=1;i2<n2;i2++)
    {
    subb(ptup1[i1],ptup2[i2],al);
    subb(ptdw1[i1],ptup2[i2],bl);
    subb(ptdw2[i2],ptup2[i2],ap);
    subb(ptdw2[i2-1],ptup2[i2],bp);
    subb(al,bl,pt);
    cross(ap,bp,norm);
    sumt=pt[X]*norm[X]+pt[Y]*norm[Y]+pt[Z]*norm[Z];
    if (sumt != 0)
     {
     t=-(bl[X]*norm[X]+bl[Y]*norm[Y]+bl[Z]*norm[Z])/sumt;
     if ((t>0) && (t<1))
      {
      pt[X]=pt[X]*t+bl[X];
      pt[Y]=pt[Y]*t+bl[Y];
      pt[Z]=pt[Z]*t+bl[Z];
      alpha = angle(ap,bp);
      alpha1= angle(ap,pt);
      alpha2= angle(bp,pt);
      if ( (dabs(pt[X])<err2 && dabs(pt[Y])<err2 && dabs(pt[Z])<err2) ||
           dabs(alpha1+alpha2-alpha) < err)
       {
       subb(pt,bp,pt1);
       subb(ap,bp,ap1);
       bp1[X]=-bp[X];
       bp1[Y]=-bp[Y];
       bp1[Z]=-bp[Z];
       beta = angle(ap1,bp1);
       beta1= angle(ap1,pt1);
       beta2= angle(bp1,pt1);
       if ( (dabs(pt1[X])<err2 && dabs(pt1[Y])<err2 && dabs(pt1[Z])<err2) ||
           dabs(beta1+beta2-beta) < err)
        {
        add(pt,ptup2[i2],ptdw1[i1]);
        }
       }
      }
     }
    }
  }
}