Esope source files

   1 : C POSPRE    SOURCE    CHAT      05/01/13    02:18:16     5004           
   2 :         SUBROUTINE POSPRE(MTRAV,NSPT,NAMRT,IG,ID,TE,NT)
   3 : C       ===============================================
   4 : C=======================================================================
   5 : C       COMPUTATION OF THE RESPONSE SPECTRA FROM A  POWER 
   6 : C                                   SPECTRA
   7 : C
   8 : C       MTRAV=F(NSPT),S(NSPT) ,ETI(NAMRT), T(NI), RES(NAMRT,NI)
   9 : C
  10 : C       F(1:NSPT)  
  11 : C       S(1:NSPT)    =FREQUENCY  AND P. SPC. VALUES
  12 : C       NSPT         =NUMBER OF P. SPC. VALUES
  13 : C
  14 : C       ETI(1:NAMRT) =DAMPING    VALUES
  15 : C       NAMRT        =NUMBER OF DAMPING VALUES
  16 : C
  17 : C       IG           =TYPE OF THE RESPONSE (1 ACCELERATION,2 VELOCITY,
  18 : C                                           3 DISPLACEMENT)
  19 : C       ID           =TYPE OF THE MAXIMUM DISTRIBUTION (1 CRAMER,
  20 : C                                                       2 NEWGUMB)
  21 : C
  22 : C       TE           =TIME DURATION
  23 : C
  24 : C       T(1:NT)      =PERIODS AVAILABLE (INPUT)
  25 : C       T(1:NT)      =PERIODS (OUTPUT)
  26 : C       NT           =NUMBER OF PERIOD AVAILABLE (INPUT)
  27 : C       NT           =NUMBER OF PERIOD (OUTPUT)
  28 : C       -->IF NT=0 T SHOULD BE COMPUTED AND NT<201
  29 : C       -->IF NT>0 T IS USED
  30 : C
  31 : C       RES(1:NAMRT,1:NT)
  32 : C                    = RESPONSE SPECTRA FOR VARIOUS DAMPING
  33 : C                      VALUES
  34 : C
  35 : C       FROM A.P. CUSTUMIZED BY P.P. 27/2/90 + 2/4/90 + 15/9/91
  36 : C       ============================
  37 : C=======================================================================
  38 :       IMPLICIT INTEGER(I-N)
  39 :         IMPLICIT REAL*8(A-H,O-Z)
  40 : C======+++++++++++ bon fontionnement erreur PSRS/PRNS +++++++++++ ======
  41 : -INC CCOPTIO
  42 :         SEGMENT,MTRAV
  43 :           IMPLIED F(NSPT),S(NSPT) ,ETI(NAMRT), T(NI), RES(NAMRT,NI)
  44 :         ENDSEGMENT
  45 :         SEGMENT,MMTRA
  46 :           IMPLIED XLTIME(NNT,3)
  47 :         ENDSEGMENT
  48 :         NNT=1
  49 :         SEGINI,MMTRA
  50 : C======+++++++++++ bon fontionnement erreur PSRS/PRNS +++++++++++ ======
  51 : C
  52 :         DPI=8*ATAN(1.D0)
  53 : C
  54 : C       SELECTION OF PERIODS
  55 : C       (IF NT >0 THE PERIODS ARE SELECTED BY THE USER)
  56 : C       ===================
  57 :             IF(NT.EQ.0)THEN
  58 :         NT=75
  59 :         TINF=.04D0
  60 :         UNPXI=EXP((LOG(TE)-LOG(TINF))/(NT-1))
  61 :         T(1)=TINF
  62 :         DO 300 I=2,NT
  63 : 300       T(I)=T(I-1)*UNPXI
  64 :         T(NT)=TE
  65 :             ENDIF
  66 : C       COMPUTING OF RESPONSE SPECTRA
  67 : C       =============================
  68 :         DO 350 J=1,NAMRT
  69 :         ETA=ETI(J)
  70 :         DO 351 I=1,NT
  71 :         FRQ=1/T(I)
  72 : C       -MOMENT
  73 : C        ------
  73 : C       CALL MOMENT(FRQ,ETA,TE,NSPT,F,S, IG, XL0,XL1,XL2)
  73 :         CALL MOMENT(FRQ,ETA,TE,NSPT,F,S, IG, XLTIME(1,1),
  74 :      >                           XLTIME(1,2),XLTIME(1,3))
  75 : C       -DISTRIBUTION
  76 : C        ------------
  76 : C       CALL DISTRI(XL0,XL1,XL2,TE, ID,DPI, VALMAX)
  76 :         CALL DISTRI(MMTRA,TE, ID,DPI, VALMAX)
  77 :         RES(J,I)=VALMAX
  78 : 351     CONTINUE
  79 : 350     CONTINUE
  80 :         SEGSUP,MMTRA
  81 :         RETURN
  82 :         END
  83 :