OBJ= modules.o main.o setupgrid.o \
constructlrstates.o \
newton2b.o \
init2.o  \
setboundaryvalues.o \
mhddriver.o \
bdryexchange.o \
inviscidfluxhll.o  \
utils.o  \
svdcmp.o  \
hqr.o  \
elmhes.o  \
dPackgmres.o

#projection.o  \
#xchange.o  \
#projection3D.o  \
#xchange3D.o  \

#phi.o \
#viscousflux.o  \


DEBUG=
LIST=
#Lahey-Fujitsu flags
FFLAGS= -c --chk[aesu] --chkglobal --fix -g --trace  --trap -DDYNAMIC -DETAPERIODIC#-DPERIODIC #-DTWO_D #-DPERIODIC# -DTWO_HALF_D
#FFLAGS= -c -O --fix -DDYNAMIC -DETAPERIODIC #-DPARALLEL -I/usr/local/lff95/mpi/include

#F90= frt
F90= lf95


all: pelletExec

pelletExec: ${OBJ}
	${F90} -o pelletExec  ${OBJ}  -lblas ##-L/usr/local/lff95/mpi/lib -lmpich -lfmpich


modules.o:modules.F
	 ${F90}  ${FFLAGS} modules.F

gammafun.o:gammafun.F
	 ${F90}  ${FFLAGS} gammafun.F

normalpred.o:normalpred.F
	 ${F90}  ${FFLAGS} normalpred.F

electronheat.o:electronheat.F
	 ${F90}  ${FFLAGS} electronheat.F

electronheatparallel2.o:electronheatparallel2.F
	 ${F90}  ${FFLAGS} electronheatparallel2.F

utils.o:utils.F
	 ${F90}  ${FFLAGS} utils.F

pelletSource.o:pelletSource.F
	 ${F90}  ${FFLAGS} pelletSource.F

pelletSourceMC.o:pelletSourceMC.F
	 ${F90}  ${FFLAGS} pelletSourceMC.F

constructlrstates.o:constructlrstates.F
	 ${F90}  ${FFLAGS} constructlrstates.F

constructlrstates2.o:constructlrstates2.F
	 ${F90}  ${FFLAGS} constructlrstates2.F

constructlrstates4.o:constructlrstates4.F
	 ${F90}  ${FFLAGS} constructlrstates4.F

inviscidfluxroe.o:inviscidfluxroe.F
	 ${F90}  ${FFLAGS} inviscidfluxroe.F

inviscidfluxcd.o:inviscidfluxcd.F
	 ${F90}  ${FFLAGS} inviscidfluxcd.F

inviscidfluxrp.o:inviscidfluxrp.F
	 ${F90}  ${FFLAGS} inviscidfluxrp.F

inviscidfluxrp2.o:inviscidfluxrp2.F
	 ${F90}  ${FFLAGS} inviscidfluxrp2.F

inviscidfluxl.o:inviscidfluxl.F
	 ${F90}  ${FFLAGS} inviscidfluxl.F

inviscidfluxhll.o:inviscidfluxhll.F
	 ${F90}  ${FFLAGS} inviscidfluxhll.F

inviscidfluxhll2.o:inviscidfluxhll2.F
	 ${F90}  ${FFLAGS} inviscidfluxhll2.F

viscousterms.o:viscousterms.F
	 ${F90}  ${FFLAGS} viscousterms.F

main.o:main.F
	 ${F90}  ${FFLAGS} main.F

newdt.o:newdt.F
	 ${F90}  ${FFLAGS} newdt.F

output.o:output.F
	 ${F90}  ${FFLAGS} output.F

setboundaryvalues.o:setboundaryvalues.F
	 ${F90}  ${FFLAGS} setboundaryvalues.F

setboundaryvaluesP.o:setboundaryvaluesP.F
	 ${F90}  ${FFLAGS} setboundaryvaluesP.F

setupgrid.o:setupgrid.F
	 ${F90}  ${FFLAGS} setupgrid.F

mhddriver.o:mhddriver.F
	 ${F90}  ${FFLAGS} mhddriver.F

ctrans.o:ctrans.F
	 ${F90}  ${FFLAGS} ctrans.F

ctrans2.o:ctrans2.F
	 ${F90}  ${FFLAGS} ctrans2.F

emsrc.o:emsrc.F
	 ${F90}  ${FFLAGS} emsrc.F

projection.o:projection.F
	 ${F90}  ${FFLAGS} projection.F

projection3D.o:projection3D.F
	${F90}  ${FFLAGS} projection3D.F

xchange.o:xchange.F
	${F90}  ${FFLAGS} xchange.F

xchange3D.o:xchange3D.F
	${F90}  ${FFLAGS} xchange3D.F

init.o:init.F
	 ${F90}  ${FFLAGS} init.F

init2.o:init2.F
	 ${F90}  ${FFLAGS} init2.F

initesc.o:initesc.F
	 ${F90}  ${FFLAGS} initesc.F

initrm.o:initrm.F
	 ${F90}  ${FFLAGS} initrm.F

initrmVW.o:initrmVW.F
	 ${F90}  ${FFLAGS} initrmVW.F

initPellet.o:initPellet.F
	 ${F90}  ${FFLAGS} initPellet.F

initwave.o:initwave.F
	 ${F90}  ${FFLAGS} initwave.F

newton.o:newton.F
	 ${F90}  ${FFLAGS} newton.F

newton2.o:newton2.F
	 ${F90}  ${FFLAGS} newton2.F

newton2b.o:newton2b.F
	 ${F90}  ${FFLAGS} newton2b.F

mhdsolvesemi.o:mhdsolvesemi.F
	 ${F90}  ${FFLAGS} mhdsolvesemi.F

mhdsolvesemiRKTVD.o:mhdsolvesemiRKTVD.F
	 ${F90}  ${FFLAGS} mhdsolvesemiRKTVD.F

mhdsolvesemiRK4.o:mhdsolvesemiRK4.F
	 ${F90}  ${FFLAGS} mhdsolvesemiRK4.F

mhdsolveunsplit.o:mhdsolveunsplit.F
	 ${F90}  ${FFLAGS} mhdsolveunsplit.F

mhdsolverk2.o:mhdsolverk2.F
	 ${F90}  ${FFLAGS} mhdsolverk2.F

mhdsolveeuler.o:mhdsolveeuler.F
	${F90}  ${FFLAGS} mhdsolveeuler.F

mhdsolveeuler2.o:mhdsolveeuler2.F
	${F90}  ${FFLAGS} mhdsolveeuler2.F


bdryexchange.o:bdryexchange.F
	 ${F90}  ${FFLAGS} bdryexchange.F

setroevariables.o:setroevariables.F
	 ${F90}  ${FFLAGS} setroevariables.F

seteigenvalues.o:seteigenvalues.F
	 ${F90}  ${FFLAGS} seteigenvalues.F

seteigenvectors.o:seteigenvectors.F
	 ${F90}  ${FFLAGS} seteigenvectors.F

seteigenvectors2.o:seteigenvectors2.F
	 ${F90}  ${FFLAGS} seteigenvectors2.F

seteigenvectors3.o:seteigenvectors3.F
	 ${F90}  ${FFLAGS} seteigenvectors3.F

seteigenvectors3b.o:seteigenvectors3b.F
	 ${F90}  ${FFLAGS} seteigenvectors3b.F

seteigenvectors4.o:seteigenvectors4.F
	 ${F90}  ${FFLAGS} seteigenvectors4.F

setsls.o:setsls.F
	 ${F90}  ${FFLAGS} setsls.F

setalphas.o:setalphas.F
	 ${F90}  ${FFLAGS} setalphas.F

setdurl.o:setdurl.F
	 ${F90}  ${FFLAGS} setdurl.F

phi2.o:phi2.F
	 ${F90}  ${FFLAGS} phi2.F

phi.o:phi.F
	 ${F90}  ${FFLAGS} phi.F

divb.o:divb.F
	 ${F90}  ${FFLAGS} divb.F

diags.o:diags.F
	 ${F90}  ${FFLAGS} diags.F

svdcmp.o:svdcmp.F
	 ${F90}  ${FFLAGS} svdcmp.F

hqr.o:hqr.F
	 ${F90}  ${FFLAGS} hqr.F

elmhes.o:elmhes.F
	 ${F90}  ${FFLAGS} elmhes.F

dPackgmres.o:dPackgmres.F
	${F90}  ${FFLAGS} dPackgmres.F




clean:
	rm ${OBJ}


renew:
	touch *