Ignore:
Timestamp:
Oct 16, 2018 8:39:37 PM (3 years ago)
Author:
knoop
Message:

Introduced global constant rd_d_rv=0.622

File:
1 edited

Legend:

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

    r3359 r3361  
    170170
    171171    USE basic_constants_and_equations_mod,                                     &
    172         ONLY:  g, kappa, lv_d_cp, lv_d_rd
     172        ONLY:  g, kappa, lv_d_cp, lv_d_rd, rd_d_rv
    173173
    174174    USE control_parameters,                                                    &
     
    30063006                         k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                  &
    30073007                                       ( q(k,j,i) - ql(k,j,i) ) *              &
    3008                               ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /       &
    3009                               ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *        &
     3008                              ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /        &
     3009                              ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *         &
    30103010                              ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    30113011                         k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
     
    30663066                               k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *            &
    30673067                                          ( q(k,j,i) - ql(k,j,i) ) *           &
    3068                                  ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /    &
    3069                                  ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *     &
     3068                                 ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /     &
     3069                                 ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *      &
    30703070                                 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    30713071                               k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
     
    31033103                            k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *               &
    31043104                                          ( q(k,j,i) - ql(k,j,i) ) *           &
    3105                                  ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /    &
    3106                                  ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *     &
     3105                                 ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /     &
     3106                                 ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *      &
    31073107                                 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    31083108                            k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
     
    31393139                            k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *               &
    31403140                                          ( q(k,j,i) - ql(k,j,i) ) *           &
    3141                                  ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /    &
    3142                                  ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *     &
     3141                                 ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /     &
     3142                                 ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *      &
    31433143                                 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    31443144                            k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
     
    31783178                            k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *               &
    31793179                                       ( q(k,j,i) - ql(k,j,i) ) *              &
    3180                               ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /       &
    3181                               ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *        &
     3180                              ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /        &
     3181                              ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *         &
    31823182                              ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    31833183                            k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
     
    37943794                   k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                        &
    37953795                                 ( q(k,j,i) - ql(k,j,i) ) *                    &
    3796                         ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /             &
    3797                         ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *              &
     3796                        ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /              &
     3797                        ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *               &
    37983798                        ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    37993799                   k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
     
    38513851                        k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                   &
    38523852                                   ( q(k,j,i) - ql(k,j,i) ) *                  &
    3853                           ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /           &
    3854                           ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *            &
     3853                          ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /            &
     3854                          ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *             &
    38553855                          ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    38563856                        k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
     
    38883888                       k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                    &
    38893889                                     ( q(k,j,i) - ql(k,j,i) ) *                &
    3890                             ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /         &
    3891                             ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *          &
     3890                            ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /          &
     3891                            ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *           &
    38923892                            ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    38933893                       k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
     
    39243924                       k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                    &
    39253925                                     ( q(k,j,i) - ql(k,j,i) ) *                &
    3926                             ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /         &
    3927                             ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *          &
     3926                            ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /          &
     3927                            ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *           &
    39283928                            ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    39293929                       k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
     
    39643964                      k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp *                     &
    39653965                                 ( q(k,j,i) - ql(k,j,i) ) *                    &
    3966                         ( 1.0_wp + 0.622_wp * lv_d_rd / temp ) ) /             &
    3967                         ( 1.0_wp + 0.622_wp * lv_d_rd * lv_d_cp *              &
     3966                        ( 1.0_wp + rd_d_rv * lv_d_rd / temp ) ) /              &
     3967                        ( 1.0_wp + rd_d_rv * lv_d_rd * lv_d_cp *               &
    39683968                        ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) )
    39693969                      k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
Note: See TracChangeset for help on using the changeset viewer.