kohnakagawa

set xlabel "q{/Symbol s}" offset 0,-0.75 font "Helvetica,30"
set ylabel "<|h(q{/Symbol s})|^{2}>/{/Symbol s}^{4}" offset -4.5,0 font "Helvetica,30"
#set xtics 0,30000,6000000
#set ytics 0,0.2,1.0
unset key
set logscale xy
set tics scale 3.0, 1.5
set tics font "Helvetica, 30"
plot "q_vs_ha.txt" lw 6 
k=4.79766

kohnakagawa

set terminal postscript eps enhanced color solid
set output "prog_chem.eps"

set lmargin 0
set rmargin 0
set tmargin 0
set bmargin 0

set xrange [0:2500]

kohnakagawa

#define variance
variance <- function(x) var(x)*(length(x)-1)/length(x)
options(digits=10)
tardir <- "./result1"
rfile <- paste(tardir,"/mlocval.txt",sep="")
dat <- read.table(rfile)
dat <- data.matrix(dat)
h <- dat[,1]
h <- unique(h)
v <- dat[,4]

kohnakagawa

set terminal postscript eps enhanced color solid
set xlabel "{/Symbol W}^{2}" offset 3.0,-0.55 font "Helvetica,30"
set ylabel "Probability" offset -1.5,0 font "Helvetica,30"
set label 1 at graph 0.45,0.2  "N_{emb}/N_{amp} = 0" font "Helvetica,22"
set label 2 at graph 0.06,0.35 "N_{emb}/N_{amp} = 0.025" font "Helvetica,22"
set label 3 at graph 0.06,0.58 "N_{emb}/N_{amp} = 0.05" font "Helvetica,22"
set label 4 at graph 0.035,0.95 "N_{emb}/N_{amp} {/Symbol \263} 0.1" font "Helvetica,22"
set xtics 0,0.5,1.0
set mxtics 5
set ytics 0,0.5,1.0

kohnakagawa

#!/usr/bin/python
import sys, math, re
argvs = sys.argv

if (len(argvs) != 4):
    print "Usage $ python %s input_filename output_filename offset_length" % argvs[0]
    quit()
print "Generate ribon bilayer. Stripe length is 2*Lx"

def ReadLineGro(line):

kohnakagawa

 #!/usr/bin/python
import sys, math, re

class MolDat:
    def __init__(self):
        self.molnum   = -1
        self.molname  = ""
        self.atomnum  = -1
        self.atomname = ""
        self.coord    = []

kohnakagawa

#include <iostream>
#include <cassert>
#include <string>
#include <vector>
#include <fstream>
#include <cstdlib>
#include <ctime>
#include <algorithm>
#include "mvector3.hpp"

kohnakagawa

#!/usr/bin/python
import sys, math

class PrtclDat:
    def __init__(self):
        self.pos = [0.0, 0.0, 0.0]
        self.vel = [0.0, 0.0, 0.0]
        self.prp = -1
        self.chem = False
        self.unit = -1

kohnakagawa

#include <iostream>
#include <cassert>
#include <string>
#include <vector>
#include <fstream>
#include <cstdlib>
#include <ctime>
#include <algorithm>
#include "mvector3.hpp"

kohnakagawa

#!/usr/bin/python
import sys, datetime, math
import numpy as np

def ReadLinePDBform(line):
    ##01234567890123456789012345678901234567890123456789012345678901234567890123456789
    ##ATOM      1  N   DOPC    1      -5.364  -3.265  20.254  1.00  0.00      MEMB
    line_type = str(line[0:6].strip())
    ret = []    
    if line_type == "ATOM" or line_type == "HETATM":