frepic.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\tomawac\frepic.f
00002 !
00049                      SUBROUTINE FREPIC
00050 !                    *****************
00051 !
00052      &( FPIC  , F     , FREQ  , NF    , NPLAN , NPOIN2, EMAX  , E     )
00053 !
00054 !***********************************************************************
00055 ! TOMAWAC   V6P1                                   15/06/2011
00056 !***********************************************************************
00057 !
00058 !
00059 !
00060 !
00061 !
00062 !
00063 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00064 !| E              |<->| WORK TABLE
00065 !| EMAX           |<->| WORK TABLE
00066 !| F              |-->| VARIANCE DENSITY DIRECTIONAL SPECTRUM
00067 !| FPIC           |<--| PEAK FREQUENCIES
00068 !| FREQ           |-->| DISCRETIZED FREQUENCIES
00069 !| NF             |-->| NUMBER OF FREQUENCIES
00070 !| NPLAN          |-->| NUMBER OF DIRECTIONS
00071 !| NPOIN2         |-->| NUMBER OF POINTS IN 2D MESH
00072 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00073 !
00074       IMPLICIT NONE
00075 !
00076 !.....VARIABLES IN ARGUMENT
00077 !     """"""""""""""""""""
00078       INTEGER  NF    , NPLAN , NPOIN2
00079       DOUBLE PRECISION F(NPOIN2,NPLAN,NF), FREQ(NF)  , FPIC(NPOIN2)
00080       DOUBLE PRECISION EMAX(NPOIN2),E(NPOIN2)
00081 !
00082 !.....LOCAL VARIABLES
00083 !     """""""""""""""""
00084       INTEGER  JP    , JF    , IP
00085 !
00086 !
00087       DO IP = 1,NPOIN2
00088         FPIC(IP) = 1.D-20
00089         EMAX(IP) = 0.D0
00090       ENDDO ! IP
00091 !
00092 !.....LOOP OVER DISCRETISED FREQUENCIES
00093 !     """""""""""""""""""""""""""""""""""""""""""""
00094       DO JF = 1,NF
00095 !
00096 !.......INTEGRATES WRT DIRECTIONS TO GET E(F)
00097 !       """""""""""""""""""""""""""""""""""""""""""""""""
00098         DO IP=1,NPOIN2
00099           E(IP) = 0.D0
00100         ENDDO ! IP
00101         DO JP = 1,NPLAN
00102           DO IP=1,NPOIN2
00103                  E(IP) = E(IP) + F(IP,JP,JF)
00104           ENDDO ! IP
00105         ENDDO ! JP
00106 !
00107 !.......KEEPS THE MAXIMUM VALUE FOR E(F) AND ASSOCIATED FREQUENCY
00108 !       """""""""""""""""""""""""""""""""""""""""""""""""""""
00109         DO IP=1,NPOIN2
00110           IF (E(IP).GT.EMAX(IP)) THEN
00111             EMAX(IP) = E(IP)
00112             FPIC(IP) = FREQ(JF)
00113           ENDIF
00114         ENDDO ! IP
00115 !
00116       ENDDO ! JF
00117 !
00118       RETURN
00119       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0