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

Timestamp:

Jun 2, 2007 4:48:38 PM (17 years ago)

Author:

raasch

Message:

further preliminary uncomplete changes for ocean version

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (1 diff)

palm/trunk/SOURCE/check_parameters.f90

@@ -1005 +1005 @@
 
 !
+!-- Boundary conditions for salinity
+    IF ( ocean )  THEN
+       IF ( bc_sa_t == 'dirichlet' )  THEN
+          ibc_sa_t = 0
+       ELSEIF ( bc_sa_t == 'neumann' )  THEN
+          ibc_sa_t = 1
+       ELSE
+          IF ( myid == 0 )  THEN
+             PRINT*, '+++ check_parameters:'
+             PRINT*, '    unknown boundary condition: bc_sa_t = ', bc_sa_t
+          ENDIF
+          CALL local_stop
+       ENDIF
+
+       IF ( top_salinityflux == 9999999.9 )  constant_top_salinityflux = .FALSE.
+
+!
+!--    A fixed salinity at the top implies Dirichlet boundary condition for
+!--    salinity. In this case specification of a constant salinity flux is
+!--    forbidden.
+       IF ( ibc_sa_t == 0  .AND.   constant_top_salinityflux  .AND. &
+            top_salinityflux /= 0.0 )  THEN
+          IF ( myid == 0 )  THEN
+             PRINT*, '+++ check_parameters:'
+             PRINT*, '    boundary_condition: bc_sa_t = ', bc_sa_t
+             PRINT*, '    is not allowed with constant_top_salinityflux = ', &
+                          '.TRUE.'
+          ENDIF
+          CALL local_stop
+       ENDIF
+
+    ENDIF
+
+!
 !-- In case of humidity or passive scalar, set boundary conditions for total 
 !-- water content / scalar