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(:,:)   :: angq

  real(dbl_kind), allocatable, dimension(:,:)   :: dx
  real(dbl_kind), allocatable, dimension(:,:)   :: dy
  real(dbl_kind), allocatable, dimension(:,:) :: xsgp1
  real(dbl_kind), allocatable, dimension(:,:) :: ysgp1

  ! 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(:,:,:) :: 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(:) :: 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

contains

  subroutine allocate_all

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

    allocate( xsgp1(0:nx,0:ny+1), ysgp1(0:nx,0:ny+1) )

    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) )

    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) )
    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