Changeset 619 for LMDZ4/trunk/libf/phylmd/cv_driver.F
- Timestamp:
- Apr 15, 2005, 2:36:17 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/cv_driver.F
r559 r619 5 5 & t1,q1,qs1,u1,v1,tra1, 6 6 & p1,ph1,iflag1,ft1,fq1,fu1,fv1,ftra1, 7 & precip1, 7 & precip1,VPrecip1, 8 8 & cbmf1,sig1,w01, 9 & delt,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1) 9 & icb1,inb1, 10 & delt,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1, 11 & da1,phi1,mp1) 10 12 C 11 13 implicit none … … 37 39 C ftra1 Real Output tracor tend 38 40 C precip1 Real Output precipitation 41 C VPrecip1 Real Output vertical profile of precipitations 39 42 C cbmf1 Real Output cloud base mass flux 40 43 C sig1 Real In/Out section adiabatic updraft … … 84 87 real precip1(len) 85 88 real cbmf1(len) 89 real VPrecip1(len,nd+1) 86 90 real Ma1(len,nd) 87 91 real upwd1(len,nd) … … 93 97 real cape1(len) 94 98 99 real da1(len,nd),phi1(len,nd,nd),mp1(len,nd) 100 real da(len,nd),phi(len,nd,nd),mp(len,nd) 95 101 real tra1(len,nd,ntra) 96 102 real ftra1(len,nd,ntra) … … 185 191 ! precip: Scalar convective precipitation rate (mm/day). 186 192 ! 193 ! VPrecip: Vertical profile of convective precipitation (kg/m2/s). 194 ! 187 195 ! wd: A convective downdraft velocity scale. For use in surface 188 196 ! flux parameterizations. See convect.ps file for details. … … 213 221 integer nk1(klon) 214 222 integer icb1(klon) 223 integer inb1(klon) 215 224 integer icbs1(klon) 216 225 … … 267 276 real ments(nloc,klev,klev), qents(nloc,klev,klev) 268 277 real sij(nloc,klev,klev), elij(nloc,klev,klev) 269 real mp(nloc,klev),qp(nloc,klev), up(nloc,klev), vp(nloc,klev)278 real qp(nloc,klev), up(nloc,klev), vp(nloc,klev) 270 279 real wt(nloc,klev), water(nloc,klev), evap(nloc,klev) 271 280 real b(nloc,klev), ft(nloc,klev), fq(nloc,klev) … … 275 284 real tps(nloc,klev), qprime(nloc), tprime(nloc) 276 285 real precip(nloc) 286 real VPrecip(nloc,klev+1) 277 287 real tra(nloc,klev,ntra), trap(nloc,klev,ntra) 278 288 real ftra(nloc,klev,ntra), traent(nloc,klev,klev,ntra) … … 324 334 clw(i,k)=0.0 325 335 gz1(i,k) = 0. 326 336 VPrecip1(i,k) = 0. 327 337 Ma1(i,k)=0.0 328 338 upwd1(i,k)=0.0 … … 346 356 wd1(i)=0.0 347 357 cape1(i)=0.0 358 VPrecip1(i,nd+1)=0.0 348 359 60 continue 349 360 … … 552 563 : ,ment,qent,uent,vent,nent,elij,traent,sig 553 564 : ,tv,tvp 554 o ,iflag,precip, ft,fq,fu,fv,ftra565 o ,iflag,precip,VPrecip,ft,fq,fu,fv,ftra 555 566 o ,upwd,dnwd,dnwd0,ma,mike,tls,tps,qcondc,wd) 556 567 endif … … 568 579 569 580 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 581 ! --- passive tracers 582 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 583 584 if (iflag_con.eq.3) then 585 CALL cv3_tracer(nloc,len,ncum,nd,nd, 586 : ment,sij,da,phi) 587 endif 588 589 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 570 590 ! --- UNCOMPRESS THE FIELDS 571 591 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 572 573 592 c set iflag1 =42 for non convective points 593 do i=1,len 594 iflag1(i)=42 595 end do 596 c 574 597 if (iflag_con.eq.3) then 575 598 CALL cv3_uncompress(nloc,len,ncum,nd,ntra,idcum 576 599 : ,iflag 577 : ,precip, sig,w0600 : ,precip,VPrecip,sig,w0 578 601 : ,ft,fq,fu,fv,ftra 602 : ,inb 579 603 : ,Ma,upwd,dnwd,dnwd0,qcondc,wd,cape 604 : ,da,phi,mp 580 605 o ,iflag1 581 o ,precip1, sig1,w01606 o ,precip1,VPrecip1,sig1,w01 582 607 o ,ft1,fq1,fu1,fv1,ftra1 583 o ,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1 ) 608 o ,inb1 609 o ,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1 610 o ,da1,phi1,mp1) 584 611 endif 585 612 … … 611 638 c -- si .TRUE., on rend la gravite plus explicite et eventuellement 612 639 c differente de 10.0 dans convect3: 613 cvflag_grav = . FALSE.640 cvflag_grav = .TRUE. 614 641 615 642 return … … 654 681 c ori t0 = RTT 655 682 t0 = 273.15 ! convect3 (RTT=273.16) 656 grav= 10. ! implicitely or explicitely used in convect3 683 c maf grav= 10. ! implicitely or explicitely used in convect3 684 grav= g ! implicitely or explicitely used in convect3 657 685 endif 658 686
Note: See TracChangeset
for help on using the changeset viewer.