latitu.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\latitu.f
00002 !
00052                      SUBROUTINE LATITU
00053 !                    *****************
00054 !
00055      &(COSLAT,SINLAT,LAMBD0,  Y,NPOIN)
00056 !
00057 !***********************************************************************
00058 ! BIEF   V6P3                                   21/08/2010
00059 !***********************************************************************
00060 !
00061 !
00062 !
00063 !
00064 !
00065 !
00066 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00067 !| COSLAT         |<--| COSINUS OF LAMBDA0
00068 !| LAMBD0         |-->| LATITUDE OF ORIGIN POINT
00069 !| NPOIN          |-->| NUMBER OF POINTS
00070 !| SINLAT         |<--| SINUS OF LAMBDA0
00071 !| Y              |-->| ORDINATES OF POINTS IN THE MESH
00072 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00073 !
00074       IMPLICIT NONE
00075       INTEGER LNG,LU
00076       COMMON/INFO/LNG,LU
00077 !
00078 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00079 !
00080       INTEGER         , INTENT(IN)    :: NPOIN
00081       DOUBLE PRECISION, INTENT(IN)    :: Y(NPOIN)
00082       DOUBLE PRECISION, INTENT(INOUT) :: COSLAT(NPOIN),SINLAT(NPOIN)
00083       DOUBLE PRECISION, INTENT(IN)    :: LAMBD0
00084 !
00085 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00086 !
00087       INTEGER I
00088 !
00089       DOUBLE PRECISION LB2RAD,SURR,PISUR4,PISUR2,XLAMB
00090 !
00091       INTRINSIC TAN,ATAN,SIN,COS,EXP
00092 !
00093 !-----------------------------------------------------------------------
00094 !
00095 ! EARTH RADIUS
00096 !
00097       SURR = 1.D0 / 6370000.D0
00098 !
00099 !-----------------------------------------------------------------------
00100 !
00101       PISUR4 = ATAN(1.D0)
00102       PISUR2 = PISUR4 + PISUR4
00103 !
00104 !     LAMBD0/2 IN RADIANS
00105 !
00106       LB2RAD = LAMBD0 * PISUR4 / 90.D0
00107 !
00108 !     1/COS(LAMBDA),COS(LAMBDA),SIN(LAMBDA)
00109 !
00110       DO I = 1 , NPOIN
00111 !
00112         XLAMB = 2.D0* ATAN(EXP(Y(I)*SURR)*TAN(LB2RAD+PISUR4))-PISUR2
00113         COSLAT(I) = COS(XLAMB)
00114         SINLAT(I) = SIN(XLAMB)
00115 !
00116       ENDDO
00117 !
00118 !-----------------------------------------------------------------------
00119 !
00120       RETURN
00121       END

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