Changeset 130 for palm/trunk/SOURCE/user_interface.f90

Timestamp:

Nov 13, 2007 2:08:40 PM (17 years ago)

Author:

letzel

Message:

DVRP output modifications:

• The user can now visualize user-defined quantities using dvrp.

data_output_dvrp calls the new user_interface subroutine
user_data_output_dvrp in case of unknown variables (CASE DEFAULT).

• Two instead of one digit are allowed to specify isosurface and slicer

variables with the parameter mode_dvrp.

File:

• palm/trunk/SOURCE/user_interface.f90 (modified) (2 diffs)

palm/trunk/SOURCE/user_interface.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! new subroutine user_data_output_dvrp
 ! +argument gls in user_init_grid
 !
@@ -833 +834 @@
 
 
+ SUBROUTINE user_data_output_dvrp( output_variable, local_pf )
+
+!------------------------------------------------------------------------------!
+!
+! Description:
+! ------------
+! Execution of user-defined dvrp output
+!------------------------------------------------------------------------------!
+
+    USE control_parameters
+    USE indices
+    USE pegrid
+    USE user
+
+    IMPLICIT NONE
+
+    CHARACTER (LEN=*) ::  output_variable
+
+    INTEGER ::  i, j, k
+
+    REAL, DIMENSION(nxl:nxr+1,nys:nyn+1,nzb:nz_do3d) ::  local_pf
+
+!
+!-- Here the user-defined DVRP output follows:
+
+!
+!-- Move original array to intermediate array
+    SELECT CASE ( output_variable )
+
+!       CASE ( 'u2', 'u2_xy', 'u2_xz', 'u2_yz',  )
+!!
+!!--       Here the user can add user_defined output quantities. 
+!!--       Uncomment and extend the following lines, if necessary.
+!          DO  i = nxl, nxr+1
+!             DO  j = nys, nyn+1
+!                DO  k = nzb, nz_do3d
+!                   local_pf(i,j,k) = u2(k,j,i)
+!                ENDDO
+!             ENDDO
+!          ENDDO
+
+
+       CASE DEFAULT
+!
+!--       The DEFAULT case is reached if output_variable contains a
+!--       wrong character string that is neither recognized in data_output_dvrp
+!--       nor here in user_data_output_dvrp.
+          IF ( myid == 0 )  THEN
+             PRINT*,'+++ (user_)data_output_dvrp: no output possible for: ', &
+                  output_variable
+          ENDIF
+
+    END SELECT
+
+
+ END SUBROUTINE user_data_output_dvrp
+
+
+
  SUBROUTINE user_data_output_2d( av, variable, found, grid, local_pf )