Changes between Version 90 and Version 91 of doc/app/userint/output

Timestamp:

Oct 17, 2019 6:59:19 AM (5 years ago)

Author:

Giersch

Comment:

--

doc/app/userint/output

@@ -46 +46 @@
 == 2. Output of user-defined time series
 This example shows the output of two time series for the maxima of the absolute values of the horizontal velocities u and v. If more than one user-defined quantity shall be output, the following steps have to be carried out in the same way for each of the quantities.\\\\
-   1. For each time series quantity you have to give a label and a unit in subroutine [../int#user_init user_init] (user_check_data_output_ts for newer versions), which will be used for the netCDF file. They must not contain more than '''13''' characters. The value of {{{dots_num}}} has to be increased by the number of new time series quantities. Its old value has to be stored in {{{dots_num_palm}}}.
-
-   {{{dots_label(dots_num+}}}[[span({{{1}}},style=color: red)]]{{{) =}}} '[[span({{{abs_umx}}},style=color: red)]]' \\
-   {{{dots_unit(dots_num+}}}[[span({{{1}}},style=color: red)]]{{{) =}}} '[[span({{{m/s}}},style=color: red)]]' \\
-   {{{dots_label(dots_num+}}}[[span({{{2}}},style=color: red)]]{{{) =}}} '[[span({{{abs_vmx}}},style=color: red)]]' \\
-   {{{dots_unit(dots_num+}}}[[span({{{2}}},style=color: red)]]{{{) =}}} '[[span({{{m/s}}},style=color: red)]]' \\
+   1. For each time series quantity you have to give a label and a unit in subroutine [../int#user_init user_init] (user_check_data_output_ts for newer versions), which will be used for the netCDF file. They must not contain more than '''13''' characters. The value of {{{dots_num}}} and {{{dots_num_user}}} have to be increased by the number of new time series quantities. The old value of {{{dots_num}}} has to be stored in {{{dots_num_palm}}}.
 
    {{{dots_num_palm = dots_num}}} \\
-   {{{dots_num = dots_num}}} [[span({{{+2}}},style=color: red)]] \\
-
+
+   {{{dots_num = dots_num + 1}}} \\
+   {{{dots_num_user = dots_num_user + 1}}} \\
+   {{{dots_label(dots_num) = 'abs_umax'}}}' \\
+   {{{dots_unit(dots_num) = 'm/s'}}} \\
+
+   {{{dots_num = dots_num + 1}}} \\
+   {{{dots_num_user = dots_num_user + 1}}} \\
+   {{{dots_label(dots_num) = 'abs_vmax'}}} \\
+   {{{dots_unit(dots_num) = 'm/s'}}} \\
 
    2. These quantities are calculated and output in subroutine [../int#user_statistics user_statistics] for every statistic region {{{sr}}} defined by the user, but at least for the region "total domain":
 
-   {{{ ts_value(dots_num_palm+1,sr) = ABS( u_max ) }}} \\
-   {{{ ts_value(dots_num_palm+2,sr) = ABS( v_max ) }}} \\\\
+   {{{ ts_value_l(1) = ABS( u_max ) }}} \\
+   {{{ ts_value_l(2) = ABS( v_max ) }}} \\
+
+Subsequently, collect / send values to PE0, because only PE0 outputs the time series. However, collection is done by taking the sum over all processors. You may have to normalize this sum, depending on the quantity that you like to calculate. For serial runs, nothing has to be done:
+
+{{{#if defined( __parallel )}}}
+       {{{IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )}}}\\
+       {{{CALL MPI_ALLREDUCE( ts_value_l(1),                            &}}}\\
+                           {{{ts_value(dots_num_palm+1,sr),             &}}}\\
+                           {{{dots_num_user, MPI_REAL, MPI_MAX, comm2d, &}}}\\
+                           {{{ierr )}}}
+{{{#else}}}\\
+       {{{ts_value(dots_num_palm+1:dots_num_palm+dots_num_user,sr) = ts_value_l}}}\\
+{{{#endif}}}\\