picaso
picaso copied to clipboard
Climate Unit Tests
Place to add inputs for climate unit tests. Current cases from the tutorials:
- WASP-39b, 1 convective zone, Guillot profile guess, CK (no on the fly)
- 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)