grdvars.F90

Go to the documentation of this file.

module grdvars

  use gengrid_kinds, only : dbl_kind, real_kind, int_kind

  implicit none

  integer :: ni
  integer :: nj
  integer :: npx

  integer :: nx
  integer :: ny

  logical :: editmask
  logical :: debug
  logical :: do_postwgts
  logical :: roottask

  integer, parameter :: nv = 4.
  integer, parameter :: ncoord = 2*4.
  integer, parameter :: nverts = 2*4.
  integer, parameter ::  nvars = ncoord + nverts                   


  real(dbl_kind)     :: sg_maxlat
  integer(int_kind)  :: ipole(2)

  integer, parameter, dimension(nv) :: ivertct = (/0, -1, -1,  0/)
  integer, parameter, dimension(nv) :: jvertct = (/0,  0, -1, -1/)
  integer, dimension(nv) :: ivertcv
  integer, dimension(nv) :: jvertcv
  integer, dimension(nv) :: ivertcu
  integer, dimension(nv) :: jvertcu
  integer, dimension(nv) :: ivertbu
  integer, dimension(nv) :: jvertbu
  ! Super-grid source grid variables
  real(dbl_kind), allocatable, dimension(:,:)   :: x
  real(dbl_kind), allocatable, dimension(:,:)   :: y
  real(dbl_kind), allocatable, dimension(:,:)   :: dx
  real(dbl_kind), allocatable, dimension(:,:)   :: dy

  ! Output grid variables
  real(dbl_kind), allocatable, dimension(:,:) :: latct
  real(dbl_kind), allocatable, dimension(:,:) :: lonct
  real(dbl_kind), allocatable, dimension(:,:) :: latcv
  real(dbl_kind), allocatable, dimension(:,:) :: loncv
  real(dbl_kind), allocatable, dimension(:,:) :: latcu
  real(dbl_kind), allocatable, dimension(:,:) :: loncu
  real(dbl_kind), allocatable, dimension(:,:) :: latbu
  real(dbl_kind), allocatable, dimension(:,:) :: lonbu
  real(dbl_kind), allocatable, dimension(:,:) :: areact
  real(dbl_kind), allocatable, dimension(:,:) :: anglet
  real(dbl_kind), allocatable, dimension(:,:) :: angle
  real(dbl_kind), allocatable, dimension(:,:) :: angchk

  real(dbl_kind), allocatable, dimension(:,:,:) :: latct_vert
  real(dbl_kind), allocatable, dimension(:,:,:) :: lonct_vert

  real(dbl_kind), allocatable, dimension(:,:,:) :: latcv_vert
  real(dbl_kind), allocatable, dimension(:,:,:) :: loncv_vert

  real(dbl_kind), allocatable, dimension(:,:,:) :: latcu_vert
  real(dbl_kind), allocatable, dimension(:,:,:) :: loncu_vert

  real(dbl_kind), allocatable, dimension(:,:,:) :: latbu_vert
  real(dbl_kind), allocatable, dimension(:,:,:) :: lonbu_vert


  real(dbl_kind), allocatable, dimension(:) :: xlonct
  real(dbl_kind), allocatable, dimension(:) :: xlatct
  real(dbl_kind), allocatable, dimension(:) :: xangct

  real(dbl_kind), allocatable, dimension(:) :: xloncu
  real(dbl_kind), allocatable, dimension(:) :: xlatcu
  real(dbl_kind), allocatable, dimension(:) :: dlatbu
  real(dbl_kind), allocatable, dimension(:) :: dlatcv
  ! MOM6 fix fields
  real(real_kind), allocatable, dimension(:,:) :: wet4
  real(dbl_kind),  allocatable, dimension(:,:) :: wet8

  real(real_kind), allocatable, dimension(:,:) :: dp4
  real(dbl_kind),  allocatable, dimension(:,:) :: dp8

  ! CICE6 fields
  real(dbl_kind), allocatable, dimension(:,:) :: ulon
  real(dbl_kind), allocatable, dimension(:,:) :: ulat
  real(dbl_kind), allocatable, dimension(:,:) ::  htn
  real(dbl_kind), allocatable, dimension(:,:) ::  hte

  real(kind=real_kind), parameter :: minimum_depth = 9.5
  real(kind=real_kind), parameter :: maximum_depth = 6500.0
  real(kind=real_kind), parameter :: masking_depth = 0.0
  real(kind=real_kind), parameter :: maximum_lat = 88.0

  integer, parameter :: nar = 6
  integer, parameter, dimension(nar) :: catm = (/48, 96, 192, 384, 768, 1152/)

contains

  subroutine allocate_all

    allocate( x(0:nx,0:ny),  y(0:nx,0:ny) )
    allocate(  dx(nx,0:ny), dy(0:nx,ny) )

    allocate( latct(ni,nj), lonct(ni,nj) )
    allocate( latcv(ni,nj), loncv(ni,nj) )
    allocate( latcu(ni,nj), loncu(ni,nj) )
    allocate( latbu(ni,nj), lonbu(ni,nj) )

    allocate( areact(ni,nj), anglet(ni,nj), angle(ni,nj), angchk(ni,nj))

    allocate( latct_vert(ni,nj,nv), lonct_vert(ni,nj,nv) )
    allocate( latcv_vert(ni,nj,nv), loncv_vert(ni,nj,nv) )
    allocate( latcu_vert(ni,nj,nv), loncu_vert(ni,nj,nv) )
    allocate( latbu_vert(ni,nj,nv), lonbu_vert(ni,nj,nv) )

    allocate( xlonct(ni), xlatct(ni), xangct(ni) )
    allocate( xloncu(ni), xlatcu(ni) )
    allocate( dlatbu(ni), dlatcv(ni) )

    allocate( wet4(ni,nj) )
    allocate( wet8(ni,nj) )

    allocate(  dp4(ni,nj) )
    allocate(  dp8(ni,nj) )

    allocate( ulon(ni,nj), ulat(ni,nj) )
    allocate(  htn(ni,nj),  hte(ni,nj) )

  end subroutine allocate_all

end module grdvars