Next: celledge_global Up: Module areaMod (File: areaMod.F90) Previous: cellarea_global   Contents   Index

## celledge_regional

INTERFACE:

```   subroutine celledge_regional (nlat   , nlon  , numlon , longxy ,  &
latixy , edgen , edgee  , edges  ,  &
edgew  , lats  , lonw   )
```
DESCRIPTION:

```   Southern and western edges of grid cells - regional grid
(can become global as special case)
Latitudes -- southern/northern edges for each latitude strip.
For grids oriented South to North, the southern
and northern edges of latitude strip [j] are:
southern = lats(j  )
northern = lats(j+1)
For grids oriented North to South: the southern
and northern edges of latitude strip [j] are:
northern = lats(j  )
southern = lats(j+1)
In both cases, [lats] must be dimensioned lats(lat+1)
Longitudes -- western edges. Longitudes for the western edge of the
cells must increase continuously and span 360 degrees. Assume that
grid starts at Dateline with western edge on Dateline Western edges
correspond to [longxy] (longitude at center of cell) and range from
-180 to 180 with negative longitudes west of Greenwich.
Partial grids that do not span 360 degrees are allowed so long as they
have the convention of Grid 1 with
western edge of grid: >= -180 and < 180
eastern edge of grid: > western edge  and <= 180
[lonw] must be dimensioned lonw(lon+1,lat) because each latitude
strip can have variable longitudinal resolution
```
ARGUMENTS:
```     implicit none
integer , intent(in) :: nlat              !dimension: number of latitude
integer , intent(in) :: nlon              !dimension: number of longitude
integer , intent(in) :: numlon(nlat)      !number of grid cells per latit
real(r8), intent(in) :: longxy(nlon,nlat) !longitude at center of grid ce
real(r8), intent(in) :: latixy(nlon,nlat) !latitude at center of grid cel
real(r8), intent(in) :: edgen             !northern edge of grid (degrees
real(r8), intent(in) :: edgee             !eastern edge of grid (degrees)
real(r8), intent(in) :: edges             !southern edge of grid (degrees
real(r8), intent(in) :: edgew             !western edge of grid (degrees)
real(r8), intent(out):: lats(nlat+1)      !grid cell latitude, southern e
real(r8), intent(out):: lonw(nlon+1,nlat) !grid cell longitude, western e
```
REVISION HISTORY:
```   Created by Mariana Vertenstein
```

Mariana Vertenstein 2004-06-21