Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- !
- ! CHECK BODY INPUTS TO MEET SPECIFICATIONS NOTED ABOVE
- ! IN READING BODY NAMELIST DATA
- !
- NNX=NX+0.5
- IF(ELLIP.EQ.0.)ELLIP=UNUSED
- IF(ELLIP.EQ.UNUSED) ELLIP=1.0
- IF(R(1).NE.UNUSED .AND. S(1).EQ.UNUSED .AND. P(1).EQ.UNUSED) &
- & GO TO 1310
- IF(S(1).NE.UNUSED .AND. R(1).EQ.UNUSED .AND. P(1).EQ.UNUSED) &
- & GO TO 1330
- IF(R(1).NE.UNUSED .AND. S(1).NE.UNUSED .AND. P(1).NE.UNUSED) &
- & GO TO 1350
- !
- ! HERE FOR RADIUS OR AREA INPUT WITH ANOTHER VARIABLE.
- ! ONLY R OR S OR R,S,P CAN BE USED. SET THE GONOGO FLAG,
- ! WRITE A MESSAGE AND CONTINUE.
- !
- GONOGO=.FALSE.
- WRITE(6,1510)
- GO TO 1370
- !
- ! HALF WIDTH, R, SPECIFIED
- !
- 1310 DO 1320 K=1,NNX
- S(K)=PI*R(K)**2*ELLIP
- P(K)=2*PI*R(K)*SQRT((ELLIP**2+1.)/2.)
- R(K)=R(K)*SQRT(ELLIP)
- 1320 END DO
- GO TO 1370
- !
- ! CROSS SECTIONAL AREA, S, SPECIFIED
- !
- 1330 DO 1340 K=1,NNX
- R(K)=SQRT(S(K)/PI)
- P(K)=2.*PI*SQRT((S(K)/ELLIP)*(ELLIP**2+1.)/2./PI)
- 1340 END DO
- GO TO 1370
- !
- ! R, S, AND P ALL SPECIFIED
- !
- 1350 DO 1360 K=1,NNX
- R(K)=R(K)*SQRT(ELLIP)
- 1360 END DO
- 1370 CONTINUE
- !
- ! ----WRITE FINAL OUTPUTS USED.
- !
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement