Advertisement
hamzaalloush

BODY_NAMELIST

Feb 17th, 2016
99
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. !
  2. !  CHECK BODY INPUTS TO MEET SPECIFICATIONS NOTED ABOVE
  3. !  IN READING BODY NAMELIST DATA
  4. !
  5.       NNX=NX+0.5
  6.       IF(ELLIP.EQ.0.)ELLIP=UNUSED
  7.       IF(ELLIP.EQ.UNUSED) ELLIP=1.0
  8.       IF(R(1).NE.UNUSED .AND. S(1).EQ.UNUSED .AND. P(1).EQ.UNUSED)      &
  9.      &   GO TO 1310
  10.       IF(S(1).NE.UNUSED .AND. R(1).EQ.UNUSED .AND. P(1).EQ.UNUSED)      &
  11.      &   GO TO 1330
  12.       IF(R(1).NE.UNUSED .AND. S(1).NE.UNUSED .AND. P(1).NE.UNUSED)      &
  13.      &   GO TO 1350
  14. !
  15. !  HERE FOR RADIUS OR AREA INPUT WITH ANOTHER VARIABLE.
  16. !  ONLY R OR S OR R,S,P CAN BE USED.  SET THE GONOGO FLAG,
  17. !  WRITE A MESSAGE AND CONTINUE.
  18. !
  19.       GONOGO=.FALSE.
  20.       WRITE(6,1510)
  21.       GO TO 1370
  22. !
  23. !     HALF WIDTH, R, SPECIFIED
  24. !
  25.  1310 DO 1320 K=1,NNX
  26.       S(K)=PI*R(K)**2*ELLIP
  27.       P(K)=2*PI*R(K)*SQRT((ELLIP**2+1.)/2.)
  28.       R(K)=R(K)*SQRT(ELLIP)
  29.  1320 END DO
  30.       GO TO 1370
  31. !
  32. !     CROSS SECTIONAL AREA, S, SPECIFIED
  33. !
  34.  1330 DO 1340 K=1,NNX
  35.       R(K)=SQRT(S(K)/PI)
  36.       P(K)=2.*PI*SQRT((S(K)/ELLIP)*(ELLIP**2+1.)/2./PI)
  37.  1340 END DO
  38.       GO TO 1370
  39. !
  40. !     R, S, AND P ALL SPECIFIED
  41. !
  42.  1350 DO 1360 K=1,NNX
  43.       R(K)=R(K)*SQRT(ELLIP)
  44.  1360 END DO
  45.  1370 CONTINUE
  46. !
  47. !     ----WRITE FINAL OUTPUTS USED.
  48. !
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement