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            
 
  ! ATM resolutions
  integer, parameter :: maxatmres = 10
  integer, allocatable, dimension(:) :: catm
 
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