prism_reducedgrid_map

prism_reducedgrid_map (grid_id, nbr_latitudes, nbr_points_per_lat, ierror)

Table 5.14: prism_reducedgrid_map arguments

For Gaussian reduced grids only - mandatory. All processes that announce a Gaussian reduced grid have to call prism_reducedgrid_map for a description of the global reduced Gaussian grid, providing the same information about the global grid. As the coordinates of a Gaussian reduced grid are expressed in a 1d arrays, this additional information is needed to speed up the reconstruction of a 2d view of the 1d arrays.

Example:

  integer, parameter         :: gnbr_lats   =     96
  integer                    :: gnbr_lons(gnbr_lats)
  integer, parameter         :: ndim = 48
  integer                    :: array(ndim)
  integer, intent(out)       :: nbr_lons(2*ndim)
  integer                    :: ierror
  data array / &
    20,    25,    36,    40,    45,    50,    60,    60,    72,    75,   &
    80,    90,    96,   100,   108,   120,   120,   120,   128,   135,   &
   144,   144,   160,   160,   160,   160,   160,   180,   180,   180,   &
   180,   180,   192,   192,   192,   192,   192,   192,   192,   192,   &
   192,   192,   192,   192,   192,   192,   192,   192/

  do i = 1, ndim
    gnbr_lons(i) = array(i)
    gnbr_lons(2*ndim+1-i) = array(i)
  enddo

  call prism_reducedgrid_map ( grid_id(1), gnbr_lats, gnbr_lons, ierror )
  if (ierror /= 0) n_errors = n_errors + 1