yom_phys_grid.F90

Go to the documentation of this file.

MODULE yom_phys_grid

USE parkind1  ,ONLY : jpim     ,jprb

IMPLICIT NONE
SAVE

PRIVATE
PUBLIC type_phys_point ,phys_grid,dyn_grid,dyn_sl,phys_sl, yphypoi, ydynpoi, &
        & jpmxnei, phys_grid_struct, sl_struct, type_dyn_point

!-------------------------------------------------------------------------
! Derived types for describing the coarse physics grid structure. 
!  The descriptors themselves
! (YGFL and YGFLC) can be found in module yom_ygfl.F90.
!-------------------------------------------------------------------------
! Modifications:

INTEGER(KIND=JPIM), PARAMETER :: jpmxnei=36 ! maximum number of neighbouring
                                            ! points for averaging or
                                            ! interpolation
                                              
TYPE phys_grid_struct
  INTEGER(KIND=JPIM) :: ngptot          ! number of physics points in the task
  INTEGER(KIND=JPIM) :: ngptotg
  INTEGER(KIND=JPIM) :: ndgsal, ndgenl
  INTEGER(KIND=JPIM) :: ndgsag, ndgeng
  INTEGER(KIND=JPIM), POINTER :: nloeng(:)
  REAL(KIND=JPRB), POINTER   :: rmu(:)
  REAL(KIND=JPRB), POINTER   :: rw(:)

  INTEGER(KIND=JPIM)   :: nresol_id
  INTEGER(KIND=JPIM)   :: ngptotmx, nspec2, nsmax
  INTEGER(KIND=JPIM)   :: nptrfloff, nump, ndlon
  INTEGER(KIND=JPIM)   ::  ndgsah, ndgenh
  INTEGER(KIND=JPIM)   ::  ndglg, ndlsur
  INTEGER(KIND=JPIM)   ::  nfrstloff, ndsur1,ndgsur
  INTEGER(KIND=JPIM)   ::  myfrstactlat, mylstactlat
  INTEGER(KIND=JPIM)   :: ngpblks
  INTEGER(KIND=JPIM)   :: nproma

  INTEGER(KIND=JPIM), POINTER, DIMENSION(:)     :: nrgri, nptrfrstlat, nfrstlat
  INTEGER(KIND=JPIM), POINTER, DIMENSION(:)     :: nlstlat, myms, nasm0
  INTEGER(KIND=JPIM), POINTER, DIMENSION(:,:)   :: nsta, nonl
  INTEGER(KIND=JPIM), POINTER, DIMENSION(:)     :: nstagp
  REAL(KIND=JPRB), POINTER   :: rsqm2(:), rlatig(:), rlati(:)
  REAL(KIND=JPRB), POINTER   :: ripi0(:), ripi1(:), ripi2(:)
  INTEGER(KIND=JPIM), POINTER, DIMENSION(:)     :: nptrlstlat

END TYPE phys_grid_struct

TYPE(phys_grid_struct) :: phys_grid, dyn_grid 

TYPE sl_struct
  INTEGER(KIND=JPIM),POINTER :: nslsta(:)
  INTEGER(KIND=JPIM),POINTER :: nslonl(:)
  INTEGER(KIND=JPIM),POINTER :: nsloff(:)
  INTEGER(KIND=JPIM),POINTER :: nslext(:,:)
  INTEGER(KIND=JPIM),POINTER :: nslsendpos(:)
  INTEGER(KIND=JPIM),POINTER :: nslrecvpos(:)
  INTEGER(KIND=JPIM),POINTER :: nslsendptr(:)
  INTEGER(KIND=JPIM),POINTER :: nslrecvptr(:)
  INTEGER(KIND=JPIM),POINTER :: nslcore(:)
  INTEGER(KIND=JPIM),POINTER :: nslcomm(:)

  INTEGER(KIND=JPIM) :: naslb1
  INTEGER(KIND=JPIM) :: nslprocs
  INTEGER(KIND=JPIM) :: nslmpbufsz
  INTEGER(KIND=JPIM) :: nslrpt
  INTEGER(KIND=JPIM) :: nslspt
  INTEGER(KIND=JPIM) :: nslwiden
  INTEGER(KIND=JPIM) :: nslwides
  INTEGER(KIND=JPIM) :: nslwidee
  INTEGER(KIND=JPIM) :: nslwidew
END TYPE sl_struct

TYPE(sl_struct) :: dyn_sl
TYPE(sl_struct) :: phys_sl

 
TYPE type_phys_point    ! Individual physics point characteristics

  REAL(KIND=JPRB) :: gelam, gelat, gemu
  REAL(KIND=JPRB) :: geclo, geslo, gm, gaw
  REAL(KIND=JPRB) :: gnordl, gnordm, gsqm2
  REAL(KIND=JPRB) :: rcolon, rsilon 
  REAL(KIND=JPRB) :: rindx, rindy 
  REAL(KIND=JPRB) :: orog
  INTEGER(KIND=JPIM) :: ngplat       ! row number in the physics grid

  INTEGER(KIND=JPIM) :: neigh       ! number of neighbours in the dynamics grid 
                                    ! for going from the dynamics to the physics grid
  INTEGER(KIND=JPIM), POINTER :: nl0(:)    ! indexes in the interpolation buffer of 
                                           ! the dynamics neighbours
  REAL(KIND=JPRB), POINTER   :: wgt(:)     ! weights for every neighbouring 
                                           ! dynamics point 
END TYPE type_phys_point

TYPE(type_phys_point),ALLOCATABLE :: yphypoi(:)

TYPE type_dyn_point     ! Individual dynamics point characteristics
  INTEGER(KIND=JPIM) :: neigh       ! number of neighbours in the physics grid 
                                    ! for going from the physics to the dynamics grid
  INTEGER(KIND=JPIM), POINTER :: nl0(:)    ! indexes in the interpolation buffer of 
                                           ! the physics neighbours
  REAL(KIND=JPRB), POINTER   :: wgt(:)     ! weights for every neighbouring 
                                           ! physics point 
END TYPE type_dyn_point

TYPE(type_dyn_point),ALLOCATABLE :: ydynpoi(:)
 
!$OMP THREADPRIVATE(dyn_grid,dyn_sl,phys_grid,phys_sl)
!$OMP THREADPRIVATE(ydynpoi,yphypoi)
END MODULE yom_phys_grid