Changeset 4272 for palm/trunk/SOURCE/ocean_mod.f90

Timestamp:

Oct 23, 2019 3:18:57 PM (5 years ago)

Author:

schwenkel

Message:

further modularization of boundary conditions: moving boundary conditions to their modules

File:

• palm/trunk/SOURCE/ocean_mod.f90 (modified) (9 diffs)

palm/trunk/SOURCE/ocean_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Further modularization of boundary conditions: moved boundary conditions to
+! respective modules
+!
+! 4196 2019-08-29 11:02:06Z gronemeier
 ! Consider rotation of model domain for calculating the Stokes drift
 ! 
@@ -81 +85 @@
         ONLY:  sums_wssas_ws_l
 
+    USE indices,                                                               &
+        ONLY:  advc_flags_s, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
+
+    USE surface_mod,                                                           &
+        ONLY:  bc_h, surf_def_v, surf_def_h, surf_lsm_h, surf_lsm_v,           &
+               surf_usm_h, surf_usm_v
 
     IMPLICIT NONE
@@ -192 +202 @@
     END INTERFACE ocean_prognostic_equations
 
+    INTERFACE ocean_boundary_conditions
+       MODULE PROCEDURE ocean_boundary_conditions
+    END INTERFACE ocean_boundary_conditions
+
     INTERFACE ocean_swap_timelevel
        MODULE PROCEDURE ocean_swap_timelevel
@@ -236 +250 @@
            ocean_swap_timelevel, ocean_rrd_global, ocean_rrd_local,            &
            ocean_wrd_global, ocean_wrd_local, ocean_3d_data_averaging,         &
-           stokes_drift_terms, wave_breaking_term
+           ocean_boundary_conditions, stokes_drift_terms, wave_breaking_term
 
 
@@ -257 +271 @@
     USE indices,                                                               &
         ONLY:  nxl, nxr, nyn, nys, nzb, nzt
-
-    USE surface_mod,                                                           &
-        ONLY :  bc_h
 
     IMPLICIT NONE
@@ -370 +381 @@
     USE indices,                                                               &
         ONLY:  nzb, nzt
-
-    USE surface_mod,                                                           &
-       ONLY :  bc_h
 
     IMPLICIT NONE
@@ -1585 +1593 @@
     USE diffusion_s_mod,                                                       &
         ONLY:  diffusion_s
-
-    USE indices,                                                               &
-        ONLY:  advc_flags_s, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
-
-    USE surface_mod,                                                           &
-        ONLY:  surf_def_v, surf_def_h, surf_lsm_h, surf_lsm_v, surf_usm_h,     &
-               surf_usm_v
 
     IMPLICIT NONE
@@ -1757 +1758 @@
     USE diffusion_s_mod,                                                       &
         ONLY:  diffusion_s
-
-    USE indices,                                                               &
-        ONLY:  advc_flags_s, nzb, nzt, wall_flags_0
-
-    USE surface_mod,                                                           &
-        ONLY:  surf_def_v, surf_def_h, surf_lsm_h, surf_lsm_v, surf_usm_h,     &
-               surf_usm_v
 
     IMPLICIT NONE
@@ -1862 +1856 @@
  END SUBROUTINE ocean_prognostic_equations_ij
 
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Boundary conditions for ocean model
+!------------------------------------------------------------------------------!
+ SUBROUTINE ocean_boundary_conditions
+
+    IMPLICIT NONE
+
+    INTEGER(iwp) ::  i                            !< grid index x direction.
+    INTEGER(iwp) ::  j                            !< grid index y direction.
+    INTEGER(iwp) ::  k                            !< grid index z direction.
+    INTEGER(iwp) ::  l                            !< running index boundary type, for up- and downward-facing walls.
+    INTEGER(iwp) ::  m                            !< running index surface elements.
+
+!
+!--    Boundary conditions for salinity
+!--    Bottom boundary: Neumann condition because salinity flux is always
+!--    given.
+       DO  l = 0, 1
+          !$OMP PARALLEL DO PRIVATE( i, j, k )
+          DO  m = 1, bc_h(l)%ns
+             i = bc_h(l)%i(m)
+             j = bc_h(l)%j(m)
+             k = bc_h(l)%k(m)
+             sa_p(k+bc_h(l)%koff,j,i) = sa_p(k,j,i)
+          ENDDO
+       ENDDO
+!
+!--    Top boundary: Dirichlet or Neumann
+       IF ( ibc_sa_t == 0 )  THEN
+           sa_p(nzt+1,:,:) = sa(nzt+1,:,:)
+       ELSEIF ( ibc_sa_t == 1 )  THEN
+           sa_p(nzt+1,:,:) = sa_p(nzt,:,:)
+       ENDIF
+
+ END SUBROUTINE ocean_boundary_conditions
 
 !------------------------------------------------------------------------------!