steppingdirection.F90

Go to the documentation of this file.

module steppingdirection_mod
  use monc_component_mod, only : component_descriptor_type
  use state_mod, only : model_state_type, centred_stepping, forward_stepping
  implicit none

#ifndef TEST_MODE
  private
#endif

  public steppingdirection_get_descriptor

contains

  type(component_descriptor_type) function steppingdirection_get_descriptor()
    steppingdirection_get_descriptor%name="stepping_direction"
    steppingdirection_get_descriptor%version=0.1
    steppingdirection_get_descriptor%initialisation=>initialisation_callback
    steppingdirection_get_descriptor%timestep=>timestep_callback
  end function steppingdirection_get_descriptor

  subroutine initialisation_callback(current_state)
    type(model_state_type), target, intent(inout) :: current_state

    ! TODO - allow configuration to use forward stepping for scalars (th and q) and momentum (u,v,w)
    current_state%scalar_stepping = centred_stepping
    current_state%momentum_stepping = centred_stepping
  end subroutine initialisation_callback  

  subroutine timestep_callback(current_state)
    type(model_state_type), target, intent(inout) :: current_state
    
    current_state%field_stepping = centred_stepping
    if (current_state%timestep .eq. current_state%start_timestep) current_state%field_stepping = forward_stepping
  end subroutine timestep_callback  
end module steppingdirection_mod