Changeset 19 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Feb 23, 2007 4:53:48 AM (17 years ago)

Author:

raasch

Message:

preliminary version of modified boundary conditions at top

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (6 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Temperature and humidity gradients at top are now calculated for nzt+1,
+! top_heatflux and respective boundary condition bc_pt_t is checked
 !
 ! Former revisions:
@@ -483 +484 @@
 !--    Store temperature gradient at the top boundary for possile Neumann
 !--    boundary condition
-       bc_pt_t_val = ( pt_init(nzt) - pt_init(nzt-1) ) / dzu(nzt)
+       bc_pt_t_val = ( pt_init(nzt+1) - pt_init(nzt) ) / dzu(nzt+1)
 
 !
@@ -539 +540 @@
 !--       Store humidity gradient at the top boundary for possile Neumann
 !--       boundary condition
-          bc_q_t_val = ( q_init(nzt) - q_init(nzt-1) ) / dzu(nzt)
+          bc_q_t_val = ( q_init(nzt+1) - q_init(nzt) ) / dzu(nzt+1)
 
        ENDIF
@@ -795 +796 @@
     ELSEIF ( bc_pt_t == 'neumann' )  THEN
        ibc_pt_t = 1
+    ELSEIF ( bc_pt_t == 'initial_gradient' )  THEN
+       ibc_pt_t = 2
     ELSE
        IF ( myid == 0 )  THEN
@@ -804 +807 @@
 
     IF ( surface_heatflux == 9999999.9 )  constant_heatflux = .FALSE.
+    IF ( top_heatflux     == 9999999.9 )  THEN
+       constant_top_heatflux = .FALSE.
+    ELSE
+       use_top_fluxes = .TRUE.    ! because this is currently the only choice
+    ENDIF
 
 !
@@ -823 +831 @@
           PRINT*, '    allowed with pt_surface_initial_change (/=0) = ', &
                   pt_surface_initial_change
+       ENDIF
+       CALL local_stop
+    ENDIF
+
+!
+!-- A given temperature at the top implies Dirichlet boundary condition for
+!-- temperature. In this case specification of a constant heat flux is
+!-- forbidden.
+    IF ( ibc_pt_t == 0  .AND.   constant_top_heatflux  .AND. &
+         top_heatflux /= 0.0 )  THEN
+       IF ( myid == 0 )  THEN
+          PRINT*, '+++ check_parameters:'
+          PRINT*, '    boundary_condition: bc_pt_t = ', bc_pt_t
+          PRINT*, '    is not allowed with constant_top_heatflux = .TRUE.'
        ENDIF
        CALL local_stop