post_initial.F

Go to the documentation of this file.

#include "cppdefs.h"
      MODULE post_initial_mod
#ifdef NONLINEAR
!
!git $Id$
!svn $Id$
!================================================== Hernan G. Arango ===
!  Copyright (c) 2002-2025 The ROMS Group                              !
!    Licensed under a MIT/X style license                              !
!    See License_ROMS.md                                               !
!=======================================================================
!                                                                      !
!  On the first timestep, it computes the initial depths and level     !
!  thicknesses from the initial free-surface field. Additionally, it   !
!  initializes the nonlinear state variables for all time levels and   !
!  applies lateral boundary conditions.                                !
!                                                                      !
!=======================================================================
!
      USE mod_param
      USE mod_parallel
# ifdef NESTING
      USE mod_scalars
# endif
!
      USE ini_fields_mod, ONLY : ini_fields, ini_zeta
# ifdef SOLVE3D
      USE set_depth_mod,  ONLY : set_depth
# endif
!
      implicit none
!
      PUBLIC  :: post_initial
      PRIVATE
!
      CONTAINS
!
!***********************************************************************
      SUBROUTINE post_initial (ng, model)
!***********************************************************************
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, model
!
!  Local variable declarations.
!
      integer :: tile
!
!-----------------------------------------------------------------------
!  Initialize free-surface and compute initial level thicknesses and
!  depths.
!-----------------------------------------------------------------------
!
      DO tile=first_tile(ng),last_tile(ng),+1
        CALL ini_zeta (ng, tile, model)
        CALL set_depth (ng, tile, model)
      END DO
!$OMP BARRIER
!
!-----------------------------------------------------------------------
!  Initialize other state variables.
!-----------------------------------------------------------------------
!
      DO tile=last_tile(ng),first_tile(ng),-1
        CALL ini_fields (ng, tile, model)
      END DO
!$OMP BARRIER
 
# ifdef NESTING
!
!-----------------------------------------------------------------------
!  Extract donor grid initial data at contact points and store it in
!  REFINED structure so it can be used for the space-time interpolation.
!-----------------------------------------------------------------------
!
      IF (refinedgrid(ng)) THEN
        CALL nesting (ng, model, ngetd)
      END IF
# endif
!
      RETURN
      END SUBROUTINE post_initial
#endif
      END MODULE post_initial_mod