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

Timestamp:

May 25, 2007 12:08:31 PM (17 years ago)

Author:

raasch

Message:

further changes concerning user-defined profiles

File:

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

palm/trunk/SOURCE/user_interface.f90

@@ -70 +70 @@
 
 
-    NAMELIST /userpar/  data_output_pr_user, data_output_user, max_pr_user, &
-                        region
+    NAMELIST /userpar/  data_output_pr_user, data_output_user, region
 
 !
@@ -89 +88 @@
     READ ( 11, userpar )
     user_defined_namelist_found = .TRUE.
-
-!
-!-- Determine the maximum number of user-profiles allowed to be calculated.
-!-- This number may be (much) larger than num_user_profiles which only gives
-!-- the number of user-defined profiles to be output.
-    IF ( data_output_pr_user(1) /= ' '  .AND.  max_pr_user == 0 )  THEN
-       max_pr_user = 100
-    ENDIF
 
 !
@@ -109 +100 @@
        DO  WHILE ( data_output_pr_user(j) /= ' '  .AND.  j <= 100 )
           data_output_pr(i) = data_output_pr_user(j)
-          num_user_profiles = num_user_profiles + 1
+          max_pr_user       = max_pr_user + 1
           i = i + 1
           j = j + 1
        ENDDO
-    ENDIF
-
-    IF ( num_user_profiles > max_pr_user )  THEN
-       IF ( myid == 0 )  THEN
-          PRINT*, '+++ user_parin: number of user-defined profiles given by '
-          PRINT*, '                data_output_pr_user (', num_user_profiles, &
-                                   ') is larger than the allowed maximum'
-          PRINT*, '                number of profiles: max_pr_user(', &
-                                   max_pr_user, ')'
-       ENDIF
-       CALL local_stop
     ENDIF
 
@@ -529 +509 @@
 !             DO  k = nzb_s_outer(j,i)+1, nzt
 !!
-!!--             Sample on how to calculate the profile for vertical velocity
-!                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +       &
-!                                                 w(k,j,i)               &
+!!--             Sample on how to calculate the profile of the resolved-scale
+!!--             horizontal momentum flux u*v*
+!                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +           &
+!                      ( 0.5 * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) * &
+!                      ( 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) ) * &
 !                                                 * rmask(j,i,sr)
 !!
@@ -770 +752 @@
 !--    the number "pr_palm+1", the second one "pr_palm+2", etc. The respective
 !--    user-profile-numbers have also to be used in routine user_statistics!
-!       CASE ( 'w_user' )                  ! quantity string as given in
-!                                          ! data_output_pr_user
+!       CASE ( 'u*v*' )                      ! quantity string as given in
+!                                            ! data_output_pr_user
 !          index = pr_palm + 1
-!          dopr_index(var_count)  = index  ! the quantities' user-profile-number
-!          dopr_unit(var_count)   = 'm/s'  ! quantity unit
-!          hom(:,2,index,:)       = SPREAD( zw, 2, statistic_regions+1 )
-!                                          ! grid on which the quantity is
-!                                          ! defined (use zu or zw)
+!          dopr_index(var_count)  = index    ! quantities' user-profile-number
+!          dopr_unit(var_count)   = 'm2/s2'  ! quantity unit
+!          hom(:,2,index,:)       = SPREAD( zu, 2, statistic_regions+1 )
+!                                            ! grid on which the quantity is
+!                                            ! defined (use zu or zw)
 
        CASE DEFAULT