#!/usr/bin/env python
# coding: utf-8
# # Plot phiinrz*
# In[ ]:
import numpy as np
import matplotlib.pyplot as plt
import glob
import f90nml
from read_f90 import read_parameters, read_time # In-house module for GKV-diag outputs
### GKV parameters from gkvp_header.f90 ###
nxw = read_parameters("../src/gkvp_header.f90", "nxw", int)
nyw = read_parameters("../src/gkvp_header.f90", "nyw", int)
### GKV parameters from gkvp_namelist ###
nml = f90nml.read("../gkvp_namelist.001")
#print(nml)
dtout_ptn = nml["times"]["dtout_ptn"]
print("nxw =", nxw)
print("dtout_ptn =", dtout_ptn)
# In[ ]:
### Load data of mominxy* ###
filelist = sorted(glob.glob("./data/phiinrz_t*dat"))
# Time
#!/usr/bin/env python
# coding: utf-8
# # Plot trninkxky*
# In[ ]:
import numpy as np
import matplotlib.pyplot as plt
import glob
import f90nml
from read_f90 import read_parameters, read_time # In-house module for GKV-diag outputs
### GKV parameters from gkvp_header.f90 ###
nx=read_parameters("../src/gkvp_header.f90", "nx", int)
global_ny=read_parameters("../src/gkvp_header.f90", "global_ny", int)
### GKV parameters from gkvp_namelist ###
nml=f90nml.read("../gkvp_namelist.001")
#print(nml)
dtout_eng=nml["times"]["dtout_eng"]
print("nx =",nx)
print("global_ny =",global_ny)
print("dtout_eng =",dtout_eng)
# In[ ]:
# coding: utf-8
# # Plot phiinz*
# Alinz* also has the same format.
# In[ ]:
import numpy as np
import matplotlib.pyplot as plt
import glob
import f90nml
from read_f90 import read_parameters, read_time # In-house module for GKV-diag outputs
### GKV parameters from gkvp_header.f90 ###
global_nz=read_parameters("../src/gkvp_header.f90", "global_nz", int)
### GKV parameters from gkvp_namelist ###
nml=f90nml.read("../gkvp_namelist.001")
#print(nml)
n_tht=nml["nperi"]["n_tht"]
dtout_ptn=nml["times"]["dtout_ptn"]
print("global_nz =",global_nz)
print("n_tht =",n_tht)
print("dtout_ptn =",dtout_ptn)
# In[ ]: