next up previous contents index
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