atm2geo.F90

Go to the documentation of this file.

!
! $Id: atm2geo.F90 2346 2015-08-21 15:13:46Z emillour $
!
SUBROUTINE atm2geo ( im, jm, pte, ptn, plon, plat, pxx, pyy, pzz )
  USE dimphy
  USE mod_phys_lmdz_para
  IMPLICIT NONE
  include 'YOMCST.h'
!
! Change wind local atmospheric coordinates to geocentric
!
  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
  
  REAL :: rad


  rad = rpi / 180.0e0
  
  pxx(:,:) = & 
       - pte(:,:) * sin(rad * plon(:,:)) &
       - ptn(:,:) * sin(rad * plat(:,:)) * cos(rad * plon(:,:))

  pyy(:,:) = &
       + pte(:,:) * cos(rad * plon(:,:)) &
       - ptn(:,:) * sin(rad * plat(:,:)) * sin(rad * plon(:,:))
  
  pzz(:,:) = &
       + ptn(:,:) * cos(rad * plat(:,:))
  
! Value at North Pole  
  IF (is_north_pole) THEN
     pxx(:, 1) = - pte(1, 1)
     pyy(:, 1) = - ptn(1, 1) 
     pzz(:, 1) = pzz(1,1)
  ENDIF

! Value at South Pole
  IF (is_south_pole) THEN
     pxx(:,jm) = pxx(1,jm)
     pyy(:,jm) = pyy(1,jm)
     pzz(:,jm) = pzz(1,jm)
  ENDIF
  
END SUBROUTINE atm2geo