source: LMDZ4/branches/IPSL-CM4_IPCC_branch/libf/phylmd/conemav.F @ 5425

Last change on this file since 5425 was 524, checked in by lmdzadmin, 21 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.3 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE conemav (dtime,paprs,pplay,t,q,u,v,tra,ntra,
5     .             work1,work2,d_t,d_q,d_u,d_v,d_tra,
6     .             rain, snow, kbas, ktop,
7     .             upwd,dnwd,dnwdbis,Ma,cape,tvp,iflag,
8     .             pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr)
9 
10c
11      IMPLICIT none
12c======================================================================
13c Auteur(s): Z.X. Li (LMD/CNRS) date: 19930818
14c Objet: schema de convection de Emanuel (1991) interface
15c======================================================================
16c Arguments:
17c dtime--input-R-pas d'integration (s)
18c s-------input-R-la valeur "s" pour chaque couche
19c sigs----input-R-la valeur "sigma" de chaque couche
20c sig-----input-R-la valeur de "sigma" pour chaque niveau
21c psolpa--input-R-la pression au sol (en Pa)
22C pskapa--input-R-exponentiel kappa de psolpa
23c h-------input-R-enthalpie potentielle (Cp*T/P**kappa)
24c q-------input-R-vapeur d'eau (en kg/kg)
25c
26c work*: input et output: deux variables de travail,
27c                            on peut les mettre a 0 au debut
28c ALE-----input-R-energie disponible pour soulevement
29c
30C d_h-----output-R-increment de l'enthalpie potentielle (h)
31c d_q-----output-R-increment de la vapeur d'eau
32c rain----output-R-la pluie (mm/s)
33c snow----output-R-la neige (mm/s)
34c upwd----output-R-saturated updraft mass flux (kg/m**2/s)
35c dnwd----output-R-saturated downdraft mass flux (kg/m**2/s)
36c dnwd0---output-R-unsaturated downdraft mass flux (kg/m**2/s)
37c Cape----output-R-CAPE (J/kg)
38c Tvp-----output-R-Temperature virtuelle d'une parcelle soulevee
39c                  adiabatiquement a partir du niveau 1 (K)
40c deltapb-output-R-distance entre LCL et base de la colonne (<0 ; Pa)
41c Ice_flag-input-L-TRUE->prise en compte de la thermodynamique de la glace
42c======================================================================
43c
44#include "dimensions.h"
45#include "dimphy.h"
46c
47      integer NTRAC
48      PARAMETER (NTRAC=nqmx-2)
49c
50       REAL dtime, paprs(klon,klev+1),pplay(klon,klev)
51       REAL t(klon,klev),q(klon,klev),u(klon,klev),v(klon,klev)
52       REAL tra(klon,klev,ntrac)
53       INTEGER ntra
54       REAL work1(klon,klev),work2(klon,klev)
55c
56       REAL d_t(klon,klev),d_q(klon,klev),d_u(klon,klev),d_v(klon,klev)
57       REAL d_tra(klon,klev,ntrac)
58       REAL rain(klon),snow(klon)
59c
60       INTEGER kbas(klon),ktop(klon)
61       REAL em_ph(klon,klev+1),em_p(klon,klev)
62       REAL upwd(klon,klev),dnwd(klon,klev),dnwdbis(klon,klev)
63       REAL Ma(klon,klev),cape(klon),tvp(klon,klev)
64       INTEGER iflag(klon)
65       REAL rflag(klon)
66       REAL pbase(klon),bbase(klon)
67       REAL dtvpdt1(klon,klev),dtvpdq1(klon,klev)
68       REAL dplcldt(klon),dplcldr(klon)
69c
70       REAL zx_t,zdelta,zx_qs,zcor
71c
72       INTEGER noff, minorig
73       INTEGER i,k,itra
74       REAL qs(klon,klev)
75       REAL cbmf(klon)
76       SAVE cbmf
77       INTEGER ifrst
78       SAVE ifrst
79       DATA ifrst /0/
80#include "YOMCST.h"
81#include "YOETHF.h"
82#include "FCTTRE.h"
83c
84c
85      IF (ifrst .EQ. 0) THEN
86         ifrst = 1
87         DO i = 1, klon
88          cbmf(i) = 0.
89         ENDDO
90      ENDIF
91
92      DO k = 1, klev+1
93         DO i=1,klon
94         em_ph(i,k) = paprs(i,k) / 100.0
95      ENDDO
96      ENDDO
97c
98      DO k = 1, klev
99         DO i=1,klon
100         em_p(i,k) = pplay(i,k) / 100.0
101      ENDDO
102      ENDDO
103
104c
105      DO k = 1, klev
106        DO i = 1, klon
107         zx_t = t(i,k)
108         zdelta=MAX(0.,SIGN(1.,rtt-zx_t))
109         zx_qs= MIN(0.5 , r2es * FOEEW(zx_t,zdelta)/em_p(i,k)/100.0)
110         zcor=1./(1.-retv*zx_qs)
111         qs(i,k)=zx_qs*zcor
112        ENDDO
113      ENDDO
114c
115      noff = 2
116      minorig = 2
117      CALL convect1(klon,klev,klev+1,noff,minorig,t,q,qs,u,v,
118     $              em_p,em_ph,iflag,
119     $              d_t,d_q,d_u,d_v,rain,cbmf,dtime,Ma)
120c
121      DO i = 1,klon
122        rain(i) = rain(i)/86400.
123        rflag(i)=iflag(i)
124      ENDDO
125c      call dump2d(iim,jjm-1,rflag(2:klon-1),'FLAG CONVECTION   ')
126c     if (klon.eq.1) then
127c        print*,'IFLAG ',iflag
128c     else
129c        write(*,'(96i1)') (iflag(i),i=2,klon-1)
130c     endif
131      DO k = 1, klev
132        DO i = 1, klon
133           d_t(i,k) = dtime*d_t(i,k)
134           d_q(i,k) = dtime*d_q(i,k)
135           d_u(i,k) = dtime*d_u(i,k)
136           d_v(i,k) = dtime*d_v(i,k)
137        ENDDO
138        DO itra = 1,ntra
139          DO i = 1, klon
140            d_tra(i,k,itra) = 0.
141          ENDDO
142        ENDDO
143      ENDDO
144 
145c
146c
147c
148      RETURN
149      END
150 
Note: See TracBrowser for help on using the repository browser.