Comments (4)
While working on this, I noticed that
p_H2O = RH * params.esat(temp)
is not valid in the stratosphere. Assuming q_H2O is uniform in the stratosphere, p_H2O has to decrease with pressure. Will fix the bug and generalize code to arbitrary dry gas mixtures.
from pyrads.
@AndrewWilliams3142 I think your suggestion might not be quite correct yet. 50000ppmv CO2 presumably means that we want the dry background gas to consist of 5% CO2, 95% rest by volume. However, once you're pushing into hot territory where the overall composition is 50% H2O-50% dry by volume, then ppv_CO2 & q_CO2 should drop due to the added H2O
from pyrads.
Hey @ddbkoll ! That makes sense to me. For context, I'm trying to reproduce the results of Seeley and Jeevanjee, 2020 with PyRADS, but PyRADS seems to be underestimating the CO2 at high temperatures compared to the line-by-line code in their paper ('RFM' in the figure). I thought that perhaps these simplifications could be the issue? I'm not entirely sure how the RFM treats these high-CO2/H2O situations...
Also, I'm currently running PyRADS with 100 log-spaced pressure levels, but maybe that needs to be improved as well (there are quite a few 'wiggles' which are probably due to the discretization).
from pyrads.
Actually I just realized it's probably because PyRADS doesn't include the CO2 self-continuum... I'll code up the analytic fits from POPC and see if that fixes the problem ^^
EDIT: It looks like it's a combination of the spectral and the vertical resolution. Upping spectral resolution to 0.01 cm-1 and using 300 levels gives you a 'bump' (see figure), but there are still a lot of wiggles, and the bump is underestimated - I'm rerunning now with 1000 vertical levels to see if that fixes it. I think this might be it, but running with hundreds of vertical levels slows the code down a lot so it might not be that practical.
from pyrads.
Related Issues (13)
- Add an open-source license HOT 1
- Stratospheric water vapor inaccurate
- Make latent heat of vaporization a function of temperature
- Use numba to accelerate OpticalThickness calculation? HOT 8
- Numba causes silent errors for mixed-component atmospheres HOT 11
- OpticalThickness.py/compute_tau_H2ON2_CO2dilute HOT 1
- Adding isotopologues
- Upload scripts to reproduce Figures in Koll & Cronin (2018).
- Change pressure broadening coefficients from "air" to "mixed" HOT 1
- conda-installed gfortran on OSX: builds fail at runtime HOT 7
- Packaging this nice package
- Python 3 compatibility
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from pyrads.