##############################################################################
#Test for multibinit bound_model = 3
##############################################################################

#print options
prt_model = 4 
dipdip_prt = 0 

# Construction of harmonic lattice part
ngqpt =  2 2 2
dipdip  = 1 
ncell = 2 2 2

# Optimize a coefficient in th effective potential 
bound_model = 3
bound_rangePower = 6 6  

#%%<BEGIN TEST_INFO>
#%% [setup]
#%% executable = multibinit
#%% system_xml = t100_model.xml
#%% md_hist = t98.nc
#%% test_chain = t98.in, t99.in, t100.in
#%% [files]
#%% files_to_test = 
#%%   t100.out, tolnlines = 0, tolabs = 0.0, tolrel = 0.0, fld_options = -easy;
#%%   t100_coeffs.xml, tolnlines = 0, tolabs = 0.0, tolrel = 0.0, fld_options = -easy;
#%% [paral_info]
#%% max_nprocs = 1
#%% [extra_info]
#%% authors = M. M. Schmitt
#%% keywords = Effective potential, multibinit
#%% description = 
#%%   Test the bound_model option 3: Create the equivalent high order terms for 
#%%   diverging low order terms. 
#%% topics = LatticeModel,DynamicsMultibinit
#%%<END TEST_INFO>