intang.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\tomawac\intang.f
00002 !
00060                      SUBROUTINE INTANG
00061 !                    *****************
00062 !
00063      &( LAVANT, LAPRES, IPLAN , NPLAN , DELTAD)
00064 !
00065 !***********************************************************************
00066 ! TOMAWAC   V6P1                                   17/06/2011
00067 !***********************************************************************
00068 !
00069 !
00070 !
00071 !reference  HASSELMANN S., HASSELMANN K. ET AL.(1985) :
00072 !+                     "COMPUTATIONS AND PARAMETERIZATIONS OF THE NONLINEAR
00073 !+                      ENERGY TRANSFER IN GRAVITY-WAVE SPECTRUM. PART2 :
00074 !+                      PARAMETERIZATIONS OF THE NONLINEAR ENERGY TRANSFER
00075 !+                      FOR APPLICATION IN WAVE MODELS". JPO, VOL 15, PP 1378-1391.
00076 !
00077 !
00078 !
00079 !
00080 !
00081 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00082 !| DELTAD         |-->| ANGULAR DEVIATION FROM THE STARTING DIRECTION
00083 !| IPLAN          |-->| STARTING DIRECTION INDEX
00084 !| LAPRES         |<--| ANGULAR INDEX FOLLOWING THE DIRECTION
00085 !| LAVANT         |<--| ANGULAR INDEX PRECEDING THE DIRECTION
00086 !| NPLAN          |-->| NUMBER OF DIRECTIONS
00087 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00088 !
00089       IMPLICIT NONE
00090 !
00091 !.....VARIABLES IN ARGUMENT
00092 !     """""""""""""""""""""
00093       INTEGER  LAVANT, LAPRES, NPLAN , IPLAN
00094       DOUBLE PRECISION DELTAD
00095 !
00096 !.....LOCAL VARIABLES
00097 !     """"""""""""""""""
00098       DOUBLE PRECISION TETA  , DTETAD
00099 !
00100 !
00101       DTETAD=360.D0/DBLE(NPLAN)
00102       TETA=DBLE(IPLAN-1)*DTETAD+DELTAD
00103 !
00104 !.....TETA IS ADJUSTED TO BE COMPRISED BETWEEN 0 AND 360 DEG.
00105 !     """""""""""""""""""""""""""""""""""""""""""""""""""""""""""
00106   100 IF (TETA.GE.360.D0) THEN
00107         TETA=TETA-360.D0
00108         GOTO 100
00109       ENDIF
00110   110 IF (TETA.LT.0.D0) THEN
00111         TETA=TETA+360.D0
00112         GOTO 110
00113       ENDIF
00114 !
00115 !.....GETS THE ANGULAR INDICES PRECEDING AND FOLLOWING TETA
00116 !     """"""""""""""""""""""""""""""""""""""""""""""""""""""""
00117       LAVANT=INT(TETA/DTETAD)+1
00118       LAPRES=LAVANT+1
00119       IF (LAPRES.GT.NPLAN) LAPRES=1
00120 !
00121       RETURN
00122       END

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