Skip to content

Commit

Permalink
Modified PBL to free troposphere transition layer.
Browse files Browse the repository at this point in the history
  • Loading branch information
lars2015 committed Dec 29, 2024
1 parent 1c31198 commit 3fd2bb6
Showing 1 changed file with 4 additions and 3 deletions.
7 changes: 4 additions & 3 deletions src/mptrac.c
Original file line number Diff line number Diff line change
Expand Up @@ -4432,13 +4432,14 @@ double pbl_weight(
const int ip) {

/* Get PBL pressure... */
double pbl;
double pbl, ps;
INTPOL_INIT;
INTPOL_2D(pbl, 1);
INTPOL_2D(ps, 1);

/* Get pressure range... */
const double p1 = pbl * 0.971832875032981;
const double p0 = pbl / 0.971832875032981;
const double p1 = pbl - 0.2 * (ps - pbl);
const double p0 = pbl;

/* Get weighting factor... */
if (atm->p[ip] > p0)
Expand Down

0 comments on commit 3fd2bb6

Please sign in to comment.