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

S. Riette 17 Nov 2022: condensation re-plugged in lima_adjust_split

NOT TESTED!
parent dbb44353
No related branches found
No related tags found
No related merge requests found
...@@ -992,7 +992,7 @@ SELECT CASE ( HCLOUD ) ...@@ -992,7 +992,7 @@ SELECT CASE ( HCLOUD )
PTHS,PRS, PSVS(:,:,:,NSV_LIMA_BEG:NSV_LIMA_END), & PTHS,PRS, PSVS(:,:,:,NSV_LIMA_BEG:NSV_LIMA_END), &
PCLDFR, PSRCS ) PCLDFR, PSRCS )
ELSE IF (LPTSPLIT) THEN ELSE IF (LPTSPLIT) THEN
CALL LIMA_ADJUST_SPLIT(KRR, KMI, TPFILE, CCONDENS, CLAMBDA3, & CALL LIMA_ADJUST_SPLIT(YLDIMPHYEX, KRR, KMI, TPFILE, CCONDENS, CLAMBDA3, &
OSUBG_COND, OSIGMAS, PTSTEP, PSIGQSAT, & OSUBG_COND, OSIGMAS, PTSTEP, PSIGQSAT, &
PRHODREF, PRHODJ, PEXNREF, PPABST, PSIGS, PMFCONV, PPABST, ZZZ, & PRHODREF, PRHODJ, PEXNREF, PPABST, PSIGS, PMFCONV, PPABST, ZZZ, &
PDTHRAD, PW_ACT, & PDTHRAD, PW_ACT, &
......
...@@ -9,7 +9,7 @@ ...@@ -9,7 +9,7 @@
! !
INTERFACE INTERFACE
! !
SUBROUTINE LIMA_ADJUST_SPLIT(KRR, KMI, TPFILE, HCONDENS, HLAMBDA3, & SUBROUTINE LIMA_ADJUST_SPLIT(D, KRR, KMI, TPFILE, HCONDENS, HLAMBDA3, &
OSUBG_COND, OSIGMAS, PTSTEP, PSIGQSAT, & OSUBG_COND, OSIGMAS, PTSTEP, PSIGQSAT, &
PRHODREF, PRHODJ, PEXNREF, PPABSM, PSIGS, PMFCONV, & PRHODREF, PRHODJ, PEXNREF, PPABSM, PSIGS, PMFCONV, &
PPABST, PZZ, PDTHRAD, PW_NU, & PPABST, PZZ, PDTHRAD, PW_NU, &
...@@ -18,7 +18,9 @@ INTERFACE ...@@ -18,7 +18,9 @@ INTERFACE
! !
USE MODD_IO, ONLY: TFILEDATA USE MODD_IO, ONLY: TFILEDATA
USE MODD_NSV, only: NSV_LIMA_BEG USE MODD_NSV, only: NSV_LIMA_BEG
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
! !
TYPE(DIMPHYEX_t), INTENT(IN) :: D
INTEGER, INTENT(IN) :: KRR ! Number of moist variables INTEGER, INTENT(IN) :: KRR ! Number of moist variables
INTEGER, INTENT(IN) :: KMI ! Model index INTEGER, INTENT(IN) :: KMI ! Model index
TYPE(TFILEDATA), INTENT(IN) :: TPFILE ! Output file TYPE(TFILEDATA), INTENT(IN) :: TPFILE ! Output file
...@@ -68,7 +70,7 @@ END INTERFACE ...@@ -68,7 +70,7 @@ END INTERFACE
END MODULE MODI_LIMA_ADJUST_SPLIT END MODULE MODI_LIMA_ADJUST_SPLIT
! !
! ########################################################################### ! ###########################################################################
SUBROUTINE LIMA_ADJUST_SPLIT(KRR, KMI, TPFILE, HCONDENS, HLAMBDA3, & SUBROUTINE LIMA_ADJUST_SPLIT(D, KRR, KMI, TPFILE, HCONDENS, HLAMBDA3, &
OSUBG_COND, OSIGMAS, PTSTEP, PSIGQSAT, & OSUBG_COND, OSIGMAS, PTSTEP, PSIGQSAT, &
PRHODREF, PRHODJ, PEXNREF, PPABSM, PSIGS, PMFCONV, & PRHODREF, PRHODJ, PEXNREF, PPABSM, PSIGS, PMFCONV, &
PPABST, PZZ, PDTHRAD, PW_NU, & PPABST, PZZ, PDTHRAD, PW_NU, &
...@@ -162,6 +164,10 @@ USE MODD_PARAM_LIMA ...@@ -162,6 +164,10 @@ USE MODD_PARAM_LIMA
USE MODD_PARAM_LIMA_COLD USE MODD_PARAM_LIMA_COLD
USE MODD_PARAM_LIMA_MIXED USE MODD_PARAM_LIMA_MIXED
USE MODD_PARAM_LIMA_WARM USE MODD_PARAM_LIMA_WARM
USE MODD_RAIN_ICE_PARAM, ONLY: RAIN_ICE_PARAM
USE MODD_NEB, ONLY: NEB
USE MODD_TURB_n, ONLY: TURBN
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
! !
use mode_budget, only: Budget_store_init, Budget_store_end use mode_budget, only: Budget_store_init, Budget_store_end
USE MODE_IO_FIELD_WRITE, only: IO_Field_write USE MODE_IO_FIELD_WRITE, only: IO_Field_write
...@@ -178,6 +184,7 @@ IMPLICIT NONE ...@@ -178,6 +184,7 @@ IMPLICIT NONE
!* 0.1 Declarations of dummy arguments : !* 0.1 Declarations of dummy arguments :
! !
! !
TYPE(DIMPHYEX_t), INTENT(IN) :: D
INTEGER, INTENT(IN) :: KRR ! Number of moist variables INTEGER, INTENT(IN) :: KRR ! Number of moist variables
INTEGER, INTENT(IN) :: KMI ! Model index INTEGER, INTENT(IN) :: KMI ! Model index
TYPE(TFILEDATA), INTENT(IN) :: TPFILE ! Output file TYPE(TFILEDATA), INTENT(IN) :: TPFILE ! Output file
...@@ -510,13 +517,13 @@ DO JITER =1,ITERMAX ...@@ -510,13 +517,13 @@ DO JITER =1,ITERMAX
ZRI=0. ZRI=0.
ZSIGS=PSIGS ZSIGS=PSIGS
ZSIGQSAT2D(:,:)=PSIGQSAT ZSIGQSAT2D(:,:)=PSIGQSAT
!CALL CONDENSATION(IIU, IJU, IKU, IIB, IIE, IJB, IJE, IKB, IKE, 1, & CALL CONDENSATION(D, CST, RAIN_ICE_PARAM, NEB, TURBN, &
! 'S', HCONDENS, HLAMBDA3, & &'S', HCONDENS, HLAMBDA3, &
! PPABST, PZZ, PRHODREF, ZT, ZRV_IN, ZRV, ZRC_IN, ZRC, ZRI_IN, ZRI,& &PPABST, PZZ, PRHODREF, ZT, ZRV_IN, ZRV, ZRC_IN, ZRC, ZRI_IN, ZRI, &
! PRRS*PTSTEP, PRSS*PTSTEP, PRGS*PTSTEP, ZSIGS, PMFCONV, PCLDFR, & &PRRS*PTSTEP, PRSS*PTSTEP, PRGS*PTSTEP, ZSIGS, .FALSE., PMFCONV, PCLDFR, PSRCS, .FALSE., &
! PSRCS, .FALSE., OSIGMAS, .FALSE., .FALSE., & &OSIGMAS, .FALSE., .FALSE., &
! ZDUM, ZDUM, ZDUM, ZDUM, ZDUM, ZSIGQSAT2D, & &ZDUM, ZDUM, ZDUM, ZDUM, ZDUM, ZSIGQSAT2D, &
! PLV=ZLV, PLS=ZLS, PCPH=ZCPH ) &ZLV, ZLS, ZCPH)
PCLDFR(:,:,:) = MIN(PCLDFR(:,:,:) + PCF_MF(:,:,:) , 1.) PCLDFR(:,:,:) = MIN(PCLDFR(:,:,:) + PCF_MF(:,:,:) , 1.)
ZRV(:,:,:) = ZRV(:,:,:) - MAX(MIN(PRC_MF(:,:,:), ZRV(:,:,:)),0.) ZRV(:,:,:) = ZRV(:,:,:) - MAX(MIN(PRC_MF(:,:,:), ZRV(:,:,:)),0.)
ZRC(:,:,:) = ZRC(:,:,:) + MAX(MIN(PRC_MF(:,:,:), ZRV(:,:,:)),0.) ZRC(:,:,:) = ZRC(:,:,:) + MAX(MIN(PRC_MF(:,:,:), ZRV(:,:,:)),0.)
......
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