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.
par_orca_r4.F90 in branches/dev_001_GM/NEMO/OPA_SRC – NEMO

source: branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r4.F90 @ 778

Last change on this file since 778 was 778, checked in by gm, 16 years ago

dev_001_GM - par_oce evolution: suppress the par_...h90 and create new module for all ocean configuration -compilation OK
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 8.2 KB
Line 
1MODULE par_orca_r4
2   !!======================================================================
3   !!                        ***  par_ORCA_R4  ***
4   !! Ocean Domain parameter : 4 degrees orca global ocean (0RCA_R4 configuration)
5   !!    including a 1D configuration case.
6   !!    NB: with AGRIF, values are not set as parameters
7   !!======================================================================
8   !! History :   9.0  !  2007-11  (G. Madec)  Original code from par_ORCA_R4.h90
9   !!----------------------------------------------------------------------
10   !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007)
11   !! $Id$
12   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
13   !!----------------------------------------------------------------------
14#if defined key_orca_r4
15   !!---------------------------------------------------------------------
16   !!   'key_orca_r4'   :                           global ocean : ORCA R4
17   !!---------------------------------------------------------------------
18   USE par_kind          ! kind parameters
19
20   IMPLICIT NONE
21   PUBLIC
22
23# if defined key_agrif
24   !------!!--------------------------------------------------------------
25   !      !!   'key_agrif'   :                                       AGRIF
26   !------!!--------------------------------------------------------------
27
28   ! Configuration, and domain size are not assigned to a parameter value
29   !    as they are different in mother and children models
30   INTEGER ::   cp_cfg, jp_cfg                          !: configuration id       (AGRIF case)
31   INTEGER ::   jpidta, jpjdta, jpkdta                  !: input data size        (AGRIF case)
32   INTEGER ::   jpiglo, jpjglo, jpk, jpizoom, jpjzoom   !: domain size and zoom   (AGRIF case)
33   INTEGER ::   jperio, jpisl, jpnisl                   !: domain characteristics (AGRIF case)
34
35# else
36   !------!!--------------------------------------------------------------
37   !      !!   Default       :                  no AGRIF: model parameters
38   !------!!--------------------------------------------------------------
39
40   ! Configuration name and resolution                 !!! *  ORCA R4  *
41   CHARACTER (len=16), PARAMETER ::   cp_cfg = "orca"   !: name of the configuration
42   INTEGER, PARAMETER            ::   jp_cfg = 4        !: resolution of the configuration (degrees)
43
44   ! input data size                                   !!! *  size of all input files  *
45   INTEGER, PARAMETER ::   jpidta  =  92                !: 1st lateral dimension ( >= jpiglo )
46   INTEGER, PARAMETER ::   jpjdta  =  76                !: 2nd    "       "      ( >= jpjglo )
47   INTEGER, PARAMETER ::   jpkdta  =  31                !: number of levels      ( >= jpk    )
48
49#  if defined key_cfg_1d
50   !------------!!--------------------------------------------------------
51   !            !!   'key_cfg_1d' :            1D configuration of ORCA R4
52   !------------!!--------------------------------------------------------
53   ! global domain size                                !!! *  1D configuration  *
54   INTEGER, PARAMETER ::   jpiglo  = 3                  !: 1st dimension of global domain --> i
55   INTEGER, PARAMETER ::   jpjglo  = 3                  !: 2nd    "                  "    --> j
56   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
57   ! starting position of the zoom
58   INTEGER, PARAMETER ::   jpizoom =   61               !: left bottom (i,j) indices of the zoom in data domain indices
59   INTEGER, PARAMETER ::   jpjzoom =   30               !: here : (jpizoom,jpjzoom) = (61,30)  for (??W,??N)
60   ! Domain characteristics
61   INTEGER, PARAMETER ::   jperio  =   0                !: lateral cond. type (between 0 and 6)
62   INTEGER, PARAMETER ::   jpisl   =   1                !: number of islands
63   INTEGER, PARAMETER ::   jpnisl  =   1                !: maximum number of points per island
64
65#  else
66   !------------!!--------------------------------------------------------
67   !            !!   Default :               Full global domain of ORCA R4
68   !------------!!--------------------------------------------------------
69   ! global domain size                                !!! *  global domain  *
70   INTEGER, PARAMETER ::   jpiglo  = jpidta             !: 1st dimension of global domain --> i
71   INTEGER, PARAMETER ::   jpjglo  = jpjdta             !: 2nd    "                  "    --> j
72   INTEGER, PARAMETER ::   jpk     = jpkdta             !: number of vertical levels
73   ! starting position of the zoom
74   INTEGER, PARAMETER ::   jpizoom =   1                !: left bottom (i,j) indices of the zoom
75   INTEGER, PARAMETER ::   jpjzoom =   1                !: in data domain indices
76   ! Domain characteristics
77   INTEGER, PARAMETER ::   jperio  =   4                !: lateral cond. type (between 0 and 6)
78   INTEGER, PARAMETER ::   jpisl   =  11                !: number of islands
79   INTEGER, PARAMETER ::   jpnisl  = 400                !: maximum number of points per island
80#  endif
81
82# endif
83
84   !!---------------------------------------------------------------------
85   !! grid parameters for domhgr and domzgr
86   !!---------------------------------------------------------------------
87   ! Values set to pp_not_used indicates that this parameter is not used in THIS config.
88   ! Values set to pp_to_be_computed  indicates that variables will be computed in domzgr
89   REAL(wp), PARAMETER ::   pp_not_used       = 999999_wp   !: default value for "not used"
90   REAL(wp), PARAMETER ::   pp_to_be_computed = 0._wp       !: default value for "to be computed"
91
92   ! horizontal mesh (domhgr)
93   !    Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 )
94   INTEGER, PARAMETER ::   jphgr_msh = 0   !: type of horizontal mesh
95   !                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc file
96
97   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T-point of the grid.
98   !   The mercator grid starts only approximately at gphi0 because
99   !   of the constraint that the equator is a T-point.
100   REAL(wp) ,PARAMETER ::   ppglam0  = pp_not_used   !: longitude of first raw and column T-point (jphgr_msh=1)
101   REAL(wp) ,PARAMETER ::   ppgphi0  = pp_not_used   !: latitude  of first raw and column T-point (jphgr_msh=1)
102   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
103   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees)
104   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees)
105   !
106   REAL(wp) ,PARAMETER ::   ppe1_m   = pp_not_used   !: zonal      grid-spacing (meters )
107   REAL(wp) ,PARAMETER ::   ppe2_m   = pp_not_used   !: meridional grid-spacing (meters )
108
109   ! Vertical mesh (domzgr)
110   !    parameter used in the analytical function defining the repartition of the level and their thickness
111   REAL(wp), PARAMETER ::   ppsur = -4762.96143546300_wp   !: ORCA r4, r2 and r05 coefficients
112   REAL(wp), PARAMETER ::   ppa0  =   255.58049070440_wp   !: (default coefficients)
113   REAL(wp), PARAMETER ::   ppa1  =   245.58132232490_wp   !:
114   REAL(wp), PARAMETER ::   ppkth =    21.43336197938_wp   !: (non dimensional): gives the approximate
115   !                                                       !: layer number above which  stretching will
116   !                                                       !: be maximum. Usually of order jpk/2.
117   REAL(wp), PARAMETER ::   ppacr =     3.00000000000_wp   !: (non dimensional): stretching factor
118   !                                                       !: for the grid. The highest zacr, the smallest
119   !                                                       !: the stretching.
120   !  If both ppa0 ppa1 and ppsur are specified to 0, then
121   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
122   REAL(wp), PARAMETER ::   ppdzmin = pp_not_used   !: (meters) vertical thickness of the top layer
123   REAL(wp), PARAMETER ::   pphmax  = pp_not_used   !: (meters) Maximum depth of the ocean gdepw(jpk)
124
125#else
126   !!---------------------------------------------------------------------
127   !!   default option  :                                     Empty module
128   !!---------------------------------------------------------------------
129#endif
130
131   !!======================================================================
132END MODULE par_orca_r4
Note: See TracBrowser for help on using the repository browser.