source: LMDZ4/branches/V3_test/libf/phylmd/cv_driver.F @ 5373

Last change on this file since 5373 was 704, checked in by Laurent Fairhead, 18 years ago

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 23.8 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE cv_driver(len,nd,ndp1,ntra,iflag_con,
5     &                   t1,q1,qs1,u1,v1,tra1,
6     &                   p1,ph1,iflag1,ft1,fq1,fu1,fv1,ftra1,
7     &                   precip1,VPrecip1,
8     &                   cbmf1,sig1,w01,
9     &                   icb1,inb1,
10     &                   delt,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1,
11     &                   da1,phi1,mp1)
12C
13      USE dimphy
14      implicit none
15C
16C.............................START PROLOGUE............................
17C
18C PARAMETERS:
19C      Name            Type         Usage            Description
20C   ----------      ----------     -------  ----------------------------
21C
22C      len           Integer        Input        first (i) dimension
23C      nd            Integer        Input        vertical (k) dimension
24C      ndp1          Integer        Input        nd + 1
25C      ntra          Integer        Input        number of tracors
26C      iflag_con     Integer        Input        version of convect (3/4)
27C      t1            Real           Input        temperature
28C      q1            Real           Input        specific hum
29C      qs1           Real           Input        sat specific hum
30C      u1            Real           Input        u-wind
31C      v1            Real           Input        v-wind
32C      tra1          Real           Input        tracors
33C      p1            Real           Input        full level pressure
34C      ph1           Real           Input        half level pressure
35C      iflag1        Integer        Output       flag for Emanuel conditions
36C      ft1           Real           Output       temp tend
37C      fq1           Real           Output       spec hum tend
38C      fu1           Real           Output       u-wind tend
39C      fv1           Real           Output       v-wind tend
40C      ftra1         Real           Output       tracor tend
41C      precip1       Real           Output       precipitation
42C      VPrecip1      Real           Output       vertical profile of precipitations
43C      cbmf1         Real           Output       cloud base mass flux
44C      sig1          Real           In/Out       section adiabatic updraft
45C      w01           Real           In/Out       vertical velocity within adiab updraft
46C      delt          Real           Input        time step
47C      Ma1           Real           Output       mass flux adiabatic updraft
48C      upwd1         Real           Output       total upward mass flux (adiab+mixed)
49C      dnwd1         Real           Output       saturated downward mass flux (mixed)
50C      dnwd01        Real           Output       unsaturated downward mass flux
51C      qcondc1       Real           Output       in-cld mixing ratio of condensed water
52C      wd1           Real           Output       downdraft velocity scale for sfc fluxes
53C      cape1         Real           Output       CAPE
54C
55C S. Bony, Mar 2002:
56C       * Several modules corresponding to different physical processes
57C       * Several versions of convect may be used:
58C               - iflag_con=3: version lmd  (previously named convect3)
59C               - iflag_con=4: version 4.3b (vect. version, previously convect1/2)
60C   + tard:     - iflag_con=5: version lmd with ice (previously named convectg)
61C S. Bony, Oct 2002:
62C       * Vectorization of convect3 (ie version lmd)
63C
64C..............................END PROLOGUE.............................
65c
66c
67cym#include "dimensions.h"
68cym#include "dimphy.h"
69
70      integer len
71      integer nd
72      integer ndp1
73      integer noff
74      integer iflag_con
75      integer ntra
76      real t1(len,nd)
77      real q1(len,nd)
78      real qs1(len,nd)
79      real u1(len,nd)
80      real v1(len,nd)
81      real p1(len,nd)
82      real ph1(len,ndp1)
83      integer iflag1(len)
84      real ft1(len,nd)
85      real fq1(len,nd)
86      real fu1(len,nd)
87      real fv1(len,nd)
88      real precip1(len)
89      real cbmf1(len)
90      real VPrecip1(len,nd+1)
91      real Ma1(len,nd)
92      real upwd1(len,nd)
93      real dnwd1(len,nd)
94      real dnwd01(len,nd)
95
96      real qcondc1(len,nd)     ! cld
97      real wd1(len)            ! gust
98      real cape1(len)     
99
100      real da1(len,nd),phi1(len,nd,nd),mp1(len,nd)
101      real da(len,nd),phi(len,nd,nd),mp(len,nd)
102      real tra1(len,nd,ntra)
103      real ftra1(len,nd,ntra)
104
105      real delt
106
107!-------------------------------------------------------------------
108! --- ARGUMENTS
109!-------------------------------------------------------------------
110! --- On input:
111!
112!  t:   Array of absolute temperature (K) of dimension ND, with first
113!       index corresponding to lowest model level. Note that this array
114!       will be altered by the subroutine if dry convective adjustment
115!       occurs and if IPBL is not equal to 0.
116!
117!  q:   Array of specific humidity (gm/gm) of dimension ND, with first
118!       index corresponding to lowest model level. Must be defined
119!       at same grid levels as T. Note that this array will be altered
120!       if dry convective adjustment occurs and if IPBL is not equal to 0.
121!
122!  qs:  Array of saturation specific humidity of dimension ND, with first
123!       index corresponding to lowest model level. Must be defined
124!       at same grid levels as T. Note that this array will be altered
125!       if dry convective adjustment occurs and if IPBL is not equal to 0.
126!
127!  u:   Array of zonal wind velocity (m/s) of dimension ND, witth first
128!       index corresponding with the lowest model level. Defined at
129!       same levels as T. Note that this array will be altered if
130!       dry convective adjustment occurs and if IPBL is not equal to 0.
131!
132!  v:   Same as u but for meridional velocity.
133!
134!  tra: Array of passive tracer mixing ratio, of dimensions (ND,NTRA),
135!       where NTRA is the number of different tracers. If no
136!       convective tracer transport is needed, define a dummy
137!       input array of dimension (ND,1). Tracers are defined at
138!       same vertical levels as T. Note that this array will be altered
139!       if dry convective adjustment occurs and if IPBL is not equal to 0.
140!
141!  p:   Array of pressure (mb) of dimension ND, with first
142!       index corresponding to lowest model level. Must be defined
143!       at same grid levels as T.
144!
145!  ph:  Array of pressure (mb) of dimension ND+1, with first index
146!       corresponding to lowest level. These pressures are defined at
147!       levels intermediate between those of P, T, Q and QS. The first
148!       value of PH should be greater than (i.e. at a lower level than)
149!       the first value of the array P.
150!
151!  nl:  The maximum number of levels to which convection can penetrate, plus 1.
152!       NL MUST be less than or equal to ND-1.
153!
154!  delt: The model time step (sec) between calls to CONVECT
155!
156!----------------------------------------------------------------------------
157! ---   On Output:
158!
159!  iflag: An output integer whose value denotes the following:
160!       VALUE   INTERPRETATION
161!       -----   --------------
162!         0     Moist convection occurs.
163!         1     Moist convection occurs, but a CFL condition
164!               on the subsidence warming is violated. This
165!               does not cause the scheme to terminate.
166!         2     Moist convection, but no precip because ep(inb) lt 0.0001
167!         3     No moist convection because new cbmf is 0 and old cbmf is 0.
168!         4     No moist convection; atmosphere is not
169!               unstable
170!         6     No moist convection because ihmin le minorig.
171!         7     No moist convection because unreasonable
172!               parcel level temperature or specific humidity.
173!         8     No moist convection: lifted condensation
174!               level is above the 200 mb level.
175!         9     No moist convection: cloud base is higher
176!               then the level NL-1.
177!
178!  ft:   Array of temperature tendency (K/s) of dimension ND, defined at same
179!        grid levels as T, Q, QS and P.
180!
181!  fq:   Array of specific humidity tendencies ((gm/gm)/s) of dimension ND,
182!        defined at same grid levels as T, Q, QS and P.
183!
184!  fu:   Array of forcing of zonal velocity (m/s^2) of dimension ND,
185!        defined at same grid levels as T.
186!
187!  fv:   Same as FU, but for forcing of meridional velocity.
188!
189!  ftra: Array of forcing of tracer content, in tracer mixing ratio per
190!        second, defined at same levels as T. Dimensioned (ND,NTRA).
191!
192!  precip: Scalar convective precipitation rate (mm/day).
193!
194!  VPrecip: Vertical profile of convective precipitation (kg/m2/s).
195!
196!  wd:   A convective downdraft velocity scale. For use in surface
197!        flux parameterizations. See convect.ps file for details.
198!
199!  tprime: A convective downdraft temperature perturbation scale (K).
200!          For use in surface flux parameterizations. See convect.ps
201!          file for details.
202!
203!  qprime: A convective downdraft specific humidity
204!          perturbation scale (gm/gm).
205!          For use in surface flux parameterizations. See convect.ps
206!          file for details.
207!
208!  cbmf: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE MUST
209!        BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT AT
210!        ITS NEXT CALL. That is, the value of CBMF must be "remembered"
211!        by the calling program between calls to CONVECT.
212!
213!  det:   Array of detrainment mass flux of dimension ND.
214!
215!-------------------------------------------------------------------
216c
217c  Local arrays
218c
219
220      integer i,k,n,il,j
221      integer icbmax
222      integer nk1(klon)
223      integer icb1(klon)
224      integer inb1(klon)
225      integer icbs1(klon)
226
227      real plcl1(klon)
228      real tnk1(klon)
229      real qnk1(klon)
230      real gznk1(klon)
231      real pnk1(klon)
232      real qsnk1(klon)
233      real pbase1(klon)
234      real buoybase1(klon)
235
236      real lv1(klon,klev)
237      real cpn1(klon,klev)
238      real tv1(klon,klev)
239      real gz1(klon,klev)
240      real hm1(klon,klev)
241      real h1(klon,klev)
242      real tp1(klon,klev)
243      real tvp1(klon,klev)
244      real clw1(klon,klev)
245      real sig1(klon,klev)
246      real w01(klon,klev)
247      real th1(klon,klev)
248c
249      integer ncum
250c
251c (local) compressed fields:
252c
253cym      integer nloc
254cym      parameter (nloc=klon) ! pour l'instant
255#define nloc klon
256      integer idcum(nloc)
257      integer iflag(nloc),nk(nloc),icb(nloc)
258      integer nent(nloc,klev)
259      integer icbs(nloc)
260      integer inb(nloc), inbis(nloc)
261
262      real cbmf(nloc),plcl(nloc),tnk(nloc),qnk(nloc),gznk(nloc)
263      real t(nloc,klev),q(nloc,klev),qs(nloc,klev)
264      real u(nloc,klev),v(nloc,klev)
265      real gz(nloc,klev),h(nloc,klev),lv(nloc,klev),cpn(nloc,klev)
266      real p(nloc,klev),ph(nloc,klev+1),tv(nloc,klev),tp(nloc,klev)
267      real clw(nloc,klev)
268      real dph(nloc,klev)
269      real pbase(nloc), buoybase(nloc), th(nloc,klev)
270      real tvp(nloc,klev)
271      real sig(nloc,klev), w0(nloc,klev)
272      real hp(nloc,klev), ep(nloc,klev), sigp(nloc,klev)
273      real frac(nloc), buoy(nloc,klev)
274      real cape(nloc)
275      real m(nloc,klev), ment(nloc,klev,klev), qent(nloc,klev,klev)
276      real uent(nloc,klev,klev), vent(nloc,klev,klev)
277      real ments(nloc,klev,klev), qents(nloc,klev,klev)
278      real sij(nloc,klev,klev), elij(nloc,klev,klev)
279      real qp(nloc,klev), up(nloc,klev), vp(nloc,klev)
280      real wt(nloc,klev), water(nloc,klev), evap(nloc,klev)
281      real b(nloc,klev), ft(nloc,klev), fq(nloc,klev)
282      real fu(nloc,klev), fv(nloc,klev)
283      real upwd(nloc,klev), dnwd(nloc,klev), dnwd0(nloc,klev)
284      real Ma(nloc,klev), mike(nloc,klev), tls(nloc,klev)
285      real tps(nloc,klev), qprime(nloc), tprime(nloc)
286      real precip(nloc)
287      real VPrecip(nloc,klev+1)
288      real tra(nloc,klev,ntra), trap(nloc,klev,ntra)
289      real ftra(nloc,klev,ntra), traent(nloc,klev,klev,ntra)
290      real qcondc(nloc,klev)  ! cld
291      real wd(nloc)           ! gust
292
293      nent(:,:)=0
294!-------------------------------------------------------------------
295! --- SET CONSTANTS AND PARAMETERS
296!-------------------------------------------------------------------
297
298c -- set simulation flags:
299c   (common cvflag)
300
301       CALL cv_flag
302
303c -- set thermodynamical constants:
304c       (common cvthermo)
305
306       CALL cv_thermo(iflag_con)
307
308c -- set convect parameters
309c
310c       includes microphysical parameters and parameters that
311c       control the rate of approach to quasi-equilibrium)
312c       (common cvparam)
313
314      if (iflag_con.eq.3) then
315       CALL cv3_param(nd,delt)
316      endif
317
318      if (iflag_con.eq.4) then
319       CALL cv_param(nd)
320      endif
321
322!---------------------------------------------------------------------
323! --- INITIALIZE OUTPUT ARRAYS AND PARAMETERS
324!---------------------------------------------------------------------
325
326      do 20 k=1,nd
327        do 10 i=1,len
328         ft1(i,k)=0.0
329         fq1(i,k)=0.0
330         fu1(i,k)=0.0
331         fv1(i,k)=0.0
332         tvp1(i,k)=0.0
333         tp1(i,k)=0.0
334         clw1(i,k)=0.0
335cym
336         clw(i,k)=0.0   
337         gz1(i,k) = 0.
338         VPrecip1(i,k) = 0.
339         Ma1(i,k)=0.0
340         upwd1(i,k)=0.0
341         dnwd1(i,k)=0.0
342         dnwd01(i,k)=0.0
343         qcondc1(i,k)=0.0
344 10     continue
345 20   continue
346
347      do 30 j=1,ntra
348       do 31 k=1,nd
349        do 32 i=1,len
350         ftra1(i,k,j)=0.0
351 32     continue   
352 31    continue   
353 30   continue   
354
355      do 60 i=1,len
356        precip1(i)=0.0
357        iflag1(i)=0
358        wd1(i)=0.0
359        cape1(i)=0.0
360        VPrecip1(i,nd+1)=0.0
361 60   continue
362
363      if (iflag_con.eq.3) then
364        do il=1,len
365         sig1(il,nd)=sig1(il,nd)+1.
366         sig1(il,nd)=amin1(sig1(il,nd),12.1)
367        enddo
368      endif
369
370!--------------------------------------------------------------------
371! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
372!--------------------------------------------------------------------
373
374      if (iflag_con.eq.3) then
375       CALL cv3_prelim(len,nd,ndp1,t1,q1,p1,ph1            ! nd->na
376     o               ,lv1,cpn1,tv1,gz1,h1,hm1,th1)
377      endif
378
379      if (iflag_con.eq.4) then
380       CALL cv_prelim(len,nd,ndp1,t1,q1,p1,ph1
381     o               ,lv1,cpn1,tv1,gz1,h1,hm1)
382      endif
383
384!--------------------------------------------------------------------
385! --- CONVECTIVE FEED
386!--------------------------------------------------------------------
387
388      if (iflag_con.eq.3) then
389       CALL cv3_feed(len,nd,t1,q1,qs1,p1,ph1,hm1,gz1           ! nd->na
390     o         ,nk1,icb1,icbmax,iflag1,tnk1,qnk1,gznk1,plcl1)
391      endif
392
393      if (iflag_con.eq.4) then
394       CALL cv_feed(len,nd,t1,q1,qs1,p1,hm1,gz1
395     o         ,nk1,icb1,icbmax,iflag1,tnk1,qnk1,gznk1,plcl1)
396      endif
397
398!--------------------------------------------------------------------
399! --- UNDILUTE (ADIABATIC) UPDRAFT / 1st part
400! (up through ICB for convect4, up through ICB+1 for convect3)
401!     Calculates the lifted parcel virtual temperature at nk, the
402!     actual temperature, and the adiabatic liquid water content.
403!--------------------------------------------------------------------
404
405      if (iflag_con.eq.3) then
406       CALL cv3_undilute1(len,nd,t1,q1,qs1,gz1,plcl1,p1,nk1,icb1  ! nd->na
407     o                        ,tp1,tvp1,clw1,icbs1)
408      endif
409
410      if (iflag_con.eq.4) then
411       CALL cv_undilute1(len,nd,t1,q1,qs1,gz1,p1,nk1,icb1,icbmax
412     :                        ,tp1,tvp1,clw1)
413      endif
414
415!-------------------------------------------------------------------
416! --- TRIGGERING
417!-------------------------------------------------------------------
418
419      if (iflag_con.eq.3) then
420       CALL cv3_trigger(len,nd,icb1,plcl1,p1,th1,tv1,tvp1      ! nd->na
421     o                 ,pbase1,buoybase1,iflag1,sig1,w01)
422      endif
423
424      if (iflag_con.eq.4) then
425       CALL cv_trigger(len,nd,icb1,cbmf1,tv1,tvp1,iflag1)
426      endif
427
428!=====================================================================
429! --- IF THIS POINT IS REACHED, MOIST CONVECTIVE ADJUSTMENT IS NECESSARY
430!=====================================================================
431
432      ncum=0
433      do 400 i=1,len
434        if(iflag1(i).eq.0)then
435           ncum=ncum+1
436           idcum(ncum)=i
437        endif
438 400  continue
439
440c       print*,'klon, ncum = ',len,ncum
441
442      IF (ncum.gt.0) THEN
443
444!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
445! --- COMPRESS THE FIELDS
446!               (-> vectorization over convective gridpoints)
447!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
448
449      if (iflag_con.eq.3) then
450       CALL cv3_compress( len,nloc,ncum,nd,ntra
451     :    ,iflag1,nk1,icb1,icbs1
452     :    ,plcl1,tnk1,qnk1,gznk1,pbase1,buoybase1
453     :    ,t1,q1,qs1,u1,v1,gz1,th1
454     :    ,tra1
455     :    ,h1,lv1,cpn1,p1,ph1,tv1,tp1,tvp1,clw1
456     :    ,sig1,w01
457     o    ,iflag,nk,icb,icbs
458     o    ,plcl,tnk,qnk,gznk,pbase,buoybase
459     o    ,t,q,qs,u,v,gz,th
460     o    ,tra
461     o    ,h,lv,cpn,p,ph,tv,tp,tvp,clw
462     o    ,sig,w0  )
463      endif
464
465      if (iflag_con.eq.4) then
466       CALL cv_compress( len,nloc,ncum,nd
467     :    ,iflag1,nk1,icb1
468     :    ,cbmf1,plcl1,tnk1,qnk1,gznk1
469     :    ,t1,q1,qs1,u1,v1,gz1
470     :    ,h1,lv1,cpn1,p1,ph1,tv1,tp1,tvp1,clw1
471     o    ,iflag,nk,icb
472     o    ,cbmf,plcl,tnk,qnk,gznk
473     o    ,t,q,qs,u,v,gz,h,lv,cpn,p,ph,tv,tp,tvp,clw
474     o    ,dph )
475      endif
476
477!-------------------------------------------------------------------
478! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :
479! ---   FIND THE REST OF THE LIFTED PARCEL TEMPERATURES
480! ---   &
481! ---   COMPUTE THE PRECIPITATION EFFICIENCIES AND THE
482! ---   FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
483! ---   &
484! ---   FIND THE LEVEL OF NEUTRAL BUOYANCY
485!-------------------------------------------------------------------
486
487      if (iflag_con.eq.3) then
488       CALL cv3_undilute2(nloc,ncum,nd,icb,icbs,nk        !na->nd
489     :                        ,tnk,qnk,gznk,t,q,qs,gz
490     :                        ,p,h,tv,lv,pbase,buoybase,plcl
491     o                        ,inb,tp,tvp,clw,hp,ep,sigp,buoy)
492      endif
493
494      if (iflag_con.eq.4) then
495       CALL cv_undilute2(nloc,ncum,nd,icb,nk
496     :                        ,tnk,qnk,gznk,t,q,qs,gz
497     :                        ,p,dph,h,tv,lv
498     o             ,inb,inbis,tp,tvp,clw,hp,ep,sigp,frac)
499      endif
500
501!-------------------------------------------------------------------
502! --- CLOSURE
503!-------------------------------------------------------------------
504
505      if (iflag_con.eq.3) then
506       CALL cv3_closure(nloc,ncum,nd,icb,inb              ! na->nd
507     :                       ,pbase,p,ph,tv,buoy
508     o                       ,sig,w0,cape,m)
509      endif
510
511      if (iflag_con.eq.4) then
512       CALL cv_closure(nloc,ncum,nd,nk,icb
513     :                ,tv,tvp,p,ph,dph,plcl,cpn
514     o                ,iflag,cbmf)
515      endif
516
517!-------------------------------------------------------------------
518! --- MIXING
519!-------------------------------------------------------------------
520
521      if (iflag_con.eq.3) then
522       CALL cv3_mixing(nloc,ncum,nd,nd,ntra,icb,nk,inb    ! na->nd
523     :                     ,ph,t,q,qs,u,v,tra,h,lv,qnk
524     :                     ,hp,tv,tvp,ep,clw,m,sig
525     o ,ment,qent,uent,vent,sij,elij,ments,qents,traent)
526      endif
527
528      if (iflag_con.eq.4) then
529       CALL cv_mixing(nloc,ncum,nd,icb,nk,inb,inbis
530     :                     ,ph,t,q,qs,u,v,h,lv,qnk
531     :                     ,hp,tv,tvp,ep,clw,cbmf
532     o                     ,m,ment,qent,uent,vent,nent,sij,elij)
533      endif
534
535!-------------------------------------------------------------------
536! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS
537!-------------------------------------------------------------------
538
539      if (iflag_con.eq.3) then
540       CALL cv3_unsat(nloc,ncum,nd,nd,ntra,icb,inb    ! na->nd
541     :               ,t,q,qs,gz,u,v,tra,p,ph
542     :               ,th,tv,lv,cpn,ep,sigp,clw
543     :               ,m,ment,elij,delt,plcl
544     o          ,mp,qp,up,vp,trap,wt,water,evap,b)
545      endif
546
547      if (iflag_con.eq.4) then
548       CALL cv_unsat(nloc,ncum,nd,inb,t,q,qs,gz,u,v,p,ph
549     :                   ,h,lv,ep,sigp,clw,m,ment,elij
550     o                   ,iflag,mp,qp,up,vp,wt,water,evap)
551      endif
552
553!-------------------------------------------------------------------
554! --- YIELD
555!     (tendencies, precipitation, variables of interface with other
556!      processes, etc)
557!-------------------------------------------------------------------
558
559      if (iflag_con.eq.3) then
560       CALL cv3_yield(nloc,ncum,nd,nd,ntra            ! na->nd
561     :                     ,icb,inb,delt
562     :                     ,t,q,u,v,tra,gz,p,ph,h,hp,lv,cpn,th
563     :                     ,ep,clw,m,tp,mp,qp,up,vp,trap
564     :                     ,wt,water,evap,b
565     :                     ,ment,qent,uent,vent,nent,elij,traent,sig
566     :                     ,tv,tvp
567     o                     ,iflag,precip,VPrecip,ft,fq,fu,fv,ftra
568     o                     ,upwd,dnwd,dnwd0,ma,mike,tls,tps,qcondc,wd)
569      endif
570
571      if (iflag_con.eq.4) then
572       CALL cv_yield(nloc,ncum,nd,nk,icb,inb,delt
573     :              ,t,q,u,v,gz,p,ph,h,hp,lv,cpn
574     :              ,ep,clw,frac,m,mp,qp,up,vp
575     :              ,wt,water,evap
576     :              ,ment,qent,uent,vent,nent,elij
577     :              ,tv,tvp
578     o              ,iflag,wd,qprime,tprime
579     o              ,precip,cbmf,ft,fq,fu,fv,Ma,qcondc)
580      endif
581
582!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
583! --- passive tracers
584!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
585
586      if (iflag_con.eq.3) then
587       CALL cv3_tracer(nloc,len,ncum,nd,nd,
588     :                  ment,sij,da,phi)
589      endif
590
591!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
592! --- UNCOMPRESS THE FIELDS
593!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
594c set iflag1 =42 for non convective points
595      do  i=1,len
596        iflag1(i)=42
597      end do
598c
599      if (iflag_con.eq.3) then
600       CALL cv3_uncompress(nloc,len,ncum,nd,ntra,idcum
601     :          ,iflag
602     :          ,precip,VPrecip,sig,w0
603     :          ,ft,fq,fu,fv,ftra
604     :          ,inb
605     :          ,Ma,upwd,dnwd,dnwd0,qcondc,wd,cape
606     :          ,da,phi,mp
607     o          ,iflag1
608     o          ,precip1,VPrecip1,sig1,w01
609     o          ,ft1,fq1,fu1,fv1,ftra1
610     o          ,inb1
611     o          ,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1
612     o          ,da1,phi1,mp1)
613      endif
614
615      if (iflag_con.eq.4) then
616       CALL cv_uncompress(nloc,len,ncum,nd,idcum
617     :          ,iflag
618     :          ,precip,cbmf
619     :          ,ft,fq,fu,fv
620     :          ,Ma,qcondc           
621     o          ,iflag1
622     o          ,precip1,cbmf1
623     o          ,ft1,fq1,fu1,fv1
624     o          ,Ma1,qcondc1 )           
625      endif
626
627      ENDIF ! ncum>0
628
6299999  continue
630
631      return
632      end
633
634!==================================================================
635      SUBROUTINE cv_flag
636      implicit none
637
638#include "cvflag.h"
639
640c -- si .TRUE., on rend la gravite plus explicite et eventuellement
641c differente de 10.0 dans convect3:
642      cvflag_grav = .TRUE.
643
644      return
645      end
646
647!==================================================================
648      SUBROUTINE cv_thermo(iflag_con)
649          implicit none
650
651c-------------------------------------------------------------
652c Set thermodynamical constants for convectL
653c-------------------------------------------------------------
654
655#include "YOMCST.h"
656#include "cvthermo.h"
657
658      integer iflag_con
659
660
661c original set from convect:
662      if (iflag_con.eq.4) then
663       cpd=1005.7
664       cpv=1870.0
665       cl=4190.0
666       rrv=461.5
667       rrd=287.04
668       lv0=2.501E6
669       g=9.8
670       t0=273.15
671       grav=g
672      endif
673
674c constants consistent with LMDZ:
675      if (iflag_con.eq.3) then
676       cpd = RCPD
677       cpv = RCPV
678       cl  = RCW
679       rrv = RV
680       rrd = RD
681       lv0 = RLVTT
682       g   = RG     ! not used in convect3
683c ori      t0  = RTT
684       t0  = 273.15 ! convect3 (RTT=273.16)
685c maf       grav= 10.    ! implicitely or explicitely used in convect3
686       grav= g    ! implicitely or explicitely used in convect3
687      endif
688
689      rowl=1000.0 !(a quelle variable de YOMCST cela correspond-il?)
690
691      clmcpv=cl-cpv
692      clmcpd=cl-cpd
693      cpdmcp=cpd-cpv
694      cpvmcpd=cpv-cpd
695      cpvmcl=cl-cpv ! for convect3
696      eps=rrd/rrv
697      epsi=1.0/eps
698      epsim1=epsi-1.0
699c      ginv=1.0/g
700      ginv=1.0/grav
701      hrd=0.5*rrd
702
703      return
704      end
705
Note: See TracBrowser for help on using the repository browser.