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

Timestamp:

Mar 6, 2007 12:28:36 PM (17 years ago)

Author:

raasch

Message:

preliminary version, several changes to be explained later

File:

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

palm/trunk/SOURCE/user_interface.f90

@@ -5 +5 @@
 ! -----------------
 ! New routine user_init_3d_model which allows the initial setting of all 3d
-! arrays under control of the user
+! arrays under control of the user,
+! routine user_statistics now has one argument (sr),
+! sample for generating time series quantities added
 !
 ! Former revisions:
@@ -174 +176 @@
 !       ENDDO
 !    ENDIF
+
+!
+!-- Sample for user-defined time series
+!-- For each time series quantity you have to give a label and a unit,
+!-- which will be used for the NetCDF file. They must not contain more than
+!-- seven characters. The value of dots_num has to be increased by the
+!-- number of new time series quantities. Its old value has to be store in
+!-- dots_num_palm. See routine user_statistics on how to output calculate
+!-- and output these quantities.
+!    dots_label(dots_num+1) = 'abs_umx'
+!    dots_unit(dots_num+1)  = 'm/s'
+!    dots_label(dots_num+2) = 'abs_vmx'
+!    dots_unit(dots_num+2)  = 'm/s'
+!
+!    dots_num_palm = dots_num
+!    dots_num = dots_num + 2
 
  END SUBROUTINE user_init
@@ -421 +439 @@
 
 
- SUBROUTINE user_statistics
+ SUBROUTINE user_statistics( sr )
 
 !------------------------------------------------------------------------------!
@@ -428 +446 @@
 ! ------------
 ! Calculation of user-defined statistics
+! This routine is called for every statistic region sr defined by the user,
+! but at least for the region "total domain" (sr=0).
 !------------------------------------------------------------------------------!
 
@@ -435 +455 @@
     IMPLICIT NONE
 
+    INTEGER ::  sr
+
+!
+!-- Sample on how to add values for the user-defined time series quantities.
+!-- These have to be defined before in routine user_init. This sample
+!-- creates two time series for the absolut values of the horizontal
+!-- velocities u and v.
+!    ts_value(dots_num_palm+1,sr) = ABS( u_max )
+!    ts_value(dots_num_palm+2,sr) = ABS( v_max )
 
  END SUBROUTINE user_statistics