- Timestamp:
- 2015-11-30T20:55:41+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r5944 r5956 219 219 & ppsur == pp_to_be_computed ) THEN 220 220 ! 221 #if defined key_agrif 222 za1 = ( ppdzmin - pphmax / FLOAT(jpkdta-1) ) & 223 & / ( TANH((1-ppkth)/ppacr) - ppacr/FLOAT(jpkdta-1) * ( LOG( COSH( (jpkdta - ppkth) / ppacr) )& 224 & - LOG( COSH( ( 1 - ppkth) / ppacr) ) ) ) 225 #else 221 226 za1 = ( ppdzmin - pphmax / FLOAT(jpkm1) ) & 222 227 & / ( TANH((1-ppkth)/ppacr) - ppacr/FLOAT(jpk-1) * ( LOG( COSH( (jpk - ppkth) / ppacr) ) & 223 228 & - LOG( COSH( ( 1 - ppkth) / ppacr) ) ) ) 229 #endif 224 230 za0 = ppdzmin - za1 * TANH( (1-ppkth) / ppacr ) 225 231 zsur = - za0 - za1 * ppacr * LOG( COSH( (1-ppkth) / ppacr ) ) … … 236 242 WRITE(numout,*) ' Uniform grid with ',jpk-1,' layers' 237 243 WRITE(numout,*) ' Total depth :', zhmax 244 #if defined key_agrif 245 WRITE(numout,*) ' Layer thickness:', zhmax/(jpkdta-1) 246 #else 238 247 WRITE(numout,*) ' Layer thickness:', zhmax/(jpk-1) 248 #endif 239 249 ELSE 240 250 IF( ppa1 == 0._wp .AND. ppa0 == 0._wp .AND. ppsur == 0._wp ) THEN … … 260 270 ! Reference z-coordinate (depth - scale factor at T- and W-points) 261 271 ! ====================== 262 IF( ppkth == 0._wp ) THEN ! uniform vertical grid 272 IF( ppkth == 0._wp ) THEN ! uniform vertical grid 273 #if defined key_agrif 274 za1 = zhmax / FLOAT(jpkdta-1) 275 #else 263 276 za1 = zhmax / FLOAT(jpk-1) 277 #endif 264 278 DO jk = 1, jpk 265 279 zw = FLOAT( jk ) … … 487 501 IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA R2 configuration 488 502 ! ! ===================== 489 IF( nn_cla == 0 ) THEN 490 ii0 = 140 ; ii1 = 140 ! Gibraltar Strait open 491 ij0 = 102 ; ij1 = 102 ! (Thomson, Ocean Modelling, 1995) 492 DO ji = mi0(ii0), mi1(ii1) 493 DO jj = mj0(ij0), mj1(ij1) 494 mbathy(ji,jj) = 15 495 END DO 503 ! 504 ii0 = 140 ; ii1 = 140 ! Gibraltar Strait open 505 ij0 = 102 ; ij1 = 102 ! (Thomson, Ocean Modelling, 1995) 506 DO ji = mi0(ii0), mi1(ii1) 507 DO jj = mj0(ij0), mj1(ij1) 508 mbathy(ji,jj) = 15 496 509 END DO 497 IF(lwp) WRITE(numout,*)498 IF(lwp) WRITE(numout,*) ' orca_r2: Gibraltar strait open at i=',ii0,' j=',ij0499 !500 ii0 = 160 ; ii1 = 160 ! Bab el mandeb Strait open501 ij0 = 88 ; ij1 = 88 ! (Thomson, Ocean Modelling, 1995)502 DO ji = mi0(ii0), mi1(ii1)503 DO jj = mj0(ij0), mj1(ij1)504 mbathy(ji,jj) = 12505 END DO510 END DO 511 IF(lwp) WRITE(numout,*) 512 IF(lwp) WRITE(numout,*) ' orca_r2: Gibraltar strait open at i=',ii0,' j=',ij0 513 ! 514 ii0 = 160 ; ii1 = 160 ! Bab el mandeb Strait open 515 ij0 = 88 ; ij1 = 88 ! (Thomson, Ocean Modelling, 1995) 516 DO ji = mi0(ii0), mi1(ii1) 517 DO jj = mj0(ij0), mj1(ij1) 518 mbathy(ji,jj) = 12 506 519 END DO 507 IF(lwp) WRITE(numout,*)508 IF(lwp) WRITE(numout,*) ' orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0509 ENDIF520 END DO 521 IF(lwp) WRITE(numout,*) 522 IF(lwp) WRITE(numout,*) ' orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0 510 523 ! 511 524 ENDIF … … 531 544 ! 532 545 IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA R2 configuration 533 ! 534 IF( nn_cla == 0 ) THEN 535 ii0 = 140 ; ii1 = 140 ! Gibraltar Strait open 536 ij0 = 102 ; ij1 = 102 ! (Thomson, Ocean Modelling, 1995) 537 DO ji = mi0(ii0), mi1(ii1) 538 DO jj = mj0(ij0), mj1(ij1) 539 bathy(ji,jj) = 284._wp 540 END DO 546 ! 547 ii0 = 140 ; ii1 = 140 ! Gibraltar Strait open 548 ij0 = 102 ; ij1 = 102 ! (Thomson, Ocean Modelling, 1995) 549 DO ji = mi0(ii0), mi1(ii1) 550 DO jj = mj0(ij0), mj1(ij1) 551 bathy(ji,jj) = 284._wp 541 552 END DO 542 IF(lwp) WRITE(numout,*)543 IF(lwp) WRITE(numout,*) ' orca_r2: Gibraltar strait open at i=',ii0,' j=',ij0544 !545 ii0 = 160 ; ii1 = 160 ! Bab el mandeb Strait open546 ij0 = 88 ; ij1 = 88 ! (Thomson, Ocean Modelling, 1995)547 DO ji = mi0(ii0), mi1(ii1)548 DO jj = mj0(ij0), mj1(ij1)549 bathy(ji,jj) = 137._wp550 END DO553 END DO 554 IF(lwp) WRITE(numout,*) 555 IF(lwp) WRITE(numout,*) ' orca_r2: Gibraltar strait open at i=',ii0,' j=',ij0 556 ! 557 ii0 = 160 ; ii1 = 160 ! Bab el mandeb Strait open 558 ij0 = 88 ; ij1 = 88 ! (Thomson, Ocean Modelling, 1995) 559 DO ji = mi0(ii0), mi1(ii1) 560 DO jj = mj0(ij0), mj1(ij1) 561 bathy(ji,jj) = 137._wp 551 562 END DO 552 IF(lwp) WRITE(numout,*)553 IF(lwp) WRITE(numout,*) ' orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0554 ENDIF563 END DO 564 IF(lwp) WRITE(numout,*) 565 IF(lwp) WRITE(numout,*) ' orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0 555 566 ! 556 567 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.