| ... | ... | @@ -96,87 +96,52 @@ The following list shows possible field names for namelist. If they are not defi |
|
|
|
`below_cloud_model.m` returns a struct "profiles" with the arrays listed below. The array sizes depends on the model inputs.
|
|
|
|
|
|
|
|
Variables that describe the profile of ambient air are of size (1xn), where n is the length of the input array "z". "z" defines the grid on which these profiles are defined (background grid).
|
|
|
|
|
|
|
|
For example:
|
|
|
|
z: (1x3001) array; p: (1x3001) array
|
|
|
|
|
|
|
|
```
|
|
|
|
z: (1x3001) array; p: (1x3001) array
|
|
|
|
z(305) = 805 m
|
|
|
|
p(305) = 93000 Pa
|
|
|
|
→ The pressure at model level 305 is 930 hPa, which corresponds to 805 m above sea level.
|
|
|
|
→ The pressure at model level 305 is 930 hPa, which corresponds to 805 m above sea level.
|
|
|
|
```
|
|
|
|
|
|
|
|
Variables that describe properties of the hydrometeor are arrays of size (1xm), where m is the number of model levels, which are at and below the starting altitude (z_form) of the hydrometeor. The corresponding grid (hydrometeor grid) is defined in "z_out" (and "height")
|
|
|
|
|
|
|
|
For example:
|
|
|
|
|
|
|
|
```
|
|
|
|
z = 0:1:3000 -> (1x3001) array
|
|
|
|
z_form = 2000
|
|
|
|
z_out contains all z<=z_form:
|
|
|
|
z_out = 0:1:2000 -> (1x2001) array
|
|
|
|
m=n if z_form=z(end) (if namelist.z_form=0)
|
|
|
|
|
|
|
|
Array name
|
|
|
|
Size
|
|
|
|
Description
|
|
|
|
z
|
|
|
|
(1xn)
|
|
|
|
Height above sea level of all model levels [m]
|
|
|
|
z_out
|
|
|
|
(1xm)
|
|
|
|
Height above sea level of all model levels at and below the formation height of the hydrometeor [m]
|
|
|
|
height
|
|
|
|
(1xm)
|
|
|
|
Same as z_out.
|
|
|
|
h
|
|
|
|
(1xn)
|
|
|
|
Relative humidity of ambient air []
|
|
|
|
T
|
|
|
|
(1xn)
|
|
|
|
Temperature of ambient air [°C]
|
|
|
|
p
|
|
|
|
(1xn)
|
|
|
|
Pressure of ambient air [Pa]
|
|
|
|
exm
|
|
|
|
(1xn)
|
|
|
|
Excess moisture of ambient air (specific humidity divided by the saturation specific humidity over ice at ambient temperature; exm>=0) [kg/kg]
|
|
|
|
m
|
|
|
|
(1xm)
|
|
|
|
Hydrometeor mass [kg]
|
|
|
|
diam
|
|
|
|
(1xm)
|
|
|
|
Hydrometeor diameter [m]
|
|
|
|
Td
|
|
|
|
(1xm)
|
|
|
|
Hydrometeor temperature [°C]
|
|
|
|
Tdelta
|
|
|
|
(1xm)
|
|
|
|
Temperature difference between hydrometeor and surrounding air [°C]
|
|
|
|
equih
|
|
|
|
(1xm)
|
|
|
|
Relative humidity of the surrounding air with respect to Td, instead of T (equivalent relative humidity) []
|
|
|
|
masstend
|
|
|
|
(1xm)
|
|
|
|
Hydrometeor mass tendency (rate of change of mass dm/dt) [kg/s]
|
|
|
|
temptend
|
|
|
|
(1xm)
|
|
|
|
Hydrometeor temperature tendency (rate of change of temperature dTd/dt) [°C/s]
|
|
|
|
f
|
|
|
|
(1xm)
|
|
|
|
Remaining fraction of mass of the hydrometeor relative to the next higher model level []
|
|
|
|
ftot
|
|
|
|
(1xm)
|
|
|
|
Remaining fraction of mass of the hydrometeor relative to the initial mass (at z_form) []
|
|
|
|
evapfrac
|
|
|
|
(1xm)
|
|
|
|
Evaporated fraction of mass (current mass relative to the initial mass of the hydrometeor; evapfrac=1-ftot) []
|
|
|
|
deltaHv
|
|
|
|
(1xn)
|
|
|
|
delta2H 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]
|
|
|
|
deltaHp
|
|
|
|
(1xm)
|
|
|
|
delta2H value of the hydrometeor [permil]
|
|
|
|
dexcp
|
|
|
|
(1xm)
|
|
|
|
```
|
|
|
|
|
|
|
|
| Array name | Size | Description |
|
|
|
|
| ------ | ------ | ------ |
|
|
|
|
| z | (1xn) | Height above sea level of all model levels [m] |
|
|
|
|
| z_out | (1xm) | Height above sea level of all model levels at and below the formation height of the hydrometeor [m] |
|
|
|
|
| height | (1xm) | Same as z_out. |
|
|
|
|
| h | (1xn) | Relative humidity of ambient air [] |
|
|
|
|
| T | (1xn) | Temperature of ambient air [°C] |
|
|
|
|
| p | (1xn) | Pressure of ambient air [Pa] |
|
|
|
|
| exm | (1xn) | Excess moisture of ambient air (specific humidity divided by the saturation specific humidity over ice at ambient temperature; exm>=0) [kg/kg] |
|
|
|
|
| m | (1xm) | Hydrometeor mass [kg] |
|
|
|
|
| diam | (1xm) | Hydrometeor diameter [m] |
|
|
|
|
| Td | (1xm) | Hydrometeor temperature [°C] |
|
|
|
|
| Tdelta | (1xm) | Temperature difference between hydrometeor and surrounding air [°C] |
|
|
|
|
| equih | (1xm) | Relative humidity of the surrounding air with respect to Td, instead of T (equivalent relative humidity) [] |
|
|
|
|
| masstend | (1xm) | Hydrometeor mass tendency (rate of change of mass dm/dt) [kg/s] |
|
|
|
|
| temptend | (1xm) | Hydrometeor temperature tendency (rate of change of temperature dTd/dt) [°C/s] |
|
|
|
|
| f | (1xm) | Remaining fraction of mass of the hydrometeor relative to the next higher model level [] |
|
|
|
|
| ftot | (1xm) | Remaining fraction of mass of the hydrometeor relative to the initial mass (at z_form) [] |
|
|
|
|
| evapfrac | (1xm) | Evaporated fraction of mass (current mass relative to the initial mass of the hydrometeor; evapfrac=1-ftot) [] |
|
|
|
|
| deltaHv | (1xn) | delta2H 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] |
|
|
|
|
| deltaHp | (1xm) | delta2H value of the hydrometeor [permil] |
|
|
|
|
| dexcp | (1xm) |
|
|
|
|
d-excess value of the hydrometeor [permil]
|
|
|
|
deltaHp_eq
|
|
|
|
(1xm)
|
| ... | ... | |
| ... | ... | |