angles.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\tomawac\angles.f
00002 !
00076                      SUBROUTINE ANGLES
00077 !                    *****************
00078 !
00079      &( XLAMD , DTPLUS, DTMOIN)
00080 !
00081 !***********************************************************************
00082 ! TOMAWAC   V6P1                                   08/06/2011
00083 !***********************************************************************
00084 !
00085 !
00086 !
00087 !reference  HASSELMANN S., HASSELMANN K. ET AL.(1985) :
00088 !+                     "COMPUTATIONS AND PARAMETERIZATIONS OF THE NONLINEAR
00089 !+                      ENERGY TRANSFER IN GRAVITY-WAVE SPECTRUM. PART2 :
00090 !+                      PARAMETERIZATIONS OF THE NONLINEAR ENERGY TRANSFER
00091 !+                      FOR APPLICATION IN WAVE MODELS". JPO, VOL 15, PP 1378-1391.
00092 !
00093 !
00094 !
00095 !
00096 !
00097 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00098 !| DTMOIN         |<--| ANGLE ASSOCIATED TO THE FREQUENCY (1-XLAMD).FREQ
00099 !| DTPLUS         |<--| ANGLE ASSOCIATED TO THE FREQUENCY (1+XLAMD).FREQ
00100 !| XLAMD          |-->| DIA STANDARD CONFIGURATION LAMBDA COEFFICIENT
00101 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00102 !
00103       USE DECLARATIONS_TOMAWAC, ONLY : GRADEG
00104 !
00105       IMPLICIT NONE
00106 !
00107       INTEGER LNG,LU
00108       COMMON/INFO/ LNG,LU
00109 !
00110 !.....VARIABLES IN ARGUMENT
00111 !     """""""""""""""""""""
00112       DOUBLE PRECISION XLAMD , DTPLUS, DTMOIN
00113 !
00114 !.....LOCAL VARIABLES
00115 !     """"""""""""""""""
00116       DOUBLE PRECISION AUX
00117 !
00118 !
00119 !.....1. CHECKS THAT LAMBDA RANGES BETWEEN 0 AND 0.5.
00120 !     """"""""""""""""""""""""""""""""""""""""""""""""""""""
00121 !
00122       IF(XLAMD.LT.0.D0.OR.XLAMD.GT.0.5D0) THEN
00123         IF(LNG.EQ.1) THEN
00124           WRITE(LU,1001) XLAMD
00125         ELSE
00126           WRITE(LU,1002) XLAMD
00127         ENDIF
00128         CALL PLANTE(1)
00129         STOP
00130       ENDIF
00131 !
00132 !.....2. COMPUTES TETA_PLUS (DTPLUS) AND TETA_MOINS (DTMOIN).
00133 !     """""""""""""""""""""""""""""""""""""""""""""""""""""""
00134       AUX=2.D0*XLAMD*(1.D0+XLAMD*XLAMD)
00135       DTPLUS=ACOS( (1.D0+AUX)/(1.D0+XLAMD)**2 )*GRADEG
00136       DTMOIN=ACOS( (1.D0-AUX)/(1.D0-XLAMD)**2 )*GRADEG
00137 !
00138 !
00139  1001 FORMAT('/!/-----------------------------------------------/!/'/
00140      &       '/!/  ARRET DU PROGRAMME DANS SUBROUTINE ANGLES    /!/'/
00141      &       '/!/-----------------------------------------------/!/'/
00142      &       '/!/  LA VALEUR CHOISIE POUR LE PARAMETRE LAMBDA   /!/'/
00143      &       '/!/  INTERVENANT DANS LA METHODE DISCRETE INTER-  /!/'/
00144      &       '/!/  ACTION APPROXIMATION DOIT ETRE INCLUSE DANS  /!/'/
00145      &       '/!/  L INTERVALLE [ 0. ; 0.5 ].                   /!/'/
00146      &       '/!/  OR LA VALEUR UTILISEE EST : ', G10.4,'       /!/'/
00147      &       '/!/-----------------------------------------------/!/')
00148 !
00149  1002 FORMAT('/!/-----------------------------------------------/!/'/
00150      &       '/!/        PROGRAM STOP IN SUBROUTINE ANGLES      /!/'/
00151      &       '/!/-----------------------------------------------/!/'/
00152      &       '/!/       THE VALUE OF THE LAMBDA PARAMETER       /!/'/
00153      &       '/!/             USED IN THE DIA METHOD            /!/'/
00154      &       '/!/  MUST BE INCLUDED IN THE INTERVAL [ 0. ; 0.5 ]/!/'/
00155      &       '/!/  THE VALUE HERE IMPOSED IS : ', G10.4,'       /!/'/
00156      &       '/!/-----------------------------------------------/!/')
00157 !
00158       RETURN
00159       END

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