Changeset 2513 for palm/trunk/SOURCE/check_parameters.f90
- Timestamp:
- Oct 4, 2017 9:24:39 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/check_parameters.f90
r2508 r2513 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Bugfix for some dopr(_initial)_index values and units connected to 28 ! passive-scalar output 29 ! 30 ! 2508 2017-10-02 08:57:09Z suehring 27 31 ! Bugfix, change default value of vertical_gradient level in order to consider 28 32 ! also ocean runs … … 2462 2466 CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 ) 2463 2467 ELSE 2464 dopr_index(i) = 11 72468 dopr_index(i) = 115 2465 2469 dopr_unit(i) = 'kg/m3' 2466 hom(:,2,11 7,:) = SPREAD( zu, 2, statistic_regions+1 )2470 hom(:,2,115,:) = SPREAD( zu, 2, statistic_regions+1 ) 2467 2471 IF ( data_output_pr(i)(1:1) == '#' ) THEN 2468 dopr_initial_index(i) = 1 172469 hom(:,2,1 17,:) = SPREAD( zu, 2, statistic_regions+1 )2470 hom(nzb,2,1 17,:) = 0.0_wp ! because zu(nzb) is negative2472 dopr_initial_index(i) = 121 2473 hom(:,2,121,:) = SPREAD( zu, 2, statistic_regions+1 ) 2474 hom(nzb,2,121,:) = 0.0_wp ! because zu(nzb) is negative 2471 2475 data_output_pr(i) = data_output_pr(i)(2:) 2472 2476 ENDIF … … 2583 2587 CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 ) 2584 2588 ELSE 2585 dopr_index(i) = 11 92589 dopr_index(i) = 117 2586 2590 dopr_unit(i) = 'kg/m3 m/s' 2587 hom(:,2,11 9,:) = SPREAD( zw, 2, statistic_regions+1 )2591 hom(:,2,117,:) = SPREAD( zw, 2, statistic_regions+1 ) 2588 2592 ENDIF 2589 2593 … … 2595 2599 CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 ) 2596 2600 ELSE 2597 dopr_index(i) = 11 62601 dopr_index(i) = 114 2598 2602 dopr_unit(i) = 'kg/m3 m/s' 2599 hom(:,2,11 6,:) = SPREAD( zw, 2, statistic_regions+1 )2603 hom(:,2,114,:) = SPREAD( zw, 2, statistic_regions+1 ) 2600 2604 ENDIF 2601 2605 … … 2607 2611 CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 ) 2608 2612 ELSE 2609 dopr_index(i) = 1 202613 dopr_index(i) = 118 2610 2614 dopr_unit(i) = 'kg/m3 m/s' 2611 hom(:,2,1 20,:) = SPREAD( zw, 2, statistic_regions+1 )2615 hom(:,2,118,:) = SPREAD( zw, 2, statistic_regions+1 ) 2612 2616 ENDIF 2613 2617 … … 2935 2939 CALL message( 'check_parameters', 'PA0185', 1, 2, 0, 6, 0 ) 2936 2940 ELSE 2937 dopr_index(i) = 11 82938 dopr_unit(i) = 'kg2/ kg2'2939 hom(:,2,11 8,:) = SPREAD( zu, 2, statistic_regions+1 )2941 dopr_index(i) = 116 2942 dopr_unit(i) = 'kg2/m6' 2943 hom(:,2,116,:) = SPREAD( zu, 2, statistic_regions+1 ) 2940 2944 ENDIF 2941 2945 … … 3162 3166 CALL message( 'check_parameters', 'PA0110', 1, 2, 0, 6, 0 ) 3163 3167 ENDIF 3164 unit = ' conc'3168 unit = 'kg/m3' 3165 3169 3166 3170 CASE ( 'sa' )
Note: See TracChangeset
for help on using the changeset viewer.