Changeset 1822 for palm/trunk/SOURCE/boundary_conds.f90

Timestamp:

Apr 7, 2016 7:49:42 AM (8 years ago)

Author:

hoffmann

Message:

changes in LPM and bulk cloud microphysics

File:

• palm/trunk/SOURCE/boundary_conds.f90 (modified) (7 diffs)

palm/trunk/SOURCE/boundary_conds.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-!
+! icloud_scheme removed. microphyisics_seifert added.
 !
 ! Former revisions:
@@ -139 +139 @@
                cloud_physics, dt_3d, humidity,                                 &
                ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_sa_t, ibc_uv_b,       &
-               ibc_uv_t, icloud_scheme, inflow_l, inflow_n, inflow_r, inflow_s,&
-               intermediate_timestep_count, large_scale_forcing, nest_domain,  &
-               nest_bound_l, nest_bound_s, nudging, ocean,                     &
-               outflow_l, outflow_n, outflow_r, outflow_s, passive_scalar,     &
-               precipitation, tsc, use_cmax
+               ibc_uv_t, inflow_l, inflow_n, inflow_r, inflow_s,               &
+               intermediate_timestep_count, large_scale_forcing,               &
+               microphysics_seifert, nest_domain, nest_bound_l, nest_bound_s,  &
+               nudging, ocean, outflow_l, outflow_n, outflow_r, outflow_s,     &
+               passive_scalar, tsc, use_cmax
 
     USE grid_variables,                                                        &
@@ -319 +319 @@
        ENDIF
 
-       IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.  precipitation )  THEN
+       IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
 !             
 !--       Surface conditions rain water (Dirichlet)
@@ -373 +373 @@
        IF ( humidity  .OR.  passive_scalar )  THEN
           q_p(:,nys-1,:) = q_p(:,nys,:)
-          IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.              &
-               precipitation)  THEN
+          IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
              qr_p(:,nys-1,:) = qr_p(:,nys,:)
              nr_p(:,nys-1,:) = nr_p(:,nys,:)
@@ -384 +383 @@
        IF ( humidity  .OR.  passive_scalar )  THEN
           q_p(:,nyn+1,:) = q_p(:,nyn,:)
-          IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.              &
-               precipitation )  THEN
+          IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
              qr_p(:,nyn+1,:) = qr_p(:,nyn,:)
              nr_p(:,nyn+1,:) = nr_p(:,nyn,:)
@@ -395 +393 @@
        IF ( humidity  .OR.  passive_scalar )  THEN
           q_p(:,:,nxl-1) = q_p(:,:,nxl)
-          IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.              &
-               precipitation )  THEN
+          IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
              qr_p(:,:,nxl-1) = qr_p(:,:,nxl)
              nr_p(:,:,nxl-1) = nr_p(:,:,nxl)
@@ -406 +403 @@
        IF ( humidity .OR. passive_scalar )  THEN
           q_p(:,:,nxr+1) = q_p(:,:,nxr)
-          IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.  precipitation )  THEN
+          IF ( cloud_physics  .AND.  microphysics_seifert )  THEN
              qr_p(:,:,nxr+1) = qr_p(:,:,nxr)
              nr_p(:,:,nxr+1) = nr_p(:,:,nxr)