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 :: mastertask

  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