- Timestamp:
- 2013-11-04T12:51:55+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OFF_SRC/domain.F90
r3632 r4147 88 88 !!---------------------------------------------------------------------- 89 89 USE ioipsl 90 INTEGER :: ios ! Local integer output status for namelist read 90 91 NAMELIST/namrun/ nn_no , cn_exp , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl, & 91 92 & nn_it000, nn_itend , nn_date0 , nn_leapy , nn_istate , nn_stock , & 92 93 & nn_write, ln_dimgnnn, ln_mskland , ln_clobber , nn_chunksz 93 NAMELIST/namdom/ nn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh , rn_hmin, & 94 & nn_acc , rn_atfp , rn_rdt , rn_rdtmin , & 95 & rn_rdtmax, rn_rdth , nn_baro , nn_closea 94 NAMELIST/namdom/ nn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh , rn_hmin, & 95 & nn_acc , rn_atfp , rn_rdt , rn_rdtmin , & 96 & rn_rdtmax, rn_rdth , nn_baro , nn_closea , & 97 & jphgr_msh, & 98 & ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & 99 & ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 100 & ppa2, ppkth2, ppacr2 96 101 NAMELIST/namcla/ nn_cla 97 102 #if defined key_netcdf4 … … 100 105 !!---------------------------------------------------------------------- 101 106 102 REWIND( numnam ) ! Namelist namrun : parameters of the run 103 READ ( numnam, namrun ) 107 REWIND( numnam_ref ) ! Namelist namrun in reference namelist : Parameters of the run 108 READ ( numnam_ref, namrun, IOSTAT = ios, ERR = 901) 109 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namrun in reference namelist', lwp ) 110 111 REWIND( numnam_cfg ) ! Namelist namrun in configuration namelist : Parameters of the run 112 READ ( numnam_cfg, namrun, IOSTAT = ios, ERR = 902 ) 113 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namrun in configuration namelist', lwp ) 114 WRITE ( numond, namrun ) 104 115 ! 105 116 IF(lwp) THEN ! control print … … 170 181 #endif 171 182 172 REWIND( numnam ) ! Domain 173 READ ( numnam, namdom ) 183 REWIND( numnam_ref ) ! Namelist namdom in reference namelist : space & time domain (bathymetry, mesh, timestep) 184 READ ( numnam_ref, namdom, IOSTAT = ios, ERR = 903) 185 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in reference namelist', lwp ) 186 187 REWIND( numnam_cfg ) ! Namelist namdom in configuration namelist : space & time domain (bathymetry, mesh, timestep) 188 READ ( numnam_cfg, namdom, IOSTAT = ios, ERR = 904 ) 189 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in configuration namelist', lwp ) 190 WRITE ( numond, namdom ) 174 191 175 192 IF(lwp) THEN … … 193 210 WRITE(numout,*) ' depth of transition rn_rdth = ', rn_rdth 194 211 WRITE(numout,*) ' suppression of closed seas (=0) nn_closea = ', nn_closea 212 WRITE(numout,*) ' type of horizontal mesh jphgr_msh = ', jphgr_msh 213 WRITE(numout,*) ' longitude of first raw and column T-point ppglam0 = ', ppglam0 214 WRITE(numout,*) ' latitude of first raw and column T-point ppgphi0 = ', ppgphi0 215 WRITE(numout,*) ' zonal grid-spacing (degrees) ppe1_deg = ', ppe1_deg 216 WRITE(numout,*) ' meridional grid-spacing (degrees) ppe2_deg = ', ppe2_deg 217 WRITE(numout,*) ' zonal grid-spacing (degrees) ppe1_m = ', ppe1_m 218 WRITE(numout,*) ' meridional grid-spacing (degrees) ppe2_m = ', ppe2_m 219 WRITE(numout,*) ' ORCA r4, r2 and r05 coefficients ppsur = ', ppsur 220 WRITE(numout,*) ' ppa0 = ', ppa0 221 WRITE(numout,*) ' ppa1 = ', ppa1 222 WRITE(numout,*) ' ppkth = ', ppkth 223 WRITE(numout,*) ' ppacr = ', ppacr 224 WRITE(numout,*) ' Minimum vertical spacing ppdzmin = ', ppdzmin 225 WRITE(numout,*) ' Maximum depth pphmax = ', pphmax 226 WRITE(numout,*) ' Use double tanf function for vertical coordinates ldbletanh = ', ldbletanh 227 WRITE(numout,*) ' Double tanh function parameters ppa2 = ', ppa2 228 WRITE(numout,*) ' ppkth2 = ', ppkth2 229 WRITE(numout,*) ' ppacr2 = ', ppacr2 195 230 ENDIF 196 231 … … 206 241 rdth = rn_rdth 207 242 208 REWIND( numnam ) ! Namelist cross land advection 209 READ ( numnam, namcla ) 243 REWIND( numnam_ref ) ! Namelist namcla in reference namelist : Cross land advection 244 READ ( numnam_ref, namcla, IOSTAT = ios, ERR = 905) 245 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in reference namelist', lwp ) 246 247 REWIND( numnam_cfg ) ! Namelist namcla in configuration namelist : Cross land advection 248 READ ( numnam_cfg, namcla, IOSTAT = ios, ERR = 906 ) 249 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in configuration namelist', lwp ) 250 WRITE( numond, namcla ) 251 210 252 IF(lwp) THEN 211 253 WRITE(numout,*) … … 216 258 #if defined key_netcdf4 217 259 ! ! NetCDF 4 case ("key_netcdf4" defined) 218 REWIND( numnam ) ! Namelist namnc4 : netcdf4 chunking parameters 219 READ ( numnam, namnc4 ) 260 REWIND( numnam_ref ) ! Namelist namnc4 in reference namelist : NETCDF 261 READ ( numnam_ref, namnc4, IOSTAT = ios, ERR = 907) 262 907 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namnc4 in reference namelist', lwp ) 263 264 REWIND( numnam_cfg ) ! Namelist namnc4 in configuration namelist : NETCDF 265 READ ( numnam_cfg, namnc4, IOSTAT = ios, ERR = 908 ) 266 908 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namnc4 in configuration namelist', lwp ) 267 WRITE( numond, namnc4 ) 220 268 IF(lwp) THEN ! control print 221 269 WRITE(numout,*) … … 257 305 !!---------------------------------------------------------------------- 258 306 INTEGER :: ioptio = 0 ! temporary integer 307 INTEGER :: ios 259 308 !! 260 309 NAMELIST/namzgr/ ln_zco, ln_zps, ln_sco 261 310 !!---------------------------------------------------------------------- 262 311 263 REWIND ( numnam ) ! Read Namelist namzgr : vertical coordinate' 264 READ ( numnam, namzgr ) 312 REWIND( numnam_ref ) ! Namelist namzgr in reference namelist : Vertical coordinate 313 READ ( numnam_ref, namzgr, IOSTAT = ios, ERR = 901 ) 314 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in reference namelist', lwp ) 315 316 REWIND( numnam_cfg ) ! Namelist namzgr in configuration namelist : Vertical coordinate 317 READ ( numnam_cfg, namzgr, IOSTAT = ios, ERR = 902 ) 318 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in configuration namelist', lwp ) 319 WRITE ( numond, namzgr ) 265 320 266 321 IF(lwp) THEN ! Control print
Note: See TracChangeset
for help on using the changeset viewer.