- 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/OPA_SRC/DOM/domain.F90
r3764 r4147 126 126 NAMELIST/namdom/ nn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh , rn_hmin, & 127 127 & nn_acc , rn_atfp , rn_rdt , rn_rdtmin , & 128 & rn_rdtmax, rn_rdth , nn_baro , nn_closea 128 & rn_rdtmax, rn_rdth , nn_baro , nn_closea , & 129 & jphgr_msh, & 130 & ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & 131 & ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 132 & ppa2, ppkth2, ppacr2 129 133 NAMELIST/namcla/ nn_cla 130 134 #if defined key_netcdf4 131 135 NAMELIST/namnc4/ nn_nchunks_i, nn_nchunks_j, nn_nchunks_k, ln_nc4zip 132 136 #endif 133 !!---------------------------------------------------------------------- 134 135 REWIND( numnam ) ! Namelist namrun : parameters of the run 136 READ ( numnam, namrun ) 137 INTEGER :: ios ! Local integer output status for namelist read 138 !!---------------------------------------------------------------------- 139 140 REWIND( numnam_ref ) ! Namelist namrun in reference namelist : Parameters of the run 141 READ ( numnam_ref, namrun, IOSTAT = ios, ERR = 901) 142 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namrun in reference namelist', lwp ) 143 144 REWIND( numnam_cfg ) ! Namelist namrun in configuration namelist : Parameters of the run 145 READ ( numnam_cfg, namrun, IOSTAT = ios, ERR = 902 ) 146 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namrun in configuration namelist', lwp ) 147 WRITE ( numond, namrun ) 137 148 ! 138 149 IF(lwp) THEN ! control print … … 143 154 WRITE(numout,*) ' job number nn_no = ', nn_no 144 155 WRITE(numout,*) ' experiment name for output cn_exp = ', cn_exp 156 WRITE(numout,*) ' file prefix restart input cn_ocerst_in= ', cn_ocerst_in 157 WRITE(numout,*) ' file prefix restart output cn_ocerst_out= ', cn_ocerst_out 145 158 WRITE(numout,*) ' restart logical ln_rstart = ', ln_rstart 146 159 WRITE(numout,*) ' control of time step nn_rstctl = ', nn_rstctl … … 200 213 #endif 201 214 202 REWIND( numnam ) ! Namelist namdom : space & time domain (bathymetry, mesh, timestep) 203 READ ( numnam, namdom ) 215 REWIND( numnam_ref ) ! Namelist namdom in reference namelist : space & time domain (bathymetry, mesh, timestep) 216 READ ( numnam_ref, namdom, IOSTAT = ios, ERR = 903) 217 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in reference namelist', lwp ) 218 219 REWIND( numnam_cfg ) ! Namelist namdom in configuration namelist : space & time domain (bathymetry, mesh, timestep) 220 READ ( numnam_cfg, namdom, IOSTAT = ios, ERR = 904 ) 221 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in configuration namelist', lwp ) 222 WRITE ( numond, namdom ) 204 223 205 224 IF(lwp) THEN … … 224 243 WRITE(numout,*) ' depth of transition rn_rdth = ', rn_rdth 225 244 WRITE(numout,*) ' suppression of closed seas (=0) nn_closea = ', nn_closea 245 WRITE(numout,*) ' type of horizontal mesh jphgr_msh = ', jphgr_msh 246 WRITE(numout,*) ' longitude of first raw and column T-point ppglam0 = ', ppglam0 247 WRITE(numout,*) ' latitude of first raw and column T-point ppgphi0 = ', ppgphi0 248 WRITE(numout,*) ' zonal grid-spacing (degrees) ppe1_deg = ', ppe1_deg 249 WRITE(numout,*) ' meridional grid-spacing (degrees) ppe2_deg = ', ppe2_deg 250 WRITE(numout,*) ' zonal grid-spacing (degrees) ppe1_m = ', ppe1_m 251 WRITE(numout,*) ' meridional grid-spacing (degrees) ppe2_m = ', ppe2_m 252 WRITE(numout,*) ' ORCA r4, r2 and r05 coefficients ppsur = ', ppsur 253 WRITE(numout,*) ' ppa0 = ', ppa0 254 WRITE(numout,*) ' ppa1 = ', ppa1 255 WRITE(numout,*) ' ppkth = ', ppkth 256 WRITE(numout,*) ' ppacr = ', ppacr 257 WRITE(numout,*) ' Minimum vertical spacing ppdzmin = ', ppdzmin 258 WRITE(numout,*) ' Maximum depth pphmax = ', pphmax 259 WRITE(numout,*) ' Use double tanf function for vertical coordinates ldbletanh = ', ldbletanh 260 WRITE(numout,*) ' Double tanh function parameters ppa2 = ', ppa2 261 WRITE(numout,*) ' ppkth2 = ', ppkth2 262 WRITE(numout,*) ' ppacr2 = ', ppacr2 226 263 ENDIF 227 264 … … 237 274 rdth = rn_rdth 238 275 239 REWIND( numnam ) ! Namelist cross land advection 240 READ ( numnam, namcla ) 276 REWIND( numnam_ref ) ! Namelist namcla in reference namelist : Cross land advection 277 READ ( numnam_ref, namcla, IOSTAT = ios, ERR = 905) 278 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in reference namelist', lwp ) 279 280 REWIND( numnam_cfg ) ! Namelist namcla in configuration namelist : Cross land advection 281 READ ( numnam_cfg, namcla, IOSTAT = ios, ERR = 906 ) 282 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in configuration namelist', lwp ) 283 WRITE( numond, namcla ) 284 241 285 IF(lwp) THEN 242 286 WRITE(numout,*) … … 244 288 WRITE(numout,*) ' cross land advection nn_cla = ', nn_cla 245 289 ENDIF 290 IF ( nn_cla .EQ. 1 ) THEN 291 IF ( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA R2 292 CONTINUE 293 ELSE 294 CALL ctl_stop( 'STOP', 'Cross land advation iplemented only for ORCA2 configuration: cp_cfg = "orca" and jp_cfg = 2 ' ) 295 ENDIF 296 ENDIF 246 297 247 298 #if defined key_netcdf4 248 299 ! ! NetCDF 4 case ("key_netcdf4" defined) 249 REWIND( numnam ) ! Namelist namnc4 : netcdf4 chunking parameters 250 READ ( numnam, namnc4 ) 300 REWIND( numnam_ref ) ! Namelist namnc4 in reference namelist : NETCDF 301 READ ( numnam_ref, namnc4, IOSTAT = ios, ERR = 907) 302 907 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namnc4 in reference namelist', lwp ) 303 304 REWIND( numnam_cfg ) ! Namelist namnc4 in configuration namelist : NETCDF 305 READ ( numnam_cfg, namnc4, IOSTAT = ios, ERR = 908 ) 306 908 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namnc4 in configuration namelist', lwp ) 307 WRITE( numond, namnc4 ) 308 251 309 IF(lwp) THEN ! control print 252 310 WRITE(numout,*)
Note: See TracChangeset
for help on using the changeset viewer.