New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4639 for branches/2013/dev_ECMWF_waves/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90 – NEMO

Ignore:
Timestamp:
2014-05-13T18:39:35+02:00 (10 years ago)
Author:
acc
Message:

Branch 2013/dev_ECMWF_waves ticket #1216. Add limiter on zphio_flux (zdftke.F90) to prevent spurious negative values in input data causing NaNs?

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_ECMWF_waves/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

    r4351 r4639  
    271271               lambda = 2.38_wp/z0 
    272272               ! The energy flux from the wave model in physical units 
    273                zphio_flux = phioc_wavepar(ji,jj) 
     273               zphio_flux = MAX(phioc_wavepar(ji,jj), 0.0_wp) 
    274274               ! TKE at surface (Mellor and Blumberg, 2004) 
    275275               e0(ji,jj) = 0.5_wp*( 15.8_wp*zphio_flux/rau0 )**0.67_wp 
Note: See TracChangeset for help on using the changeset viewer.