| name |
type |
meaning |
note |
| nx0 |
integer(4) |
grid number of x-direction (in full domain) |
grid number in full domain |
| ny0 |
integer(4) |
grid number of y-direction (in full domain) |
grid number in full domain |
| nz0 |
integer(4) |
grid number of z-direction (in full domain) |
grid number in full domain |
| nx |
integer(4) |
grid number of x-direction (in divided domain for each MPI processes) |
grid number in divided domain (includes halo-region) |
| ny |
integer(4) |
grid number of y-direction (in divided domain for each MPI processes) |
grid number in divided domain (includes halo-region) |
| nz |
integer(4) |
grid number of z-direction (in divided domain for each MPI processes) |
grid number in divided domain (includes halo-region). At present, same as nz0. |
| ngm |
integer(4) |
layer number used for calculating of ground temperature's forecast |
It usually set 4. |
| nx_hivi |
integer(4) |
grid number of x-direction (in divided domain) |
grid number of x-direction after changing line in the HI-VI
|
| ixstlc |
integer(4) |
start index of x-direction in divided domain (scalar point)
|
definded for the region excluding halo-region or physical boundary. it's usually 2.
|
| ixenlc |
integer(4) |
end index of x-direction in divided domain (scalar point)
|
definded for the region excluding halo-region or physical boundary. it's usually nx-1.
|
| jystlc |
integer(4) |
start index of y-direction in divided domain (scalar point)
|
definded for the region excluding halo-region or physical boundary. it's usually 2.
|
| jyenlc |
integer(4) |
end index of y-direction in divided domain (scalar point)
|
definded for the region excluding halo-region or physical boundary. it's usually ny-1.
|
| ixstlc_u |
integer(4) |
start index of x-direction in divided domain (vector(u) point)
|
definded for the region excluding halo-region or physical boundary. it's usually 2.
|
| ixenlc_u |
integer(4) |
end index of x-direction in divided domain (vector(u) point)
|
definded for the region excluding halo-region or physical boundary. it's usually nx.
|
| jystlc_v |
integer(4) |
start index of y-direction in divided domain (vector(v) point)
|
definded for the region excluding halo-region or physical boundary. it's usually 2.
|
| jyenlc |
integer(4) |
end index of y-direction in divided domain (vector(v) point)
|
definded for the region excluding halo-region or physical boundary. it's usually ny.
|
| ixstgb |
integer(4) |
start index of x-direction in full domain (scalar point)
|
definded for the region excluding physical boundary. it's usually 2.
|
| ixengb |
integer(4) |
end index of x-direction in full domain (scalar point)
|
definded for the region excluding physical boundary. it's usually nx-1.
|
| jystgb |
integer(4) |
start index of y-direction in full domain (scalar point)
|
definded for the region excluding physical boundary. it's usually 2.
|
| jyengb |
integer(4) |
end index of y-direction in full domain (scalar point)
|
definded for the region excluding physical boundary. it's usually ny-1.
|
| kzst |
integer(4) |
start index of z-direction in full domain (scalar point)
|
definded for the region excluding lowest boundary. it's usually 2.
|
| kzen |
integer(4) |
end index of z-direction in full domain (scalar point)
|
definded for the region excluding top boundary. it's usually nz-1.
|
| vdx(nx0) |
real(r_size) |
grid interval of x-direction (between vector points)
|
the index is defined in full-domain even in parallel execution case. |
| vdy(ny0) |
real(r_size) |
grid interval of y-direction (between vector points)
|
the index is defined in full-domain even in parallel execution case. |
| vdz(nz0) |
real(r_size) |
grid interval of z-direction (between half levels)
|
the index is defined in full-domain even in parallel execution case. |
| vdz2(nz0) |
real(r_size) |
grid interval of z-direction (between full levels)
|
the index is defined in full-domain even in parallel execution case. |
| vrdx(nx0) |
real(r_size) |
a reciprocal of grid interval of x-direction (between vector points)
|
= 1 / vdx. The index is defined in full-domain even in parallel execution case. |
| vrdy(ny0) |
real(r_size) |
a reciprocal of grid interval of y-direction (between vector points)
|
= 1 / vdy. The index is defined in full-domain even in parallel execution case. |
| vrdz(nz0) |
real(r_size) |
a reciprocal of grid interval of z-direction (between half levels)
|
= 1 / vdz. The index is defined in full-domain even in parallel execution case. |
| vrdx2(nx0) |
real(r_size) |
a reciprocal of grid interval of x-direction (between scalar points)
|
vdx2 is not variable. The index is defined in full-domain even in parallel execution case. |
| vrdy2(ny0) |
real(r_size) |
a reciprocal of grid interval of y-direction (between scalar points)
|
vdy2 is not variable. The index is defined in full-domain even in parallel execution case. |
| vrdz2(nz0) |
integer(4) |
real(r_size) |
= 1 / vdz2. The index is defined in full-domain even in parallel execution case. |
| zrp(nz0) |
real(r_size) |
height of full level |
|
| zrw(nz0) |
real(r_size) |
height of half level |
|
| dz_g(ngm) |
real(r_size) |
grid interval of layer in under ground
|
|
| name |
type |
meaning |
note |
| dt2 |
real(r_size) |
dt * 2 |
|
| swcmpt |
integer(4) |
switch parameter |
0: inelastic, 1: elastic(HI-VI), 2: elastic(HE-VI) |
| prcmpt |
integer(4) |
switch parameter |
0: inelastic, 1: elastic |
| swbuoy |
integer(4) |
switch parameter(buoyancy) |
0: calculate from perturbation of density, 1: calculate from perturbation of potential temperature |
| nsound |
integer(4) |
interval numbers of short time step |
|
| u(nx, ny, nz, 2) |
real(r_size) |
momentum of x-direction (U) |
= ρ * G^{1/2} * u / m_2 |
| v(nx, ny, nz, 2) |
real(r_size) |
momentum of y-direction (V) |
= ρ * G^{1/2} * v / m_1 |
| w(nx, ny, nz, 2) |
real(r_size) |
momentum of z-direction (W) |
= ρ * G^{1/2} * w / m_3 |
| pt(nx, ny, nz, 2) |
real(r_size) |
deviation of potential temperature (\theta) from ptrf |
pt = θ - ptrf |
| qv(nx, ny, nz, 2) |
real(r_size) |
mixing ratio of water vapor |
unit(kg/kg) |
| qc(nx, ny, nz, 2) |
real(r_size) |
mixing ratio of cloud water |
unit(kg/kg) |
| qr(nx, ny, nz, 2) |
real(r_size) |
mixing ratio of rain |
unit(kg/kg) |
| qci(nx, ny, nz, 2) |
real(r_size) |
mixing ratio of cloud ice |
unit(kg/kg) |
| qs(nx, ny, nz, 2) |
real(r_size) |
mixing ratio of snow |
unit(kg/kg) |
| qg(nx, ny, nz, 2) |
real(r_size) |
mixing ratio of graupel |
unit(kg/kg) |
| qh(nx, ny, nz, 2) |
real(r_size) |
mixing ratio of hail |
unit(kg/kg) |
| qncw(nx, ny, nz, 2) |
real(r_size) |
number concentration of cloud water |
unit(1 / m^3) |
| qnr(nx, ny, nz, 2) |
real(r_size) |
number concentration of rain |
unit(1 / m^3) |
| qnci(nx, ny, nz, 2) |
real(r_size) |
number concentration of cloud ice |
unit(1 / m^3) |
| qns(nx, ny, nz, 2) |
real(r_size) |
number concentration of snow |
unit(1 / m^3) |
| qng(nx, ny, nz, 2) |
real(r_size) |
number concentration of graupel |
unit(1 / m^3) |
| qnh(nx, ny, nz, 2) |
real(r_size) |
number concentration of hail |
unit(1 / m^3) |
| qcn(nx, ny, nz, 2) |
real(r_size) |
mixing ratio of condensation nucleus |
unit(kg/kg) |
| qncn(nx, ny, nz, 2) |
real(r_size) |
number concentration of condensation nucleus |
unit(1 / m^3) |
| dns(nx, ny, nz) |
real(r_size) |
air density |
|
| dnsg2(nx, ny, nz) |
real(r_size) |
dns * G^{1/2} |
|
| rdnsg2(nx, ny, nz) |
real(r_size) |
1 / (dns * G^{1/2}) |
|
| omw(nx, ny, nz, 2) |
real(r_size) |
momentum of z* -direction (W*) |
= 1 / G^{1/2} * (W + m_1 * m_2 / m_3 * (G^{1/2}G^{13} * U + G^{1/2}G^{23} * V)) |
| advu(nx, ny, nz) |
real(r_size) |
advection term of x-direction |
|
| advv(nx, ny, nz) |
real(r_size) |
advection term of y-direction |
|
| advw(nx, ny, nz) |
real(r_size) |
advection term of z-direction |
|
| ptm(nx, ny, nz, 2) |
real(r_size) |
mass virtual potential temperature |
(pt + ptrf) * (1. + 0.608 qv) * (1 - qc - qci - qr - qs - qh - qg) |
| pft(nx, ny, nz) |
real(r_size) |
heat expansion of air |
(ρ G^{1/2} / θ) * (∂θ_m / ∂t) |
| buoyd(nx, ny, nz) |
real(r_size) |
buoyancy |
|
| prs(nx, ny, nz) |
real(r_size) |
deviation of pressure |
prs = (prs_loc _ prsrft) * g2p |
| div(nx, ny, nz) |
real(r_size) |
divergence of z^* system |
= m_1 * m_2 * (∂U/∂x^ + ∂V/∂y^ ) + m_3 * ∂W^ /∂z^ |
| prc(nx, ny, nz) |
real(r_size) |
fall of precipitation |
= m_1 * m_2 * (∂U/∂x^ + ∂V/∂y^ ) + m_3 * ∂W^ /∂z^ |
| cs2(nx, ny, nz) |
real(r_size) |
square of the speed of sound |
= (C_p / C_v) * R_d * T_v
|
| cs2mn(nz) |
real(r_size) |
horizontal mean of cs2 |
|
| fgrtmn(nz) |
real(r_size) |
horizontal mean of g / m_3 / cs2 ^2 |
|
| prs_loc(nx, ny, nz) |
real(r_size) |
pressure |
= prsrft + prs * g2invp |
| pi_loc(nx, ny, nz) |
real(r_size) |
Exner function |
= (prs_loc / presrf) ** rdvcp |
| dudtbc(ny, nz, 2) |
real(r_size) |
production term of u on the lateral boundary of the y-direction |
|
| rstptw(nx, ny, 2) |
real(r_size) |
the earth's surface flux by turbulence |
rstptw(:, :, 1) : θ, rstptw(:, :, 1) : u |
| rstqvw(nx, ny, 2) |
real(r_size) |
the earth's surface flux by turbulence |
rstqvw(:, :, 1) : qv, rstqvw(:, :, 1) : v |
| radpt(nx, ny, nz) |
real(r_size) |
production term of potential temperature by radiation |
unit(K/s) |
| solar(nx, ny) |
real(r_size) |
Net downward shortwave radiation flux at surface(dswb - uswb) |
unit(W/m^2) |
| dswb(nx, ny) |
real(r_size) |
Downward shortwave radiation flux at surface |
unit(W/m^2) |
| dlwb(nx, ny) |
real(r_size) |
Downward longwave radiation flux at surface |
unit(W/m^2) |
| zmean(nx, ny) |
real(r_size) |
Cosine of solar zenith angle |
|
| tradpt(nx, ny, nz) |
real(r_size) |
Time integrated radpt from initial |
unit(K), for MRI format |
| tsolar(nx, ny) |
real(r_size) |
Time integrated solar from initial |
unit(J/m^2), for MRI format |
| tdswb(nx, ny) |
real(r_size) |
Time integrated dswb from initial |
unit(J/m^2), for MRI format |
| tdlwb(nx, ny) |
real(r_size) |
Time integrated dlwb from initial |
unit(J/m^2), for MRI format |
| tin(nx, ny, ngm) |
real(r_size) |
underground temperature |
tin(:, :, 1) is the 1st layer of underground.
As the dimension is bigger, it becomes deeper. |
| ptgrd(nx, ny) |
real(r_size) |
ground potential temperature |
|
| qvgrd(nx, ny) |
real(r_size) |
ground qv |
= wet * qsatg + (1. - wet) * qv(:, :, 2, kt) |
| wet(nx, ny) |
real(r_size) |
wetness (the efficient of evapotranspiration on the surface) |
|
| w_g(nx, ny, 2) |
real(r_size) |
volume water content |
w_g(:,:,1) = volume water content at the 1st ground layer, w_g(:,:,2) = mean volume water content from surface to 50cm in depth |
| fktg(nx, ny, ngm) |
real(r_size) |
heat diffusion's coefficient of ground or underground |
|
| roctg(nx, ny) |
real(r_size) |
ground heat capacity |
|
| albed(nx, ny) |
real(r_size) |
albedo |
|
| urban(nx, ny) |
real(r_size) |
urban heat |
|
| qsatg(nx, ny) |
real(r_size) |
ground saturation water vapor volume |
|
| snheat(nx, ny) |
real(r_size) |
sensible flux |
unit(W/m^2) |
| evheat(nx, ny) |
real(r_size) |
latent heat flux |
unit(W/m^2) |
| skind(nx, ny) |
integer(4) |
classification on the land's surface |
(1: land, 2: sea, 3: snow, 4: ice) |
| ssti(nx, ny) |
real(r_size) |
sea surface temperature |
|
| smqr(nx, ny) |
real(r_size) |
accumulated precipitation of rain |
unit(mm) |
| smqs(nx, ny) |
real(r_size) |
accumulated precipitation of snow |
unit(mm) |
| smqi(nx, ny) |
real(r_size) |
accumulated precipitation of (cloud) ice |
unit(mm) |
| smqg(nx, ny) |
real(r_size) |
accumulated precipitation of graupel |
unit(mm) |
| smqh(nx, ny) |
real(r_size) |
accumulated precipitation of hail |
unit(mm) |
| ppt(nx, ny, nz) |
real(r_size) |
production term of potential temperature |
|
| pqv(nx, ny, nz) |
real(r_size) |
production term of qv |
|
| pqr(nx, ny, nz) |
real(r_size) |
production term of qr |
|
| pqs(nx, ny, nz) |
real(r_size) |
production term of qs |
|
| pqg(nx, ny, nz) |
real(r_size) |
production term of qg |
|
| pqh(nx, ny, nz) |
real(r_size) |
production term of qh |
|
| pqci(nx, ny, nz) |
real(r_size) |
production term of qci |
|
| pqcw(nx, ny, nz) |
real(r_size) |
production term of qc |
|
| uf(nx, ny, nz) |
real(r_size) |
momentum of x-direction (U) |
use during short time step |
| vf(nx, ny, nz) |
real(r_size) |
momentum of y-direction (V) |
use during short time step |
| wf(nx, ny, nz) |
real(r_size) |
momentum of z-direction (W) |
use during short time step |
| ptf(nx, ny, nz, nn) |
real(r_size) |
potential temperature (θ) |
Use during short time step.
In case of integrating potential temperature with leap-frog during short time step, nn is 2.
The other case, nn is 1.
|
| omwf(nx, ny, nz) |
real(r_size) |
momentum of z^*-direction (W^*) |
use during short time step |
| prsf(nx, ny, nz) |
real(r_size) |
perturbation of pressure |
use during short time step |
| dts |
real(r_size) |
size of short time step |
|
| beta |
real(r_size) |
implicit rate |
It's usually used 1.0 (backward difference). |
| dtb |
real(r_size) |
dts * β |
|
| udiag_vc(nx, ny, nz, 2) |
real(r_size) |
wind (u) of x-direction on a vector point |
udiag_vc = u * m_2 / dnsg2 |
| vdiag_vc(nx, ny, nz, 2) |
real(r_size) |
wind (v) of y-direction on a vector point |
vdiag_vc = v * m_1 / dnsg2 |
| wdiag_vc(nx, ny, nz, 2) |
real(r_size) |
wind (w) of z-direction on a half level |
wdiag_vc = w * m_3 / dnsg2 |
| udiag_sc(nx, ny, nz, 2) |
real(r_size) |
wind (u) of x-direction on a scalar point |
udiag_vc = u * m_2 / dnsg2 |
| vdiag_vc(nx, ny, nz, 2) |
real(r_size) |
wind (v) of y-direction on a scalar point |
vdiag_vc = v * m_1 / dnsg2 |
| wdiag_vc(nx, ny, nz, 2) |
real(r_size) |
wind (w) of z-direction on a full level |
wdiag_vc = w * m_3 / dnsg2 |
| z_0m(nx, ny) |
real(r_size) |
momentum roughness index |
|
| z_0h(nx, ny) |
real(r_size) |
heat roughness index |
|
| z_0q(nx, ny) |
real(r_size) |
water vapor roughness index |
|
| z_0q(nx, ny) |
real(r_size) |
water vapor roughness index |
|
| sf_u_w(nx, ny) |
real(r_size) |
|
|
| sf_w_u(nx, ny) |
real(r_size) |
|
|
| sf_pt_w(nx, ny) |
real(r_size) |
|
|
| sf_qv_w(nx, ny) |
real(r_size) |
|
|
| c_m(nx, ny) |
real(r_size) |
|
|
| c_h(nx, ny) |
real(r_size) |
|
|
| c_q(nx, ny) |
real(r_size) |
|
|
| v_abs_a(nx, ny) |
real(r_size) |
|
|
| v_abs_a_u(nx, ny) |
real(r_size) |
|
|
| v_abs_a_v(nx, ny) |
real(r_size) |
|
|
| l_mo_inv(nx, ny) |
real(r_size) |
|
|
| ri_b(nx, ny) |
real(r_size) |
|
|
| eturb(nx, ny, nz, 2) |
real(r_size) |
turbulent energy |
|
| h_pbl(nx, ny) |
real(r_size) |
height of boundary layer |
|
| l_mix_q(nx, ny) |
real(r_size) |
mixed length |
|
| eddykm_x(nx, ny, nz, 2) |
real(r_size) |
momentum diffusion's coefficient of x-direction |
|
| eddykm_y(nx, ny, nz, 2) |
real(r_size) |
momentum diffusion's coefficient of y-direction |
|
| eddykm_z(nx, ny, nz, 2) |
real(r_size) |
momentum diffusion's coefficient of y-direction |
|
| eddykh_x(nx, ny, nz, 2) |
real(r_size) |
heat diffusion's coefficient of x-direction |
|
| eddykh_y(nx, ny, nz, 2) |
real(r_size) |
heat diffusion's coefficient of y-direction |
|
| eddykh_z(nx, ny, nz, 2) |
real(r_size) |
heat diffusion's coefficient of y-direction |
|
| l_mix_x(nx, ny, nz) |
real(r_size) |
|
|
| l_mix_y(nx, ny, nz) |
real(r_size) |
|
|
| l_mix_z(nx, ny, nz) |
real(r_size) |
|
|
| prandtl_t_x(nx, ny, nz) |
real(r_size) |
|
|
| prandtl_t_y(nx, ny, nz) |
real(r_size) |
|
|
| prandtl_t_z(nx, ny, nz) |
real(r_size) |
|
|
| qke(nx, ny, nz, 2) |
real(r_size) |
eturb * 2.0 |
|
| tsq(nx, ny, nz, 2) |
real(r_size) |
(θ') ** 2 |
|
| qsq(nx, ny, nz, 2) |
real(r_size) |
(qw') ** 2 |
|
| cov(nx, ny, nz, 2) |
real(r_size) |
(θ') * (qw') |
|
| cu(nx, ny, nz, 2) |
real(r_size) |
|
|
| cv(nx, ny, nz, 2) |
real(r_size) |
|
|
| qm(nx, ny, nz, 2) |
real(r_size) |
a part of the degree of cloudiness |
variable has the range between 0.0 and 1.0. |
| name |
type |
meaning |
note |
| rshrt(nx, ny, nz) |
real(r_size) |
PT time change rate by shortwave radiation heating(before solar zenith angle correction) |
unit(K/s) |
| rlong(nx, ny, nz) |
real(r_size) |
PT time change rate by longwave raiation heating(before surface skin temperature correction) |
unit(K/s) |
| coefc(nx, ny, nz + 1) |
real(r_size) |
Transmissivity estimated from surface |
no dim |
| ttrans(nx, ny, nz + 1) |
real(r_size) |
Transmissivity estimated from surface |
no dim |
| t1sav(nx, ny) |
real(r_size) |
Longwave radiation flux emitted from bottom of the atmosphere |
unit(W/m^2) |
| dswbm(nx, ny) |
real(r_size) |
Downward shortwave radiation flux at surface(before solar zenith angle correction) |
unit(W/m^2) |
| uswbm(nx, ny) |
real(r_size) |
Upward shortwave radiation flux at surface(before solar zenith angle correction) |
unit(W/m^2) |
| ozon_3dim(nx, ny, nz) |
real(r_size) |
Monthly climate of 3D ozone amount |
unit(cm/STP) |
| aod_2d(nx, ny) |
real(r_size) |
Monthly climate of aerosol optical depth(vertical integrated) |
no dim |
| dswtm(nx, ny) |
real(r_size) |
Downward shortwave radiation flux at TOA(before solar zenith angle correction) |
unit(W/m^2) |
| uswtm(nx, ny) |
real(r_size) |
Upwnward shortwave radiation flux at TOA(before solar zenith angle correction) |
unit(W/m^2) |
| ulwtm(nx, ny) |
real(r_size) |
Upwnward longwave radiation flux at TOA(before surface skin temperature correction) |
unit(W/m^2) |
| rvisbm(nx, ny) |
real(r_size) |
Downward direct radiation flux at surface(uv and visible, before solar zenith angle correction) |
unit(W/m^2) |
| rvisdm(nx, ny) |
real(r_size) |
Downward diffuse radiation flux at surface(uv and visible, before solar zenith angle correction) |
unit(W/m^2) |
| rnirbm(nx, ny) |
real(r_size) |
Downward direct radiation flux at surface(near-infrared, before solar zenith angle correction) |
unit(W/m^2) |
| rnirdm(nx, ny) |
real(r_size) |
Downward diffuse radiation flux at surface(near-infrared, before solar zenith angle correction) |
unit(W/m^2) |
| dswbcsm (nx, ny) |
real(r_size) |
Downward shortwave radiation flux at surface(clear sky, before solar zenith angle correction) |
unit(W/m^2) |
| uswbcsm (nx, ny) |
real(r_size) |
Upwward shortwave radiation flux at surface(clear sky, before solar zenith angle correction) |
unit(W/m^2) |
| dlwbcs (nx, ny) |
real(r_size) |
Upward longwave radiation flux at surface(clear sky) |
unit(W/m^2) |
| uswtcsm (nx, ny) |
real(r_size) |
Upward shortwave radiation flux at TOA(clear sky, before solar zenith angle correction) |
unit(W/m^2) |
| ulwtcsm (nx, ny) |
real(r_size) |
Upward longwave radiation flux at TOA(clear sky, before surface skin temperature correction) |
unit(W/m^2) |
| coefccs (nx, ny) |
real(r_size) |
Transmissivity estimated from surface(clear sky) |
no dim |
| ttranscs(nx, ny) |
real(r_size) |
Transmissiviy estimated from surface(clear sky) |
no dim |
| cllrad(nx, ny) |
real(r_size) |
Low level cloud fraction in radiation calculation |
no dim |
| clmrad(nx, ny) |
real(r_size) |
Middle level cloud fraction in radiation calculation |
no dim |
| clhrad(nx, ny) |
real(r_size) |
High level cloud fraction in radiation calculation |
no dim |
| clarad(nx, ny) |
real(r_size) |
Total cloud fraction in radiation calculationt |
no dim |
| name |
type |
meaning |
note |
| qvrf1(nz) |
real(r_size) |
horizontal mean of qv (full level)
|
mean in the level. |
| qvsrf1(nz) |
real(r_size) |
horizontal mean of qvsat(saturation water vapor volume) (full level)
|
mean in the level. |
| ptrf1(nz) |
real(r_size) |
horizontal mean of pt (full level)
|
mean in the level. |
| vptrf1(nz) |
real(r_size) |
horizontal mean of imaginary potential temperature(Į_v) (full level)
|
vptrf1(kz) = ptrf1(kz) * (1.0 + 0.608 * qvrf1(kz)), mean in the level. |
| vptrf2(nz) |
real(r_size) |
horizontal mean of imaginary potential temperature(θ_v) (half level)
|
vptrf1(kz) = 0.5 * (vptrf2(kz) + vptrf1(kz+1))
|
| pairf1(nz) |
real(r_size) |
horizontal mean of Exner function(Π) (full level)
|
get by piling vptrf2, mean in the level. |
| pairf2(nz) |
real(r_size) |
horizontal mean of Exner function(Π) (half level)
|
get by piling vptrf1, mean in the level. |
| prsrf1(nz) |
real(r_size) |
horizontal mean of pressure (full level)
|
get by pairf1, mean in the level. |
| prsrf2(nz) |
real(r_size) |
horizontal mean of pressure (half level)
|
get by pairf2, mean in the level. |
| trf2(nz) |
real(r_size) |
horizontal mean of T (half level)
|
trf1(kz) = pairf1(kz) * ptrf1(kz), mean in the level. |
| dnsrf1(nz) |
real(r_size) |
horizontal mean of ρ (full level)
|
mean in the level. |
| dnsrf2(nz) |
real(r_size) |
horizontal mean of ρ (half level)
|
mean in the level. |
| ptrft(nx, ny, nz) |
real(r_size) |
basic field of potential temperature(θ)
|
calculated from ptrf1 |
| vptrft(nx, ny, nz) |
real(r_size) |
basic field of imaginary potential temperature(θ)
|
calculated from vptrf1 |
| pairft(nx, ny, nz) |
real(r_size) |
basic field of Exner function(Π)
|
calculated from pairf1 |
| prsrft(nx, ny, nz) |
real(r_size) |
basic field of pressure(p)
|
get by calculated pairft (It's changed from the old one to the new.) |
| dnsrft(nx, ny, nz) |
real(r_size) |
basic field of density(ρ) |
not use in the case of fully compressible system |