- Timestamp:
- 2019-01-08T11:15:29+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/dev_10448_WAD_SBC_BUGFIX/src/OCE/SBC/sbcmod.F90
r10456 r10466 466 466 !!$!clem: it looks like it is necessary for the north fold (in certain circumstances). Don't know why. 467 467 !!$ CALL lbc_lnk( 'sbcmod', emp, 'T', 1. ) 468 IF ( ll_wd ) THEN ! If near WAD point limit the flux for now469 zthscl = atanh(rn_wd_sbcfra) ! taper frac default is .999470 zwdht(:,:) = sshn(:,:) + ht_0(:,:) - rn_wdmin1 ! do this calc of water468 IF ( ll_wd ) THEN ! If near WAD point limit the flux for now 469 zthscl = atanh(rn_wd_sbcfra) ! taper frac default is .999 470 zwdht(:,:) = sshn(:,:) + ht_0(:,:) - rn_wdmin1 ! do this calc of water 471 471 ! depth above wd limit once 472 WHERE( zwdht(:,:) <= 0.0 )472 WHERE( zwdht(:,:) <= 0.0 ) 473 473 taum(:,:) = 0.0 474 474 utau(:,:) = 0.0 … … 478 478 emp (:,:) = min(emp(:,:),0.0) !can allow puddles to grow but not shrink 479 479 sfx (:,:) = 0.0 480 END WHERE481 zwght(:,:) = tanh(zthscl*zwdht(:,:))482 WHERE( zwdht(:,:) > 0.0 .and. zwdht(:,:) < rn_wd_sbcdep ) ! 5 m hard limit here is arbitrary480 END WHERE 481 zwght(:,:) = tanh(zthscl*zwdht(:,:)) 482 WHERE( zwdht(:,:) > 0.0 .and. zwdht(:,:) < rn_wd_sbcdep ) ! 5 m hard limit here is arbitrary 483 483 qsr (:,:) = qsr(:,:) * zwght(:,:) 484 484 qns (:,:) = qns(:,:) * zwght(:,:) … … 486 486 utau (:,:) = utau(:,:) * zwght(:,:) 487 487 vtau (:,:) = vtau(:,:) * zwght(:,:) 488 sfx (:,:)= sfx(:,:) * zwght(:,:)488 sfx (:,:) = sfx(:,:) * zwght(:,:) 489 489 emp (:,:) = emp(:,:) * zwght(:,:) 490 END WHERE491 ENDIF490 END WHERE 491 ENDIF 492 492 ! 493 493 IF( kt == nit000 ) THEN ! set the forcing field at nit000 - 1 !
Note: See TracChangeset
for help on using the changeset viewer.