- Timestamp:
- 2013-10-28T14:50:08+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3987_UKMO6_C1D/NEMOGCM/NEMO/OPA_SRC/par_GYRE.h90
r2715 r4144 8 8 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 9 9 !!---------------------------------------------------------------------- 10 CHARACTER (len=16) &10 CHARACTER (len=16) & 11 11 #if !defined key_agrif 12 , PARAMETER &12 , PARAMETER & 13 13 #endif 14 :: &15 cp_cfg = "gyre" !: name of the configuration16 INTEGER &14 :: & 15 cp_cfg = "gyre" !: Name of the configuration 16 INTEGER & 17 17 #if !defined key_agrif 18 , PARAMETER &18 , PARAMETER & 19 19 #endif 20 :: &21 jp_cfg = 1 , & !:20 :: & 21 jp_cfg = 1 , & !: Resolution of the configuration (degrees) 22 22 23 ! data size !!! * size of all the input files*24 jpidta = 30*jp_cfg+2, & !: 1st horizontal dimension ( >= jpi)25 jpjdta = 20*jp_cfg+2, & !: 2nd " " ( >= jpj)26 jpkdta = 31, & !: number of levels ( >= jpk)23 ! Data domain size !!! * Size of all input files * 24 jpidta = 30*jp_cfg+2, & !: 1st lateral dimension ( >= jpiglo ) 25 jpjdta = 20*jp_cfg+2, & !: 2nd lateral dimension ( >= jpjglo ) 26 jpkdta = 31 , & !: Number of levels ( >= jpk ) 27 27 28 ! global domain size !!! * full domain * 29 jpiglo = jpidta, & !: 1st dimension of global domain --> i30 jp jglo = jpjdta, & !: 2nd " " --> j31 ! zoom starting position32 jpizoom = 1 , & !: left bottom (i,j) indices of the zoom33 jp jzoom = 1 , & !: in data indices28 #if defined key_c1d 29 ! Zoom domain size !!! * C1D zoom * 30 jpiglo = 3 , & !: 1st dimension of global domain --> i 31 jpjglo = 3 , & !: 2nd dimension of global domain --> j 32 ! Domain characteristics 33 jperio = 0 !: Lateral cond. type (between 0 and 6) 34 34 35 INTEGER & 36 :: & 37 ! Starting position of the zoom 38 jpizoom = 1 , & !: Left bottom (i,j) indices of the zoom 39 jpjzoom = 1 !: in data domain indices 40 #else 41 ! Global domain size !!! * Global domain * 42 jpiglo = jpidta, & !: 1st dimension of global domain --> i 43 jpjglo = jpjdta, & !: 2nd dimension of global domain --> j 44 ! Starting position of the zoom 45 jpizoom = 1 , & !: Left bottom (i,j) indices of the zoom 46 jpjzoom = 1 , & !: in data domain indices 35 47 ! Domain characteristics 36 jperio = 0 !: lateral cond. type (between 0 and 6) 48 jperio = 0 !: Lateral cond. type (between 0 and 6) 49 #endif 50 37 51 38 52 !! Values set to pp_not_used indicates that this parameter is not used in THIS config. 39 53 !! Values set to pp_to_be_computed indicates that variables will be computed in domzgr 40 REAL(wp), PARAMETER :: & !: 41 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 = 5 !: 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 ! ! =5 beta-plane with regular grid-spacing and rotated domain (GYRE configuration) 54 55 REAL(wp), PARAMETER :: & 56 pp_not_used = 999999._wp , & !: 57 pp_to_be_computed = 0._wp !: 58 59 60 !! Coefficients associated with the horizontal coordinate system 61 62 INTEGER, PARAMETER :: & 63 jphgr_msh = 5 !: type of horizontal mesh 64 ! !: = 0 curvilinear coordinate on the sphere 65 ! !: read in coordinate.nc file 66 ! !: = 1 geographical mesh on the sphere 67 ! !: with regular grid-spacing 68 ! !: = 2 f-plane with regular grid-spacing 69 ! !: = 3 beta-plane with regular grid-spacing 70 ! !: = 4 Mercator grid with T/U point at the equator with 71 ! !: isotropic resolution (e1_deg) 72 ! !: = 5 beta-plane with regular grid-spacing and rotated domain (GYRE configuration) 57 73 58 74 ! ppglam0 , ppgphi0: coordinates of the lower leftmost T point of the grid. 59 75 ! The mercator grid starts only approximately at gphi0 because 60 76 ! of the constraint that the equator be a T point. 61 REAL(wp) & 77 78 REAL(wp) & 62 79 #if !defined key_agrif 63 , PARAMETER &80 , PARAMETER & 64 81 #endif 65 :: & !:66 ppglam0 = 0.0_wp, & !: longitude of first raw and column T-point(jphgr_msh = 1)67 ppgphi0 = 29.0_wp, & !: latitude of first raw and column T-point(jphgr_msh = 1)68 ! ! latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3)69 ppe1_deg = pp_not_used ,& !: zonal grid-spacing (degrees)70 ppe2_deg = pp_not_used ,& !: meridional grid-spacing (degrees)82 :: & 83 ppglam0 = 0.0_wp , & !: longitude of first raw and column T-point (jphgr_msh = 1) 84 ppgphi0 = 29.0_wp , & !: latitude of first raw and column T-point (jphgr_msh = 1) 85 ! !: latitude for the Coriolis or Beta parameter (jphgr_msh = 2 or 3) 86 ppe1_deg = pp_not_used, & !: zonal grid-spacing (degrees) 87 ppe2_deg = pp_not_used, & !: meridional grid-spacing (degrees) 71 88 ! 72 ppe1_m = pp_not_used, & !: zonal grid-spacing (meters ) 73 ppe2_m = pp_not_used !: meridional grid-spacing (meters ) 74 !! 89 ppe1_m = pp_not_used, & !: zonal grid-spacing (meters ) 90 ppe2_m = pp_not_used !: meridional grid-spacing (meters ) 91 92 75 93 !! Coefficients associated with the vertical coordinate system 76 !!77 94 78 REAL(wp), PARAMETER :: & !: 79 & ppsur = -2033.194295283385_wp , & !: 80 & ppa0 = 155.8325369664153_wp , & !: 81 & ppa1 = 146.3615918601890_wp , & !: 82 ! 83 & ppkth = 17.28520372419791_wp , & !: (non dimensional): gives the approximate 84 ! ! layer number above which stretching will 85 ! ! be maximum. Usually of order jpk/2. 86 & ppacr = 5.000000000000000_wp !: (non dimensional): stretching factor 87 ! ! for the grid. The highest zacr, the smallest 88 ! ! the stretching. 95 REAL(wp), PARAMETER :: & 96 & ppsur = -2033.194295283385_wp , & !: 97 & ppa0 = 155.8325369664153_wp , & !: 98 & ppa1 = 146.3615918601890_wp , & !: 99 & ppkth = 17.28520372419791_wp , & !: (non dimensional): gives the approximate 100 ! !: layer number above which stretching will 101 ! !: be maximum. Usually of order jpk/2. 102 & ppacr = 5.000000000000000_wp !: (non dimensional): stretching factor 103 ! !: for the grid. The higher zacr, the smaller 104 ! !: the stretching. 89 105 90 !! 91 !! If all ppa0 ppa1 and ppsur are specified to 0, then 92 !! they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 93 !! 94 REAL(wp), PARAMETER :: & !: 95 & ppdzmin = pp_not_used , & !: (meters): depth of the top (first) model layer 96 ! ! depth of second "w" level 97 & pphmax = pp_not_used !: (meters): maximum depth of the ocean 98 ! ! depth of the last "w" level 99 LOGICAL, PARAMETER :: & 106 ! If both ppa0 ppa1 and ppsur are specified to pp_to_be_computed, then 107 ! they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr 108 109 REAL(wp), PARAMETER :: & 110 & ppdzmin = pp_not_used , & !: (meters) vertical thickness of the top layer 111 & pphmax = pp_not_used !: (meters) Maximum depth of the ocean gdepw(jpk) 112 LOGICAL, PARAMETER :: & 100 113 & ldbletanh = .FALSE. !: Use/do not use double tanf function for vertical coordinates 101 REAL(wp), PARAMETER :: &114 REAL(wp), PARAMETER :: & 102 115 & ppa2 = pp_not_used , & !: Double tanh function parameters 103 116 & ppkth2 = pp_not_used , & !: 104 117 & ppacr2 = pp_not_used !: 105 106 118 !!---------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.