ginit.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	ginit (msave, first, nvar, naux, start_time)
	Initializes solution on all grids at level of grid msave, by calling qinit(). More...

Function/Subroutine Documentation

ginit()

subroutine ginit	(		msave,
		logical	first,
			nvar,
			naux,
			start_time
	)

Initializes solution on all grids at level of grid msave, by calling qinit().

If first = true, (first call to init), then allocate the soln storage area too, else was already allocated.

Parameters

msave	First grid in the grid list that stores all grids at this level
first	Is this the first call to init?
nvar	number of equations for the system
naux	number of auxiliary variables
start_time	start time of current simulation

Definition at line 18 of file ginit.f.

References amr_module::alloc, amr_module::cornxlo, amr_module::cornylo, amr_module::hxposs, amr_module::hyposs, igetsp(), amr_module::levelptr, amr_module::mxnest, amr_module::ndihi, amr_module::ndilo, amr_module::ndjhi, amr_module::ndjlo, amr_module::needs_to_be_set, amr_module::nestlevel, amr_module::nghost, amr_module::node, qinit(), amr_module::rnode, setaux(), amr_module::store1, amr_module::store2, amr_module::storeaux, and amr_module::timemult.

Referenced by domain(), and setgrd().

c
      use amr_module
      implicit double precision (a-h,o-z)
      logical first
      

c ::::::::::::::::::::::::::::: GINIT ::::::::::::::::::::::::
c
c  initializes soln on all grids at 'level'  by calling qinit
c  if first = true, (first call to init), then allocate the
c  soln storage area too, else was already allocated.
c
c :::::::::::::::::::::::::::::::::::::::;::::::::::::::::::::

      if (msave .eq. 0) go to 99

      level = node(nestlevel,msave)
      hx    = hxposs(level)
      hy    = hyposs(level)
      mptr  = msave
 
 10       nx      = node(ndihi,mptr) - node(ndilo,mptr) + 1
          ny      = node(ndjhi,mptr) - node(ndjlo,mptr) + 1
          mitot   = nx + 2*nghost
          mjtot   = ny + 2*nghost
          corn1   = rnode(cornxlo,mptr)
          corn2   = rnode(cornylo,mptr)
          if(.not. (first)) go to 20
              loc                 = igetsp(mitot*mjtot*nvar)
              node(store1,mptr)   = loc
              if (naux .gt. 0) then
                locaux              = igetsp(mitot*mjtot*naux)
                do k = 1, mitot*mjtot*naux,naux  ! set first component of aux to signal that it
                   alloc(locaux+k-1) = needs_to_be_set ! needs val, wasnt copied from other grids
                end do
                
                call setaux(nghost,nx,ny,corn1,corn2,hx,hy,
     &                    naux,alloc(locaux))
              else 
                locaux = 1
              endif
              node(storeaux,mptr) = locaux
              if (level .lt. mxnest) then
                loc2              = igetsp(mitot*mjtot*nvar)
                node(store2,mptr) = loc2
              endif
              rnode(timemult, mptr) = start_time
              go to 30
 20       continue
c
c  if 2nd time through, put initial values in store2 so finer grids
c  can be advanced with interpolation of their boundary values.
c  new time soln should still be in location store1.
c
          loc     = node(store2,mptr)
          locaux  = node(storeaux,mptr)
c
   30     continue
          call qinit(nvar,nghost,nx,ny,corn1,corn2,hx,hy,
     &               alloc(loc),naux,alloc(locaux))

c
          mptr  = node(levelptr, mptr)
      if (mptr .ne. 0) go to 10
c
c
 99   continue
      return

Here is the call graph for this function:

Here is the caller graph for this function: