C -*- Mode: Fortran; -*-
c-----------------------------------------------------------------
c     Ravi Samtaney
c     Copyright 2001
c     Princeton Plasma Physics Laboratory
c     All Rights Reserved
c-----------------------------------------------------------------
c     $Log: modules.F,v $
c     Revision 1.4  2005/10/11 22:01:54  samtaney
c     Added static Bflag.
c
c     Revision 1.3  2005/10/02 19:09:32  samtaney
c     Put in several options.
c
c     Revision 1.2  2005/09/23 18:50:42  samtaney
c     BoundaryPressTemp module added.
c
c     Revision 1.1.1.1  2005/09/12 17:28:10  samtaney
c     Original source.
c
c-----------------------------------------------------------------
c
#ifdef DYNAMIC
	include 'mesh_uparms.h'
	include 'mesh_parms.h'
#else
	include 'mesh_uparms_static.h'
	include 'mesh_parms_static.h'
#endif
	include 'mesh_common.h'
	include 'mesh.h'
	include 'properties.h'

#ifdef PARALLEL
	include 'mpistuff.h'
#endif
c
c        module time
c        save
c        double precision:: dt, ttot
c        end module time
c
	include 'profiles.h'
c
	module iounits
        integer, parameter:: istderr=0
        integer, parameter:: istdin=5
        integer, parameter:: istdout=6
        integer, parameter:: ihis=8
        integer, parameter:: imesh=16
        integer, parameter:: ifluid=21
        integer, parameter:: idump=71
        integer, parameter:: ibinoutput=41
        integer, parameter:: itecoutput=61
        integer, parameter:: itime=81
        integer, parameter:: idiag=91
	integer, parameter:: icons=92
	integer, parameter:: ieng=94
	integer, parameter:: iflux=97
	integer, parameter:: itflux=98
	integer, parameter:: ipflux=99
	integer, parameter:: imass=95
	integer, parameter:: irad=88
	end module iounits
	module BoundaryPressTemp
	save
c	double precision,allocatable:: pbdry(:)
c	double precision,allocatable:: pbdry(:)
	double precision:: pbdry, tempBdry
	end module BoundaryPressTemp
c
        module MagAxisToroidalField
        save
        double precision:: Btaxis
        end module MagAxisToroidalField
c
c
c       Options module
	module OptionFlags
	save
c       Which inviscid flux function to use
c        0 = Rusanov
c        1 = HLL
	integer:: fluxType
c       Most flags 0=no 1=yes
c       Use non conservative source term or not
	integer:: nonConFlag
c       Use constrained transport or not
	integer:: cTransFlag
c       Use EMCT scheme - no upwinding on B. - not
c       compatible with use of CTransFlag
	integer:: emCTFlag
c
c       1 = ideal MHD 0 = Resistive
	integer:: idealMHDFlag
c       eqsrcFlag: 0 = no eqsrc
c                  1 = use as if Ideal MHD
c                  2 = use resistive
	integer:: eqsrcFlag
c
c       Heating model flag
	integer:: qModelFlag
c       qModelFlag: 0 = No model
c                   1 = Parks Model
c                   2 = Instantaneous temp equilibration model
c                   3 = Energy source based on pressure
c                   4 = Electron heat flux source where mass source is active

c       Temp boundary condition flag
c        0 = adiabatic wall
c        1 = isothermal wall
	integer:: tempBCFlag
c       Pellet launch type -1: LFS, 1: HFS
	integer::  launchType
c       Static B field
	integer::  staticBFlag
	end module OptionFlags