| ... | @@ -141,29 +141,25 @@ For example: |
... | @@ -141,29 +141,25 @@ For example: |
|
|
| deltaOv | (1xn) | delta18O value of the ambient vapour [permil] |
|
|
| deltaOv | (1xn) | delta18O value of the ambient vapour [permil] |
|
|
|
| dexcv | (1xn) | d-excess value of the ambient vapour [permil] |
|
|
| dexcv | (1xn) | d-excess value of the ambient vapour [permil] |
|
|
|
| deltaHp | (1xm) | delta2H value of the hydrometeor [permil] |
|
|
| deltaHp | (1xm) | delta2H value of the hydrometeor [permil] |
|
|
|
| dexcp | (1xm) |
|
|
| dexcp | (1xm) | d-excess value of the hydrometeor [permil] |
|
|
|
d-excess value of the hydrometeor [permil]
|
|
| deltaHp_eq | (1xm) | delta2H value of equilibrium vapour from the hydrometeor at ambient temperature [permil] |
|
|
|
deltaHp_eq
|
|
| dexcp_eq | (1xm) | d-excess value of equilibrium vapour from the hydrometeor at ambient temperature [permil] |
|
|
|
(1xm)
|
|
| z_intsteps | (1xk) | Height above sea level of the individual integration steps [m] (z_intsteps defines the integration grid and is different from z if time integration was chosen). k = number of integration steps. |
|
|
|
delta2H value of equilibrium vapour from the hydrometeor at ambient temperature [permil]
|
|
|
|
|
dexcp_eq
|
|
|
|
|
(1xm)
|
|
|
|
|
d-excess value of equilibrium vapour from the hydrometeor at ambient temperature [permil]
|
|
|
|
|
z_intsteps
|
|
|
|
|
(1xk)
|
|
|
|
|
Height above sea level of the individual integration steps [m] (z_intsteps defines the integration grid and is different from z if time integration was chosen). k = number of integration steps.
|
|
|
|
|
|
|
|
|
|
|
**Examples:**
|
|
|
|
|
|
|
|
|
A hydrometeor is falling from 3500 meters above sea level to the ground at 500 m a.s.l. (with surface pressure of 950 hPa). The hydrometeor falls through background air, which is defined by a (moist-)adiabatically ascending air parcel from the surface with starting values of T=12°C., h=75%, delta2H=-150 permil and d-excess=10 permil. Default settings are used if not defined by the namelist (see above), i.e., integration after time, mass increase allowed, explicit treatment of all isotope mass equations, wegener-bergeron-findeisen formation mechanism, liquid terminal velocity, drop temperature different from that of the environment, which can be modified by both evaporation and heat diffusion and which is used to calculate the fractionation factor.
|
|
|
|
|
|
|
|
Examples:
|
|
```
|
|
|
A hydrometeor falling from 3500 meters above sea level to the ground at 500 m a.s.l. (with surface pressure of 950 hPa). The hydrometeor falls through background air, which is defined by a (moist-)adiabatically ascending air parcel from the surface with starting values of T=12°C., h=75%, delta2H=-150 permil and d-excess=10 permil. Default settings are used if not defined by the namelist (see above), i.e., integration after time, mass increase allowed, explicit treatment of all isotope mass equations, wegener-bergeron-findeisen formation mechanism, liquid terminal velocity, drop temperature different from that of the environment, which can be modified by both evaporation and heat diffusion and which is used to calculate the fractionation factor.
|
|
|
|
|
|
|
|
|
|
namelist.p0=95000 % Surface pressure 950hPa
|
|
namelist.p0=95000 % Surface pressure 950hPa
|
|
|
namelist.z_form=3500; % Formation height at 3500 m above sea level
|
|
namelist.z_form=3500; % Formation height at 3500 m above sea level
|
|
|
z=linspace(500,3500,3000) % 3000 height steps from 500 to 3500 m
|
|
z=linspace(500,3500,3000) % 3000 height steps from 500 to 3500 m
|
|
|
[T,h,~,~,~]=rayleigh_ascent(12,0.75,namelist.p0,-150,10,z);
|
|
[T,h,~,~,~]=rayleigh_ascent(12,0.75,namelist.p0,-150,10,z);
|
|
|
% Profiles of T and h are defined by an air parcel that ascends % (moist-)adiabatically from 950hPa with starting temperature of 12°C.
|
|
% Profiles of T and h are defined by an air parcel that ascends
|
|
|
% and starting relative humidity of 75%. Isotopic values (-150 and 10) % have to be given to rayleigh_ascent() as input, but are not of
|
|
% (moist-)adiabatically from 950hPa with starting temperature of 12°C.
|
|
|
|
% and starting relative humidity of 75%. Isotopic values (-150 and 10)
|
|
|
|
% have to be given to rayleigh_ascent() as input, but are not of
|
|
|
% relevance here.
|
|
% relevance here.
|
|
|
diam=1; % Hydrometeor diameter is 1mm
|
|
diam=1; % Hydrometeor diameter is 1mm
|
|
|
[profiles]=below_cloud_model({diam,’end’},[-150 0],[10 0],h,T,z,namelist)
|
|
[profiles]=below_cloud_model({diam,’end’},[-150 0],[10 0],h,T,z,namelist)
|
|
|
|
``` |