# L. Fita, LMD. October 2016
# Compilation of all required Fortran modules in `camelot' with gcc compilers
# A following python 'library' creation is required

# FOLLOWING: http://stackoverflow.com/questions/17156125/f2py-exposing-parameters-from-used-modules
# FOLLOWING: http://osdir.com/ml/python.f2py.user/2006-08/msg00019.html
# module load gnu/4.4.7
# module load python/2.7-anaconda
# export LD_LIBRARY_PATH=/opt/anaconda/lib/:$LD_LIBRARY_PATH
# Diagnostics compilation _______
# f2py -m module_ForDiag -h module_ForDiag.pyf module_generic.f90 module_ForDiagnosticsVars.f90 module_ForDiagnostics.f90 >& run_f2py_ForDiag_wrapper.log
# make
# f2py -L/opt/anaconda/lib/ -m module_ForDiag -c --fcompiler=gnu95 module_ForDiag.pyf module_generic.f90 module_ForDiagnosticsVars.f90 module_ForDiagnostics.f90 >& run_f2py_ForDiag_compile.log
# Interpolation compilation _______
# f2py -m module_ForInt -h module_ForInt.pyf module_generic.f90 module_ForInterpolate.f90  >& run_f2py_ForInt_wrapper.log
# f2py -L/opt/anaconda/lib/ -m module_ForInt -c --fcompiler=gnu95 module_ForInt.pyf module_generic.f90 module_ForInterpolate.f90  >& run_f2py_ForInt_compile.log
#
# using netcdf4.3.3.1gfortran
## module load netcdf4/4.3.3.1-gfortran

FC 		= /usr/bin/gfortran
FCF 	= -c -fPIC
LIB_INC         = 
RM 		= rm -f
F2PY    = /opt/canopy-1.3.0/Canopy_64bit/User/bin/f2py -L/opt/anaconda/lib/ --fcompiler=gnu95
#DBGFLAGS = -g -Wall -Wextra -Warray-temporaries -Wconversion -fimplicit-none -fbacktrace -ffree-line-length-0 -fcheck=all -ffpe-trap=zero,overflow,underflow -finit-real=nan
NCLIBFOLD	= /opt/netcdf43/gfortran/bin/
NCINCFOLD	= /opt/netcdf43/gfortran/include/
LIB_NETCDF	= -L$(NCLIBFOLD) -lnetcdff -lnetcdf -I$(NCINCFOLD) -static-libgfortran

FCFLAGS = $(FCF) $(DBGFLAGS)

# Sources for f2py
srcs = module_definitions.f90 module_basic.f90 module_generic.f90 module_scientific.f90
distrisrcs = $(srcs) module_DistriCorrection.f90
diagsrcs = $(srcs) module_ForDiagnosticsVars.f90 module_ForDiagnostics.f90
intsrcs = $(srcs) module_ForInterpolate.f90
scisrcs = $(srcs)
defsrcs = module_definitions.f90

####### ###### ##### #### ### ## #

MODULES = \
	module_definitions.o \
	module_basic.o \
	module_generic.o \
	module_NCgeneric.o \
	module_scientific.o \
	module_ForInterpolate.o \
	module_ForDiagnosticsVars.o \
	module_ForDiagnostics.o \
	module_DistriCorrection.o

FINTMODULES = \
	module_definitions.o \
	module_basic.o \
	module_generic.o \
	module_scientific.o \
	module_ForInterpolate.o

all : \
	module_ForDiagnostics.o \
	interpolate.o \
	DistriCorrection.o \
	pydistrimods.o \
	pydiagmods.o \
	pyintmods.o \
	pyscimods.o \
	pydefmods.o \
	trajectories_overlap.o

diags : \
	module_ForDiagnostics.o

clean :
	$(RM) *.mod *.o interpolate module_ForDiag.so module_ForInt.so module_ForDistriCorrect.so module_ForSci.so module_ForDef.so

########    #######
##
    ########

module_definitions.o:
	$(FC) $(FCFLAGS) module_definitions.f90

module_basic.o: module_definitions.o
	$(FC) $(FCFLAGS) module_basic.f90

module_generic.o: module_definitions.o module_basic.o
	$(FC) $(FCFLAGS) $(LIB_INC) module_generic.f90

module_NCgeneric.o: module_definitions.o module_basic.o module_generic.o
	$(FC) $(FCFLAGS) $(LIB_NETCDF) $(LIB_INC) module_NCgeneric.f90

module_scientific.o: module_definitions.o module_basic.o module_generic.o
	$(FC) $(FCFLAGS) module_scientific.f90

module_ForInterpolate.o: module_definitions.o module_generic.o
	$(FC) $(FCFLAGS) $(LIB_INC) module_ForInterpolate.f90

module_ForDiagnosticsVars.o: module_definitions.o module_basic.o module_generic.o
	$(FC) $(FCFLAGS) $(LIB_INC) module_ForDiagnosticsVars.f90

module_ForDiagnostics.o: module_definitions.o module_basic.o module_generic.o module_ForDiagnosticsVars.o
	$(FC) $(FCFLAGS) $(LIB_INC) module_ForDiagnostics.f90

module_DistriCorrection.o: module_definitions.o module_basic.o module_generic.o module_scientific.o
	$(FC) $(FCFLAGS) module_DistriCorrection.f90

FCEXEFINTCFLAGS = $(FINTMODULES) $(LIB_NETCDF) $(LIB_INC)
FCEXECFLAGS = $(MODULES) $(LIB_NETCDF) $(LIB_INC)

interpolate.o: module_definitions.o module_basic.o module_generic.o module_scientific.o module_ForInterpolate.o
	$(FC) $(FCEXEFINTCFLAGS) interpolate.f90 -o interpolate

DistriCorrection.o: module_definitions.o module_basic.o module_generic.o module_NCgeneric.o module_scientific.o module_DistriCorrection.o
	$(FC) $(FCEXECFLAGS) DistriCorrection.f90 -o DistriCorrection

pydistrimods.o:
	$(F2PY) -c -I$(NCINCFOLD) -m module_ForDistriCorrect $(distrisrcs) -L$(NCLIBFOLD)

pydiagmods.o:
	$(F2PY) -c -I$(NCINCFOLD) -m module_ForDiag $(diagsrcs) -L$(NCLIBFOLD)

pyintmods.o:
	$(F2PY) -c -I$(NCINCFOLD) -m module_ForInt $(intsrcs) -L$(NCLIBFOLD)

pyscimods.o:
	$(F2PY) -c -I$(NCINCFOLD) -m module_ForSci $(scisrcs) -L$(NCLIBFOLD)

pydefmods.o:
	$(F2PY) -c -I$(NCINCFOLD) -m module_ForDef $(defsrcs) -L$(NCLIBFOLD)

trajectories_overlap.o: module_definitions.o module_basic.o module_generic.o module_NCgeneric.o module_scientific.o
	$(FC) $(FCEXECFLAGS) trajectories_overlap.f90 -o trajectories_overlap