# -*-shell-script-*-
# Ising 2D Monte Carlo path sampling
#setnolog
setoverwrite
#dirname = runs/ising2d-path      #path sampling
dirname = runs/ising2d-path-H0055-T15 #path sampling
#dirname = runs/ising2d-path-anl  #analysis
#--------------------------------------------
# Initializtion
#
NX = 100  NY = 100 input = [-1] initspin
#incnfile = ../ising2d/s100.cn readcn
#--------------------------------------------
# Physical parameters
J = 1     # interaction strength
H = 0.055  # applied field
#H = 0.06  # applied field
kBT = 1.5 # temperature
#kBT = 1.8 # temperature
srand48bytime #randomize number generator
#--------------------------------------------
# Plot Configuration
#
atomradius = 1.0 win_width = 480 win_height = 480 rotateangles = [ 0 0 0 1.7]
color00 = white color01 = blue backgroundcolor = gray80
openwin alloccolors rotate saverot plot
#openwin reversergb alloccolors rotate saverot plot
#sleep quit
#------------------------------------------------------------
# Path Sampling (manual test)
#input = [ 1 ] initpath finalcnfile = path0.pcn writepath
#input = [ 1 3500 ] copyPathtoS 
#input = [ 0.03 0.5 10000000 1 ] plotfreq = 10000 printfreq = 10000 SampleMCpath 
#finalcnfile = path1.pcn input = [ 1 ] writepath quit
#sleep
#------------------------------------------------------------
# Path Sampling (automatic)
input = [ 0.03 0.5 100000000 ] plotfreq = 1 printfreq = 10000
totalsteps = 30000   savecn = 1  savecnfreq = 100  openintercnfile
WalkonChain

finalcnfile = pathA.pcn input = [ 1 ] writepath
finalcnfile = pathB.pcn input = [ 2 ] writepath
finalcnfile = pathC.pcn input = [ 3 ] writepath
quit
#sleep
#------------------------------------------------------------
# Saddle points analysis
#input = [ 1 108 100 0.03 0.5 100000000 ] plotfreq = 1 printfreq = 10
#AnalyzeConfigs
#quit
#sleep
#------------------------------------------------------------
# Monte Carlo simulation
#totalsteps = 10000 finalcnfile = s100.cn writecn #equilibration
#
#totalsteps = 1000001   plotfreq = 1000 #1 
#saveprop = 1 savepropfreq = 100 printfreq = 1000
#savecn = 1   savecnfreq   = 1000 openintercnfile
#srand48bytime
#openpropfile MCrun
#finalcnfile = s100-a.cn writecn
#sleep
#quit
quit

