Changeset 1159 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

May 21, 2013 11:58:22 AM (11 years ago)

Author:

fricke

Message:

Bugfix: In case of non-cyclic lateral boundary conditions, Neumann boundary conditions for the velocity components at the outflow are in fact radiation boundary conditions using the maximum phase velocity that ensures numerical stability (CFL-condition).
Logical operator use_cmax is now used instead of bc_lr_dirneu/_neudir.

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (5 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -23 +23 @@
 ! Current revisions:
 ! ------------------
+! -bc_lr_dirneu, bc_lr_neudir, bc_ns_dirneu, bc_ns_neudir
 !
 !
@@ -1695 +1696 @@
     ptdf_x = 0.0
     ptdf_y = 0.0
-    IF ( bc_lr_dirrad  .OR.  bc_lr_dirneu )  THEN
+    IF ( bc_lr_dirrad )  THEN
        DO  i = nxl, nxr
           IF ( ( i * dx ) < pt_damping_width )  THEN
@@ -1703 +1704 @@
           ENDIF
        ENDDO
-    ELSEIF ( bc_lr_raddir  .OR.  bc_lr_neudir )  THEN
+    ELSEIF ( bc_lr_raddir )  THEN
        DO  i = nxl, nxr
           IF ( ( i * dx ) > ( nx * dx - pt_damping_width ) )  THEN
@@ -1711 +1712 @@
           ENDIF
        ENDDO 
-    ELSEIF ( bc_ns_dirrad  .OR.  bc_ns_dirneu )  THEN
+    ELSEIF ( bc_ns_dirrad )  THEN
        DO  j = nys, nyn
           IF ( ( j * dy ) > ( ny * dy - pt_damping_width ) )  THEN
@@ -1719 +1720 @@
           ENDIF
        ENDDO 
-    ELSEIF ( bc_ns_raddir  .OR.  bc_ns_neudir )  THEN
+    ELSEIF ( bc_ns_raddir )  THEN
        DO  j = nys, nyn
           IF ( ( j * dy ) < pt_damping_width )  THEN