CC---PROGRAM TO GENERATE NODE FILE SUBROUTINE ABQMAIN C IMPLICIT REAL*8(A-H,O-Z) PARAMETER (ZERO=0.0d0,ONE=1.0d0,TWO=2.0d0,TEN=10.0d0) OPEN(UNIT=16,STATUS='NEW',FILE='RINGSHL.NOD') FI=ZERO NDIVS=4 TMP=ONE DFI=ASIN(TMP)/DBLE(NDIVS) R=TEN*TEN A=ONE/TEN NNODE=NDIVS+1 DO 10 K1=1,NNODE RR=R+A*COS(TWO*FI) X=RR*COS(FI) Y=RR*SIN(FI) KNODE=100+K1 Z=ZERO WRITE(16,1)KNODE,X,Y,Z KNODE=200+K1 Z=ONE/TWO WRITE(16,1)KNODE,X,Y,Z KNODE=300+K1 Z=ONE WRITE(16,1)KNODE,X,Y,Z FI=FI+DFI 10 CONTINUE 1 FORMAT(I5,1H,,F15.8,1H,,F15.8,1H,,F15.8) REWIND 16 STOP END