Changes between Version 76 and Version 77 of doc/app/maskedoutput

Timestamp:

Oct 26, 2018 8:03:09 PM (7 years ago)

Author:

gronemeier

Comment:

--

doc/app/maskedoutput

@@ -6 +6 @@
 NetCDF4 file format is already available for masked data output.\\
 Default quantities (instantaneous and time averaged) are declared with the parameter [../d3par#data_output_masks data_output_masks], user-defined quantities can be output with [../userpar#data_output_masks_user data_output_masks_user] (see also [[../userint/output|User-defined output quantities]]).\\\\
+Since r3436, terrain-following masked output is implemented. This allows to output masked quantities at a specified height above the surface.
 == Output steering ==
 Different [../d3par runtime parameters] steer the output of the quantities assigned with [../d3par#data_output_masks data_output_masks] and/or [../userpar#data_output_masks_user data_output_masks_user]:[=#list]
@@ -15 +16 @@
 ||[../d3par#mask_y_loop mask_y_loop] ||loop begin, end and stride for y-coordinates of mask locations for masks (in multiples of mask scale) ||
 ||[../d3par#mask_z_loop mask_z_loop] ||loop begin, end and stride for z-coordinates of mask locations for masks (in multiples of mask scale) ||
+||[../d3par#mask_k_over_surface mask_k_over_surface] ||vertical grid index above a surface to use for terrain-following masked output ||
 These parameters define the masks. Each mask can be steered with a separate set of x-, y- and z-parameters. By default all gridpoints along the respective direction are used.
 \\\\
@@ -36 +38 @@
 By default 50 different masks can be assigned. The amount of masks is steered with '''max_masks''' (=50 by default), which can only be changed in {{{modules.f90}}}. If more than 50 masks are desired, '''max_masks''' has to be changed to the corresponding value.
 \\\\
-== Further assignments ==
-Besides assigning the parameters listed [#list above], the required output files have to be commented out in the '''mrun'''-[../configexample configuration file]:
-{{{
-#----------------------------------------------------------------------------
-# List of output-files
-#----------------------------------------------------------------------------
-DATA_MASK_01_NETCDF       out:loc      ma#:maf         $base_data/$fname/OUTPUT      _m01     nc
-DATA_MASK_01_AV_NETCDF    out:loc      ma#:maf         $base_data/$fname/OUTPUT      _m01_av  nc
-DATA_MASK_02_NETCDF       out:loc      ma#:maf         $base_data/$fname/OUTPUT      _m02     nc
-DATA_MASK_02_AV_NETCDF    out:loc      ma#:maf         $base_data/$fname/OUTPUT      _m02_av  nc
-#DATA_MASK_03_NETCDF      out:loc      ma#:maf         $base_data/$fname/OUTPUT      _m03     nc
-#DATA_MASK_03_AV_NETCDF   out:loc      ma#:maf         $base_data/$fname/OUTPUT      _m03_av  nc
-...
-}}}
-Here output of '''masks 1''' and '''2''' is done. '''Masks 3''' will not be output even if the parameters are assigned.\\\\\\
-Within the '''mrun''' call ''ma#'' has to be added to the '''mrun'''-option -r:
-{{{
-mrun -d example_cbl -h... -r "d3# ma#"
-}}} 
-\\
 == Examples ==
 The following examples are given for the [../examples/cbl example_cbl run] with a model domain size of 2000m x 2000m x 2000m and a resolution of 50.0m.
@@ -213 +195 @@
 }}}
 Output starts from the beginning for '''mask 1''' and after 1h simulation time (see [../d3par#skip_time_domask skip_time_domask]) for '''mask 2, 3''' and '''4'''. It is done every 30 minutes for '''mask 1''' and at time intervals of [../d3par#dt_data_output dt_data_output] for '''mask 2, 3''' and '''4''' (see [../d3par#dt_domask dt_domask]).
+=== 6. Output of terrain-following mask ===
+'''Mask 1:''' xy cross section from x = 500m to x = 1000m and from y = 1000m to 2000m at the second, fith and tenth grid point above the surface. Output of instantaneous and time averaged data of the three wind components.\\
+{{{
+&initialization_parameters
+
+    nx = 39, ny = 39, nz = 40,
+    dx = 50.0, dy = 50.0, dz = 50.0,
+    .../
+
+&runtime_parameters
+
+    ...
+    data_output_masks(1,:) = 'u','v','w','u_av','v_av','w_av',
+   
+    mask_x_loop(1,:) = 500.,1000.,50. ,
+    mask_y_loop(1,:) = 1000.,2000.,50. ,
+    mask_k_over_surface(1,:) = 2, 5, 10,
+
+    dt_domask = 600.,
+    skip_time_domask = 3600.,/
+}}}
+Output starts after 1h simulation time (see [../d3par#skip_time_domask skip_time_domask]) for both masks and is done every 10 minutes for '''mask 1''' and every 30 minutes for '''mask 2''' (see [../d3par#dt_domask dt_domask]).