New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
usrdef_zgr.F90 in NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/MY_SRC – NEMO

source: NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/MY_SRC/usrdef_zgr.F90 @ 11637

Last change on this file since 11637 was 11637, checked in by laurent, 5 years ago

LB: preliminary inclusion of "STATION_ASF" test-case!
File size: 4.4 KB
Line 
1MODULE usrdef_zgr
2   !!======================================================================
3   !!                     ***  MODULE usrdef_zgr  ***
4   !!
5   !!                       ===  STATION_ASF case  ===
6   !!
7   !! user defined :  vertical coordinate system of a user configuration
8   !!======================================================================
9   !! History :  4.0  ! 2019-10  (L. Brodeau)  Original code
10   !!----------------------------------------------------------------------
11
12   !!----------------------------------------------------------------------
13   !!   usr_def_zgr   : user defined vertical coordinate system (required)
14   !!---------------------------------------------------------------------
15   USE oce            ! ocean variables
16   !USE dom_oce        ! ocean domain
17   !USE depth_e3       ! depth <=> e3
18   USE usrdef_nam     ! User defined : namelist variables
19   !
20   USE in_out_manager ! I/O manager
21   USE lbclnk         ! ocean lateral boundary conditions (or mpp link)
22   USE lib_mpp        ! distributed memory computing library
23   USE timing         ! Timing
24
25   IMPLICIT NONE
26   PRIVATE
27
28   PUBLIC   usr_def_zgr   ! called by domzgr.F90
29
30   !! * Substitutions
31#  include "vectopt_loop_substitute.h90"
32   !!----------------------------------------------------------------------
33   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
34   !! $Id: usrdef_zgr.F90 10072 2018-08-28 15:21:50Z nicolasmartin $
35   !! Software governed by the CeCILL license (see ./LICENSE)
36   !!----------------------------------------------------------------------
37CONTAINS
38
39   SUBROUTINE usr_def_zgr( ld_zco  , ld_zps  , ld_sco  , ld_isfcav,    &   ! type of vertical coordinate
40      &                    pdept_1d, pdepw_1d, pe3t_1d , pe3w_1d  ,    &   ! 1D reference vertical coordinate
41      &                    pdept , pdepw ,                             &   ! 3D t & w-points depth
42      &                    pe3t  , pe3u  , pe3v   , pe3f ,             &   ! vertical scale factors
43      &                    pe3w  , pe3uw , pe3vw         ,             &   !     -      -      -
44      &                    k_top  , k_bot    )                             ! top & bottom ocean level
45      !!---------------------------------------------------------------------
46      !!              ***  ROUTINE usr_def_zgr  ***
47      !!
48      !! ** Purpose :   User defined the vertical coordinates
49      !!
50      !!----------------------------------------------------------------------
51      LOGICAL                   , INTENT(  out) ::   ld_zco, ld_zps, ld_sco      ! vertical coordinate flags ( read in namusr_def )
52      LOGICAL                   , INTENT(  out) ::   ld_isfcav                   ! under iceshelf cavity flag
53      REAL(wp), DIMENSION(:)    , INTENT(  out) ::   pdept_1d, pdepw_1d          ! 1D grid-point depth     [m]
54      REAL(wp), DIMENSION(:)    , INTENT(  out) ::   pe3t_1d , pe3w_1d           ! 1D grid-point depth     [m]
55      REAL(wp), DIMENSION(:,:,:), INTENT(  out) ::   pdept, pdepw                ! grid-point depth        [m]
56      REAL(wp), DIMENSION(:,:,:), INTENT(  out) ::   pe3t , pe3u , pe3v , pe3f   ! vertical scale factors  [m]
57      REAL(wp), DIMENSION(:,:,:), INTENT(  out) ::   pe3w , pe3uw, pe3vw         ! i-scale factors
58      INTEGER , DIMENSION(:,:)  , INTENT(  out) ::   k_top, k_bot                ! first & last ocean level
59      !!----------------------------------------------------------------------
60      !
61      IF(lwp) WRITE(numout,*)
62      IF(lwp) WRITE(numout,*) 'usr_def_zgr : LOLO STATION_ASF configuration, setting first level properties.'
63      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~'
64      !
65
66      ld_zco    = .TRUE.         ! z-coordinate without ocean cavities
67      ld_zps    = .FALSE.
68      ld_sco    = .FALSE.
69      ld_isfcav = .FALSE.
70     
71      pdept_1d(1) = rn_dept1 ! depth (m) at which the SST is taken/measured == depth of first T point!
72      pdepw_1d(1) = 0._wp
73      pe3t_1d(1)  = 2._wp*rn_dept1
74      pe3w_1d(1)  = rn_dept1 ! LB???
75
76      pdept(:,:,:) = rn_dept1
77      pdepw(:,:,:) = 0._wp
78      pe3t(:,:,:) = 2._wp*rn_dept1
79      pe3u(:,:,:) = 2._wp*rn_dept1
80      pe3v(:,:,:) = 2._wp*rn_dept1
81      pe3f(:,:,:) = 2._wp*rn_dept1
82      pe3w(:,:,:)  = rn_dept1  ! LB???
83      pe3uw(:,:,:) = rn_dept1  ! LB???
84      pe3vw(:,:,:) = rn_dept1  ! LB???
85      k_top = 1
86      k_bot = 1
87      !
88   END SUBROUTINE usr_def_zgr
89   !!======================================================================
90END MODULE usrdef_zgr
Note: See TracBrowser for help on using the repository browser.