Changeset 778 for branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r025.F90
- Timestamp:
- 2007-12-22T11:30:56+01:00 (16 years ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_GM/NEMO/OPA_SRC/par_orca_r025.F90
r772 r778 1 MODULE 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 1 15 !!--------------------------------------------------------------------- 2 !! *** par_ORCA_R025.h90 *** 3 !! Ocean Domain : 0.25 degrees resolution global ocean 4 !! (0RCA_R025 configuration) 16 !! 'key_orca_r025' : global ocean : ORCA R025 5 17 !!--------------------------------------------------------------------- 6 !!---------------------------------------------------------------------- 7 !! OPA 9.0 , LOCEAN-IPSL (2005) 8 !! $Header$ 9 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 10 !!---------------------------------------------------------------------- 11 CHARACTER (len=16) & 12 #if !defined key_agrif 13 , PARAMETER & 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 !!--------------------------------------------------------------------- 14 133 #endif 15 :: &16 cp_cfg = "orca" !: name of the configuration17 INTEGER &18 #if !defined key_agrif19 , PARAMETER &20 #endif21 :: &22 jp_cfg = 025 , & !: resolution of the configuration (degrees)23 ! Original data size24 jpidta = 1442, & !: first horizontal dimension > or = to jpi25 jpjdta = 1021, & !: second > or = to jpj26 jpkdta = 46 , & !: number of levels > or = to jpk27 ! total domain matrix size28 jpiglo = jpidta, & !: first dimension of global domain --> i29 jpjglo = jpjdta, & !: second dimension of global domain --> j30 jpk = jpkdta, & !: number of vertical levels31 ! starting position of the zoom32 jpizoom = 1 , & !: left bottom (i,j) indices of the zoom33 jpjzoom = 1 , & !: in data indices34 ! Domain characteristics35 jperio = 4 , & !: lateral cond. type (between 0 and 6)36 jpisl = 150 , & !: number of islands37 jpnisl = 3000 !: maximum number of points per island38 134 39 !! Values set to pp_not_used indicates that this parameter is not used in THIS config. 40 !! Values set to pp_to_be_computed indicates that variables will be computed in domzgr 41 REAL,PARAMETER :: pp_not_used = 999999_wp , & 42 & pp_to_be_computed = 0._wp 43 !! 44 !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 45 !! 46 INTEGER, PARAMETER :: & ! 47 jphgr_msh = 0 !: type of horizontal mesh 48 ! ! = 0 curvilinear coordinate on the sphere 49 ! ! read in coordinate.nc file 50 ! ! = 1 geographical mesh on the sphere 51 ! ! with regular grid-spacing 52 ! ! = 2 f-plane with regular grid-spacing 53 ! ! = 3 beta-plane with regular grid-spacing 54 ! ! = 4 Mercator grid with T/U point at the equator with 55 ! ! isotropic resolution (e1_deg) 56 57 ! ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 58 ! The mercator grid starts only approximately at gphi0 because 59 ! of the constraint that the equator be a T point. 60 REAL(wp), PARAMETER :: & ! 61 ppglam0 = pp_not_used, & !: longitude of first raw and column T-point (jphgr_msh = 1) 62 ppgphi0 = pp_not_used, & !: latitude of first raw and column T-point (jphgr_msh = 1) 63 ! ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 64 ppe1_deg = pp_not_used, & !: zonal grid-spacing (degrees) 65 ppe2_deg = pp_not_used, & !: meridional grid-spacing (degrees) 66 ! 67 ppe1_m = pp_not_used, & !: zonal grid-spacing (meters ) 68 ppe2_m = pp_not_used !: meridional grid-spacing (meters ) 69 70 !! Coefficients associated with the vertical coordinate system 71 72 REAL(wp), PARAMETER :: & 73 & ppsur = pp_to_be_computed , & !: Computed in domzgr, set ppdzmin and pphmax below 74 & ppa0 = pp_to_be_computed , & !: " " 75 & ppa1 = pp_to_be_computed , & !: " " 76 ! 77 & ppkth = 23.563_wp , & !: (non dimensional): gives the approximate 78 ! !: layer number above which stretching will 79 ! !: be maximum. Usually of order jpk/2. 80 & ppacr = 9.00000000000_wp !: (non dimensional): stretching factor 81 ! !: for the grid. The highest zacr, the smallest 82 ! !: the stretching. 83 84 !! 85 !! If both ppa0 ppa1 and ppsur are specified to 0, then 86 !! they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 87 !! 88 REAL(wp), PARAMETER :: & 89 & ppdzmin = 6._wp , & !: (meters) vertical thickness of the top layer 90 & pphmax = 5750._wp !: (meters) Maximum depth of the ocean gdepw(jpk) 91 !!--------------------------------------------------------------------- 135 !!====================================================================== 136 END MODULE par_orca_r025
Note: See TracChangeset
for help on using the changeset viewer.
