source: LMDZ5/trunk/libf/phy_common/regular_lonlat_mod.F90 @ 2489

Last change on this file since 2489 was 2346, checked in by Ehouarn Millour, 9 years ago

Physics/dynamics separation:

  • remove all references to dimensions.h from physics. nbp_lon (==iim) , nbp_lat (==jjm+1) and nbp_lev (==llm) from mod_grid_phy_lmdz should be used instead.
  • added module regular_lonlat_mod in phy_common to store information about the global (lon-lat) grid cell boundaries and centers.

EM

File size: 1.9 KB
Line 
1MODULE regular_lonlat_mod
2
3  ! Store information on the global physics grid
4  ! for a regular (longitude-latitude) grid
5 
6  INTEGER, PARAMETER :: north_east=1     ! boundaries of regular lontlat
7  INTEGER, PARAMETER :: north_west=2     ! boundaries of regular lontlat
8  INTEGER, PARAMETER :: south_west=3     ! boundaries of regular lontlat
9  INTEGER, PARAMETER :: south_east=4     ! boundaries of regular lontlat
10  INTEGER, PARAMETER :: east=1           ! boundaries of regular lontlat
11  INTEGER, PARAMETER :: west=2           ! boundaries of regular lontlat
12  INTEGER, PARAMETER :: north=1          ! boundaries of regular lontlat
13  INTEGER, PARAMETER :: south=2          ! boundaries of regular lontlat
14
15! Global definition, shared by all threads
16! Do not set threadprivate directives
17
18  REAL,SAVE,ALLOCATABLE :: lon_reg(:)      ! value of longitude cell (rad)
19
20  REAL,SAVE,ALLOCATABLE :: lat_reg(:)      ! value of longitude cell (rad)
21
22  REAL,SAVE,ALLOCATABLE :: boundslon_reg(:,:)      ! value of boundaries cell (1=>east, 2=>west)(rad)
23
24  REAL,SAVE,ALLOCATABLE :: boundslat_reg(:,:)      ! value of longitude cell (1=>north, 2=>south)(rad)
25
26
27CONTAINS
28 
29  SUBROUTINE init_regular_lonlat(nbp_lon, nbp_lat, lon_reg_, lat_reg_, boundslon_reg_, boundslat_reg_)
30    IMPLICIT NONE
31    INTEGER,INTENT(IN) :: nbp_lon
32    INTEGER,INTENT(IN) :: nbp_lat
33    REAL,   INTENT(IN) :: lon_reg_(nbp_lon) 
34    REAL,   INTENT(IN) :: lat_reg_(nbp_lat)
35    REAL,   INTENT(IN) :: boundslon_reg_(nbp_lon,2)
36    REAL,   INTENT(IN) :: boundslat_reg_(nbp_lat,2) 
37
38   
39    ALLOCATE(lon_reg(nbp_lon))
40    lon_reg(:)=lon_reg_(:)
41
42    ALLOCATE(lat_reg(nbp_lat))
43    lat_reg(:)=lat_reg_(:)
44
45    ALLOCATE(boundslon_reg(nbp_lon,2))
46    boundslon_reg(:,:)=boundslon_reg_(:,:)
47
48    ALLOCATE(boundslat_reg(nbp_lat,2)) 
49    boundslat_reg(:,:)=boundslat_reg_(:,:)
50
51 
52  END SUBROUTINE init_regular_lonlat   
53
54END MODULE regular_lonlat_mod
55
Note: See TracBrowser for help on using the repository browser.