libf/phylmd/atm2geo.f

!
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/atm2geo.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
!
C
      SUBROUTINE atm2geo ( im, jm, pte, ptn, plon, plat, pxx, pyy, pzz )
cc
cc Change wind local atmospheric coordinates to
cc geocentric
cc
c$$$      INCLUDE 'param.h'
c
      INTEGER, INTENT (in)              :: im, jm
      REAL, DIMENSION (im,jm), INTENT (in) :: pte, ptn
      REAL, DIMENSION (im,jm), INTENT (in) :: plon, plat
      REAL, DIMENSION (im,jm), INTENT(out) :: pxx, pyy, pzz
c
      REAL, PARAMETER :: rpi = 3.141592653E0
      REAL, PARAMETER :: rad = rpi / 180.0E0
c
      REAL, DIMENSION (im,jm) :: zsinlon, zcoslon
      REAL, DIMENSION (im,jm) :: zsinlat, zcoslat
c
      LOGICAL, SAVE :: linit = .FALSE.
c
c$$$      IF ( .NOT. linit ) THEN 
          zsinlon = SIN (rad * plon)
          zcoslon = COS (rad * plon)
          zsinlat = SIN (rad * plat)
          zcoslat = COS (rad * plat)
          linit = .TRUE.
c$$$      ENDIF 
c
      pxx = - zsinlon * pte - zsinlat * zcoslon * ptn
      pyy =   zcoslon * pte - zsinlat * zsinlon * ptn
      pzz =   zcoslat * ptn
c
c Value at North Pole
      pxx ( :,  1) = - ptn ( 1, 1)
      pyy ( :,  1) = - pte ( 1, 1)
      pzz ( :,  1) = 0.0
c Value at South Pole
      pxx ( :, jm) = + ptn ( 1, jm)
      pyy ( :, jm) = + pte ( 1, jm)
      pzz ( :, jm) = 0.0
c
      RETURN 
      END SUBROUTINE atm2geo
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/phylmd/atm2geo.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
3			!
4			C
5			SUBROUTINE atm2geo ( im, jm, pte, ptn, plon, plat, pxx, pyy, pzz )
6			cc
7			cc Change wind local atmospheric coordinates to
8			cc geocentric
9			cc
10			c$$$ INCLUDE 'param.h'
11			c
12			INTEGER, INTENT (in) :: im, jm
13			REAL, DIMENSION (im,jm), INTENT (in) :: pte, ptn
14			REAL, DIMENSION (im,jm), INTENT (in) :: plon, plat
15			REAL, DIMENSION (im,jm), INTENT(out) :: pxx, pyy, pzz
16			c
17			REAL, PARAMETER :: rpi = 3.141592653E0
18			REAL, PARAMETER :: rad = rpi / 180.0E0
19			c
20			REAL, DIMENSION (im,jm) :: zsinlon, zcoslon
21			REAL, DIMENSION (im,jm) :: zsinlat, zcoslat
22			c
23			LOGICAL, SAVE :: linit = .FALSE.
24			c
25			c$$$ IF ( .NOT. linit ) THEN
26			zsinlon = SIN (rad * plon)
27			zcoslon = COS (rad * plon)
28			zsinlat = SIN (rad * plat)
29			zcoslat = COS (rad * plat)
30			linit = .TRUE.
31			c$$$ ENDIF
32			c
33			pxx = - zsinlon * pte - zsinlat * zcoslon * ptn
34			pyy = zcoslon * pte - zsinlat * zsinlon * ptn
35			pzz = zcoslat * ptn
36			c
37			c Value at North Pole
38			pxx ( :, 1) = - ptn ( 1, 1)
39			pyy ( :, 1) = - pte ( 1, 1)
40			pzz ( :, 1) = 0.0
41			c Value at South Pole
42			pxx ( :, jm) = + ptn ( 1, jm)
43			pyy ( :, jm) = + pte ( 1, jm)
44			pzz ( :, jm) = 0.0
45			c
46			RETURN
47			END SUBROUTINE atm2geo