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

WIP -- generate valid value masks from surface pressure #625

Open
wants to merge 1 commit into
base: develop
Choose a base branch
from

Conversation

burlen
Copy link
Collaborator

@burlen burlen commented Jun 3, 2021

resolves #526

This algorithm computes surface pressure from surface temperature, mean
sea level pressure and surface elevation using the dry adiabatic lapse
rate.
void compute_surface_pressure(size_t n_elem, const num_t *tas,
const num_t *psl, coord_t *z, num_t *ps)
{
num_t cp = 1.0; // kJ/kg.K
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

cp and cv should be in units of J/Kg/K. This would matter for gam since the units cancel, but it does matter for the GAM * zi parts.

I suggest cp=1004. and cv = 718. as the values used

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

Successfully merging this pull request may close these issues.

pressure from sea level pressure
2 participants