Changeset 2513 for palm/trunk
 Timestamp:
 Oct 4, 2017 9:24:39 AM (7 years ago)
 Location:
 palm/trunk/SOURCE
 Files:

 2 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 ! passivescalar output 29 ! 30 ! 2508 20171002 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' ) 
palm/trunk/SOURCE/init_3d_model.f90
r2350 r2513 25 25 !  26 26 ! $Id$ 27 ! Bugfix in storing initial scalar profile (wrong index) 28 ! 29 ! 2350 20170815 11:48:26Z kanani 27 30 ! Bugfix in nopointer version 28 31 ! … … 1301 1304 ! 1302 1305 ! Store initial scalar profile 1303 hom(:,1,1 15,:) = SPREAD( s(:,nys,nxl), 2, statistic_regions+1 )1306 hom(:,1,121,:) = SPREAD( s(:,nys,nxl), 2, statistic_regions+1 ) 1304 1307 ENDIF 1305 1308
Note: See TracChangeset
for help on using the changeset viewer.