Changeset 3274 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
- Timestamp:
- Sep 24, 2018 3:42:55 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
r3271 r3274 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Modularization of all bulk cloud physics code components 29 ! 30 ! 3271 2018-09-24 08:20:34Z suehring 28 31 ! Comment revised 29 32 ! … … 231 234 USE arrays_3d, & 232 235 ONLY: e, kh, nc, nr, pt, q, ql, qc, qr, s, u, v, vpt, w, zu, zw, & 233 drho_air_zw, rho_air_zw 236 drho_air_zw, rho_air_zw, d_exner 237 238 USE basic_constants_and_equations_mod, & 239 ONLY: g, kappa, lv_d_cp, pi 234 240 235 241 USE chem_modules, & 236 242 ONLY: constant_csflux, nvar 237 243 238 USE cloud_parameters, &239 ONLY: l_d_cp, pt_d_t240 241 USE constants, &242 ONLY: pi243 244 244 USE cpulog 245 245 246 246 USE control_parameters, & 247 ONLY: air_chemistry, cloud_droplets, cloud_physics,&248 constant_heatflux, constant_scalarflux, & 249 constant_waterflux, coupling_mode, g, humidity, ibc_e_b,&250 ibc_pt_b, initializing_actions, kappa,&247 ONLY: air_chemistry, cloud_droplets, & 248 constant_heatflux, constant_scalarflux, & 249 constant_waterflux, coupling_mode, humidity, ibc_e_b, & 250 ibc_pt_b, initializing_actions, & 251 251 intermediate_timestep_count, intermediate_timestep_count_max, & 252 land_surface, large_scale_forcing, lsf_surf, & 253 message_string, microphysics_morrison, microphysics_seifert, & 252 land_surface, large_scale_forcing, lsf_surf, message_string, & 254 253 most_method, neutral, passive_scalar, pt_surface, q_surface, & 255 254 run_coupled, surface_pressure, simulated_time, terminate_run, & … … 264 263 265 264 USE kinds 265 266 USE bulk_cloud_model_mod, & 267 ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert 266 268 267 269 USE pegrid … … 1467 1469 1468 1470 ENDDO 1469 ELSEIF ( cloud_physics.OR. cloud_droplets ) THEN1471 ELSEIF ( bulk_cloud_model .OR. cloud_droplets ) THEN 1470 1472 !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo ) 1471 1473 DO m = 1, surf%ns … … 1596 1598 k = surf%k(m) 1597 1599 1598 IF ( cloud_physics) THEN1599 surf%pt1(m) = pt(k,j,i) + l _d_cp * pt_d_t(k) * ql(k,j,i)1600 IF ( bulk_cloud_model ) THEN 1601 surf%pt1(m) = pt(k,j,i) + lv_d_cp * d_exner(k) * ql(k,j,i) 1600 1602 surf%qv1(m) = q(k,j,i) - ql(k,j,i) 1601 1603 ELSEIF( cloud_droplets ) THEN 1602 surf%pt1(m) = pt(k,j,i) + l _d_cp * pt_d_t(k) * ql(k,j,i)1604 surf%pt1(m) = pt(k,j,i) + lv_d_cp * d_exner(k) * ql(k,j,i) 1603 1605 surf%qv1(m) = q(k,j,i) 1604 1606 ELSE … … 1820 1822 ENDIF 1821 1823 1822 IF ( cloud_physics.OR. cloud_droplets ) THEN1824 IF ( bulk_cloud_model .OR. cloud_droplets ) THEN 1823 1825 !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo ) 1824 1826 DO m = 1, surf%ns … … 1927 1929 ! 1928 1930 !-- If required compute qc* and nc* 1929 IF ( cloud_physics .AND. microphysics_morrison .AND.&1931 IF ( bulk_cloud_model .AND. microphysics_morrison .AND. & 1930 1932 .NOT. surf_vertical ) THEN 1931 1933 !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo ) … … 1953 1955 ! 1954 1956 !-- If required compute qr* and nr* 1955 IF ( cloud_physics .AND. microphysics_seifert .AND.&1957 IF ( bulk_cloud_model .AND. microphysics_seifert .AND. & 1956 1958 .NOT. surf_vertical ) THEN 1957 1959 !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo ) … … 2145 2147 ! 2146 2148 !-- Compute (turbulent) fluxes of cloud water content and cloud drop conc. 2147 IF ( cloud_physics .AND. microphysics_morrison .AND.&2149 IF ( bulk_cloud_model .AND. microphysics_morrison .AND. & 2148 2150 .NOT. downward) THEN 2149 2151 !$OMP PARALLEL DO PRIVATE( i, j ) … … 2159 2161 ! 2160 2162 !-- Compute (turbulent) fluxes of rain water content and rain drop conc. 2161 IF ( cloud_physics .AND. microphysics_seifert .AND.&2163 IF ( bulk_cloud_model .AND. microphysics_seifert .AND. & 2162 2164 .NOT. downward) THEN 2163 2165 !$OMP PARALLEL DO PRIVATE( i, j )
Note: See TracChangeset
for help on using the changeset viewer.