| |
| | [No title] |
 | | do k = kts,klpbl do i = its,ite if(pblflg(i).and.k.ge.kpbl(i))then entfac(i,k) = ((zq(i,k)-hpbl(i))/delta(i))**2. |
 | | do k = kts,klpbl do i = its,ite if(k.lt.kpbl(i)) then zfac(i,k) = min(max((1.-(zq(i,k+1)-zl1(i)) & /(hpbl(i)-zl1(i))),zfmin),1.) xkzo = ckz*dza(i,k+1) zfacent(i,k) = (1.-zfac(i,k))**3. |
 | | do k = kts,kte-1 do i = its,ite xkzo = ckz*dza(i,k+1) if(k.ge.kpbl(i)) then ss = ((ux(i,k+1)-ux(i,k))*(ux(i,k+1)-ux(i,k))+(vx(i,k+1)- & vx(i,k))*(vx(i,k+1)-vx(i,k)))/(dza(i,k+1)*dza(i,k+1))+ & 1.e-9 govrthv = g/(0.5*(thvx(i,k+1)+thvx(i,k))) ri = govrthv*(thvx(i,k+1)-thvx(i,k))/(ss*dza(i,k+1)) if(imvdif.eq.1)then if((qcx(i,k)+qix(i,k)).gt.0.01e-3.and.(qcx(i,k+1)+ & qix(i,k+1)).gt.0.01e-3)then ! |
| www.mmm.ucar.edu /wrf/WG2/WRFV2/phys/module_bl_ysu.F (1992 words) |
|