/prep7 ET,1,45 csys,1 pi=4*atan(1) rad1=20 rad2=25 delta=rad2-rad1 dicke=1 amplitude=2 anz=10 teiler=360 elsize=360/teiler *do,i,1,teiler,1 ys=i zs=amplitude*sin(anz*pi*i/180) k,i,rad1,ys,zs k,i+teiler,rad2,ys,zs *enddo *do,i,1,teiler,1 *if,i,eq,teiler,then CSKP,100+i,0,i,1,i+teiler *else CSKP,100+i,0,i,i+1,i+teiler *endif CSYS,100+i k,2*teiler+i,0,0,-dicke k,3*teiler+i,0,delta,-dicke *enddo *do,i,1,teiler,1 *if,i,eq,teiler,then V,teiler,teiler*2,teiler+1,1,3*teiler,4*teiler,3*teiler+1,2*teiler+1 *else V,i,i+teiler,i+teiler+1,i+1,2*teiler+i,3*teiler+i,3*teiler+i+1,2*teiler+i+1 *endif *enddo !esize,elsize MSHAPE,0,3D MSHKEY,1 vmesh,all