next up previous contents
Next: celledge_global Up: Fortran: Module Interface areaMod Previous: cellarea_global   Contents

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 points
     integer , intent(in) :: nlon              !dimension: number of longitude points
     integer , intent(in) :: numlon(nlat)      !number of grid cells per latitude strip
     real(r8), intent(in) :: longxy(nlon,nlat) !longitude at center of grid cell
     real(r8), intent(in) :: latixy(nlon,nlat) !latitude at center of grid cell
     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 edge (degrees)
     real(r8), intent(out):: lonw(nlon+1,nlat) !grid cell longitude, western edge (degrees)
REVISION HISTORY:
   Created by Mariana Vertenstein



Mariana Vertenstein 2003-01-14