clm_varmap.f90

来自「CCSM Research Tools: Community Atmospher」· F90 代码 · 共 91 行

#include <misc.h>
#include <preproc.h>

module clm_varmap

!----------------------------------------------------------------------- 
! 
! Purpose: 
! Module of mapping arrays
! 
! Method: 
! The land surface model works by gathering all the land points on a
! [lon]x[lat] grid into a vector of [numland] land points. This
! is then expanded into a vector of [numpatch] subgrid patches, allowing
! for up to [maxpatch] subgrid patches per land point. 
! [ixy], [jxy], [patch], and [land] are indices for the mapping: 
! [lsmlon]x[lsmlat] grid <-> [numland] vector of land points <-> [numpatch] vector 
! of subgrid points. [landvec%wtxy] are the weights to obtain the 
! grid average from the subgrid patches.
! 
! Author: Mariana Vertenstein
! 
!-----------------------------------------------------------------------
! $Id: clm_varmap.F90,v 1.1.22.3 2002/04/27 15:38:46 erik Exp $
!-----------------------------------------------------------------------

  use precision
  implicit none

  integer :: numland                 !number of land points
  integer :: begland                 !beginning land index (minimum value is 1)
  integer :: endland                 !ending land index (maximum value is  numpatch)
  integer :: numpatch                !total number of patches allowing for subgrid patches
  integer :: begpatch                !beginning patch index (minimum value is 1)
  integer :: endpatch                !ending patch index (maximum value is  numpatch)

  type land1d
     integer ,pointer :: patch(:,:)  !patch vector index: 1 to numpatch
     integer, pointer :: ixy(:)      !longitude index for each land point: 1 to lsmlon
     integer, pointer :: jxy(:)      !latitude index for each land point: 1 to lsmlat
     real(r8),pointer :: wtxy(:,:)   !subgrid weights 
  end type land1d

  type subgrid1d
     integer , pointer :: ixy(:)     !longitude index for each patch point (1 to lsmlon)
     integer , pointer :: jxy(:)     !latitude index for each patch point (1 to lsmlat)
     integer , pointer :: mxy(:)     !subgrid patch type index (1 to maxpatch)
     real(r8), pointer :: wtxy(:)    !subgrid weight for each patch point 
     integer , pointer :: land(:)    !land index for each patch point
  end type subgrid1d

  type (land1d)    :: landvec
  type (subgrid1d) :: patchvec

  SAVE

!=======================================================================
CONTAINS
!=======================================================================

  subroutine mapvar_ini

    use precision
    use infnan
    use clm_varpar, only : maxpatch
    implicit none

! Initializes mapping vectors

    allocate (landvec%ixy(numland))
    allocate (landvec%jxy(numland))
    allocate (landvec%wtxy(numland,maxpatch))
    allocate (landvec%patch(numland,maxpatch))

    allocate (patchvec%ixy(numpatch))
    allocate (patchvec%jxy(numpatch))
    allocate (patchvec%mxy(numpatch))
    allocate (patchvec%wtxy(numpatch))
    allocate (patchvec%land(numpatch))

! Initialize appropriate components as infinity

    landvec%wtxy(:,:) = inf
    patchvec%wtxy(:)  = inf

  end subroutine mapvar_ini

!=======================================================================

end module clm_varmap