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 を読み込み、そのコピー内の移動度を編集する
           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()