206 | | Note that {{{two_d = .TRUE.}}} is necessary for output of a 2d data slice. |
207 | | |
208 | | |
209 | | 9. For output of 3d volume data, the gridpoint data of the quantity has to be resorted to array local_pf in subroutine user_data_output_3d.: |
210 | | |
211 | | CASE ( 'u2' ) |
212 | | IF ( av == 0 ) THEN |
213 | | DO i = nxl-1, nxr+1 |
214 | | DO j = nys-1, nyn+1 |
215 | | DO k = nzb, nz_do |
216 | | local_pf(i,j,k) = u2(k,j,i) |
217 | | ENDDO |
218 | | ENDDO |
219 | | ENDDO |
220 | | ELSE |
221 | | DO i = nxl-1, nxr+1 |
222 | | DO j = nys-1, nyn+1 |
223 | | DO k = nzb, nz_do |
224 | | local_pf(i,j,k) = u2_av(k,j,i) |
225 | | ENDDO |
226 | | ENDDO |
227 | | ENDDO |
228 | | ENDIF |
229 | | |
230 | | The ELSE case is only needed in case that output of time-averaged data is requested. |
| 206 | Note that {{{two_d = .TRUE.}}} is necessary for output of a 2d data slice. \\ |
| 207 | |
| 208 | 9. For output of 3d volume data, the gridpoint data of the quantity has to be resorted to {{{array local_pf}}} in subroutine [../int#user_data_output_3d user_data_output_3d].: |
| 209 | |
| 210 | {{{ CASE ( 'u2' ) }}} \\ |
| 211 | {{{ IF ( av == 0 ) THEN }}} \\ |
| 212 | {{{ DO i = nxld, nxrd }}} \\ |
| 213 | {{{ DO j = nysd, nynd }}} \\ |
| 214 | {{{ DO k = nzb, nz_do }}} \\ |
| 215 | {{{ local_pf(i,j,k) = u2(k,j,i) }}} \\ |
| 216 | {{{ ENDDO }}} \\ |
| 217 | {{{ ENDDO }}} \\ |
| 218 | {{{ ENDDO }}} \\ |
| 219 | {{{ ELSE }}} \\ |
| 220 | {{{ DO i = nxld, nxrd }}} \\ |
| 221 | {{{ DO j = nysd, nynd }}} \\ |
| 222 | {{{ DO k = nzb, nz_do }}} \\ |
| 223 | {{{ local_pf(i,j,k) = u2_av(k,j,i) }}} \\ |
| 224 | {{{ ENDDO }}} \\ |
| 225 | {{{ ENDDO }}} \\ |
| 226 | {{{ ENDDO }}} \\ |
| 227 | {{{ ENDIF }}} \\ |
| 228 | |
| 229 | The {{{ELSE}}} case is only needed in case that output of time-averaged data is requested. |