picaso icon indicating copy to clipboard operation
picaso copied to clipboard

Climate Unit Tests

Open natashabatalha opened this issue 5 months ago • 0 comments

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)

natashabatalha avatar Jan 17 '24 19:01 natashabatalha