#!/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[ ]: