Climate Analysis Branch
In order to help ensure that this web analysis page remains available, we would greatly appreciate feedback on its use, particularily in the clasroom, for presentations or for research. Mail to Cathy Smith at (
Background Information
Dataset Information
Related Plot/Analysis
The ESRL website will be unavailable for 24 hours starting Friday, March 6th at 5:00pm MT due to building maintenance.

Variable Calculation Details

Calculated as T(p+1)-T(p)/dp where p is the level above the one selected (a one sided difference). Results per displayed as per 100mb to match d/dz values better.
Calculated as T(p+1)-T(p-1)/dz where dz is from theNCEP/NCAR Reanalysis geopotential height. Results displated as K/km. Positive values indicate a decrease of temperature upward.
500-700mb lapse rate:
(T(500)-T(700)(/(500mb-700mb). Positive values indicate a decrease of temperature upward.
LCL surface (.995sigma):
Based on Fortran program calcsound of Kerry Emanuel
Dew Point:
X = 1.-0.01*RH
DPD =(14.55+0.114*A)*X+ & ((2.5+0.007*A)*X)**3+(15.9+0.117*A)*X**14
Mixing Ratio
Calculated AT a level.
rh= decimel, p=mb, T=K
w=(rh*wsb)/(1.+((1.-rh)*wsb)/e) e=.622
Equivalent Potential Temperature
Calculated AT a level.
I use formula supplied by Paul Neiman which uses entropy. Will post later.
Freezing Level
Level (mb for now) which is less than freezing (273.15K). Pressure level is interpolated using log pressure weighting. If surface is freezing, level is "missing".
Moisture Flux:
specific humidity times the vector wind (u,v) at a level. units are g/kg * m/s
Integrated Moisture Flux:
specific humidity times vector wind integrated from the surface to 300mb.Each layer is multiplied by the pressure of that level and the final result is divided by g. Values are computed between levels. For example, the 925 to 850 layer as the average of 1 and vector wind for that layer multipled by 12500 Pascals. Units are xxx.