Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Climate Unit Tests #181

Open
natashabatalha opened this issue Jan 17, 2024 · 0 comments
Open

Climate Unit Tests #181

natashabatalha opened this issue Jan 17, 2024 · 0 comments

Comments

@natashabatalha
Copy link
Owner

Place to add inputs for climate unit tests. Current cases from the tutorials:

  1. WASP-39b, 1 convective zone, Guillot profile guess, CK (no on the fly)
  2. Teff~900 BD, 500K guess, 1 convective zone guess, CK (no on the fly)

For additional suggestions post the full snippet that includes a code block e.g.:

opacity_ck = jdi.opannection(ck_db=ck_db)

cl_run = jdi.inputs(calculation="browndwarf", climate = True)
teff= 1000 # Effective Temperature of your Brown Dwarf in K
grav = 1000 # Gravity of your brown dwarf in m/s/s

cl_run.gravity(gravity=grav, gravity_unit=u.Unit('m/(s**2)')) # input gravity
cl_run.effective_temp(teff) # input effective temperature

nofczns = 1 # number of convective zones initially. Let's not play with this for now.

nstr_upper = 83 # top most level of guessed convective zone
nstr_deep = nlevel -2 # this is always the case. Dont change this
nstr = np.array([0,nstr_upper,nstr_deep,0,0,0]) # initial guess of convective zones

# Here are some other parameters needed for the code.
rfacv = 0.0 #we are focused on a brown dwarf so let's keep this as is

cl_run.inputs_climate(temp_guess= temp_guess, pressure= pressure,
                      nstr = nstr, nofczns = nofczns , rfacv = rfacv)
                      
out = cl_run.climate(opacity_ck, save_all_profiles=True,with_spec=True)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant