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

Christine 26/02/2020: missing condition for budget on RC and SV in drag_veg;...

Christine 26/02/2020: missing condition for budget on RC and SV in drag_veg; missing allocation of INPRC and ACPRC with deposition; correction of typo on budget SV CEVA_BU_RSV instead of REVA_BU_RSV
parent a832dc9e
No related branches found
No related tags found
No related merge requests found
...@@ -70,6 +70,7 @@ SUBROUTINE DRAG_VEG(PTSTEP,PUT,PVT,PTKET,ODEPOTREE, PVDEPOTREE, & ...@@ -70,6 +70,7 @@ SUBROUTINE DRAG_VEG(PTSTEP,PUT,PVT,PTKET,ODEPOTREE, PVDEPOTREE, &
!! C.Lac 07/2016 : Add droplet deposition !! C.Lac 07/2016 : Add droplet deposition
!! C.Lac 10/2017 : Correction on deposition !! C.Lac 10/2017 : Correction on deposition
!! C.Lac 11/2019 : Correction in the drag formula and application to building in addition to tree !! C.Lac 11/2019 : Correction in the drag formula and application to building in addition to tree
!! C.Lac 02/2020 : Correction missing condition for budget on RC and SV
!!--------------------------------------------------------------- !!---------------------------------------------------------------
! !
! !
...@@ -285,8 +286,10 @@ SUBROUTINE DRAG_VEG(PTSTEP,PUT,PVT,PTKET,ODEPOTREE, PVDEPOTREE, & ...@@ -285,8 +286,10 @@ SUBROUTINE DRAG_VEG(PTSTEP,PUT,PVT,PTKET,ODEPOTREE, PVDEPOTREE, &
! !
IF (LBUDGET_U) CALL BUDGET (PRUS,1,'DRAG_BU_RU') IF (LBUDGET_U) CALL BUDGET (PRUS,1,'DRAG_BU_RU')
IF (LBUDGET_V) CALL BUDGET (PRVS,2,'DRAG_BU_RV') IF (LBUDGET_V) CALL BUDGET (PRVS,2,'DRAG_BU_RV')
IF (LBUDGET_RC) CALL BUDGET (PRRS(:,:,:,2),7,'DEPOTR_BU_RRC') IF (ODEPOTREE) THEN
IF (LBUDGET_SV) CALL BUDGET (PSVS(:,:,:,NSV_C2R2BEG+1),14+(NSV_C2R2BEG-1),'DEPOTR_BU_RSV') IF (LBUDGET_RC) CALL BUDGET (PRRS(:,:,:,2),7,'DEPOTR_BU_RRC')
IF (LBUDGET_SV) CALL BUDGET (PSVS(:,:,:,NSV_C2R2BEG+1),14+(NSV_C2R2BEG-1),'DEPOTR_BU_RSV')
END IF
! !
!* 3. Computations of TKE tendency due to canopy drag !* 3. Computations of TKE tendency due to canopy drag
! ------------------------------------------------ ! ------------------------------------------------
......
...@@ -52,7 +52,7 @@ END MODULE MODI_INI_MICRO_n ...@@ -52,7 +52,7 @@ END MODULE MODI_INI_MICRO_n
!! C.LAc 10/2016 Add budget for droplet deposition !! C.LAc 10/2016 Add budget for droplet deposition
!! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O !! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
!! P.Wautelet 01/2019: bug: add missing allocations !! P.Wautelet 01/2019: bug: add missing allocations
!! !! C.Lac 02/2020: add missing allocation of INPRC and ACPRC with deposition
!! -------------------------------------------------------------------------- !! --------------------------------------------------------------------------
! !
!* 0. DECLARATIONS !* 0. DECLARATIONS
...@@ -161,9 +161,9 @@ ELSE ...@@ -161,9 +161,9 @@ ELSE
ALLOCATE(XACPRR(0,0)) ALLOCATE(XACPRR(0,0))
END IF END IF
! !
IF (( CCLOUD(1:3) == 'ICE' .AND.LSEDIC) .OR. & IF (( CCLOUD(1:3) == 'ICE' .AND.(LSEDIC .OR. LDEPOSC)) .OR. &
((CCLOUD=='C2R2' .OR. CCLOUD=='C3R5' .OR. CCLOUD=='KHKO').AND.LSEDC) .OR. & ((CCLOUD=='C2R2' .OR. CCLOUD=='C3R5' .OR. CCLOUD=='KHKO').AND.(LSEDC .OR. LDEPOC)) .OR. &
( CCLOUD=='LIMA' .AND.MSEDC)) THEN ( CCLOUD=='LIMA' .AND.(MSEDC .OR. MDEPOC))) THEN
ALLOCATE(XINPRC(IIU,IJU)) ALLOCATE(XINPRC(IIU,IJU))
ALLOCATE(XACPRC(IIU,IJU)) ALLOCATE(XACPRC(IIU,IJU))
XINPRC(:,:)=0.0 XINPRC(:,:)=0.0
......
...@@ -1720,7 +1720,7 @@ ENDIF ...@@ -1720,7 +1720,7 @@ ENDIF
IF (LBUDGET_RV) CALL BUDGET (PRVS(:,:,:)*PRHODJ(:,:,:),6 ,'REVA_BU_RRV') IF (LBUDGET_RV) CALL BUDGET (PRVS(:,:,:)*PRHODJ(:,:,:),6 ,'REVA_BU_RRV')
IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8 ,'REVA_BU_RRR') IF (LBUDGET_RR) CALL BUDGET (PRRS(:,:,:)*PRHODJ(:,:,:),8 ,'REVA_BU_RRR')
IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4 ,'REVA_BU_RTH') IF (LBUDGET_TH) CALL BUDGET (PTHS(:,:,:)*PRHODJ(:,:,:),4 ,'REVA_BU_RTH')
IF (LBUDGET_SV) CALL BUDGET (PCRS(:,:,:)*PRHODJ(:,:,:),15+(NSV_C2R2BEG-1),'CEVA_BU_RSV') IF (LBUDGET_SV) CALL BUDGET (PCRS(:,:,:)*PRHODJ(:,:,:),15+(NSV_C2R2BEG-1),'REVA_BU_RSV')
! !
END SUBROUTINE C2R2_KHKO_EVAPORATION END SUBROUTINE C2R2_KHKO_EVAPORATION
! !
......
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