Changes between Version 80 and Version 81 of doc/app/userint/output

Timestamp:

Sep 6, 2018 9:00:23 AM (6 years ago)

Author:

Giersch

Comment:

--

doc/app/userint/output

@@ -75 +75 @@
    3. In order to store the quantities' data within PALM, a 3d data array has to be declared in module [../int#user user] (user_module.f90): \\
 
-   {{{ REAL, DIMENSION(:,:,:), ALLOCATABLE ::  u2, u2_av }}} //
+   {{{ REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  u2, u2_av }}} //
 
    The second array {{{u2_av}}} is needed in case that output of time averaged data is requested. It is used to store the sum of the data of the respective time levels over which the average has to be carried out. \\
@@ -124 +124 @@
       {{{    ... }}} \\
       {{{       CASE ( 'u2' ) }}} \\
-         {{{       DO  i = nxlg, nxrg }}} \\
-            {{{       DO  j = nysg, nyng }}} \\
-               {{{       DO  k = nzb, nzt+1 }}} \\
-                  {{{       u2_av(k,j,i) = u2_av(k,j,i) + u2(k,j,i) }}} \\
+         {{{       IF ( ALLOCATED( u2_av ) ) THEN }}} \\
+            {{{       DO  i = nxlg, nxrg }}} \\
+               {{{       DO  j = nysg, nyng }}} \\
+                  {{{       DO  k = nzb, nzt+1 }}} \\
+                     {{{       u2_av(k,j,i) = u2_av(k,j,i) + u2(k,j,i) }}} \\
+                  {{{ ENDDO }}} \\
                {{{ ENDDO }}} \\
             {{{ ENDDO }}} \\
-         {{{ ENDDO }}} \\
+         {{{ ENDIF}}} \\
          {{{    ...  }}} \\
    {{{ ELSEIF ( mode == 'average' )  THEN }}} \\
       {{{ ... }}} \\
       {{{    CASE ( 'u2' ) }}} \\
-         {{{    DO  i = nxlg, nxrg }}} \\
-            {{{    DO  j = nysg, nyng }}} \\
-               {{{    DO  k = nzb, nzt+1 }}} \\
-                  {{{    u2_av(k,j,i) = u2_av(k,j,i) / REAL( average_count_3d ) }}} \\
+         {{{       IF ( ALLOCATED( u2_av ) ) THEN }}} \\
+            {{{    DO  i = nxlg, nxrg }}} \\
+               {{{    DO  j = nysg, nyng }}} \\
+                  {{{    DO  k = nzb, nzt+1 }}} \\
+                     {{{    u2_av(k,j,i) = u2_av(k,j,i) / REAL( average_count_3d, KIND=wp ) }}} \\
+                  {{{ ENDDO }}} \\
                {{{ ENDDO }}} \\
             {{{ ENDDO }}} \\
-         {{{ ENDDO }}} \\
+         {{{ ENDIF}}} \\
+         {{{    ...  }}} \\
        
    9. For output of '''2d cross sections''', the data of the quantity has to be resorted to array {{{local_pf}}} in subroutine [../int#user_data_output_2d user_data_output_2d]. Also the vertical grid, on which the quantity is defined, has to be set again: \\
@@ -147 +152 @@
    {{{ CASE ( 'u2_xy', 'u2_xz', 'u2_yz' ) }}} \\
       {{{    IF ( av == 0 )  THEN }}} \\
-         {{{    DO  i = nxlg, nxrg }}} \\
-            {{{    DO  j = nysg, nyng }}} \\
-               {{{    DO  k = nzb, nzt+1 }}} \\
+         {{{    DO  i = nxl, nxr }}} \\
+            {{{    DO  j = nys, nyn }}} \\
+               {{{    DO  k = nzb_do, nzt_do }}} \\
                   {{{    local_pf(i,j,k) = u2(k,j,i) }}} \\
                {{{    ENDDO }}} \\
@@ -155 +160 @@
          {{{    ENDDO }}} \\
       {{{    ELSE }}} \\
-         {{{    DO  i = nxlg, nxrg }}} \\
-            {{{    DO  j = nysg, nyng }}} \\
-               {{{    DO  k = nzb, nzt+1 }}} \\
+         {{{    IF ( .NOT. ALLOCATED( u2_av ) ) THEN  }}} \\
+            {{{    ALLOCATE( u2_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) }}} \\
+            {{{    u2_av = REAL( fill_value, KIND = wp ) }}} \\
+         {{{    ENDIF }}} \\
+         {{{    DO  i = nxl, nxr }}} \\
+            {{{    DO  j = nys, nyn }}} \\
+               {{{    DO  k = nzb_do, nzt_do }}} \\
                   {{{    local_pf(i,j,k) = u2_av(k,j,i) }}} \\
                {{{    ENDDO }}} \\