Skip to content
Snippets Groups Projects
Commit 31d4d0c5 authored by RIETTE Sébastien's avatar RIETTE Sébastien
Browse files

S. Riette 4 oct 2022: bugfix for OCEAN key in mode_tke_eps_sources

parent 7872f1ae
No related branches found
No related tags found
No related merge requests found
...@@ -266,7 +266,7 @@ END IF ...@@ -266,7 +266,7 @@ END IF
!* 2.2 Explicit TKE sources except horizontal turbulent transport !* 2.2 Explicit TKE sources except horizontal turbulent transport
! !
! extrapolate the dynamic production with a 1/Z law from its value at the ! extrapolate the dynamic production with a 1/Z law from its value at the
IF (LOCEAN) THEN IF (OOCEAN) THEN
! W(IKE) value stored in PDP(IKE) to the mass localization of tke(IKE) ! W(IKE) value stored in PDP(IKE) to the mass localization of tke(IKE)
PDP(:,:,IKE) = PDP(:,:,IKE) * (1. + PDZZ(:,:,IKE)/PDZZ(:,:,IKE+1)) PDP(:,:,IKE) = PDP(:,:,IKE) * (1. + PDZZ(:,:,IKE)/PDZZ(:,:,IKE+1))
ELSE ELSE
...@@ -286,7 +286,7 @@ ZSOURCE(:,:,:) = ( PRTKES(:,:,:) + PRTKEMS(:,:,:) ) / PRHODJ(:,:,:) & ...@@ -286,7 +286,7 @@ ZSOURCE(:,:,:) = ( PRTKES(:,:,:) + PRTKEMS(:,:,:) ) / PRHODJ(:,:,:) &
! !
! To add here in ZSOURCE surface flux of TKE ! To add here in ZSOURCE surface flux of TKE
!(assumed to be 0 for ATM, !(assumed to be 0 for ATM,
IF (LOCEAN) THEN IF (OOCEAN) THEN
!for ocean:wave breaking simple/very rough param wE = 100 Ustar**3 where ustar is the Tau_atmi/rhocea !for ocean:wave breaking simple/very rough param wE = 100 Ustar**3 where ustar is the Tau_atmi/rhocea
ZSOURCE (:,:,IKE)=ZSOURCE(:,:,IKE)-1.E2*((PSFUM(:,:)**2 + PSFVM(:,:)**2)**1.5) /PDZZ(:,:,IKE) ZSOURCE (:,:,IKE)=ZSOURCE(:,:,IKE)-1.E2*((PSFUM(:,:)**2 + PSFVM(:,:)**2)**1.5) /PDZZ(:,:,IKE)
END IF END IF
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment