import json
       import os
       import sys

       f=open('sim.json')              #sim.json 파일 열기
       lines=f.readlines()
       f.close()
       lines="".join(lines)    #텍스트를 Python JSON 객체로 변환
     data = json.loads(lines)

     #값 편집하기(Firefox를 JSON 뷰어로 사용하면
     # 어떤 값을 편집해야 하는지 파악하는 데 도움이 됩니다)
     # 여기서는 layer 1의 이동도를 편집합니다
     data['epitaxy']['segment1']['shape_dos']['mue_y']=1.0


     #JSON 객체를 다시 문자열로 변환
     jstr = json.dumps(data, sort_keys=False, indent='\t')

     #다시 디스크에 쓰기
     f=open('sim.json',"w")
     f.write(jstr)
     f.close()

     #oghma_core를 사용해 시뮬레이션 실행
     os.system("oghma_core.exe")

f=open('sim_info.dat')
lines=f.readlines()
f.close()
lines="".join(lines)
data = json.loads(lines)

f=open('out.dat',"a")
f.write(str(data["Voc"])+"\n");
f.close()

       import json
       import os
       import shutil

       # 스크립트는 기준 시뮬레이션 디렉터리에서 시작되어야 합니다
       # (즉, 기준 sim.json이 들어 있는 디렉터리).
       base_dir = os.getcwd()

       mobilities = [1e-5, 1e-6, 1e-7, 1e-8]

       for mu in mobilities:
           # 읽기 쉽고 정렬도 잘 되는 디렉터리 이름 생성
           run_dir = os.path.join(base_dir, f"mu_{mu:.0e}")
           os.makedirs(run_dir, exist_ok=True)

           # 기준 sim.json을 실행 디렉터리로 복사
           src_sim = os.path.join(base_dir, "sim.json")
           dst_sim = os.path.join(run_dir, "sim.json")
           shutil.copyfile(src_sim, dst_sim)

           # 복사된 sim.json을 불러와 THAT copy에서 이동도 편집
           f = open(dst_sim)
           lines = f.readlines()
           f.close()
           lines = "".join(lines)
           data = json.loads(lines)

           # 값 편집: 첫 번째 소자 층의 이동도 설정
           data['epitaxy']['segment1']['shape_dos']['mue_y'] = mu

           jstr = json.dumps(data, sort_keys=False, indent='\t')
           f = open(dst_sim, "w")
           f.write(jstr)
           f.close()

           # 실행 디렉터리로 이동하여 այնտեղ서 해석기 실행
           os.chdir(run_dir)

           # 명령이 "augment.nano"인 경우 아래 문자열을 바꾸십시오.
           os.system("oghma_core.exe")

           # 다음 실행을 위해 기준 디렉터리로 복귀
           os.chdir(base_dir)

       import json
       import os

       # 선택 사항: 플로팅(matplotlib 필요)
       import matplotlib.pyplot as plt

       base_dir = os.getcwd()

       # 이동도 값은 앞서 생성한 실행 디렉터리와 일치해야 합니다
       mobilities = [1e-5, 1e-6, 1e-7, 1e-8]

       inv_mu = []
       vocs = []

       for mu in mobilities:
           run_dir = os.path.join(base_dir, f"mu_{mu:.0e}")
           info_path = os.path.join(run_dir, "sim_info.dat")

           # JSON 출력 파일을 읽고 Voc 추출
           f = open(info_path)
           lines = f.readlines()
           f.close()
           lines = "".join(lines)
           data = json.loads(lines)

           voc = data["Voc"]

           inv_mu.append(1.0 / mu)
           vocs.append(voc)

       # 작은 요약 표를 디스크에 기록
       out = open("voc_vs_inv_mobility.dat", "w")
       out.write("# inv_mobility(1/mu)    Voc(V)\n")
       for x, y in zip(inv_mu, vocs):
           out.write(f"{x:.6e}    {y}\n")
       out.close()

       # 역이동도에 대해 Voc 플로팅
       plt.plot(inv_mu, vocs, "o-")
       plt.xlabel("Inverse mobility (1/μ)")
       plt.ylabel("Open-circuit voltage Voc (V)")
       plt.title("Voc vs inverse mobility")
       plt.grid(True)
       plt.tight_layout()
       plt.show()