grids.F90

Go to the documentation of this file.

 
module grids_mod
  use datadefn_mod, only : default_precision
  implicit none
 
#ifndef TEST_MODE
  private
#endif
 
  integer, parameter, public :: z_index = 1, y_index = 2, x_index = 3
  integer, parameter, public :: primal_grid=1, dual_grid=2
 
  type, public :: horizontal_grid_configuration_type
     real(kind=default_precision) dx, dy     
     real(kind=default_precision) cx, cy     
     real(kind=default_precision) cx2, cy2   
     real(kind=default_precision) cxy        
     real(kind=default_precision) tcx, tcy
  end type horizontal_grid_configuration_type
 
  type, public :: vertical_grid_configuration_type
     real(kind=default_precision), dimension(:), allocatable :: &
          z,&      
          zn,&     
          dz,&     
          dzn,&    
          rdz,&    
          rdzn,&   
          rho,&    
          rhon,&   
          thref,&  
          dthref,& 
          tref,&   
          theta_init,& 
          temp_init,&  
          rh_init, &   
          u_init,&     
          v_init,&     
          prefn,& 
          pdiff,& 
          prefrcp,&
          rprefrcp,&
          czb,&   
          cza,&   
          czg,&   
          czh,&   
          tzc1,&  
          tzc2,&  
          tzd1,&  
          tzd2,&  
          w_subs,& 
          olubar,& 
          savolubar,&          
          olvbar,& 
          savolvbar,&
          olthbar,& 
          olzubar,& 
          olzvbar,& 
          olzthbar,& 
          dmpco,& 
          dmpcoz,& 
          tstarpr,& ! Temperature about which Taylor Expansion
          qsat,&
          dqsatdt,&
          qsatfac,&
          rneutml,&
          rneutml_sq,&
          buoy_co, &
          theta_rand, &
          theta_force, & 
          u_force, & 
          v_force, & 
          w_up,    & 
          w_dwn,   & 
          w_rand     
 
 
 
     real(kind=default_precision), dimension(:,:), allocatable :: q_init 
     real(kind=default_precision), dimension(:,:), allocatable :: q_rand 
     real(kind=default_precision), dimension(:,:), allocatable :: q_force 
 
     real(kind=default_precision), dimension(:,:), allocatable :: olqbar,olzqbar
    ! time varying forcing terms
     real(kind=default_precision), dimension(:,:), allocatable :: wsubs_time_vary 
     
     real(kind=default_precision), dimension(:), allocatable :: hgd
     real(kind=default_precision) :: czn, zlogm, zlogth, vk_on_zlogm
     integer, dimension(:), allocatable :: kgd
     integer :: kdmpmin
  end type vertical_grid_configuration_type
 
  type, public :: grid_configuration_type
     type(horizontal_grid_configuration_type) :: horizontal
     type(vertical_grid_configuration_type) :: vertical
  end type grid_configuration_type
 
  type, public :: global_grid_type
     type(grid_configuration_type) :: configuration
     real(kind=default_precision), dimension(3) :: bottom,&     
          top,&        
          resolution   
     integer, dimension(3) ::    size
     logical, dimension(3) :: active = (/ .false., .false., .false. /) 
     integer :: dimensions = 0
  end type global_grid_type
 
  type, public :: local_grid_type
     integer, dimension(3) :: start,& !< Global start coordinates of local grid
          end,& !< Global end coordinates of local grid
          size,& !< Grid points held locally
          halo_size,& !< Grid point size of the halo (halo_depth)
          local_domain_start_index,& !< The start index of computation data (local data is halo->data->halo so this precomputes the data start)
          local_domain_end_index
     logical, dimension(3) :: active = (/ .false., .false., .false. /) 
     integer, dimension(:,:), allocatable :: neighbours , corner_neighbours
     integer :: dimensions = 0
  end type local_grid_type
end module grids_mod