state.F90

Go to the documentation of this file.

 
module state_mod
  use collections_mod, only : hashmap_type
  use grids_mod, only : global_grid_type, local_grid_type
  use prognostics_mod, only : prognostic_field_type
  use communication_types_mod, only : halo_communication_type
  use datadefn_mod, only : default_precision
  implicit none
 
#ifndef TEST_MODE
  private
#endif
 
  integer, parameter, public :: centred_stepping=0, forward_stepping=1, prescribed_surface_fluxes=0, prescribed_surface_values=1
  integer, parameter, public :: time_termination_reason=0, timestep_termination_reason=1, message_termination_reason=2, &
       walltime_termination_reason=3
 
  type, public :: parallel_state_type
     integer :: processes, & !> Total number of processes
          my_rank,&           !> My process rank in the system
          neighbour_comm,&     !> Neighbour communicator
          monc_communicator=-1, io_communicator=-1, corresponding_io_server_process
     integer, dimension(3) :: &
          my_coords,&         !> My process coordinates in each dimension
          dim_sizes
     logical, dimension(3,2) :: wrapped_around
     procedure(), nopass, pointer :: decomposition_procedure => null()
  end type parallel_state_type
  
  type, public :: model_state_type
    logical :: continue_timestep=.true., initialised=.false., continuation_run=.false.
    logical :: use_viscosity_and_diffusion=.true., &
       use_surface_boundary_conditions=.true., backscatter=.true.
 
    type(hashmap_type) :: options_database
    type(global_grid_type) :: global_grid
    type(local_grid_type) :: local_grid
    type(parallel_state_type) :: parallel
    type(prognostic_field_type) :: u, w, v, th, p, zu, zw, zv, zth, su, sw, sv, sth, savu, savv, & 
         savw, vis_coefficient, &
         diff_coefficient, dis, dis_th, &
         ! Heating rates from socrates
         sth_lw, sth_sw
    type(prognostic_field_type), dimension(:), allocatable :: q, zq, sq, disq
    ! longwave and shortwave downwelling flux at the surface
    real(kind=default_precision), dimension(:,:), allocatable :: sw_down_surf, lw_down_surf
    type(halo_communication_type) :: viscosity_halo_swap_state, diffusion_halo_swap_state
    real(kind=default_precision) :: time=.0_default_precision,& ! Model time in seconds
            dtm,& ! Modeltimestep (s)
            absolute_new_dtm, &
            thref0,&
            rhobous,&
            tsmth=1e-2_default_precision,&
            timestep_runtime,&
            local_divmax, global_divmax, cvis=0.0_default_precision, surface_temperature_flux, &
            surface_vapour_flux, theta_surf, surface_vapour_mixing_ratio, fbuoy, &
            fbuoynew, theta_virtual_surf, cmbc, rcmbc, ellmocon, velmax, velmin, aloginv, cneut, cfc, &
            surface_pressure=100000.0_default_precision, surface_reference_pressure = 100000.0_default_precision, &
            cvel, cvel_x, cvel_y, cvel_z, dtm_new, rmlmax, geostrophic_wind_rate_of_change_in_x, &
            geostrophic_wind_rate_of_change_in_y, surface_geostrophic_wind_x, surface_geostrophic_wind_y, &
            local_zumin, local_zumax, local_zvmin, local_zvmax, local_cvel_z
    real(kind=default_precision), dimension(:), allocatable :: lookup_table_velocity, &
         lookup_table_ustr, cq, abswmax
    real(kind=default_precision), dimension(:,:), allocatable :: psrce_recv_buffer_x, psrce_recv_buffer_y
    integer :: timestep=1, column_global_x, column_global_y, column_local_x, column_local_y, field_stepping, scalar_stepping, &
         momentum_stepping, number_q_fields=0, start_timestep=1, type_of_surface_boundary_conditions, lookup_table_entries, &
         cfl_frequency, termination_reason
    integer :: water_vapour_mixing_ratio_index=0, liquid_water_mixing_ratio_index=0, &
         rain_water_mixing_ratio_index=0, ice_water_mixing_ratio_index=0, &
         snow_water_mixing_ratio_index=0, graupel_water_mixing_ratio_index=0, & 
         psrce_x_hs_send_request, psrce_y_hs_send_request, psrce_x_hs_recv_request, psrce_y_hs_recv_request
    logical :: first_timestep_column, last_timestep_column, halo_column, first_nonhalo_timestep_column, &
         passive_q=.false., passive_th=.false., &
         use_time_varying_surface_values, use_anelastic_equations, & ! use_anelastic_equations or use Boussinesq
         saturated_surface, update_dtm=.false., calculate_th_and_q_init, origional_vertical_grid_setup=.true., &
         io_server_enabled
    logical, allocatable :: l_forceq(:)
    double precision :: model_start_wtime
 
    logical :: galilean_transformation=.true., fix_ugal=.false., fix_vgal=.false.
    real(kind=default_precision) :: ugal=0.,vgal=0.
    ! SOCRATES time variables are included in state since they need to be dumped
    real(kind=default_precision) :: rad_last_time=0.0
 
  end type model_state_type
end module state_mod