AtHeartEngineer/kaiju
1
0
Fork 0
mirror of https://github.com/JHUAPL/kaiju.git synced 2026-01-09 15:17:56 -05:00
Code Issues Packages Projects Releases Wiki Activity

last changes for OH run

Browse Source
This commit is contained in:
Elena Provornikova
2022-01-17 19:12:11 -07:00
parent 994e394018
commit b5b271df99
7 changed files with 38 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
kaipy/gamhelio/ConfigScripts/startup.config
View File

@@ -21,4 +21,14 @@ NO2 = 4
B0 = 1.e-3 ; 100 nT
n0 = 200. ; 200/cc
[Grid]
gridfile = /glade/u/home/elenap/gameraHelio/grid/heliogrid.h5
tMin = 0.1
tMax = 0.9
Rin = 1.
Rout = 10.
Ni = 256
Nj = 128
Nk = 256

2
kaipy/gamhelio/ConfigScripts/startupOH.config
View File

@@ -7,7 +7,7 @@ gameraIbcFile = innerbc.h5
IbcDir = /glade/u/home/elenap/gameraHelio/OHelio/ibc
[WSA]
wsafile = /glade/u/home/elenap/gameraHelio/OHelio/1AUfiles/wsa.h5
wsafile = /glade/u/home/elenap/gameraHelio/OHelio/1AUfiles/CR2095/wsa.h5
density_temperature_infile = no
gauss_smooth_width = 0 ; 8
plots = yes

3
kaipy/gamhelio/helioViz.py
View File

@@ -13,7 +13,8 @@ Tsolar = 1.e6
VMax = 800.
VMin = 300.
MagVCM = "inferno"
#MagVCM = "inferno"
MagVCM = "rainbow"
#inner helio
#DMax = 150.

4
kaipy/gamhelio/heliosphere.py
View File

@@ -205,9 +205,9 @@ class GamsphPipe(GameraPipe):
Q = self.GetVar(vID,sID,vScl,doVerb)
#cell centered values from the last cell
Qi = Q[-1,:,:]
#Qi = Q[-1,:,:]
#cell centered values from the first cell
#Qi = Q[0,:,:]
Qi = Q[0,:,:]
#jd_c = self.MJDs[sID]
#print ('jd_c = ', jd_c)
return Qi

11
kaipy/gamhelio/wsa2gamera/params.py
View File

@@ -21,3 +21,14 @@ class params():
self.B0 = config.getfloat('Normalization','B0')
self.n0 = config.getfloat('Normalization','n0')
self.gridfile = config['Grid']['gridfile']
self.tMin = config.getfloat('Grid','tMin')
self.tMax = config.getfloat('Grid','tMax')
self.Rin = config.getfloat('Grid','Rin')
self.Rout = config.getfloat('Grid','Rout')
self.Ni = config.getint('Grid','Ni')
self.Nj = config.getint('Grid','Nj')
self.Nk = config.getint('Grid','Nk')

18
scripts/helio/heliopic.py
View File

@@ -121,15 +121,15 @@ if __name__ == "__main__":
else:
print ("Pic is empty. Choose pic1 or pic2 or pic3")
#add time (upper left)
if (pic == "pic1" or pic == "pic2"):
gsph.AddTime(nStp,AxL0,xy=[0.025,0.875],fs="x-large")
elif (pic == "pic3"):
gsph.AddTime(nStp,AxL0,xy=[0.015,0.82],fs="small")
elif (pic == "pic4"):
gsph.AddTime(nStp,Ax,xy=[0.015,0.92],fs="small")
else:
print ("Pic is empty. Choose pic1 or pic2 or pic3")
# add time (upper left)
# if (pic == "pic1" or pic == "pic2"):
# gsph.AddTime(nStp,AxL0,xy=[0.025,0.875],fs="x-large")
# elif (pic == "pic3"):
# gsph.AddTime(nStp,AxL0,xy=[0.015,0.82],fs="small")
# elif (pic == "pic4"):
# gsph.AddTime(nStp,Ax,xy=[0.015,0.92],fs="small")
# else:
# print ("Pic is empty. Choose pic1 or pic2 or pic3")
#TO DO: Add wsa info (lower left) instead of Solar wind params
#gsph.AddSW(nStp,AxL,xy=[0.625,0.025],fs="small")

5
scripts/helio/ih2oh.py
View File

@@ -47,7 +47,7 @@ print (B0OH, n0OH, V0OH, T0OH)
############### READ GAMERA solution at 1 AU #####################
f = h5py.File(prm.wsaFile,'r')
#the latest Step saved in inner helio solution wsa.h5
step = 'Step#4'
step = 'Step#2'
f[step].attrs.keys()
Nphi, Nth, Nr = np.shape(f[step]['Vx'])
@@ -76,7 +76,8 @@ print (theta_wsa_c.shape, phi_wsa_c.shape)
#these are normilized according to inner helio normalization
Vr = (f[step]['Vx'][:]*x[:] + f[step]['Vy'][:]*y[:] + f[step]['Vz'][:]*z[:])/r[:]
Br = f[step]['Br'][:]
#Br = f[step]['Br'][:]
Br = (f[step]['Bx'][:]*x[:] + f[step]['By'][:]*y[:] + f[step]['Bz'][:]*z[:])/r[:]
Rho = f[step]['D'][:]
T = f[step]['P'][:]/f[step]['D'][:]