Changeset 2528 for trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90
- Timestamp:
- 2010-12-27T18:33:53+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90
- Property svn:eol-style deleted
r1792 r2528 4 4 !! Ocean initialization : domain initialization 5 5 !!============================================================================== 6 !! History : ! 88-03 (G. Madec) 7 !! ! 91-11 (G. Madec) 8 !! ! 92-06 (M. Imbard) 9 !! ! 96-01 (G. Madec) terrain following coordinates 10 !! ! 97-02 (G. Madec) print mesh informations 11 !! ! 99-11 (M. Imbard) NetCDF format with IO-IPSL 12 !! ! 00-08 (D. Ludicone) Reduced section at Bab el Mandeb 13 !! ! 01-09 (M. Levy) eel config: grid in km, beta-plane 14 !! 8.5 ! 02-08 (G. Madec) F90: Free form and module, namelist 15 !! 9.0 ! 04-01 (A.M. Treguier, J.M. Molines) Case 4 (Mercator mesh) 16 !! use of parameters in par_CONFIG-Rxx.h90, not in namelist 17 !! ! 04-05 (A. Koch-Larrouy) Add Gyre configuration 6 !! History : OPA ! 1988-03 (G. Madec) Original code 7 !! 7.0 ! 1996-01 (G. Madec) terrain following coordinates 8 !! 8.0 ! 1997-02 (G. Madec) print mesh informations 9 !! 8.1 ! 1999-11 (M. Imbard) NetCDF format with IO-IPSL 10 !! 8.2 ! 2000-08 (D. Ludicone) Reduced section at Bab el Mandeb 11 !! - ! 2001-09 (M. Levy) eel config: grid in km, beta-plane 12 !! NEMO 1.0 ! 2002-08 (G. Madec) F90: Free form and module, namelist 13 !! - ! 2004-01 (A.M. Treguier, J.M. Molines) Case 4 (Mercator mesh) 14 !! use of parameters in par_CONFIG-Rxx.h90, not in namelist 15 !! - ! 2004-05 (A. Koch-Larrouy) Add Gyre configuration 18 16 !!---------------------------------------------------------------------- 19 17 20 18 !!---------------------------------------------------------------------- 21 !! dom_hgr 22 !! hgr_read 19 !! dom_hgr : initialize the horizontal mesh 20 !! hgr_read : read "coordinate" NetCDF file 23 21 !!---------------------------------------------------------------------- 24 !! * Modules used 25 USE dom_oce ! ocean space and time domain 26 USE phycst ! physical constants 27 USE in_out_manager ! I/O manager 28 USE lib_mpp 22 USE dom_oce ! ocean space and time domain 23 USE phycst ! physical constants 24 USE in_out_manager ! I/O manager 25 USE lib_mpp ! MPP library 29 26 30 27 IMPLICIT NONE 31 28 PRIVATE 32 29 33 !! * Module variables 34 REAL(wp) :: glam0, gphi0 ! variables corresponding to parameters 35 ! ! ppglam0 ppgphi0 set in par_oce 36 37 !! * Routine accessibility 38 PUBLIC dom_hgr ! called by domain.F90 30 REAL(wp) :: glam0, gphi0 ! variables corresponding to parameters ppglam0 ppgphi0 set in par_oce 31 32 PUBLIC dom_hgr ! called by domain.F90 33 39 34 !!---------------------------------------------------------------------- 40 !! OPA 9.0 , LOCEAN-IPSL (2005)35 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 41 36 !! $Id$ 42 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)37 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 43 38 !!---------------------------------------------------------------------- 44 45 39 CONTAINS 46 40 … … 100 94 !! Madec, Imbard, 1996, Clim. Dyn. 101 95 !!---------------------------------------------------------------------- 102 INTEGER :: ji, jj ! dummy loop indices 103 INTEGER :: ii0, ii1, ij0, ij1 ! temporary integers 104 INTEGER :: ijeq ! index of equator T point (used in case 4) 105 REAL(wp) :: & 106 zti, zui, zvi, zfi, & ! temporary scalars 107 ztj, zuj, zvj, zfj, & ! 108 zphi0, zbeta, znorme, & ! 109 zarg, zf0, zminff, zmaxff 110 REAL(wp) :: & 111 zlam1, zcos_alpha, zim1 , zjm1 , ze1, ze1deg, & 112 zphi1, zsin_alpha, zim05, zjm05 96 INTEGER :: ji, jj ! dummy loop indices 97 INTEGER :: ii0, ii1, ij0, ij1 ! temporary integers 98 INTEGER :: ijeq ! index of equator T point (used in case 4) 99 REAL(wp) :: zti, zui, zvi, zfi ! local scalars 100 REAL(wp) :: ztj, zuj, zvj, zfj ! - - 101 REAL(wp) :: zphi0, zbeta, znorme ! 102 REAL(wp) :: zarg, zf0, zminff, zmaxff 103 REAL(wp) :: zlam1, zcos_alpha, zim1 , zjm1 , ze1, ze1deg 104 REAL(wp) :: zphi1, zsin_alpha, zim05, zjm05 113 105 !!---------------------------------------------------------------------- 114 106 … … 138 130 IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA R2 configuration 139 131 ! ! ===================== 140 IF( n _cla == 0 ) THEN132 IF( nn_cla == 0 ) THEN 141 133 ! 142 134 ii0 = 139 ; ii1 = 140 ! Gibraltar Strait (e2u = 20 km) … … 157 149 IF(lwp) WRITE(numout,*) 158 150 IF(lwp) WRITE(numout,*) ' orca_r2: Danish Straits : e2u reduced to 10 km' 151 ! 152 ENDIF 153 154 ! ! ===================== 155 IF( cp_cfg == "orca" .AND. jp_cfg == 1 ) THEN ! ORCA R1 configuration 156 ! ! ===================== 157 158 ii0 = 281 ; ii1 = 282 ! Gibraltar Strait (e2u = 20 km) 159 ij0 = 200 ; ij1 = 200 ; e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 20.e3 160 IF(lwp) WRITE(numout,*) 161 IF(lwp) WRITE(numout,*) ' orca_r1: Gibraltar : e2u reduced to 20 km' 162 163 ii0 = 314 ; ii1 = 315 ! Bhosporus Strait (e2u = 10 km) 164 ij0 = 208 ; ij1 = 208 ; e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 10.e3 165 IF(lwp) WRITE(numout,*) 166 IF(lwp) WRITE(numout,*) ' orca_r1: Bhosporus : e2u reduced to 10 km' 167 168 ii0 = 44 ; ii1 = 44 ! Lombok Strait (e1v = 13 km) 169 ij0 = 124 ; ij1 = 125 ; e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 13.e3 170 IF(lwp) WRITE(numout,*) 171 IF(lwp) WRITE(numout,*) ' orca_r1: Lombok : e1v reduced to 10 km' 172 173 ii0 = 48 ; ii1 = 48 ! Sumba Strait (e1v = 8 km) [closed from bathy_11 on] 174 ij0 = 124 ; ij1 = 125 ; e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 8.e3 175 IF(lwp) WRITE(numout,*) 176 IF(lwp) WRITE(numout,*) ' orca_r1: Sumba : e1v reduced to 8 km' 177 178 ii0 = 53 ; ii1 = 53 ! Ombai Strait (e1v = 13 km) 179 ij0 = 124 ; ij1 = 125 ; e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 13.e3 180 IF(lwp) WRITE(numout,*) 181 IF(lwp) WRITE(numout,*) ' orca_r1: Ombai : e1v reduced to 13 km' 182 183 ii0 = 56 ; ii1 = 56 ! Timor Passage (e1v = 20 km) 184 ij0 = 124 ; ij1 = 125 ; e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 20.e3 185 IF(lwp) WRITE(numout,*) 186 IF(lwp) WRITE(numout,*) ' orca_r1: Timor Passage : e1v reduced to 20 km' 187 188 ii0 = 55 ; ii1 = 55 ! West Halmahera Strait (e1v = 30 km) 189 ij0 = 141 ; ij1 = 142 ; e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 30.e3 190 IF(lwp) WRITE(numout,*) 191 IF(lwp) WRITE(numout,*) ' orca_r1: W Halmahera : e1v reduced to 30 km' 192 193 ii0 = 58 ; ii1 = 58 ! East Halmahera Strait (e1v = 50 km) 194 ij0 = 141 ; ij1 = 142 ; e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 50.e3 195 IF(lwp) WRITE(numout,*) 196 IF(lwp) WRITE(numout,*) ' orca_r1: E Halmahera : e1v reduced to 50 km' 197 198 ! 199 200 ! 201 ! 202 ! 159 203 ! 160 204 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.