Changeset 89 for palm/trunk/SOURCE
- Timestamp:
- May 25, 2007 12:08:31 PM (17 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/check_parameters.f90
r87 r89 1405 1405 CASE ( 'Sw' ) 1406 1406 dopr_index(i) = 39 1407 dopr_unit(i) = 'none' 1407 1408 hom(:,2,39,:) = SPREAD( zw, 2, statistic_regions+1 ) 1408 1409 -
palm/trunk/SOURCE/modules.f90
r87 r89 5 5 ! Actual revisions: 6 6 ! ----------------- 7 ! +data_output_pr_user, max_pr_user, num_user_profiles, size of data_output_pr,8 ! dopr_in dex, dopr_initial_index and dopr_unit enlarged,7 ! +data_output_pr_user, max_pr_user, size of data_output_pr, dopr_index, 8 ! dopr_initial_index and dopr_unit enlarged, 9 9 ! var_hom and var_sum renamed pr_palm 10 10 ! … … 275 275 mg_switch_to_pe0_level = 0, ngsrb = 2, nsor = 20, & 276 276 nsor_ini = 100, n_sor, normalizing_region = 0, & 277 n um_user_profiles = 0, nz_do1d, nz_do3d = -9999, &278 outflow_damping_width = -1, prt_time_count = 0, runnr= 0, &279 skip_do_avs = 0,timestep_count = 0277 nz_do1d, nz_do3d = -9999, outflow_damping_width = -1, & 278 prt_time_count = 0, runnr = 0, skip_do_avs = 0, & 279 timestep_count = 0 280 280 281 281 INTEGER :: dist_nxl(0:1), dist_nxr(0:1), dist_nyn(0:1), dist_nys(0:1), & -
palm/trunk/SOURCE/user_interface.f90
r87 r89 70 70 71 71 72 NAMELIST /userpar/ data_output_pr_user, data_output_user, max_pr_user, & 73 region 72 NAMELIST /userpar/ data_output_pr_user, data_output_user, region 74 73 75 74 ! … … 89 88 READ ( 11, userpar ) 90 89 user_defined_namelist_found = .TRUE. 91 92 !93 !-- Determine the maximum number of user-profiles allowed to be calculated.94 !-- This number may be (much) larger than num_user_profiles which only gives95 !-- the number of user-defined profiles to be output.96 IF ( data_output_pr_user(1) /= ' ' .AND. max_pr_user == 0 ) THEN97 max_pr_user = 10098 ENDIF99 90 100 91 ! … … 109 100 DO WHILE ( data_output_pr_user(j) /= ' ' .AND. j <= 100 ) 110 101 data_output_pr(i) = data_output_pr_user(j) 111 num_user_profiles = num_user_profiles+ 1102 max_pr_user = max_pr_user + 1 112 103 i = i + 1 113 104 j = j + 1 114 105 ENDDO 115 ENDIF116 117 IF ( num_user_profiles > max_pr_user ) THEN118 IF ( myid == 0 ) THEN119 PRINT*, '+++ user_parin: number of user-defined profiles given by '120 PRINT*, ' data_output_pr_user (', num_user_profiles, &121 ') is larger than the allowed maximum'122 PRINT*, ' number of profiles: max_pr_user(', &123 max_pr_user, ')'124 ENDIF125 CALL local_stop126 106 ENDIF 127 107 … … 529 509 ! DO k = nzb_s_outer(j,i)+1, nzt 530 510 !! 531 !!-- Sample on how to calculate the profile for vertical velocity 532 ! sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) + & 533 ! w(k,j,i) & 511 !!-- Sample on how to calculate the profile of the resolved-scale 512 !!-- horizontal momentum flux u*v* 513 ! sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) + & 514 ! ( 0.5 * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) * & 515 ! ( 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) ) * & 534 516 ! * rmask(j,i,sr) 535 517 !! … … 770 752 !-- the number "pr_palm+1", the second one "pr_palm+2", etc. The respective 771 753 !-- user-profile-numbers have also to be used in routine user_statistics! 772 ! CASE ( ' w_user' )! quantity string as given in773 ! ! data_output_pr_user754 ! CASE ( 'u*v*' ) ! quantity string as given in 755 ! ! data_output_pr_user 774 756 ! index = pr_palm + 1 775 ! dopr_index(var_count) = index ! thequantities' user-profile-number776 ! dopr_unit(var_count) = 'm /s' ! quantity unit777 ! hom(:,2,index,:) = SPREAD( z w, 2, statistic_regions+1 )778 ! ! grid on which the quantity is779 ! ! defined (use zu or zw)757 ! dopr_index(var_count) = index ! quantities' user-profile-number 758 ! dopr_unit(var_count) = 'm2/s2' ! quantity unit 759 ! hom(:,2,index,:) = SPREAD( zu, 2, statistic_regions+1 ) 760 ! ! grid on which the quantity is 761 ! ! defined (use zu or zw) 780 762 781 763 CASE DEFAULT
Note: See TracChangeset
for help on using the changeset viewer.