Ignore:
Timestamp:
Oct 4, 2017 9:24:39 AM (7 years ago)
Author:
kanani
Message:

Bugfixes related to profile output of passive scalar

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/check_parameters.f90

    r2508 r2513  
    2525! -----------------
    2626! $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
    2731! Bugfix, change default value of vertical_gradient level in order to consider
    2832! also ocean runs
     
    24622466                CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 )
    24632467             ELSE
    2464                 dopr_index(i) = 117
     2468                dopr_index(i) = 115
    24652469                dopr_unit(i)  = 'kg/m3'
    2466                 hom(:,2,117,:) = SPREAD( zu, 2, statistic_regions+1 )
     2470                hom(:,2,115,:) = SPREAD( zu, 2, statistic_regions+1 )
    24672471                IF ( data_output_pr(i)(1:1) == '#' )  THEN
    2468                    dopr_initial_index(i) = 117
    2469                    hom(:,2,117,:)        = SPREAD( zu, 2, statistic_regions+1 )
    2470                    hom(nzb,2,117,:)      = 0.0_wp    ! because zu(nzb) is negative
     2472                   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
    24712475                   data_output_pr(i)     = data_output_pr(i)(2:)
    24722476                ENDIF
     
    25832587                CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 )
    25842588             ELSE
    2585                 dopr_index(i) = 119
     2589                dopr_index(i) = 117
    25862590                dopr_unit(i)  = 'kg/m3 m/s'
    2587                 hom(:,2,119,:) = SPREAD( zw, 2, statistic_regions+1 )
     2591                hom(:,2,117,:) = SPREAD( zw, 2, statistic_regions+1 )
    25882592             ENDIF
    25892593
     
    25952599                CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 )
    25962600             ELSE
    2597                 dopr_index(i) = 116
     2601                dopr_index(i) = 114
    25982602                dopr_unit(i)  = 'kg/m3 m/s'
    2599                 hom(:,2,116,:) = SPREAD( zw, 2, statistic_regions+1 )
     2603                hom(:,2,114,:) = SPREAD( zw, 2, statistic_regions+1 )
    26002604             ENDIF
    26012605
     
    26072611                CALL message( 'check_parameters', 'PA0093', 1, 2, 0, 6, 0 )
    26082612             ELSE
    2609                 dopr_index(i) = 120
     2613                dopr_index(i) = 118
    26102614                dopr_unit(i)  = 'kg/m3 m/s'
    2611                 hom(:,2,120,:) = SPREAD( zw, 2, statistic_regions+1 )
     2615                hom(:,2,118,:) = SPREAD( zw, 2, statistic_regions+1 )
    26122616             ENDIF
    26132617
     
    29352939                CALL message( 'check_parameters', 'PA0185', 1, 2, 0, 6, 0 )
    29362940             ELSE
    2937                 dopr_index(i) = 118
    2938                 dopr_unit(i)  = 'kg2/kg2'
    2939                 hom(:,2,118,:) = 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 )
    29402944             ENDIF
    29412945
     
    31623166                CALL message( 'check_parameters', 'PA0110', 1, 2, 0, 6, 0 )
    31633167             ENDIF
    3164              unit = 'conc'
     3168             unit = 'kg/m3'
    31653169
    31663170          CASE ( 'sa' )
Note: See TracChangeset for help on using the changeset viewer.