source: LMDZ6/branches/SETHET_DECOUPLE/libf/phylmd/cv3_mixscale.F90 @ 5446

Last change on this file since 5446 was 3496, checked in by jyg, 6 years ago

Implementation of the ejection of liquid precipitation from the adiabatic ascents.
New flags:
+cvflag_prec_eject: logical

n -> old code, y -> new code

+ejectliq: real; possible values 0. & 1.

  1. -> no liquid precipitation is ejected
  2. -> all liquid precipitation is ejected

+ejectice: real; any value between 0. and 1.

fraction of solid precipitation ejected at each level

Note that the adiabatic ascent mass flux decrease due to precipitation ejection is not taken into account.

Attempts to do it led to water conservation violation.
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 884 bytes
Line 
1SUBROUTINE cv3_mixscale(nloc, ncum, na, ment, m)
2  ! **************************************************************
3  ! *
4  ! CV3_MIXSCALE                                                *
5  ! *
6  ! *
7  ! written by   : Jean-Yves Grandpeix, 30/05/2003, 16.34.37    *
8  ! modified by :                                               *
9  ! **************************************************************
10
11  IMPLICIT NONE
12
13  include "cv3param.h"
14
15!inputs:
16  INTEGER, INTENT (IN)                               :: ncum, na, nloc
17  REAL, DIMENSION (nloc, na), INTENT (IN)            :: m
18!input/outputs:
19  REAL, DIMENSION (nloc, na, na), INTENT (INOUT)     :: ment
20
21!local variables:
22  INTEGER i, j, il
23
24    DO j = 1, nl
25      DO i = 1, nl
26        DO il = 1, ncum
27          ment(il, i, j) = m(il, i)*ment(il, i, j)
28        END DO
29      END DO
30    END DO
31
32
33  RETURN
34END SUBROUTINE cv3_mixscale
Note: See TracBrowser for help on using the repository browser.