Home

Context Navigation

runtime_parameters

Timestamp:: Oct 25, 2018 12:49:49 PM (7 years ago)
Author:: gronemeier
Comment:: --

Legend:

: Unmodified
: Added
: Removed
: Modified

doc/app/runtime_parameters

-                      v290
+                      v291
 The default assignment is:\\\\
 '' ' u v ' ''\\
 '' ' pt ', ''\\
 '' ' w"pt" w*pt* w*pt*BC wpt wptBC ', ''\\
+'' ' theta ', ''\\
+'' ' w"theta" w*theta* w*theta*BC wtheta wthetaBC ', ''\\
 '' ' w"u" w*u* wu w"v"w*v* wv ', ''\\
 '' ' km kh ', ''\\
 …
 If output of vertical profiles is produced (see  [#data_output_pr data_output_pr]) the appropriate data are written to a NetCDF file. Simultaneously, the model produces an attribute in the header of the NetCDF file which describes the layout for a plot to be generated with the plot software '''palmplot'''. The parameter '''cross_profiles''' determines how many coordinate systems (panels) the plot contains and which profiles are supposed to be drawn into which panel. (Currently, '''palmplot''' is limited to three profiles per panel.) '''cross_profiles''' expects a character string (up to 100 characters long) for each coordinate system, which consists of the names of the profiles to be drawn into this system (all available profiles and their respective names are described at parameter [#data_output_pr data_output_pr]). The single names have to be separated by one blank (' ') and a blank must be spent also at the beginning and at the end of the string.\\\\
 Example:
    '''cross_profiles''' = '' ' u v ', ' pt ' ''
 In this case, the plot consists of two coordinate systems (panels) with the first panel containing the profiles of the horizontal velocity components ('' 'u' '' and '' 'v' '') of all output times (see [#dt_dopr dt_dopr]) and the second one containing the profiles of the potential temperature ('' 'pt' '').\\\\
+   '''cross_profiles''' = '' ' u v ', ' theta ' ''
+In this case, the plot consists of two coordinate systems (panels) with the first panel containing the profiles of the horizontal velocity components ('' 'u' '' and '' 'v' '') of all output times (see [#dt_dopr dt_dopr]) and the second one containing the profiles of the potential temperature ('' 'theta' '').\\\\
 Whether the coordinate systems are actually drawn, depends on whether data of the appropriate profiles were output during the run (profiles to be output have to be selected with the parameter [#data_output_pr data_output_pr]). For example, if '''data_output_pr''' = '' 'u', 'v' '' was assigned, then the plot only consists of one panel, since no profiles of the potential temperature were output. On the other hand, if profiles were assigned to '''data_output_pr''' whose names do not appear in '''cross_profiles''', this profiles will be plotted separately behind the profiles defined in [#cross_profiles cross_profiles].\\\\
 The arrangement of the panels in the plot can be controlled with the parameters [#profile_columns profile_columns] and [#profile_rows profile_rows]. Up to 100 panels systems are allowed in a plot (however, they may be distributed on several pages).
 …
 ||='''Quantity name''' =||='''Meaning''' =||='''Unit''' =||='''Remarks''' =||
 ||e  ||SGS turbulence kinetic energy||m^2^/s^2^ || ||
 ||lpt ||liquid water potential temperature ||K ||requires [[../inipar#cloud_physics|cloud_physics]] = ''.T.''||
+||thetal ||liquid water potential temperature ||K ||requires [[../inipar#cloud_physics|cloud_physics]] = ''.T.''||
 ||lwp* ||liquid water path ||kg/m^2^ ||only horizontal cross section is allowed, requires [[../inipar#cloud_physics|cloud_physics]] = ''.T.'' ||
 ||nr ||rain drop number density ||1/m^3^ ||requires [[../inipar#cloud_physics|cloud_physics]] = ''.T.'' and [[../inipar#cloud_scheme|cloud_scheme]] = ''seifert_beheng'' and [[../inipar#precipitation|precipitation]] = ''.T.''||
 …
 ||prr* ||precipitation rate ||mm/s ||only horizontal cross section is allowed, requires [[../inipar#cloud_scheme|cloud_scheme]] = ''kessler'' and [[../inipar#precipitation|precipitation]] = ''.T.'' ||
 ||prr ||precipitation rate ||kg/kg m/s ||requires [[../inipar#cloud_scheme|cloud_scheme]] = ''seifert_beheng'' and [[../inipar#precipitation|precipitation]] = ''.T.''||
 ||pt ||potential temperature ||K || ||
+||theta ||potential temperature ||K || ||
 ||q ||water vapor mixing ratio (or total water mixing ratio if cloud physics is switched on) ||kg/kg ||requires [[../inipar#humidity|humidity]] = ''.T.'' ||
 ||qc ||cloud water mixing ratio ||kg/kg ||requires [[../inipar#cloud_physics|cloud_physics]] = ''.T.'' and [[../inipar#cloud_scheme|cloud_scheme]] = ''seifert_beheng'' ||
 …
 ||t* ||(near surface) characteristic temperature ||K ||only horizontal cross section is allowed ||
 ||u ||u-component of the velocity ||m/s || ||
 ||u* ||(near surface) friction velocity ||m/s ||only horizontal cross section is allowed ||
+||us* ||(near surface) friction velocity ||m/s ||only horizontal cross section is allowed ||
 ||v ||v-component of the velocity ||m/s || ||
 ||vpt ||virtual potential temperature ||K ||requires [[../inipar#humidity|humidity]] = ''.T.'' ||
+||thetav ||virtual potential temperature ||K ||requires [[../inipar#humidity|humidity]] = ''.T.'' ||
 ||w ||w-component of the velocity ||m/s || ||
 ||z0* ||roughness length ||m || ||
 …
 By assigning the pure strings from the above table, 3d volume data is output. Cross section data can be output by appending the string '_xy', '_xz', or '_yz' to the respective quantities. Time-averaged output is created by appending the string '_av' (for cross section data, this string must be appended to the cross section string). Cross section data can also be (additionally) averaged along the direction normal to the respective section (see below). Assignments of quantities can be given in arbitrary order:\\\\
 Example:
    '''data_output''' = '' 'u', 'pt_xz_av', 'w_xy', 'u_av' ''.
+   '''data_output''' = '' 'u', 'theta_xz_av', 'w_xy', 'u_av' ''.
 This example will create the following output: instantaneous 3d volume data of u-velocity component (by default on file [[../iofiles#DATA_3D_NETCDF|DATA_3D_NETCDF]]), temporally averaged 3d volume data of u-velocity component (by default on file [[../iofiles#DATA_3D_AV_NETCDF|DATA_3D_AV_NETCDF]]), instantaneous horizontal cross section data of w-velocity component (by default on file [[../iofiles#DATA_2D_XY_NETCDF|DATA_2D_XY_NETCDF]]), and temporally averaged vertical cross section data of potential temperature (by default on file [[../iofiles#DATA_2D_XZ_AV_NETCDF|DATA_2D_XZ_AV_NETCDF]]).\\\\
 The user is allowed to extend the above list of quantities by defining his own output quantities (see the user-parameter [[../userpar#data_output_user|data_output_user]]).\\\\
 …
 ||='''Quantity name''' =||='''Meaning''' =||='''Unit''' =||='''Remarks''' =||
 ||e  ||SGS turbulence kinetic energy|| m^2^/s^2^ || ||
 ||lpt ||liquid water potential temperature ||K ||requires [[../inipar#cloud_physics|cloud_physics]] = ''.T.''||
+||thetal ||liquid water potential temperature ||K ||requires [[../inipar#cloud_physics|cloud_physics]] = ''.T.''||
 ||nr ||rain drop number density ||1/m^3^ ||requires [[../inipar#cloud_physics|cloud_physics]] = ''.T.'' and [[../inipar#cloud_scheme|cloud_scheme]] = ''seifert_beheng'' and [[../inipar#precipitation|precipitation]] = ''.T.''||
 ||p ||perturpation pressure ||N/m^2^, Pa || ||
 ||pc ||particle/droplet concentration ||#/gridbox || ||
 ||pr ||mean particle/droplet radius ||m || ||
 ||pt ||potential temperature ||K || ||
+||theta ||potential temperature ||K || ||
 ||q ||water vapor mixing ratio (or total water mixing ratio if cloud physics is switched on) ||kg/kg ||requires [[../inipar#humidity|humidity]] = ''.T.'' ||
 ||qc ||cloud water mixing ratio ||kg/kg ||requires [[../inipar#cloud_physics|cloud_physics]] = ''.T.'' and [[../inipar#cloud_scheme|cloud_scheme]] = ''seifert_beheng'' ||
 …
 ||u ||u-component of the velocity ||m/s || ||
 ||v ||v-component of the velocity ||m/s || ||
 ||vpt ||virtual potential temperature ||K ||requires [[../inipar#humidity|humidity]] = ''.T.'' ||
+||thetav ||virtual potential temperature ||K ||requires [[../inipar#humidity|humidity]] = ''.T.'' ||
 ||w ||w-component of the velocity ||m/s || ||
 '''Example:'''
   '''data_output_masks(1,:)''' = '' 'u', 'v', 'u_av', 'v_av', ''\\
   '''data_output_masks(2,:)''' = '' 'w', 'w_av', ''\\
   '''data_output_masks(3,:)''' = '' 'pt', 'u', 'v', 'w', ''\\\\
+  '''data_output_masks(3,:)''' = '' 'theta', 'u', 'v', 'w', ''\\\\
   This example will create the following output:\\
   '''Mask 1''': Instantaneous data of u- and v-velocity component (by default on file [[../iofiles#DATA_MASK_01_NETCDF|DATA_MASK_01_NETCDF]]) and temporally averaged
 …
 ||[[span(w ,style=color: green)]] ||w-component of the velocity ||m/s ||
 ||[[span(w_subs ,style=color: red)]] ||large-scale vertical velocity ||m/s ||
 ||[[span(pt ,style=color: red)]] ||Potential temperature ||K ||
 ||[[span(vpt ,style=color: red)]] ||Virtual potential temperature ||K ||
 ||[[span(lpt ,style=color: red)]] ||Liquid water potential temperature ||K ||
+||[[span(theta ,style=color: red)]] ||Potential temperature ||K ||
+||[[span(thetav ,style=color: red)]] ||Virtual potential temperature ||K ||
+||[[span(thetal ,style=color: red)]] ||Liquid water potential temperature ||K ||
 ||[[span(q ,style=color: red)]] ||Total water mixing ratio ||kg/kg ||
 ||[[span(qv ,style=color: red)]] ||Water vapor mixing ratio ||kg/kg ||
 ||[[span(ql ,style=color: red)]] ||Liquid water mixing ratio ||kg/kg ||
 ||[[span(hyp ,style=color: red)]] ||Hydrostatic Pressure ||hPa ||
 ||[[span(rho_air ,style=color: red)]] ||Air Density ||kg/m^3^ ||
 ||[[span(rho_ocean ,style=color: red)]] ||Ocean Density ||kg/m^3^ ||
+||[[span(rho ,style=color: red)]] ||Air Density ||kg/m^3^ ||
+||[[span(rho_sea_water ,style=color: red)]] ||Ocean Density ||kg/m^3^ ||
 ||[[span(prho ,style=color: red)]] ||Potential density ||kg/m^3^ ||
 ||[[span(s ,style=color: red)]] ||Scalar concentration (requires [[../inipar#passive_scalar|passive_scalar]] = ''.T.'')  ||kg m^-3^ or ppm ||
 …
 ||[[span(w*v* ,style=color: green)]] ||Covariance or v-component of the resolved vertical momentum flux (in case of [../inipar#momentum_advec momentum_advec] = 'ws-scheme', see [../../tec/advection#statistical_evaluation details]) ||m^2^/s^2^ ||
 ||[[span(wv ,style=color: green)]] ||v-component of the total vertical momentum flux (w"v" + w*v*) ||m^2^/s^2^ ||
 ||[[span(w"pt" ,style=color: green)]] ||Subgrid-scale vertical sensible heat flux ||K m/s ||
 ||[[span(w*pt* ,style=color: green)]] ||Covariance or resolved vertical sensible heat flux (in case of [../inipar#scalar_advec scalar_advec] = 'ws-scheme', see [../../tec/advection#statistical_evaluation details]) ||K m/s ||
 ||[[span(wpt ,style=color: green)]] ||Total vertical sensible heat flux (w"pt" + w*pt*) ||K m/s ||
 ||[[span(w*pt*BC ,style=color: green)]] ||Subgrid-scale vertical sensible heat flux using the Bott-Chlond scheme ||K m/s ||
 ||[[span(wptBC ,style=color: green)]] ||Total vertical sensible heat flux using the Bott-Chlond scheme (w"pt" + w*pt*BC) ||K m/s ||
 ||[[span(w"vpt" ,style=color: green)]] ||Subgrid-scale vertical buoyancy flux ||K m/s ||
 ||[[span(w*vpt* ,style=color: green)]] ||Resolved vertical buoyancy flux ||K m/s ||
 ||[[span(wvpt ,style=color: green)]] ||Total vertical buoyancy flux (w"vpt" + w*vpt*) ||K m/s ||
+||[[span(w"theta" ,style=color: green)]] ||Subgrid-scale vertical sensible heat flux ||K m/s ||
+||[[span(w*theta* ,style=color: green)]] ||Covariance or resolved vertical sensible heat flux (in case of [../inipar#scalar_advec scalar_advec] = 'ws-scheme', see [../../tec/advection#statistical_evaluation details]) ||K m/s ||
+||[[span(wtheta ,style=color: green)]] ||Total vertical sensible heat flux (w"theta" + w*theta*) ||K m/s ||
+||[[span(w*theta*BC ,style=color: green)]] ||Subgrid-scale vertical sensible heat flux using the Bott-Chlond scheme ||K m/s ||
+||[[span(wthetaBC ,style=color: green)]] ||Total vertical sensible heat flux using the Bott-Chlond scheme (w"theta" + w*theta*BC) ||K m/s ||
+||[[span(w"thetav" ,style=color: green)]] ||Subgrid-scale vertical buoyancy flux ||K m/s ||
+||[[span(w*thetav* ,style=color: green)]] ||Resolved vertical buoyancy flux ||K m/s ||
+||[[span(wthetav ,style=color: green)]] ||Total vertical buoyancy flux (w"thetav" + w*thetav*) ||K m/s ||
 ||[[span(w"q" ,style=color: green)]] ||Subgrid-scale vertical water flux ||kg/kg m/s ||
 ||[[span(w*q* ,style=color: green)]] ||Covariance or resolved vertical water flux (in case of [../inipar#scalar_advec scalar_advec] = 'ws-scheme', see [../../tec/advection#statistical_evaluation details]) ||kg/kg m/s ||
 …
 ||[[span(v*2 ,style=color: red)]] ||Variance or horizontal momentum flux (in case of [../inipar#momentum_advec momentum_advec] = 'ws-scheme', see [../../tec/advection#statistical_evaluation details]) of the v-velocity component (resolved) ||m^2^/s^2^ ||
 ||[[span(w*2 ,style=color: green)]] ||Variance or vertical momentum flux (in case of [../inipar#momentum_advec momentum_advec] = 'ws-scheme', see [../../tec/advection#statistical_evaluation details]) of the w-velocity component (resolved) ||m^2^/s^2^ ||
 ||[[span(pt*2 ,style=color: red)]] ||Variance of the potential temperature (resolved) ||K^2^ ||
+||[[span(theta*2 ,style=color: red)]] ||Variance of the potential temperature (resolved) ||K^2^ ||
 ||[[span(q*2 ,style=color: red)]] ||Variance of the water vapor / total water mixing ratio (in case of  [../inipar#cloud_physics cloud_physics] = .T.) (resolved) ||kg^2^/kg^2^ ||
 ||[[span(s*2 ,style=color: red)]] ||Variance of the passive scalar (resolved, requires [../inipar#passive_scalar passive_scalar] = .T.) ||(kg m^-3^)^2^ ||
 ||[[span(w*3 ,style=color: green)]] ||Third moment of the w-velocity component (resolved) ||m^3^/s^3^ ||
 ||[[span(Sw ,style=color: green)]] ||Skewness of the w-velocity component (resolved, Sw = w^3^/(w^2^)^1.5^) ||m^3^/s^2^ / (m^2^/s^2^)^1.5^ ||
 ||[[span(w*2pt* ,style=color: green)]] ||Third moment (resolved) ||K m^2^/s^2^ ||
 ||[[span(w*pt*2 ,style=color: green)]] ||Third moment (resolved) ||K^2^ m/s ||
+||[[span(w*2theta* ,style=color: green)]] ||Third moment (resolved) ||K m^2^/s^2^ ||
+||[[span(w*theta*2 ,style=color: green)]] ||Third moment (resolved) ||K^2^ m/s ||
 ||[[span(w*u*u*:dz ,style=color: red)]] ||Transport of resolved-scale TKE due to turbulence (term in resolved TKE budget)||m^2^/s^3^ ||
 ||[[span(w*p*:dz ,style=color: red)]] ||Transport of resolved-scale TKE due to pressure fluctuations (term in resolved TKE budget)||Pa m/s^2^ ||
 ||[[span(hyp ,style=color: red)]] ||Hydrostatic pressure ||dbar ||
 ||[[span(td_lsa_lpt ,style=color: red)]] ||horizontal large scale advection tendency for temperature||K/s ||
+||[[span(td_lsa_thetal ,style=color: red)]] ||horizontal large scale advection tendency for temperature||K/s ||
 ||[[span(td_lsa_q ,style=color: red)]] ||horizontal large scale advection tendency for humidity||kg/(kg s) ||
 ||[[span(td_sub_lpt ,style=color: red)]] ||horizontal large scale subsidence tendency for temperature||K/s ||
+||[[span(td_sub_thetal ,style=color: red)]] ||horizontal large scale subsidence tendency for temperature||K/s ||
 ||[[span(td_sub_q ,style=color: red)]] ||horizontal large scale subsidence tendency for humidity||kg/(kg s) ||
 ||[[span(td_nud_lpt ,style=color: red)]] ||nudging tendency for temperature||K/s ||
+||[[span(td_nud_thetal ,style=color: red)]] ||nudging tendency for temperature||K/s ||
 ||[[span(td_nud_q ,style=color: red)]] ||nudging tendency for humidity||kg/(kg s) ||
 ||[[span(td_nud_u ,style=color: red)]] ||nudging tendency for u-velocity component||m/s^2^ ||
 …
 \\\\
 Beyond that, initial profiles (t=0) of some variables can additionally be output (this output is only done once with the first plot output and not repeated with the profile output at later times). The names of these profiles result from the ones specified above led by a hash "#".  Allowed values are:
    #u, #v, #pt, #km, #kh, #l, #lpt, #q, #qv, #s, #sa, #vpt (, #t_soil, #m_soil)
+   #u, #v, #theta, #km, #kh, #l, #thetal, #q, #qv, #s, #sa, #thetav (, #t_soil, #m_soil)
 Profile names preceded by a hash automatically imply that profiles for these variables are also output at later times. It is not necessary and not allowed to specify the same profile name with and without hash simultaneously(this would lead to a netCDF error).\\\\
 These initial profiles have been either set by the user or have been calculated by a 1d-model prerun.\\\\
 …
 ||E* ||Resolved-scale turbulence kinetic energy of the flow (3D domain average) ||m^2^/s^2^ ||
 ||dt ||Time step size ||s ||
 ||u* ||Friction velocity (horizontal average) ||m/s ||
+||us* ||Friction velocity (horizontal average) ||m/s ||
 ||w* ||Vertical velocity scale of the CBL (horizontal average) ||m/s ||
 ||th* ||Temperature scale (constant flux layer), defined as w"pt"0 / u* (horizontal average)||K ||
 ||q* ||Humidity scale (constant flux layer), defined as w"q"0 / u* (horizontal average)||kg/kg ||
+||th* ||Temperature scale (constant flux layer), defined as w"theta"0 / us* (horizontal average)||K ||
+||q* ||Humidity scale (constant flux layer), defined as w"q"0 / us* (horizontal average)||kg/kg ||
 ||umax ||Maximum u-component of the velocity ||m/s ||
 ||vmax ||Maximum v-component of the velocity ||m/s ||
 …
 ||div_old ||Divergence of the velocity field before the pressure solver has been called (normalized with respect to the total number of grid points) ||1/s ||
 ||div_new ||Divergence of the velocity field after the pressure solver has been called (normalized with respect to the total number of grid points) ||1/s ||
 ||z_i_wpt ||Height of the convective boundary layer (horizontal average) determined by the height of the minimum sensible heat flux ||m ||
 ||z_i_pt ||Height of the convective boundary layer (horizontal average) determined by the temperature profile, following the criterion of Sullivan et al. (1998) ||m ||
 ||w"pt"0 ||Subgrid-scale sensible heat flux at k=0 (horizontal average) within the constant flux layer ||K m/s ||
 ||w"pt" ||Subgrid-scale heat flux (horizontal average) for z = zw(1) ||K m/s ||
 ||wpt ||Total heat flux (horizontal average) for z = zw(1) ||K m/s ||
+||zi_wtheta ||Height of the convective boundary layer (horizontal average) determined by the height of the minimum sensible heat flux ||m ||
+||zi_theta ||Height of the convective boundary layer (horizontal average) determined by the temperature profile, following the criterion of Sullivan et al. (1998) ||m ||
+||w"theta"0 ||Subgrid-scale sensible heat flux at k=0 (horizontal average) within the constant flux layer ||K m/s ||
+||w"theta" ||Subgrid-scale heat flux (horizontal average) for z = zw(1) ||K m/s ||
+||wtheta ||Total heat flux (horizontal average) for z = zw(1) ||K m/s ||
 ||w"u"0 ||Subgrid-scale momentum flux (u-component) at k=0 (horizontal average) within the constant flux layer||m^2^/s^2^ ||
 ||w"v"0 ||Subgrid-scale momentum flux (v-component) at k=0 (horizontal average) within the constant flux layer ||m^2^/s^2^ ||
 ||w"q"0 ||Subgrid-scale humidity flux at k=0 (horizontal average) within the constant flux layer, zero values are output if humidity is not used ||kg/kg m/s ||
 ||pt(0) ||Potential temperature at the surface (horizontal average) ||K ||
 ||pt(z_mo) ||Potential temperature for z = zu(1) (horizontal average) ||K ||
+||theta(0) ||Potential temperature at the surface (horizontal average) ||K ||
+||theta(z_mo) ||Potential temperature for z = zu(1) (horizontal average) ||K ||
 ||ol ||Monin-Obukhov length ||m ||
 \\
 …
 If output of horizontal cross sections is selected (see [#data_output data_output]), this parameter can be used to define the position(s) of the cross section(s). Up to 100 positions of cross sections can be selected by assigning '''section_xy''' the corresponding vertical grid point index/indices k of the requested cross section(s). The exact location (height level) of the cross section depends on the variable for which the output is made: zu(k) for scalars and horizontal velocities, zw(k) for the vertical velocity. Information about the exact location of the cross section is contained in the netCDF output file.\\\\
 Assigning '''section_xy''' = ''-1'' creates the output of horizontal cross sections averaged along z. In the netCDF output file these (averaged) cross sections are given the z-coordinate -1.0.\\\\
 Assignments to '''section_xy''' does not effect the output of horizontal cross sections of variable u* and theta* and the liquid water path lwp*. For these quantities always only one cross section (for z=zu(1)) is output.
+Assignments to '''section_xy''' does not effect the output of horizontal cross sections of variable us* and theta* and the liquid water path lwp*. For these quantities always only one cross section (for z=zu(1)) is output.
 }}}
 |----------------