Skip to content
Snippets Groups Projects
Commit 4a79e573 authored by RODIER Quentin's avatar RODIER Quentin
Browse files

P.Tulet 30/06/2021: bugfix diag with chemistry variables : no time splitted...

P.Tulet 30/06/2021: bugfix diag with chemistry variables : no time splitted needed at DIAG so SVT must be computed with XSVT and GSPLIT=F
parent c67c4b14
No related branches found
No related tags found
No related merge requests found
......@@ -609,10 +609,15 @@ ZDTSOLVER = PTSTEP / NCH_SUBSTEPS
!
IF (LORILAM) THEN
ALLOCATE( ZSVT(SIZE(XSVT,1),SIZE(XSVT,2),SIZE(XSVT,3),SIZE(XSVT,4)) )
DO JSV = 1, SIZE(XSVT,4)
ZSVT(:,:,:,JSV) = XRSVS(:,:,:,JSV) *PTSTEP / XRHODJ(:,:,:)
END DO
IF (CPROGRAM /='DIAG ') THEN
DO JSV = 1, SIZE(XSVT,4)
ZSVT(:,:,:,JSV) = XRSVS(:,:,:,JSV) *PTSTEP / XRHODJ(:,:,:)
END DO
ELSE
DO JSV = 1, SIZE(XSVT,4)
ZSVT(:,:,:,JSV) = XSVT(:,:,:,JSV)
END DO
END IF
ZSVT(:,:,:,NSV_CHEMBEG:NSV_CHEMEND) = MAX(ZSVT(:,:,:,NSV_CHEMBEG:NSV_CHEMEND), XMNH_TINY)
ZSVT(:,:,:,NSV_AERBEG:NSV_AEREND) = MAX(ZSVT(:,:,:,NSV_AERBEG:NSV_AEREND), XMNH_TINY)
!
......@@ -640,9 +645,7 @@ SELECT CASE (CCH_TDISCRETIZATION)
END SELECT
!
!
IF (CPROGRAM=='DIAG ') THEN
IF (LEN_TRIM(CSPEC_BU_DIAG)/=0.OR.LEN_TRIM(CSPEC_DIAG)/=0) GSPLIT=.FALSE. ! Modif. for DIAG
END IF
IF (CPROGRAM=='DIAG ') GSPLIT=.FALSE. ! Modif. for DIAG
!
!
!* 1.6 allocate tables
......
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