Changeset 1795 for LMDZ5/branches/testing/libf/phylmd/cva_driver.F
- Timestamp:
- Jul 18, 2013, 10:20:28 AM (11 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
-
. (modified) (1 prop)
-
libf/phylmd/cva_driver.F (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk merged: 1747-1749,1751,1753-1767,1769,1771-1772,1774-1776,1778-1794
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/cva_driver.F
r1750 r1795 22 22 & ,lalim_conv, 23 23 & da1,phi1,mp1,phi21,d1a1,dam1,sigij1,clw1, ! RomP 24 & elij1,evap1,ep1, 24 & elij1,evap1,ep1,epmlmMm1,eplaMm1, ! RomP 25 25 & wdtrainA1,wdtrainM1) ! RomP 26 26 *************************************************************** … … 39 39 C.............................START PROLOGUE............................ 40 40 C 41 ! 42 ! All argument names (except len,nd,ntra,nloc,delt and the flags) have a "1" appended. 43 ! The "1" is removed for the corresponding compressed variables. 41 44 C PARAMETERS: 42 45 C Name Type Usage Description … … 99 102 C supmax01 Real Output 100 103 C asupmaxmin1 Real Output 104 ! 105 ! ftd1 Real Output Array of temperature tendency due to precipitations (K/s) of dimension ND, 106 ! defined at same grid levels as T, Q, QS and P. 107 ! 108 ! fqd1 Real Output Array of specific humidity tendencies due to precipitations ((gm/gm)/s) 109 ! of dimension ND, defined at same grid levels as T, Q, QS and P. 110 ! 111 ! wdtrainA1 Real Output precipitation detrained from adiabatic draught; 112 ! used in tracer transport (cvltr) 113 ! wdtrainM1 Real Output precipitation detrained from mixed draughts; 114 ! used in tracer transport (cvltr) 115 ! da1 Real Output used in tracer transport (cvltr) 116 ! phi1 Real Output used in tracer transport (cvltr) 117 ! mp1 Real Output used in tracer transport (cvltr) 118 ! 119 ! phi21 Real Output used in tracer transport (cvltr) 120 ! 121 ! d1a1 Real Output used in tracer transport (cvltr) 122 ! dam1 Real Output used in tracer transport (cvltr) 123 ! 124 ! epmlmMm1 Real Output used in tracer transport (cvltr) 125 ! eplaMm1 Real Output used in tracer transport (cvltr) 126 ! 127 ! evap1 Real Output 128 ! ep1 Real Output 129 ! sigij1 Real Output 130 ! elij1 Real Output 131 132 C 101 133 C S. Bony, Mar 2002: 102 134 C * Several modules corresponding to different physical processes … … 188 220 ! RomP >>> 189 221 real wdtrainA1(len,nd), wdtrainM1(len,nd) 190 real wdtrainA(nloc,klev),wdtrainM(nloc,klev)191 222 real da1(len,nd),phi1(len,nd,nd),mp1(len,nd) 192 real da(len,nd),phi(len,nd,nd)223 real epmlmMm1(len,nd,nd),eplaMm1(len,nd) 193 224 real evap1(len,nd),ep1(len,nd) 194 225 real sigij1(len,nd,nd),elij1(len,nd,nd) 195 real phi2(len,nd,nd)196 real d1a(len,nd), dam(len,nd)197 226 real phi21(len,nd,nd) 198 227 real d1a1(len,nd), dam1(len,nd) 199 228 ! RomP <<< 229 ! 230 !------------------------------------------------------------------- 231 ! Prolog by Kerry Emanuel. 200 232 !------------------------------------------------------------------- 201 233 ! --- ARGUMENTS … … 320 352 ! 321 353 ! det: Array of detrainment mass flux of dimension ND. 322 !323 ! ftd: Array of temperature tendency due to precipitations (K/s) of dimension ND,324 ! defined at same grid levels as T, Q, QS and P.325 !326 ! fqd: Array of specific humidity tendencies due to precipitations ((gm/gm)/s)327 ! of dimension ND, defined at same grid levels as T, Q, QS and P.328 !329 354 !------------------------------------------------------------------- 330 355 c … … 446 471 real wghti(nloc,nd) 447 472 real hnk(nloc),unk(nloc),vnk(nloc) 473 ! 474 ! RomP >>> 475 real wdtrainA(nloc,klev),wdtrainM(nloc,klev) 476 real da(len,nd),phi(len,nd,nd) 477 real epmlmMm(nloc,klev,klev),eplaMm(nloc,klev) 478 real phi2(len,nd,nd) 479 real d1a(len,nd), dam(len,nd) 480 ! RomP <<< 481 ! 448 482 logical, save :: first=.true. 449 483 c$OMP THREADPRIVATE(first) … … 498 532 nword4=len*nd*nd 499 533 500 ! call izilch(iflag1 ,nword1) 501 ! call zilch(iflag1 ,nword1) 502 do i=1,len 503 iflag1(i)=0 504 ktop1(i)=0 505 kbas1(i)=0 506 enddo 507 call zilch(ft1 ,nword2) 508 call zilch(fq1 ,nword2) 509 call zilch(fu1 ,nword2) 510 call zilch(fv1 ,nword2) 511 call zilch(ftra1 ,nword3) 512 call zilch(precip1 ,nword1) 513 ! call izilch(kbas1 ,nword1) 514 ! call zilch(kbas1 ,nword1) 515 ! call izilch(ktop1 ,nword1) 516 ! call zilch(ktop1 ,nword1) 517 call zilch(cbmf1 ,nword1) 518 call zilch(ptop21 ,nword1) 519 sigd1(:)=0. 520 call zilch(Ma1 ,nword2) 521 call zilch(mip1 ,nword2) 522 ! call zilch(Vprecip1,nword2) 523 Vprecip1=0. 524 call zilch(upwd1 ,nword2) 525 call zilch(dnwd1 ,nword2) 526 call zilch(dnwd01 ,nword2) 527 call zilch(qcondc1 ,nword2) 528 !test 529 ! call zilch(qcondc ,nword2) 530 call zilch(wd1 ,nword1) 531 call zilch(cape1 ,nword1) 532 call zilch(cin1 ,nword1) 533 call zilch(tvp1 ,nword2) 534 call zilch(ftd1 ,nword2) 535 call zilch(fqd1 ,nword2) 536 call zilch(Plim11 ,nword1) 537 call zilch(Plim21 ,nword1) 538 call zilch(asupmax1,nword2) 539 call zilch(supmax01,nword1) 540 call zilch(asupmaxmin1,nword1) 534 iflag1(:) = 0 535 ktop1(:) = 0 536 kbas1(:) = 0 537 ft1(:,:) = 0.0 538 fq1(:,:) = 0.0 539 fu1(:,:) = 0.0 540 fv1(:,:) = 0.0 541 ftra1(:,:,:) = 0. 542 precip1(:) = 0. 543 cbmf1(:) = 0. 544 ptop21(:) = 0. 545 sigd1(:) = 0. 546 Ma1(:,:) = 0. 547 mip1(:,:) = 0. 548 Vprecip1(:,:) = 0. 549 upwd1 (:,:) = 0. 550 dnwd1 (:,:) = 0. 551 dnwd01 (:,:) = 0. 552 qcondc1 (:,:) = 0. 553 wd1 (:) = 0. 554 cape1 (:) = 0. 555 cin1 (:) = 0. 556 tvp1 (:,:) = 0. 557 ftd1 (:,:) = 0. 558 fqd1 (:,:) = 0. 559 Plim11 (:) = 0. 560 Plim21 (:) = 0. 561 asupmax1(:,:) = 0. 562 supmax01(:) = 0. 563 asupmaxmin1(:)= 0. 541 564 c 542 565 DO il = 1,len … … 552 575 endif 553 576 577 ! RomP >>> 578 wdtrainA1(:,:) = 0. 579 wdtrainM1(:,:) = 0. 580 da1(:,:) = 0. 581 phi1(:,:,:) = 0. 582 epmlmMm1(:,:,:) = 0. 583 eplaMm1(:,:) = 0. 584 mp1(:,:) = 0. 585 evap1(:,:) = 0. 586 ep1(:,:) = 0. 587 sigij1(:,:,:) = 0. 588 elij1(:,:,:) = 0. 589 phi21(:,:,:) = 0. 590 d1a1(:,:) = 0. 591 dam1(:,:) = 0. 592 ! RomP <<< 554 593 !--------------------------------------------------------------------- 555 594 ! --- INITIALIZE LOCAL ARRAYS AND PARAMETERS … … 941 980 CALL cv3_tracer(nloc,len,ncum,nd,nd, 942 981 : ment,sigij,da,phi,phi2,d1a,dam, 943 : ep,Vprecip,elij,clw,icb,inb) 982 : ep,Vprecip,elij,clw,epmlmMm,eplaMm, 983 : icb,inb) 944 984 !RomP <<< 945 985 endif … … 964 1004 : ,asupmaxmin 965 1005 : ,da,phi,mp,phi2,d1a,dam,sigij ! RomP 966 : , wdtrainA,wdtrainM,elij,clw! RomP967 : , evap,ep! RomP1006 : ,clw,elij,evap,ep,epmlmMm,eplaMm ! RomP 1007 : ,wdtrainA,wdtrainM ! RomP 968 1008 o ,iflag1,kbas1,ktop1 969 1009 o ,precip1,cbmf1,plcl1,plfc1,wbeff1,sig1,w01,ptop21 … … 976 1016 o ,asupmaxmin1 977 1017 o ,da1,phi1,mp1,phi21,d1a1,dam1,sigij1 ! RomP 978 o , wdtrainA1,wdtrainM1,elij1,clw1! RomP979 o , evap1,ep1)! RomP1018 o ,clw1,elij1,evap1,ep1,epmlmMm1,eplaMm1! RomP 1019 o ,wdtrainA1,wdtrainM1) ! RomP 980 1020 endif 981 1021
Note: See TracChangeset
for help on using the changeset viewer.
