From 9c3a6a9cf6242acbabf6293b33a3d583f1f535b1 Mon Sep 17 00:00:00 2001 From: Aron Roland Date: Wed, 22 Jan 2025 11:23:40 +0100 Subject: [PATCH] ww3_gse2: latest work ... --- model/src/w3profsmd_pdlib.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/model/src/w3profsmd_pdlib.F90 b/model/src/w3profsmd_pdlib.F90 index 3d5990f7b..0b9f01f2a 100644 --- a/model/src/w3profsmd_pdlib.F90 +++ b/model/src/w3profsmd_pdlib.F90 @@ -1767,7 +1767,7 @@ SUBROUTINE PDLIB_W3XYPFSFCT2 ( ISP, C, LCALC, RD10, RD20, DT, ACIN, ACOUT) END IF END IF ST(NI) = ST(NI) + THETA_L(:,IE) ! the 2nd term are the theta values of each node ... - THETA_H = (1./3.+DT/(2.*PDLIB_TRIA(IE)) * KELEM(:,IE) ) * FT ! LAX + THETA_H = (1./3.+DT/ITER(IK,ITH)/(2*PDLIB_TRIA(IE)) * KELEM(:,IE) ) * FT ! LAX ! THETA_H = (1./3.+2./3.*KELEM(:,IE)/SUM(MAX(ZERO,KELEM(:,IE))))*FT ! CENTRAL ... can be tested as well a bit more dispersive then LAX THETA_ACE(:,IE) = THETA_H-THETA_L(:,IE) PP(NI) = PP(NI) + MAX(ZERO, -THETA_ACE(:,IE)) * DTSI(NI)