Changeset 778 for branches/dev_001_GM/NEMO/OPA_SRC/par_eel_r2.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_eel_r2.F90
r772 r778 1 MODULE par_eel_r2 2 !!====================================================================== 3 !! *** par_eel_r2 *** 4 !! Ocean Domain parameter : 2 km resolution cyclic Channel (EEL_R2 configuration) 5 !!====================================================================== 6 !! History : 9.0 ! 2007-11 (G. Madec) Original code from par_EEL_R2.h90 7 !!---------------------------------------------------------------------- 8 !! NEMO/OPA 2.4 , LOCEAN-IPSL (2007) 9 !! $Header: $ 10 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 11 !!---------------------------------------------------------------------- 12 #if defined key_eel_r2 1 13 !!--------------------------------------------------------------------- 2 !! *** par_EEL_R2.h90 *** 3 !! Ocean Domain : 2 km resolution Channel (EEL_R2 configuration) 14 !! 'key_eel_r2' : EEL configuration : R2 4 15 !!--------------------------------------------------------------------- 5 !!---------------------------------------------------------------------- 6 !! OPA 9.0 , LOCEAN-IPSL (2005) 7 !! $Header$ 8 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 9 !!---------------------------------------------------------------------- 10 CHARACTER (len=16) & 11 #if !defined key_agrif 12 , PARAMETER & 16 USE par_kind ! kind parameters 17 18 IMPLICIT NONE 19 PUBLIC 20 21 # if defined key_agrif 22 !------!!-------------------------------------------------------------- 23 ! !! 'key_agrif' : AGRIF 24 !------!!-------------------------------------------------------------- 25 26 ! Configuration, and domain size are not assigned to a parameter value 27 ! as they are different in mother and children models 28 INTEGER :: cp_cfg, jp_cfg !: configuration id (AGRIF case) 29 INTEGER :: jpidta, jpjdta, jpkdta !: input data size (AGRIF case) 30 INTEGER :: jpiglo, jpjglo, jpk, jpizoom, jpjzoom !: domain size and zoom (AGRIF case) 31 INTEGER :: jperio, jpisl, jpnisl !: domain characteristics (AGRIF case) 32 33 # else 34 !------!!-------------------------------------------------------------- 35 ! !! Default : no AGRIF: model parameters 36 !------!!-------------------------------------------------------------- 37 38 ! Configuration name and resolution !!! * GYRE R2 * 39 CHARACTER (len=16) :: cp_cfg = "eel" !: name of the configuration 40 INTEGER :: jp_cfg = 2 !: resolution of the configuration [km] 41 42 ! input data size !!! * size of all input files * 43 INTEGER, PARAMETER :: jpidta = 83 !: 1st lateral dimension ( >= jpiglo ) 44 INTEGER, PARAMETER :: jpjdta = 242 !: 2nd " " ( >= jpjglo ) 45 INTEGER, PARAMETER :: jpkdta = 30 !: number of levels ( >= jpk ) 46 47 48 # if defined key_cfg_1d 49 !------------!!-------------------------------------------------------- 50 ! !! 'key_cfg_1d' : 1D configuration of GYRE_R2 51 !------------!!-------------------------------------------------------- 52 ! global domain size !!! * 1D configuration * 53 INTEGER, PARAMETER :: jpiglo = 3 !: 1st dimension of global domain --> i 54 INTEGER, PARAMETER :: jpjglo = 3 !: 2nd " " --> j 55 INTEGER, PARAMETER :: jpk = jpkdta !: number of vertical levels 56 ! starting position of the zoom 57 INTEGER, PARAMETER :: jpizoom = 16 !: left bottom (i,j) indices of the zoom 58 INTEGER, PARAMETER :: jpjzoom = 11 !: in data domain indices 59 ! Domain characteristics 60 INTEGER, PARAMETER :: jperio = 0 !: lateral cond. type (between 0 and 6) 61 INTEGER, PARAMETER :: jpisl = 0 !: number of islands 62 INTEGER, PARAMETER :: jpnisl = 0 !: maximum number of points per island 63 64 # else 65 !------------!!-------------------------------------------------------- 66 ! !! Default : Full global domain of GYRE_R2 67 !------------!!-------------------------------------------------------- 68 ! global domain size !!! * global domain * 69 INTEGER, PARAMETER :: jpiglo = jpidta !: 1st dimension of global domain --> i 70 INTEGER, PARAMETER :: jpjglo = jpjdta !: 2nd " " --> j 71 INTEGER, PARAMETER :: jpk = jpkdta !: number of vertical levels 72 ! starting position of the zoom 73 INTEGER, PARAMETER :: jpizoom = 1 !: left bottom (i,j) indices of the zoom 74 INTEGER, PARAMETER :: jpjzoom = 1 !: in data domain indices 75 ! Domain characteristics 76 INTEGER, PARAMETER :: jperio = 1 !: lateral cond. type (between 0 and 6) 77 INTEGER, PARAMETER :: jpisl = 1 !: number of islands 78 INTEGER, PARAMETER :: jpnisl = jpiglo !: maximum number of points per island 79 # endif 80 81 # endif 82 83 !!--------------------------------------------------------------------- 84 !! grid parameters for domhgr and domzgr 85 !!--------------------------------------------------------------------- 86 ! Values set to pp_not_used indicates that this parameter is not used in THIS config. 87 ! Values set to pp_to_be_computed indicates that variables will be computed in domzgr 88 REAL(wp), PARAMETER :: pp_not_used = 999999_wp !: default value for "not used" 89 REAL(wp), PARAMETER :: pp_to_be_computed = 0._wp !: default value for "to be computed" 90 91 ! horizontal mesh (domhgr) 92 ! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 93 INTEGER, PARAMETER :: jphgr_msh = 3 !: type of horizontal mesh 94 ! ! = 3 beta-plane with regular grid-spacing 95 96 ! ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 97 ! The mercator grid starts only approximately at gphi0 because 98 ! of the constraint that the equator be a T point. 99 REAL(wp) ,PARAMETER :: ppglam0 = 0.0_wp !: longitude of first raw and column T-point (jphgr_msh=1) 100 REAL(wp) ,PARAMETER :: ppgphi0 = 35.0_wp !: latitude of first raw and column T-point (jphgr_msh=1) 101 ! ! latitude of the Coriolis/Beta parameter (jphgr_msh=2 or 3) 102 REAL(wp) ,PARAMETER :: ppe1_deg = pp_not_used !: zonal grid-spacing (degrees) 103 REAL(wp) ,PARAMETER :: ppe2_deg = pp_not_used !: meridional grid-spacing (degrees) 104 ! 105 REAL(wp) ,PARAMETER :: ppe1_m = 2000.0_wp !: zonal grid-spacing (meters ) 106 REAL(wp) ,PARAMETER :: ppe2_m = 2000.0_wp !: meridional grid-spacing (meters ) 107 108 ! Vertical mesh (domzgr) 109 ! parameter used in the analytical function defining the repartition of the level and their thickness 110 REAL(wp), PARAMETER :: ppsur = -2033.194295283385_wp !: ORCA r4, r2 and r05 coefficients 111 REAL(wp), PARAMETER :: ppa0 = 155.8325369664153_wp !: (default coefficients) 112 REAL(wp), PARAMETER :: ppa1 = 146.3615918601890_wp !: 113 REAL(wp), PARAMETER :: ppkth = 17.28520372419791_wp !: (non dimensional): gives the approximate 114 ! !: layer number above which stretching will 115 ! !: be maximum. Usually of order jpk/2. 116 REAL(wp), PARAMETER :: ppacr = 5.00000000000_wp !: (non dimensional): stretching factor 117 ! !: for the grid. The highest zacr, the smallest 118 ! !: the stretching. 119 ! If both ppa0 ppa1 and ppsur are specified to 0, then 120 ! they are computed from ppdzmin, pphmax, ppkth, ppacr in dom_zgr 121 REAL(wp), PARAMETER :: ppdzmin = pp_not_used !: (meters) vertical thickness of the top layer 122 REAL(wp), PARAMETER :: pphmax = pp_not_used !: (meters) Maximum depth of the ocean gdepw(jpk) 123 124 #else 125 !!--------------------------------------------------------------------- 126 !! default option : Empty module 127 !!--------------------------------------------------------------------- 13 128 #endif 14 :: &15 cp_cfg = "eel" !: name of the configuration16 INTEGER &17 #if !defined key_agrif18 , PARAMETER &19 #endif20 :: &21 jp_cfg = 2 , & !: resolution of the configuration (km)22 129 23 ! data size !!! * size of all the input files * 24 jpidta = 83, & !: 1st horizontal dimension ( >= jpi ) 25 jpjdta = 242, & !: 2nd " " ( >= jpj ) 26 jpkdta = 30, & !: number of levels ( >= jpk ) 27 28 ! global domain size !!! * full domain * 29 jpiglo = jpidta, & !: 1st dimension of global domain --> i 30 jpjglo = jpjdta, & !: 2nd " " --> j 31 jpk = jpkdta, & !: number of vertical levels 32 ! zoom starting position 33 jpizoom = 1 , & !: left bottom (i,j) indices of the zoom 34 jpjzoom = 1 , & !: in data indices 35 36 ! Domain characteristics 37 jperio = 1 , & !: lateral cond. type (between 0 and 6) 38 jpisl = 1 , & !: number of islands 39 jpnisl = jpiglo !: maximum number of points per island 40 41 !! Values set to pp_not_used indicates that this parameter is not used in THIS config. 42 !! Values set to pp_to_be_computed indicates that variables will be computed in domzgr 43 REAL(wp), PARAMETER :: & !: 44 pp_not_used = 999999._wp , & !: ??? 45 pp_to_be_computed = 0._wp !: ??? 46 !! 47 !! Coefficients associated with the horizontal coordinate system (jphgr_msh /= 0 ) 48 !! 49 INTEGER,PARAMETER :: & !: 50 jphgr_msh = 3 !: type of horizontal mesh 51 ! ! = 0 curvilinear coordinate on the sphere 52 ! ! read in coordinate.nc file 53 ! ! = 1 geographical mesh on the sphere 54 ! ! with regular grid-spacing 55 ! ! = 2 f-plane with regular grid-spacing 56 ! ! = 3 beta-plane with regular grid-spacing 57 ! ! = 4 Mercator grid with T/U point at the equator with 58 ! ! isotropic resolution (e1_deg) 59 60 ! ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 61 ! The mercator grid starts only approximately at gphi0 because 62 ! of the constraint that the equator be a T point. 63 REAL(wp) ,PARAMETER :: & !: 64 ppglam0 = 0.0_wp, & !: longitude of first raw and column T-point (jphgr_msh = 1) 65 ppgphi0 = 35.0_wp, & !: latitude of first raw and column T-point (jphgr_msh = 1) 66 ! ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 67 ppe1_deg = pp_not_used , & !: zonal grid-spacing (degrees) 68 ppe2_deg = pp_not_used , & !: meridional grid-spacing (degrees) 69 ! 70 ppe1_m = 2000.0_wp, & !: zonal grid-spacing (meters ) 71 ppe2_m = 2000.0_wp !: meridional grid-spacing (meters ) 72 !! 73 !! Coefficients associated with the vertical coordinate system 74 !! 75 REAL(wp) & 76 #if !defined key_agrif 77 , PARAMETER & 78 #endif 79 :: & !: 80 & ppsur = -2033.194295283385_wp , & !: Computed in domzgr 81 & ppa0 = 155.8325369664153_wp , & !: 82 & ppa1 = 146.3615918601890_wp , & !: 83 ! 84 & ppkth = 17.28520372419791_wp , & !: (non dimensional): gives the approximate 85 ! !: layer number above which stretching will 86 ! !: be maximum. Usually of order jpk/2. 87 & ppacr = 5.000000000000000_wp !: (non dimensional): stretching factor 88 ! !: for the grid. The highest zacr, the smallest 89 ! !: the stretching. 90 91 !! 92 !! If both ppa0 ppa1 and ppsur are specified to 0, then 93 !! they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 94 !! 95 REAL(wp), PARAMETER :: & 96 & ppdzmin = pp_not_used , & !: (meters) vertical thickness of the top layer 97 & pphmax = pp_not_used !: (meters) Maximum depth of the ocean gdepw(jpk) 98 !!--------------------------------------------------------------------- 130 !!====================================================================== 131 END MODULE par_eel_r2
Note: See TracChangeset
for help on using the changeset viewer.