- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/TOOLS/OBSTOOLS/src/par_oce.F90
r3000 r6808 25 25 INTEGER, PUBLIC, PARAMETER :: jprecj = 1 !: number of rows for overlap 26 26 27 !! Ocean Domain sizes 28 !! ------------------ 29 !! data domain (jpidta,jpjdta) 30 !! global or zoom domain (jpiglo,jpjglo) 31 !! local domain ( jpi , jpj ) 32 33 !#if defined key_orca_r4 34 ! !!--------------------------------------------------------------------- 35 ! !! 'key_orca_r4' : global ocean : ORCA R4 36 ! !!--------------------------------------------------------------------- 37 !# include "par_ORCA_R4.h90" 38 !#elif defined key_orca_r2 39 ! !!--------------------------------------------------------------------- 40 ! !! 'key_orca_r2' : global ocean : ORCA R4 41 ! !!--------------------------------------------------------------------- 42 !# include "par_ORCA_R2.h90" 43 !#elif defined key_orca_r1 44 ! !!--------------------------------------------------------------------- 45 ! !! 'key_orca_r1' : global ocean : ORCA R1 46 ! !!--------------------------------------------------------------------- 47 !# include "par_ORCA_R1.h90" 48 !#elif defined key_orca_r05 49 ! !!--------------------------------------------------------------------- 50 ! !! 'key_orca_r05' : global ocean : ORCA R05 51 ! !!--------------------------------------------------------------------- 52 !# include "par_ORCA_R05.h90" 53 !#elif defined key_orca_r025 54 ! !!--------------------------------------------------------------------- 55 ! !! 'key_orca_r025' : global ocean : ORCA R025 56 ! !!--------------------------------------------------------------------- 57 !# include "par_ORCA_R025.h90" 58 !#elif defined key_eel_r2 59 ! !!--------------------------------------------------------------------- 60 ! !! 'key_eel_r2' : channel : EEL R2 61 ! !!--------------------------------------------------------------------- 62 !# include "par_EEL_R2.h90" 63 !#elif defined key_eel_r5 64 ! !!--------------------------------------------------------------------- 65 ! !! 'key_eel_r5' : channel : EEL R5 66 ! !!--------------------------------------------------------------------- 67 !# include "par_EEL_R5.h90" 68 !#elif defined key_eel_r6 69 ! !!--------------------------------------------------------------------- 70 ! !! 'key_eel_r6' : channel : EEL R6 71 ! !!--------------------------------------------------------------------- 72 !# include "par_EEL_R6.h90" 73 !#elif defined key_gyre 74 ! !!--------------------------------------------------------------------- 75 ! !! 'key_gyre' : mid-latitude basin : GYRE 76 ! !!--------------------------------------------------------------------- 77 !# include "par_GYRE.h90" 78 !#elif defined key_pomme_r025 79 ! !!--------------------------------------------------------------------- 80 ! !! 'key_pomme_r025': regional basin : POMME025 81 ! !!--------------------------------------------------------------------- 82 !# include "par_POMME_R025.h90" 83 !#else 84 !!--------------------------------------------------------------------- 85 !! default option : small closed basin 86 !!--------------------------------------------------------------------- 87 CHARACTER(len=16), PUBLIC, PARAMETER :: cp_cfg = "default" !: name of the configuration 88 INTEGER , PUBLIC, PARAMETER :: jp_cfg = 0 !: resolution of the configuration 27 !!---------------------------------------------------------------------- 28 !! namcfg namelist parameters 29 !!---------------------------------------------------------------------- 30 CHARACTER(lc) :: cp_cfg !: name of the configuration 31 CHARACTER(lc) :: cp_cfz !: name of the zoom of configuration 32 INTEGER :: jp_cfg !: resolution of the configuration 89 33 90 34 ! data size !!! * size of all input files * 91 INTEGER , PUBLIC, PARAMETER :: jpidta = 10!: 1st lateral dimension ( >= jpi )92 INTEGER , PUBLIC, PARAMETER :: jpjdta = 12!: 2nd " " ( >= jpj )93 INTEGER , PUBLIC, PARAMETER :: jpkdta = 31!: number of levels ( >= jpk )35 INTEGER :: jpidta !: 1st lateral dimension ( >= jpi ) 36 INTEGER :: jpjdta !: 2nd " " ( >= jpj ) 37 INTEGER :: jpkdta !: number of levels ( >= jpk ) 94 38 95 39 ! global or zoom domain size !!! * computational domain * 96 INTEGER , PUBLIC, PARAMETER :: jpiglo = jpidta!: 1st dimension of global domain --> i97 INTEGER , PUBLIC, PARAMETER :: jpjglo = jpjdta!: 2nd - - --> j98 INTEGER, PUBLIC :: jpk = jpkdta !: number of vertical levels 40 INTEGER :: jpiglo !: 1st dimension of global domain --> i 41 INTEGER :: jpjglo !: 2nd - - --> j 42 99 43 ! zoom starting position 100 INTEGER , PUBLIC, PARAMETER :: jpizoom = 1!: left bottom (i,j) indices of the zoom101 INTEGER , PUBLIC, PARAMETER :: jpjzoom = 1!: in data domain indices44 INTEGER :: jpizoom !: left bottom (i,j) indices of the zoom 45 INTEGER :: jpjzoom !: in data domain indices 102 46 103 47 ! Domain characteristics 104 INTEGER, PUBLIC, PARAMETER :: jperio = 0 !: lateral cond. type (between 0 and 6) 105 ! ! = 0 closed ; = 1 cyclic East-West 106 ! ! = 2 equatorial symmetric ; = 3 North fold T-point pivot 107 ! ! = 4 cyclic East-West AND North fold T-point pivot 108 ! ! = 5 North fold F-point pivot 109 ! ! = 6 cyclic East-West AND North fold F-point pivot 48 INTEGER :: jperio !: lateral cond. type (between 0 and 6) 49 ! ! = 0 closed ; = 1 cyclic East-West 50 ! ! = 2 equatorial symmetric ; = 3 North fold T-point pivot 51 ! ! = 4 cyclic East-West AND North fold T-point pivot 52 ! ! = 5 North fold F-point pivot 53 ! ! = 6 cyclic East-West AND North fold F-point pivot 54 55 ! Input file read offset 56 LOGICAL :: ln_use_jattr !: Use file global attribute: open_ocean_jstart to determine start j-row 57 ! when reading input from those netcdf files that have the 58 ! attribute defined. This is designed to enable input files associated 59 ! with the extended grids used in the under ice shelf configurations to 60 ! be used without redundant rows when the ice shelves are not in use. 110 61 111 62 !! Values set to pp_not_used indicates that this parameter is not used in THIS config. 112 63 !! Values set to pp_to_be_computed indicates that variables will be computed in domzgr 113 REAL(wp) , PUBLIC, PARAMETER:: pp_not_used = 999999._wp !: vertical grid parameter114 REAL(wp) , PUBLIC, PARAMETER:: pp_to_be_computed = 999999._wp !: - - -64 REAL(wp) :: pp_not_used = 999999._wp !: vertical grid parameter 65 REAL(wp) :: pp_to_be_computed = 999999._wp !: - - - 115 66 116 67 117 !! Horizontal grid parameters for domhgr118 !! =====================================119 INTEGER, PUBLIC, PARAMETER :: jphgr_msh = 0 !: type of horizontal mesh120 ! ! = 0 curvilinear coordinate on the sphere read in coordinate.nc121 ! ! = 1 geographical mesh on the sphere with regular grid-spacing122 ! ! = 2 f-plane with regular grid-spacing123 ! ! = 3 beta-plane with regular grid-spacing124 ! ! = 4 Mercator grid with T/U point at the equator125 126 REAL(wp) , PUBLIC, PARAMETER :: ppglam0 = 0.0_wp !: longitude of first raw and column T-point (jphgr_msh = 1)127 REAL(wp) , PUBLIC, PARAMETER :: ppgphi0 = -35.0_wp !: latitude of first raw and column T-point (jphgr_msh = 1)128 ! ! used for Coriolis & Beta parameters (jphgr_msh = 2 or 3)129 REAL(wp) , PUBLIC, PARAMETER :: ppe1_deg = 1.0_wp !: zonal grid-spacing (degrees)130 REAL(wp) , PUBLIC, PARAMETER :: ppe2_deg = 0.5_wp !: meridional grid-spacing (degrees)131 REAL(wp) , PUBLIC, PARAMETER :: ppe1_m = 5000.0_wp !: zonal grid-spacing (degrees)132 REAL(wp) , PUBLIC, PARAMETER :: ppe2_m = 5000.0_wp !: meridional grid-spacing (degrees)133 134 !! Vertical grid parameter for domzgr135 !! ==================================136 REAL(wp), PUBLIC, PARAMETER :: ppsur = -4762.96143546300_wp !: ORCA r4, r2 and r05 coefficients137 REAL(wp), PUBLIC, PARAMETER :: ppa0 = 255.58049070440_wp !: (default coefficients)138 REAL(wp), PUBLIC, PARAMETER :: ppa1 = 245.58132232490_wp !:139 REAL(wp), PUBLIC, PARAMETER :: ppkth = 21.43336197938_wp !:140 REAL(wp), PUBLIC, PARAMETER :: ppacr = 3.00000000000_wp !:141 !142 ! If both ppa0 ppa1 and ppsur are specified to 0, then143 ! they are computed from ppdzmin, pphmax , ppkth, ppacr in dom_zgr144 REAL(wp), PUBLIC, PARAMETER :: ppdzmin = 10._wp !: Minimum vertical spacing145 REAL(wp), PUBLIC, PARAMETER :: pphmax = 5000._wp !: Maximum depth146 !147 LOGICAL , PUBLIC, PARAMETER :: ldbletanh = .TRUE. !: Use/do not use double tanf function for vertical coordinates148 REAL(wp), PUBLIC, PARAMETER :: ppa2 = 100.760928500000_wp !: Double tanh function parameters149 REAL(wp), PUBLIC, PARAMETER :: ppkth2= 48.029893720000_wp !:150 REAL(wp), PUBLIC, PARAMETER :: ppacr2= 13.000000000000_wp !:151 !152 !#endif153 68 154 69 … … 171 86 INTEGER, PUBLIC :: jpi ! = ( jpiglo-2*jpreci + (jpni-1) ) / jpni + 2*jpreci !: first dimension 172 87 INTEGER, PUBLIC :: jpj ! = ( jpjglo-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj !: second dimension 173 INTEGER, PUBLIC :: jpk ! = jpkdta !: third dimension88 INTEGER, PUBLIC :: jpk ! = jpkdta 174 89 INTEGER, PUBLIC :: jpim1 ! = jpi-1 !: inner domain indices 175 90 INTEGER, PUBLIC :: jpjm1 ! = jpj-1 !: - - - … … 177 92 INTEGER, PUBLIC :: jpij ! = jpi*jpj !: jpi x jpj 178 93 179 !!---------------------------------------------------------------------180 !! Optimization/control flags181 !!---------------------------------------------------------------------182 #if defined key_esopa183 LOGICAL, PUBLIC, PARAMETER :: lk_esopa = .TRUE. !: flag to activate the all options184 #else185 LOGICAL, PUBLIC, PARAMETER :: lk_esopa = .FALSE. !: flag to activate the all options186 #endif187 188 #if defined key_vectopt_loop189 LOGICAL, PUBLIC, PARAMETER :: lk_vopt_loop = .TRUE. !: vector optimization flag190 #else191 LOGICAL, PUBLIC, PARAMETER :: lk_vopt_loop = .FALSE. !: vector optimization flag192 #endif193 194 94 !!---------------------------------------------------------------------- 195 95 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 196 !! $Id: par_oce.F90 2715 2011-03-30 15:58:35Z rblod$96 !! $Id: par_oce.F90 5974 2015-12-02 10:52:05Z timgraham $ 197 97 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 198 98 !!======================================================================
Note: See TracChangeset
for help on using the changeset viewer.