- Timestamp:
- 2014-12-10T11:58:55+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4765_CNRS_agrif/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r4687 r4982 216 216 & ppsur == pp_to_be_computed ) THEN 217 217 ! 218 #if defined key_agrif 219 za1 = ( ppdzmin - pphmax / FLOAT(jpkdta-1) ) & 220 & / ( TANH((1-ppkth)/ppacr) - ppacr/FLOAT(jpkdta-1) * ( LOG( COSH( (jpkdta - ppkth) / ppacr) )& 221 & - LOG( COSH( ( 1 - ppkth) / ppacr) ) ) ) 222 #else 218 223 za1 = ( ppdzmin - pphmax / FLOAT(jpkm1) ) & 219 224 & / ( TANH((1-ppkth)/ppacr) - ppacr/FLOAT(jpk-1) * ( LOG( COSH( (jpk - ppkth) / ppacr) ) & 220 225 & - LOG( COSH( ( 1 - ppkth) / ppacr) ) ) ) 226 #endif 221 227 za0 = ppdzmin - za1 * TANH( (1-ppkth) / ppacr ) 222 228 zsur = - za0 - za1 * ppacr * LOG( COSH( (1-ppkth) / ppacr ) ) … … 233 239 WRITE(numout,*) ' Uniform grid with ',jpk-1,' layers' 234 240 WRITE(numout,*) ' Total depth :', zhmax 241 #if defined key_agrif 242 WRITE(numout,*) ' Layer thickness:', zhmax/(jpkdta-1) 243 #else 235 244 WRITE(numout,*) ' Layer thickness:', zhmax/(jpk-1) 245 #endif 236 246 ELSE 237 247 IF( ppa1 == 0._wp .AND. ppa0 == 0._wp .AND. ppsur == 0._wp ) THEN … … 257 267 ! Reference z-coordinate (depth - scale factor at T- and W-points) 258 268 ! ====================== 259 IF( ppkth == 0._wp ) THEN ! uniform vertical grid 269 IF( ppkth == 0._wp ) THEN ! uniform vertical grid 270 #if defined key_agrif 271 za1 = zhmax / FLOAT(jpkdta-1) 272 #else 260 273 za1 = zhmax / FLOAT(jpk-1) 274 #endif 261 275 DO jk = 1, jpk 262 276 zw = FLOAT( jk )
Note: See TracChangeset
for help on using the changeset viewer.