Changeset 3361 for palm/trunk/SOURCE/turbulence_closure_mod.f90
- Timestamp:
- Oct 16, 2018 8:39:37 PM (6 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/turbulence_closure_mod.f90 (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/turbulence_closure_mod.f90
r3359 r3361 170 170 171 171 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 173 173 174 174 USE control_parameters, & … … 3006 3006 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3007 3007 ( 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 * & 3010 3010 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3011 3011 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp ) … … 3066 3066 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3067 3067 ( 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 * & 3070 3070 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3071 3071 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp ) … … 3103 3103 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3104 3104 ( 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 * & 3107 3107 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3108 3108 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp ) … … 3139 3139 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3140 3140 ( 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 * & 3143 3143 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3144 3144 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp ) … … 3178 3178 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3179 3179 ( 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 * & 3182 3182 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3183 3183 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp ) … … 3794 3794 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3795 3795 ( 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 * & 3798 3798 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3799 3799 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp ) … … 3851 3851 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3852 3852 ( 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 * & 3855 3855 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3856 3856 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp ) … … 3888 3888 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3889 3889 ( 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 * & 3892 3892 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3893 3893 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp ) … … 3924 3924 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3925 3925 ( 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 * & 3928 3928 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3929 3929 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp ) … … 3964 3964 k1 = ( 1.0_wp - q(k,j,i) + 1.61_wp * & 3965 3965 ( 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 * & 3968 3968 ( q(k,j,i) - ql(k,j,i) ) / ( temp * temp ) ) 3969 3969 k2 = theta * ( lv_d_cp / temp * k1 - 1.0_wp )
Note: See TracChangeset
for help on using the changeset viewer.
