Changeset 7030
- Timestamp:
- 2016-10-14T14:24:48+02:00 (8 years ago)
- Location:
- branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r7022 r7030 529 529 zdta(:,jpjdta-2:jpjdta) = -4._wp 530 530 ! ! =========================== 531 CASE ( 6 ) ! WAD 6 configuration 532 ! ! ==================== 533 ! 534 IF(lwp) WRITE(numout,*) 535 IF(lwp) WRITE(numout,*) 'zgr_bat : Parabolic channel with gaussian ridge' 536 IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 537 ! 538 DO ji = 1, jpidta 539 zi = MAX(1.0-FLOAT((ji-25)**2)/484.0, -2.0 ) 540 zj = 0.95*MAX(EXP(-1.0*FLOAT((ji-25)**2)/32.0) , 0.0 ) 541 zdta(ji,:) = MAX(rn_bathy*(zi-zj), -2.0) 542 IF(lwp)write(numout,*) 'ZDTA ',ji,zi,zdta(ji,1) 543 END DO 544 zdta(1:2,:) = -4._wp 545 zdta(jpidta-1:jpidta,:) = -4._wp 546 zdta(:,1) = -4._wp 547 zdta(:,jpjdta) = -4._wp 548 zdta(:,1:3) = -4._wp 549 zdta(:,jpjdta-2:jpjdta) = -4._wp 550 ! ! =========================== 531 551 CASE DEFAULT 532 552 ! ! =========================== -
branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90
r6152 r7030 461 461 IF(ll_tmp1.AND.ll_tmp2) THEN 462 462 zcpx(ji,jj) = 1.0_wp 463 wduflt(ji,jj) = 1.0_wp464 463 ELSE IF(ll_tmp3) THEN 465 464 ! no worries about sshn(ji+1,jj)-sshn(ji,jj) = 0, it won't happen ! here 466 465 zcpx(ji,jj) = ABS((sshn(ji+1,jj) + bathy(ji+1,jj) - sshn(ji,jj) - bathy(ji,jj)) / & 467 466 & (sshn(ji+1,jj) - sshn(ji,jj))) 468 wduflt(ji,jj) = 1.0_wp469 467 ELSE 470 468 zcpx(ji,jj) = 0._wp 471 wduflt(ji,jj) = 0.0_wp472 469 END IF 473 470 … … 479 476 IF(ll_tmp1.AND.ll_tmp2) THEN 480 477 zcpy(ji,jj) = 1.0_wp 481 wdvflt(ji,jj) = 1.0_wp482 478 ELSE IF(ll_tmp3) THEN 483 479 ! no worries about sshn(ji,jj+1)-sshn(ji,jj) = 0, it won't happen ! here 484 480 zcpy(ji,jj) = ABS((sshn(ji,jj+1) + bathy(ji,jj+1) - sshn(ji,jj) - bathy(ji,jj)) / & 485 481 & (sshn(ji,jj+1) - sshn(ji,jj))) 486 wdvflt(ji,jj) = 1.0_wp487 482 ELSE 488 483 zcpy(ji,jj) = 0._wp 489 wdvflt(ji,jj) = 0.0_wp490 484 END IF 491 485 END DO -
branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DYN/wet_dry.F90
r7022 r7030 33 33 !! --------------------------------------------------------------------- 34 34 35 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: wduflt, wdvflt !: u- and v- filter36 35 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: wdmask !: u- and v- limiter 37 36 … … 88 87 89 88 IF(ln_wd) THEN 90 ALLOCATE( wd uflt(jpi,jpj), wdvflt(jpi,jpj), wdmask(jpi,jpj), STAT=ierr )89 ALLOCATE( wdmask(jpi,jpj), STAT=ierr ) 91 90 IF( ierr /= 0 ) CALL ctl_stop('STOP', 'wad_init : Array allocation error') 92 91 ENDIF … … 474 473 475 474 ! 475 ! ! =========================== 476 CASE ( 6 ) ! WAD 6 configuration 477 ! ! ==================== 478 ! 479 IF(lwp) WRITE(numout,*) 480 IF(lwp) WRITE(numout,*) 'istate_wad : Parobolic EW channel with gaussian ridge' 481 IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 482 ! 483 do ji = 1,jpi 484 !6a 485 sshn(ji,:) = ( -5.5_wp + 9.0_wp*FLOAT(jpidta - mig(ji))/FLOAT(jpidta-1))*tmask(ji,:,1) 486 !Some variations in initial slope that have been tested 487 !6b 488 !sshn(ji,:) = ( -5.5_wp + 6.5_wp*FLOAT(jpidta - mig(ji))/FLOAT(jpidta-1))*tmask(ji,:,1) 489 !6c 490 !sshn(ji,:) = ( -5.5_wp + 7.5_wp*FLOAT(jpidta - mig(ji))/FLOAT(jpidta-1))*tmask(ji,:,1) 491 !6d 492 !sshn(ji,:) = ( -4.5_wp + 8.0_wp*FLOAT(jpidta - mig(ji))/FLOAT(jpidta-1))*tmask(ji,:,1) 493 end do 494 495 ! 476 496 ! ! =========================== 477 497 CASE DEFAULT ! NONE existing configuration … … 487 507 do jj = 1,jpj 488 508 do ji = 1,jpi 489 IF( bathy(ji,jj) + sshn(ji,jj) < 0.4_wp) THEN490 sshn(ji,jj) = tmask(ji,jj,1)*( 0.4_wp - bathy(ji,jj))509 IF( bathy(ji,jj) + sshn(ji,jj) < rn_wdmin1 ) THEN 510 sshn(ji,jj) = tmask(ji,jj,1)*( rn_wdmin1 - bathy(ji,jj) ) 491 511 ENDIF 492 512 end do
Note: See TracChangeset
for help on using the changeset viewer.