var_set_momentum Function

public function var_set_momentum(ndir) result(iw)

Set momentum variables

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: ndir

Return Value integer, (ndir)


Calls

proc~~var_set_momentum~~CallsGraph proc~var_set_momentum var_set_momentum proc~errormsg errormsg proc~var_set_momentum->proc~errormsg

Called by

proc~~var_set_momentum~~CalledByGraph proc~var_set_momentum var_set_momentum proc~phys_init phys_init proc~phys_init->proc~var_set_momentum proc~phys_init~2 phys_init proc~phys_init~2->proc~var_set_momentum proc~phys_init~3 phys_init proc~phys_init~3->proc~var_set_momentum proc~phys_init~4 phys_init proc~phys_init~4->proc~var_set_momentum proc~phys_activate phys_activate proc~phys_activate->proc~phys_init proc~phys_activate~2 phys_activate proc~phys_activate~2->proc~phys_init~2 proc~phys_activate~3 phys_activate proc~phys_activate~3->proc~phys_init~3 proc~phys_activate~4 phys_activate proc~phys_activate~4->proc~phys_init~4

Source Code

  function var_set_momentum(ndir) result(iw)
    integer, intent(in) :: ndir
    integer             :: iw(ndir), idir

    if (allocated(iw_mom)) call errormsg("Error: set_mom was already called")
    allocate(iw_mom(ndir))

    do idir = 1, ndir
      nwflux       = nwflux + 1
      nwfluxbc     = nwfluxbc + 1
      nw           = nw + 1
      iw_mom(idir) = nwflux
      iw(idir)     = nwflux
      write(cons_wnames(nwflux),"(A1,I1)") "m", idir
      write(prim_wnames(nwflux),"(A1,I1)") "v", idir
    end do
    
  end function var_set_momentum