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

source: branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r025.F90 @ 4373

Last change on this file since 4373 was 778, checked in by gm, 17 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.6 KB
Line 
1MODULE par_orca_r025
2   !!======================================================================
3   !!                        ***  par_ORCA_R025  ***
4   !! Ocean Domain parameter : 0.25 degrees ORCA global ocean (0RCA_R025 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_R025.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_r025
15   !!---------------------------------------------------------------------
16   !!   'key_orca_r025'   :                       global ocean : ORCA R025
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 R025  *
41   CHARACTER (len=16), PARAMETER ::   cp_cfg = "orca"   !: name of the configuration
42   INTEGER           , PARAMETER ::   jp_cfg = 025      !: resolution of the configuration (degrees)
43
44   ! input data size                                   !!! *  size of all input files  *
45   INTEGER, PARAMETER ::   jpidta  = 1442               !: 1st lateral dimension ( >= jpiglo )
46   INTEGER, PARAMETER ::   jpjdta  = 1021               !: 2nd    "       "      ( >= jpjglo )
47   INTEGER, PARAMETER ::   jpkdta  =   46               !: number of levels      ( >= jpk    )
48   
49#  if defined key_cfg_1d
50   !------------!!--------------------------------------------------------
51   !            !!   'key_cfg_1d' :          1D configuration of ORCA R025
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 =  110               !: here : (jpizoom,jpjzoom) = (61,110)  for (??W,??N)
60   ! Domain characteristics
61   INTEGER, PARAMETER ::   jperio  =   0                !: lateral cond. type (between 0 and 6)
62   INTEGER, PARAMETER ::   jpisl   =  18                !: number of islands
63   INTEGER, PARAMETER ::   jpnisl  = 800                !: maximum number of points per island
64   
65#  else
66   !------------!!--------------------------------------------------------
67   !            !!   Default :             Full global domain of ORCA R025
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   =  150              !: number of islands
79   INTEGER, PARAMETER ::   jpnisl  = 3000              !: 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   !                                       !  = 1 geographical mesh on the sphere with regular grid-spacing
97   !                                       !  = 2 f-plane with regular grid-spacing
98   !                                       !  = 3 beta-plane with regular grid-spacing
99   !                                       !  = 4 Mercator grid with T/U point at the equator
100
101   !   ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid.
102   !   The mercator grid starts only approximately at gphi0 because
103   !   of the constraint that the equator be a T point.
104   REAL(wp) ,PARAMETER ::   ppglam0  = pp_not_used   !: longitude of first raw and column T-point (jphgr_msh=1)
105   REAL(wp) ,PARAMETER ::   ppgphi0  = pp_not_used   !: latitude  of first raw and column T-point (jphgr_msh=1)
106   !                                                 !  latitude  of the Coriolis/Beta parameter  (jphgr_msh=2 or 3)
107   REAL(wp) ,PARAMETER ::   ppe1_deg = pp_not_used   !: zonal      grid-spacing (degrees)
108   REAL(wp) ,PARAMETER ::   ppe2_deg = pp_not_used   !: meridional grid-spacing (degrees)
109   !
110   REAL(wp) ,PARAMETER ::   ppe1_m   = pp_not_used   !: zonal      grid-spacing (meters )
111   REAL(wp) ,PARAMETER ::   ppe2_m   = pp_not_used   !: meridional grid-spacing (meters )
112
113   ! Vertical mesh (domzgr)
114   !    parameter used in the analytical function defining the repartition of the level and their thickness
115   REAL(wp), PARAMETER ::   ppsur = pp_to_be_computed      !: ORCA r4, r2 and r05 coefficients
116   REAL(wp), PARAMETER ::   ppa0  = pp_to_be_computed      !: (default coefficients)
117   REAL(wp), PARAMETER ::   ppa1  = pp_to_be_computed      !:
118   REAL(wp), PARAMETER ::   ppkth =    23.563_wp           !: (non dimensional): gives the approximate
119   !                                                       !: layer number above which  stretching will
120   !                                                       !: be maximum. Usually of order jpk/2.
121   REAL(wp), PARAMETER ::   ppacr =     9.00000000000_wp   !: (non dimensional): stretching factor
122   !                                                       !: for the grid. The highest zacr, the smallest
123   !                                                       !: the stretching.
124   !  If both ppa0 ppa1 and ppsur are specified to 0, then
125   !  they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr
126   REAL(wp), PARAMETER ::   ppdzmin =    6._wp             !: (meters) vertical thickness of the top layer
127   REAL(wp), PARAMETER ::   pphmax  = 5750._wp             !: (meters) Maximum depth of the ocean gdepw(jpk)
128
129#else
130   !!---------------------------------------------------------------------
131   !!   default option  :                                     Empty module
132   !!---------------------------------------------------------------------
133#endif
134
135   !!======================================================================
136END MODULE par_orca_r025
Note: See TracBrowser for help on using the repository browser.