source: LMDZ6/trunk/libf/phylmd/geo2atm.F90 @ 3982

Last change on this file since 3982 was 3465, checked in by Laurent Fairhead, 6 years ago

Further modifications for DYNAMICO/LMDZ convergence. These are based
on Yann's LMDZ6_V2 sources. Compiles on irene and converges with revision 3459
in a bucket configuration
YM/LF

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 2.2 KB
RevLine 
[1072]1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/geo2atm.F90,v 1.1 2008-12-05 17:56:40 lsce Exp $
3!
4SUBROUTINE geo2atm(im, jm, px, py, pz, plon, plat, pu, pv, pr)
5  USE dimphy
6  USE mod_phys_lmdz_para
[3465]7  USE mod_grid_phy_lmdz, only: grid_type, unstructured, regular_lonlat
[1072]8  IMPLICIT NONE
[1146]9  INCLUDE 'YOMCST.h'
[3465]10  CHARACTER (len = 6)                :: clmodnam
11  CHARACTER (len = 20)               :: modname = 'geo2atm'
12  CHARACTER (len = 80)               :: abort_message
[1072]13
[1146]14! Change wind coordinates from cartesian geocentric to local spherical
[1072]15! NB! Fonctionne probablement uniquement en MPI seul (sans OpenMP)
16!
[3465]17! Geocentric :
18  ! axe x is eastward : crosses (0N,90E) point.
19  ! axe y  crosses (0N,180E) point.
20  ! axe z is 'up' : crosses north pole.
21!
22! NB! Aux poles, fonctionne probablement uniquement en MPI seul (sans OpenMP)
23
[1072]24  INTEGER, INTENT (IN)                 :: im, jm
25  REAL, DIMENSION (im,jm), INTENT(IN)  :: px, py, pz
26  REAL, DIMENSION (im,jm), INTENT(IN)  :: plon, plat
27  REAL, DIMENSION (im,jm), INTENT(OUT) :: pu, pv, pr
28
[3465]29  REAL :: rad,reps
[1146]30
31
32  rad = rpi / 180.0E0
[3465]33  reps = 1.0e-5
34
[1146]35  pu(:,:) = &
36       - px(:,:) * SIN(rad * plon(:,:)) &
37       + py(:,:) * COS(rad * plon(:,:))
[1072]38
[1146]39  pv(:,:) = &
40       - px(:,:) * SIN(rad * plat(:,:)) * COS(rad * plon(:,:)) &
41       - py(:,:) * SIN(rad * plat(:,:)) * SIN(rad * plon(:,:)) &
42       + pz(:,:) * COS(rad * plat(:,:)) 
[1072]43
[1146]44  pr(:,:) = &
45       + px(:,:) * COS(rad * plat(:,:)) * COS(rad * plon(:,:)) &
46       + py(:,:) * COS(rad * plat(:,:)) * SIN(rad * plon(:,:)) &
47       + pz(:,:) * SIN(rad * plat(:,:))
[1072]48
[3465]49  IF (grid_type==regular_lonlat) THEN
50    ! Value at North Pole
51    IF (is_north_pole_dyn) THEN
52       pu(:, 1) = -px (1,1)
53       pv(:, 1) = -py (1,1)
54       pr(:, 1) = 0.0
55    ENDIF
[1146]56 
[3465]57    ! Value at South Pole     
58    IF (is_south_pole_dyn) THEN
59       pu(:,jm) = -px (1,jm)
60       pv(:,jm) = -py (1,jm)
61       pr(:,jm) = 0.0
62    ENDIF
63
64  ELSE IF (grid_type==unstructured) THEN
65     ! Pole nord pour Dynamico
66     WHERE ( plat(:,:) >= 90.0-reps )
67        pu(:,:) =  py(:,:)
68        pv(:,:) = -px(:,:)
69        pr(:,:) = 0.0e0
70     END WHERE
71
72  ELSE
73     abort_message='Problem: unknown grid type'
74     CALL abort_physic(modname,abort_message,1)
75  END IF
76
[1146]77 
[3465]78 
79 
[1072]80END SUBROUTINE geo2atm
Note: See TracBrowser for help on using the repository browser.