diff --git a/src/MNH/resolved_cloud.f90 b/src/MNH/resolved_cloud.f90 index 4aebcdb849ba9534050bec97b7a9ceadca2e0992..a5437038beb7651763b8d4e7d85c47ec3e209cd6 100644 --- a/src/MNH/resolved_cloud.f90 +++ b/src/MNH/resolved_cloud.f90 @@ -449,7 +449,7 @@ INTEGER :: JK,JI,JL ! ! REAL, DIMENSION(SIZE(PZZ,1),SIZE(PZZ,2),SIZE(PZZ,3)):: ZDZZ -REAL, DIMENSION(SIZE(PZZ,1),SIZE(PZZ,2),SIZE(PZZ,3)):: ZT,ZEXN,ZLV,ZLS,ZCPH +real, dimension(:,:,:), allocatable :: ZEXN REAL, DIMENSION(SIZE(PZZ,1),SIZE(PZZ,2),SIZE(PZZ,3)):: ZCOR ! for the correction of negative rv REAL, DIMENSION(SIZE(PZZ,1),SIZE(PZZ,2),SIZE(PZZ,3)):: ZZZ @@ -575,6 +575,7 @@ IF (HCLOUD == 'C2R2' .OR. HCLOUD == 'C3R5' .OR. HCLOUD == 'KHKO' & PSVT(:,:,IKE+1,ISVBEG:ISVEND) = PSVT(:,:,IKE,ISVBEG:ISVEND) ENDIF ! +! !* 3. REMOVE NEGATIVE VALUES ! ---------------------- ! @@ -721,6 +722,8 @@ SELECT CASE ( HCLOUD ) !* 9. MIXED-PHASE MICROPHYSICAL SCHEME (WITH 3 ICE SPECIES) ! ----------------------------------------------------- ! + allocate( zexn( size( pzz, 1 ), size( pzz, 2 ), size( pzz, 3 ) ) ) + ZEXN(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD) ! !* 9.1 Compute the explicit microphysical sources ! @@ -794,12 +797,16 @@ SELECT CASE ( HCLOUD ) PRS=PRS(:,:,:,5)*PTSTEP, & PRG=PRS(:,:,:,6)*PTSTEP ) END IF + + deallocate( zexn ) ! CASE ('ICE4') ! !* 10. MIXED-PHASE MICROPHYSICAL SCHEME (WITH 4 ICE SPECIES) ! ----------------------------------------------------- ! + allocate( zexn( size( pzz, 1 ), size( pzz, 2 ), size( pzz, 3 ) ) ) + ZEXN(:,:,:)= (PPABST(:,:,:)/XP00)**(XRD/XCPD) ! !* 10.1 Compute the explicit microphysical sources ! @@ -880,6 +887,8 @@ SELECT CASE ( HCLOUD ) PRG=PRS(:,:,:,6)*PTSTEP, & PRH=PRS(:,:,:,7)*PTSTEP ) END IF + + deallocate( zexn ) ! ! !* 12. 2-MOMENT MIXED-PHASE MICROPHYSICAL SCHEME LIMA