Mie scattering of clouds, and Refelection Spectrum (Forward modelling)

[HajimeKawahara]

Finally, I implemented the radiative transfer of the reflection spectrum with Ackerman and Marley cloud model + Mie scattering, but currently only for the forward modelling. #431 #432

Check the following jupyter notebook samples:

• documents/tutorials/Jupiter_cloud_model_using_amp.ipynb (see here) for the NH3 cloud modeling of Jupiter
• documents/tutorials/Jupiter_Hires_Modeling.ipynb (see here) for the forward modeling of the Subaru/HDS real Jupiter reflection spectrum, i.e. NH3 cloud+CH4 absorption.

Some important changes

• removed resolution argument in sop (SopInstProfile and SopRotation etc)
• numpy<2 added in setup.py because of installation error on pandas, numba etc
• from jax.config import config is removed -> from jax import config

New features we focus on

We added new art, opa, pdb and also amp (Atmospheric Micro Physics class) for this new enhancement.

spec.atmrt.ArtReflectPure

art class for the pure (i.e. no emission) reflection spectrum (ArtReflectEmis, reflection+emis has not been verified)

The RT algorithm is basically the flux-adding treatment by Robinson & Crisp 2018; Robinson & Salvador 2023,

We added the opacity computation of clouds in art

dtau_cld = art.opacity_profile_cloud_lognormal(sigma_extinction, rhoc, MMRc, rg, sigmag, gravity)

spec.pardb.PdbCloud

pdb (particulate database) for clouds/condensates.

from exojax.spec.pardb import PdbCloud
pdb_nh3 = PdbCloud("NH3")
#condensate substance density
rhoc = pdb_nh3.condensate_substance_density #g/cc

spec.opacont.OpaMie

opa for Mie scattering. The following code computes the Mie parameters using PyMieScatt. So, OpaMie is a kind of a wrapper of PyMieScatt currently. I will update this for retrieval.

from exojax.spec.opacont import OpaMie
opa = OpaMie(pdb_nh3, nus)
sigma_extinction, sigma_scattering, asymmetric_factor = opa.mieparams_vector_direct_from_pymiescatt(rg, sigmag)

spec.amp.AmpAmcloud

amp for the Akerman and Marley cloud model

For instance, amp.calc_ammodel computes rg (condensate typical radius) and mass mixing ratio from the cloud parameters

rg_layer, MMRc = amp.calc_ammodel(Parr, Tarr, mu, molmass_nh3, gravity, fsed, sigmag, Kzz, MMRbase_nh3)

AM cloud model in ExoJAX

The tutorial of the implementation of the AM cloud model in ExoJAX can be found in documents/tutorials/Ackerman_and_Marley_cloud_model.ipynb, see here.

We also compared the code with VIRGA:

• documents/userguide/amclouds_comparison_virga.ipynb (see here) if you are interested in the difference between the ExoJAX/amcloud and VIRGA.

[github-actions[bot]]

Unit Test Results

154 tests  +15   154 :heavy_check_mark: +15   35s :stopwatch: +3s     1 suites ±  0       0 :zzz: ±  0      1 files   ±  0       0 :x: ±  0 

Results for commit ed0d5001. ± Comparison against base commit 4de023c0.

:recycle: This comment has been updated with latest results.

[ykawashima]

Thank you very much!

I got a bit different result for the final figure of Ackerman and Marley cloud model tutorial. Is the figure outdated or did I do something wrong?

-> thanks! No, your results are correct. done it. (H.K.)

As far as I understand, it seems that vmap is not used here. https://github.com/HajimeKawahara/exojax/blob/74fb1dd4387d08037b375713ee90249a1b89ac1f/documents/tutorials/Ackerman_and_Marley_cloud_model.rst?plain=1#L102 https://github.com/HajimeKawahara/exojax/blob/74fb1dd4387d08037b375713ee90249a1b89ac1f/documents/tutorials/Ackerman_and_Marley_cloud_model.ipynb#L213

[HajimeKawahara]

@ykawashima Thanks a lot!