Skip to content

Commit

Permalink
bugfix in conductivity calculations
Browse files Browse the repository at this point in the history
  • Loading branch information
zettergm authored and zettergm committed Feb 20, 2024
1 parent 1bbb24c commit c280fb7
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions src/gemini3d/conductivity.py
Original file line number Diff line number Diff line change
Expand Up @@ -7,19 +7,19 @@
"""

import numpy as np
from .phys_const import ms,qs,kB
from .phys_const import ms,qs
import gemini3d.msis
import scipy as sp
#import scipy as sp

def conductivity_reconstruct(time,dat,cfg,xg):
# neutral atmospheric information
atmos=gemini3d.msis.msis_setup(cfg,xg)

# composition and extraction (or creation) of full state variables
if "ns" in dat.keys(): # proxy for full state output file
ns=dat["ns"]
Ts=dat["Ts"]
vs1=dat["vs1"]
ns=np.array(dat["ns"]).transpose((1,2,3,0))
Ts=np.array(dat["Ts"]).transpose((1,2,3,0))
vs1=np.array(dat["vs1"]).transpose((1,2,3,0))
else:
p=1/2+1/2*np.tanh((xg["alt"]-220e3)/20e3)
shapevar=np.concatenate( (xg["lx"][:],[7]) )
Expand Down

0 comments on commit c280fb7

Please sign in to comment.