#Job performance
 NPAR = 8       #Use 4 or 8 with 32 cores
 KPAR = 2       #1 or 2 or 4

#Method
 GGA = PE
 IVDW = 12              #12 is d3bj dispersion correction
 #LHFCALC = .TRUE.       #True is default for b3lyp
 LDIPOL=.TRUE.          #True should be used with polarized slabs
 IDIPOL = 3             #direction of the normal of slab surface
 DIPOL=0.5 0.5 0.5     #do not set it, VASP default is good for orthorombic
 #KGAMMA                #Best not to include it! 
 #ADDGRID               #Do not include it

 #Type of job
 IBRION = -1            #-1 is SP, 2 is most reliable, 1 is best when close to the minimum

#File writing
 LWAVE = FALSE          #Write WAVECAR? It can be hefty
 LCHARG = FALSE          #Write CHGCAR. Nice to have when restarting (through ICHARG) with small changes

#Calculation accuracy
 ENCUT = 900.000     	#Energy cut-off
 EDIFF = 1.0e-06       #SC-loop condition. 1E-5 is enough for optimizations
 SIGMA = 0.02           #Smearing width in eV
 PREC = Accurate        #Just use accurate
 ISMEAR = 0             #Type of smearing. Use 0 but idk
 NSW = 0              #Max number of ionic steps 0 if IBRION=-1, maybe 200 otherwise

#Calculation performance
 ALGO = Normal        #Fast or VeryFast are ok. Damped when orbitals are read from WAVECAR. But doesn't really matter maybe
 ISTART = 1            #read the wavecar? 0|1|2|3. Default is 1 if WAVECAR exists
 ICHARG = 11            #read CHGCAR? 0|1|2|4. I only understand 1
 NELM = 100             #Max number of SC steps 
 AMIN = 0.01            #I don't understand it. Use default (usually 0.1) for bulk and 0.01 for slabs
 NSIM = 4
 NELMIN = 5

# DOS
 LORBIT = 11
 NEDOS = 2001
 EMIN = -4.5
 EMAX = -3.5
 #NBANDS = 120


