Skip to content
Snippets Groups Projects
Commit 33e53c4f authored by RIETTE Sébastien's avatar RIETTE Sébastien Committed by RODIER Quentin
Browse files

S. Riette 11 Jan 2024: BF (fixes lima crash on belenos)

parent 92f3b703
Branches
Tags
No related merge requests found
......@@ -424,10 +424,10 @@ IF( INUCT >= 1 ) THEN
ZZW1(:) = MIN(XCSTDCRIT*ZZW6(:)/(((ZZT(:)*ZSMAX(:))**3)*ZRHODREF(:)),1.E-5)
END WHERE
!
PTOT_RV_HENU(:,:,:) = 0.
IF(PRESENT(PTOT_RV_HENU)) PTOT_RV_HENU(:,:,:) = 0.
IF (.NOT.LSUBG_COND) THEN
ZW(:,:,:) = MIN( UNPACK( ZZW1(:),MASK=GNUCT(:,:,:),FIELD=0.0 ),PRVT(:,:,:) )
PTOT_RV_HENU(:,:,:) = ZW(:,:,:)
IF(PRESENT(PTOT_RV_HENU)) PTOT_RV_HENU(:,:,:) = ZW(:,:,:)
PTHT(:,:,:) = PTHT(:,:,:) + ZW(:,:,:) * (CST%XLVTT+(CST%XCPV-CST%XCL)*(PT(:,:,:)-CST%XTT))/ &
(PEXNREF(:,:,:)*(CST%XCPD+CST%XCPV*PRVT(:,:,:)+CST%XCL*(PRCT(:,:,:)+PRRT(:,:,:))))
PRVT(:,:,:) = PRVT(:,:,:) - ZW(:,:,:)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment