Symbol Description Uses
L Temps for string locations 2618



File Line number Procedure Code
./pkg/fizhi/do_fizhi.F 130 DO_FIZHI
      integer i,j,L
./pkg/fizhi/do_fizhi.F 202 DO_FIZHI
      do L = 1,Nrphys
./pkg/fizhi/do_fizhi.F 205 DO_FIZHI
       u(i,j,L) = uphy(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 206 DO_FIZHI
       v(i,j,L) = vphy(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 207 DO_FIZHI
       t(i,j,L) = thphy(i,j,L,bi,bj)/p0kappa
./pkg/fizhi/do_fizhi.F 208 DO_FIZHI
       q(i,j,L,1) = sphy(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 209 DO_FIZHI
       pl(i,j,L,bi,bj) = (pephy(i,j,L,bi,bj)+pephy(i,j,L+1,bi,bj))/200.
./pkg/fizhi/do_fizhi.F 210 DO_FIZHI
       dpres(i,j,L,bi,bj)=(pephy(i,j,L+1,bi,bj)-pephy(i,j,L,bi,bj))/100.
./pkg/fizhi/do_fizhi.F 211 DO_FIZHI
       ple(i,j,L,bi,bj) = pephy(i,j,L,bi,bj)/100.
./pkg/fizhi/do_fizhi.F 212 DO_FIZHI
       if (ple(i,j,L,bi,bj).gt.0.) then
./pkg/fizhi/do_fizhi.F 213 DO_FIZHI
        pkle(i,j,L,bi,bj) = ple(i,j,L,bi,bj) **kappa
./pkg/fizhi/do_fizhi.F 215 DO_FIZHI
        pkle(i,j,L,bi,bj) = 0.
./pkg/fizhi/do_fizhi.F 217 DO_FIZHI
       cldtot_lw(i,j,L) = cldtotlwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 218 DO_FIZHI
       cldlsp_lw(i,j,L) = cldlsplwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 219 DO_FIZHI
       cldras_lw(i,j,L) = cldraslwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 220 DO_FIZHI
       lwlz(i,j,L) = lwlzin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 221 DO_FIZHI
       qliqavelw(i,j,L) = qliqavelwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 222 DO_FIZHI
       fccavelw(i,j,L) = fccavelwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 223 DO_FIZHI
       cldtot_sw(i,j,L) = cldtotswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 224 DO_FIZHI
       cldlsp_sw(i,j,L) = cldlspswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 225 DO_FIZHI
       cldras_sw(i,j,L) = cldrasswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 226 DO_FIZHI
       swlz(i,j,L) = swlzin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 227 DO_FIZHI
       qliqavesw(i,j,L) = qliqaveswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 228 DO_FIZHI
       fccavesw(i,j,L) = fccaveswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 266 DO_FIZHI
      do L = 1,Nrphys
./pkg/fizhi/do_fizhi.F 269 DO_FIZHI
       cldtotlwin(i,j,L,bi,bj) = cldtot_lw(i,j,L)
./pkg/fizhi/do_fizhi.F 270 DO_FIZHI
       cldlsplwin(i,j,L,bi,bj) = cldlsp_lw(i,j,L)
./pkg/fizhi/do_fizhi.F 271 DO_FIZHI
       cldraslwin(i,j,L,bi,bj) = cldras_lw(i,j,L)
./pkg/fizhi/do_fizhi.F 272 DO_FIZHI
       lwlzin(i,j,L,bi,bj) = lwlz(i,j,L)
./pkg/fizhi/do_fizhi.F 273 DO_FIZHI
       qliqavelwin(i,j,L,bi,bj) = qliqavelw(i,j,L)
./pkg/fizhi/do_fizhi.F 274 DO_FIZHI
       fccavelwin(i,j,L,bi,bj) = fccavelw(i,j,L)
./pkg/fizhi/do_fizhi.F 275 DO_FIZHI
       cldtotswin(i,j,L,bi,bj) = cldtot_sw(i,j,L)
./pkg/fizhi/do_fizhi.F 276 DO_FIZHI
       cldlspswin(i,j,L,bi,bj) = cldlsp_sw(i,j,L)
./pkg/fizhi/do_fizhi.F 277 DO_FIZHI
       cldrasswin(i,j,L,bi,bj) = cldras_sw(i,j,L)
./pkg/fizhi/do_fizhi.F 278 DO_FIZHI
       swlzin(i,j,L,bi,bj) = swlz(i,j,L)
./pkg/fizhi/do_fizhi.F 279 DO_FIZHI
       qliqaveswin(i,j,L,bi,bj) = qliqavesw(i,j,L)
./pkg/fizhi/do_fizhi.F 280 DO_FIZHI
       fccaveswin(i,j,L,bi,bj) = fccavesw(i,j,L)
./pkg/fizhi/do_fizhi.F 295 DO_FIZHI
      do L = 1,Nrphys
./pkg/fizhi/do_fizhi.F 298 DO_FIZHI
       duphy(i,j,L,bi,bj) = moistu(i,j,L,bi,bj) + turbu(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 299 DO_FIZHI
       dvphy(i,j,L,bi,bj) = moistv(i,j,L,bi,bj) + turbv(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 300 DO_FIZHI
       dthphy(i,j,L,bi,bj) = ((moistt(i,j,L,bi,bj)+turbt(i,j,L,bi,bj)+
./pkg/fizhi/do_fizhi.F 301 DO_FIZHI
     .   lwdt(i,j,L,bi,bj) +
./pkg/fizhi/do_fizhi.F 302 DO_FIZHI
     .   dlwdtg(i,j,L,bi,bj) * (tgz(i,j,bi,bj)-tg0(i,j,bi,bj)) +
./pkg/fizhi/do_fizhi.F 303 DO_FIZHI
     .   swdt(i,j,L,bi,bj)*radswt(i,j,bi,bj) )*p0kappa ) / p(i,j,bi,bj)
./pkg/fizhi/do_fizhi.F 304 DO_FIZHI
       dsphy(i,j,L,bi,bj) = (moistq(i,j,L,1,bi,bj)+turbq(i,j,L,1,bi,bj))
./pkg/fizhi/fizhi_fillnegs.F 196 TRACER_FILL
      if( pq(i,j,L).lt.0.0 ) then
./pkg/fizhi/fizhi_fillnegs.F 198 TRACER_FILL
      qval     = pq(i  ,j,L)*d2p(i  ,j)*dp(i,j,L)
./pkg/fizhi/fizhi_fillnegs.F 199 TRACER_FILL
      array(1) = pq(ip1,j,L)*d2p(ip1,j)*dp(i,j,L)
./pkg/fizhi/fizhi_fillnegs.F 200 TRACER_FILL
      array(2) = pq(im1,j,L)*d2p(im1,j)*dp(i,j,L)
./pkg/fizhi/fizhi_fillnegs.F 205 TRACER_FILL
      array(3) = pq(i,j+1,L)*d2p(i,j+1)*dp(i,j,L)
./pkg/fizhi/fizhi_fillnegs.F 210 TRACER_FILL
      array(4) = pq(i,j-1,L)*d2p(i,j-1)*dp(i,j,L)
./pkg/fizhi/fizhi_fillnegs.F 212 TRACER_FILL
      if( L.eq.1    ) then
./pkg/fizhi/fizhi_fillnegs.F 215 TRACER_FILL
      array(5) = pq(i,j,L-1)*d2p(i,j)*dp(i,j,L)
./pkg/fizhi/fizhi_fillnegs.F 217 TRACER_FILL
      if( L.eq.lm   ) then
./pkg/fizhi/fizhi_fillnegs.F 220 TRACER_FILL
      array(6) = pq(i,j,L+1)*d2p(i,j)*dp(i,j,L)
./pkg/fizhi/fizhi_fillnegs.F 232 TRACER_FILL
               if( array(1).gt.0 ) pq(ip1,j,L) = pq(ip1,j,L) * fact
./pkg/fizhi/fizhi_fillnegs.F 233 TRACER_FILL
               if( array(2).gt.0 ) pq(im1,j,L) = pq(im1,j,L) * fact
./pkg/fizhi/fizhi_fillnegs.F 234 TRACER_FILL
               if( array(3).gt.0 ) pq(i,j+1,L) = pq(i,j+1,L) * fact
./pkg/fizhi/fizhi_fillnegs.F 235 TRACER_FILL
               if( array(4).gt.0 ) pq(i,j-1,L) = pq(i,j-1,L) * fact
./pkg/fizhi/fizhi_fillnegs.F 236 TRACER_FILL
               if( array(5).gt.0 ) pq(i,j,L-1) = pq(i,j,L-1) * fact
./pkg/fizhi/fizhi_fillnegs.F 237 TRACER_FILL
               if( array(6).gt.0 ) pq(i,j,L+1) = pq(i,j,L+1) * fact
./pkg/fizhi/fizhi_fillnegs.F 238 TRACER_FILL
                                   pq(i,j,L)   = 0.0
./pkg/fizhi/fizhi_fillnegs.F 240 TRACER_FILL
               pq(i,j,L) = 0.0
./pkg/fizhi/fizhi_fillnegs.F 243 TRACER_FILL
          if( imax.eq.1 ) pq(ip1,j,L) = pq(ip1,j,L) +
./pkg/fizhi/fizhi_fillnegs.F 244 TRACER_FILL
     .                                 pq(i,j,L)*d2p(i,j)/d2p(ip1,j)
./pkg/fizhi/fizhi_fillnegs.F 245 TRACER_FILL
          if( imax.eq.2 ) pq(im1,j,L) = pq(im1,j,L) +
./pkg/fizhi/fizhi_fillnegs.F 246 TRACER_FILL
     .                                 pq(i,j,L)*d2p(i,j)/d2p(im1,j)
./pkg/fizhi/fizhi_fillnegs.F 247 TRACER_FILL
          if( imax.eq.3 ) pq(i,j+1,L) = pq(i,j+1,L) +
./pkg/fizhi/fizhi_fillnegs.F 248 TRACER_FILL
     .                                 pq(i,j,L)*d2p(i,j)/d2p(i,j+1)
./pkg/fizhi/fizhi_fillnegs.F 249 TRACER_FILL
          if( imax.eq.4 ) pq(i,j-1,L) = pq(i,j-1,L) +
./pkg/fizhi/fizhi_fillnegs.F 250 TRACER_FILL
     .                                 pq(i,j,L)*d2p(i,j)/d2p(i,j-1)
./pkg/fizhi/fizhi_fillnegs.F 251 TRACER_FILL
          if( imax.eq.5 ) pq(i,j,L-1) = pq(i,j,L-1) +
./pkg/fizhi/fizhi_fillnegs.F 252 TRACER_FILL
     .                                 pq(i,j,L)*dp(i,j,L) /dp(i,j,L-1)
./pkg/fizhi/fizhi_fillnegs.F 253 TRACER_FILL
          if( imax.eq.6 ) pq(i,j,L+1) = pq(i,j,L+1) +
./pkg/fizhi/fizhi_fillnegs.F 254 TRACER_FILL
     .                                 pq(i,j,L)*dp(i,j,L) /dp(i,j,L+1)
./pkg/fizhi/fizhi_fillnegs.F 255 TRACER_FILL
                          pq(i,j,L)   = 0.0
./pkg/fizhi/fizhi_fillnegs.F 31 QCHECK
      integer i,j,L,LM1
./pkg/fizhi/fizhi_fillnegs.F 36 QCHECK
      DO L=2,Ldim
./pkg/fizhi/fizhi_fillnegs.F 37 QCHECK
      LM1 = L-1
./pkg/fizhi/fizhi_fillnegs.F 40 QCHECK
       ddsig = dp(i,j,LM1,bi,bj)/dp(i,j,L,bi,bj)
./pkg/fizhi/fizhi_fillnegs.F 42 QCHECK
        qz(i,j,L,bi,bj  ) = qz(i,j,L,bi,bj) + qz(i,j,LM1,bi,bj)*ddsig
./pkg/fizhi/fizhi_gwdrag.F 135 GWDRAG
      do L = 1,Lm
./pkg/fizhi/fizhi_gwdrag.F 138 GWDRAG
      tv(i,j,L) = tz(i,j,L)*pkz(i,j,L)*(1.+.609*qz(i,j,L))
./pkg/fizhi/fizhi_gwdrag.F 143 GWDRAG
      do L = 1,Lm
./pkg/fizhi/fizhi_gwdrag.F 146 GWDRAG
       dragu(i,j,L) = 0.
./pkg/fizhi/fizhi_gwdrag.F 147 GWDRAG
       dragv(i,j,L) = 0.
./pkg/fizhi/fizhi_gwdrag.F 148 GWDRAG
       dragt(i,j,L) = 0.
./pkg/fizhi/fizhi_gwdrag.F 184 GWDRAG
      do L = 1,Lm
./pkg/fizhi/fizhi_gwdrag.F 187 GWDRAG
         dragu(i,j,L) = sign( min(0.006 _d 0,abs(dragu(i,j,L))), dragu(i
./pkg/fizhi/fizhi_gwdrag.F 188 GWDRAG
     $        ,j,L) ) 
./pkg/fizhi/fizhi_gwdrag.F 189 GWDRAG
         dragv(i,j,L) = sign( min(0.006 _d 0,abs(dragv(i,j,L))), dragv(i
./pkg/fizhi/fizhi_gwdrag.F 190 GWDRAG
     $        ,j,L) ) 
./pkg/fizhi/fizhi_gwdrag.F 191 GWDRAG
      dragt(i,j,L) = -( uz(i,j,L)*dragu(i,j,L)+vz(i,j,L)*dragv(i,j,L) )
./pkg/fizhi/fizhi_gwdrag.F 193 GWDRAG
       dudt(i,j,L) = dudt(i,j,L) + dragu(i,j,L)
./pkg/fizhi/fizhi_gwdrag.F 194 GWDRAG
       dvdt(i,j,L) = dvdt(i,j,L) + dragv(i,j,L)
./pkg/fizhi/fizhi_gwdrag.F 195 GWDRAG
       dtdt(i,j,L) = dtdt(i,j,L) + dragt(i,j,L)*pz(i,j)/pkz(i,j,L)
./pkg/fizhi/fizhi_gwdrag.F 203 GWDRAG
      do L = 1,Lm
./pkg/fizhi/fizhi_gwdrag.F 208 GWDRAG
        tmpdiag(i,j) = dragu(i,j,L)*86400
./pkg/fizhi/fizhi_gwdrag.F 211 GWDRAG
       call DIAGNOSTICS_FILL(tmpdiag,'GWDU    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_gwdrag.F 217 GWDRAG
        tmpdiag(i,j) = dragv(i,j,L)*86400
./pkg/fizhi/fizhi_gwdrag.F 220 GWDRAG
       call DIAGNOSTICS_FILL(tmpdiag,'GWDV    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_gwdrag.F 226 GWDRAG
        tmpdiag(i,j) = dragt(i,j,L)*86400
./pkg/fizhi/fizhi_gwdrag.F 229 GWDRAG
       call DIAGNOSTICS_FILL(tmpdiag,'GWDT    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_gwdrag.F 254 GWDRAG
      do L = 1,Lm
./pkg/fizhi/fizhi_gwdrag.F 257 GWDRAG
      sumu(i,j) = sumu(i,j) + dragu(i,j,L)*dpres(i,j,L)/pz(i,j)
./pkg/fizhi/fizhi_gwdrag.F 277 GWDRAG
      do L = 1,Lm
./pkg/fizhi/fizhi_gwdrag.F 280 GWDRAG
      sumu(i,j) = sumu(i,j) + dragv(i,j,L)*dpres(i,j,L)/pz(i,j)
./pkg/fizhi/fizhi_gwdrag.F 367 GWDD
      integer  i,L
./pkg/fizhi/fizhi_gwdrag.F 385 GWDD
       ro(i,L) = pl(i,Lm+1-L)/(rgas*t(i,Lm+1-L))
./pkg/fizhi/fizhi_gwdrag.F 393 GWDD
       roiave  = ( 1./ro(i,L-1) + 1./ro(i,L) )*0.5
./pkg/fizhi/fizhi_gwdrag.F 394 GWDD
       dz(i,L) = agrav*roiave*( pl(i,Lm+2-L)-pl(i,Lm+1-L) )
./pkg/fizhi/fizhi_gwdrag.F 412 GWDD
      do  L = 1,nbase(i)-1
./pkg/fizhi/fizhi_gwdrag.F 413 GWDD
       robar(i) = robar(i) + ro(i,L) * (ple(i,Lm+2-L)-ple(i,Lm+1-L))
./pkg/fizhi/fizhi_gwdrag.F 414 GWDD
       ubar(i) = ubar(i) + u(i,Lm+1-L) * (ple(i,Lm+2-L)-ple(i,Lm+1-L))
./pkg/fizhi/fizhi_gwdrag.F 415 GWDD
       vbar(i) = vbar(i) + v(i,Lm+1-L) * (ple(i,Lm+2-L)-ple(i,Lm+1-L))
./pkg/fizhi/fizhi_gwdrag.F 432 GWDD
        vai1 = (t(i,Lm+1-L)-t(i,Lm+2-L))/dz(i,L)+gocp
./pkg/fizhi/fizhi_gwdrag.F 436 GWDD
        vai2    = 2.0*grav/( t(i,Lm+1-L)+t(i,Lm+2-L) )
./pkg/fizhi/fizhi_gwdrag.F 438 GWDD
        bv(i,L) = sqrt(vsqua)
./pkg/fizhi/fizhi_gwdrag.F 449 GWDD
       nbar(i) = nbar(i) + bv(i,L)*(pl(i,Lm+2-L)-pl(i,Lm+1-L))
./pkg/fizhi/fizhi_gwdrag.F 477 GWDD
         tensio(i,L) = tensio(i,1)
./pkg/fizhi/fizhi_gwdrag.F 478 GWDD
         xtens(i,L) = xtens(i,1)
./pkg/fizhi/fizhi_gwdrag.F 479 GWDD
         ytens(i,L) = ytens(i,1)
./pkg/fizhi/fizhi_gwdrag.F 491 GWDD
        velco = 0.5*( (u(i,Lm+1-L)*ubar(i) + v(i,Lm+1-L)*vbar(i))
./pkg/fizhi/fizhi_gwdrag.F 492 GWDD
     .            + (u(i,Lm+2-L)*ubar(i) + v(i,Lm+2-L)*vbar(i))  )
./pkg/fizhi/fizhi_gwdrag.F 496 GWDD
        roave = 0.5*(ro(i,L-1)+ro(i,L)) * 100.0     
./pkg/fizhi/fizhi_gwdrag.F 499 GWDD
         tensio(i,L) = tensio(i,L-1)
./pkg/fizhi/fizhi_gwdrag.F 505 GWDD
        fro2 = bv(i,L)/(akwnmb*roave*velco*velco*velco)*tensio(i,L-1)
./pkg/fizhi/fizhi_gwdrag.F 506 GWDD
        deluu = u(i,Lm+1-L)-u(i,Lm+2-L)
./pkg/fizhi/fizhi_gwdrag.F 507 GWDD
        delvv = v(i,Lm+1-L)-v(i,Lm+2-L)
./pkg/fizhi/fizhi_gwdrag.F 513 GWDD
         delz = dz(i,L)
./pkg/fizhi/fizhi_gwdrag.F 514 GWDD
         vsqua = bv(i,L)*bv(i,L)
./pkg/fizhi/fizhi_gwdrag.F 521 GWDD
         tensio(i,L) = tensio(i,L-1)
./pkg/fizhi/fizhi_gwdrag.F 533 GWDD
         tensio(i,L) = crif2/fro2*tensio(i,L-1)
./pkg/fizhi/fizhi_gwdrag.F 535 GWDD
         tensio(i,L) = tensio(i,L-1)
./pkg/fizhi/fizhi_gwdrag.F 539 GWDD
        xtens(i,L) = tensio(i,L)*cos(ang(i))
./pkg/fizhi/fizhi_gwdrag.F 540 GWDD
        ytens(i,L) = tensio(i,L)*sin(ang(i))
./pkg/fizhi/fizhi_gwdrag.F 556 GWDD
        tensio(i,L) = 0.0
./pkg/fizhi/fizhi_gwdrag.F 565 GWDD
         tensio(i,L) = tensio(i,L-1)
./pkg/fizhi/fizhi_gwdrag.F 569 GWDD
        roave = 0.5*(ro(i,L-1)+ro(i,L)) * 100.0
./pkg/fizhi/fizhi_gwdrag.F 570 GWDD
        vai1  = (t(i,Lm+1-L)-t(i,Lm+2-L))/dz(i,L)+gocp
./pkg/fizhi/fizhi_gwdrag.F 574 GWDD
         tensio(i,L) = 0.0
./pkg/fizhi/fizhi_gwdrag.F 578 GWDD
        vai2  = 2.0*grav/(t(i,Lm+1-L)+t(i,Lm+2-L))
./pkg/fizhi/fizhi_gwdrag.F 582 GWDD
        velco = 0.5*( (u(i,Lm+1-L)*ubar(i) + v(i,Lm+1-L)*vbar(i))
./pkg/fizhi/fizhi_gwdrag.F 583 GWDD
     .            + (u(i,Lm+2-L)*ubar(i) + v(i,Lm+2-L)*vbar(i))  )
./pkg/fizhi/fizhi_gwdrag.F 588 GWDD
         tensio(i,L) = 0.0
./pkg/fizhi/fizhi_gwdrag.F 594 GWDD
        fro2 = vaisd/(akwnmb*roave*velco*velco*velco)*tensio(i,L-1)
./pkg/fizhi/fizhi_gwdrag.F 595 GWDD
        deluu = u(i,Lm+1-L)-u(i,Lm+2-L)
./pkg/fizhi/fizhi_gwdrag.F 596 GWDD
        delvv = v(i,Lm+1-L)-v(i,Lm+2-L)
./pkg/fizhi/fizhi_gwdrag.F 602 GWDD
         delz   = dz(i,L)
./pkg/fizhi/fizhi_gwdrag.F 609 GWDD
         tensio(i,L) = 0.0
./pkg/fizhi/fizhi_gwdrag.F 621 GWDD
         tensio(i,L) = crif2/fro2*tensio(i,L-1)
./pkg/fizhi/fizhi_gwdrag.F 623 GWDD
         tensio(i,L) = tensio(i,L-1)
./pkg/fizhi/fizhi_gwdrag.F 627 GWDD
        xtens(i,L) = tensio(i,L)*cos(ang(i))
./pkg/fizhi/fizhi_gwdrag.F 628 GWDD
        ytens(i,L) = tensio(i,L)*sin(ang(i))
./pkg/fizhi/fizhi_gwdrag.F 638 GWDD
       coef = -grav*ps(i)/dpres(i,Lm+1-L)
./pkg/fizhi/fizhi_gwdrag.F 639 GWDD
       dudt(i,Lm+1-L) = coef*(xtens(i,L+1)-xtens(i,L))
./pkg/fizhi/fizhi_gwdrag.F 640 GWDD
       dvdt(i,Lm+1-L) = coef*(ytens(i,L+1)-ytens(i,L))
./pkg/fizhi/fizhi_gwdrag.F 657 GWDD
         dudt(i,Lm+1-L) = dudt(i,Lm)
./pkg/fizhi/fizhi_gwdrag.F 658 GWDD
         dvdt(i,Lm+1-L) = dvdt(i,Lm)
./pkg/fizhi/fizhi_gwdrag.F 667 GWDD
       dudt(i,L) = - dudt(i,L)/ps(i)*0.01
./pkg/fizhi/fizhi_gwdrag.F 668 GWDD
       dvdt(i,L) = - dvdt(i,L)/ps(i)*0.01
./pkg/fizhi/fizhi_gwdrag.F 78 GWDRAG
      integer n,i,j,L
./pkg/fizhi/fizhi_init_chem.F 203 GET_METHANE_N2O
      integer L,L1,L2,lup,ldn
./pkg/fizhi/fizhi_init_chem.F 230 GET_METHANE_N2O
      do L=1,46
./pkg/fizhi/fizhi_init_chem.F 231 GET_METHANE_N2O
      plevz(L) = plevc(47-L)
./pkg/fizhi/fizhi_init_chem.F 232 GET_METHANE_N2O
      hghtz(L) = hghtc(47-L)
./pkg/fizhi/fizhi_init_chem.F 271 GET_METHANE_N2O
      do L = 1,Nrphys
./pkg/fizhi/fizhi_init_chem.F 272 GET_METHANE_N2O
       if( hght(L).gt.26. ) then
./pkg/fizhi/fizhi_init_chem.F 273 GET_METHANE_N2O
        n2o(L) = 120.* exp( (26.- hght(L)) / 5.69 ) * 1.e-9
./pkg/fizhi/fizhi_init_chem.F 274 GET_METHANE_N2O
       else if( hght(L).gt.16. ) then
./pkg/fizhi/fizhi_init_chem.F 275 GET_METHANE_N2O
        n2o(L) = 307.* exp( (16.- hght(L)) /10.47 ) * 1.e-9
./pkg/fizhi/fizhi_init_chem.F 277 GET_METHANE_N2O
        n2o(L) = 307.e-9
./pkg/fizhi/fizhi_init_chem.F 281 GET_METHANE_N2O
      do L = 1,Nrphys
./pkg/fizhi/fizhi_init_chem.F 282 GET_METHANE_N2O
       if( hght(L).gt.55. ) then
./pkg/fizhi/fizhi_init_chem.F 283 GET_METHANE_N2O
        methane(L) = 0.2e-6
./pkg/fizhi/fizhi_init_chem.F 284 GET_METHANE_N2O
       else if( hght(L).gt.14. ) then
./pkg/fizhi/fizhi_init_chem.F 285 GET_METHANE_N2O
        methane(L) = 1.7* exp( (14.- hght(L)) /19.16 ) * 1.e-6
./pkg/fizhi/fizhi_init_chem.F 287 GET_METHANE_N2O
        methane(L) = 1.7e-6
./pkg/fizhi/fizhi_init_fixed.F 128 FIZHI_INIT_FIXED
       do L = 2,Nrphys+1
./pkg/fizhi/fizhi_init_fixed.F 129 FIZHI_INIT_FIXED
        pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100.
./pkg/fizhi/fizhi_init_fixed.F 133 FIZHI_INIT_FIXED
       do L = 1,Nrphys
./pkg/fizhi/fizhi_init_fixed.F 134 FIZHI_INIT_FIXED
        pressure(L)=(pressure0(Nrphys+2-L)+pressure0(Nrphys+1-L))/2.
./pkg/fizhi/fizhi_init_fixed.F 44 FIZHI_INIT_FIXED
       integer i,j,L,bi,bj
./pkg/fizhi/fizhi_init_vars.F 102 FIZHI_INIT_VARS
         do L = 1,Nr
./pkg/fizhi/fizhi_init_vars.F 103 FIZHI_INIT_VARS
          pedyn(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 117 FIZHI_INIT_VARS
         do L = Lbotij+1,Nr+1
./pkg/fizhi/fizhi_init_vars.F 118 FIZHI_INIT_VARS
          pedyn(i,j,L,bi,bj) = pedyn(i,j,L-1,bi,bj) -
./pkg/fizhi/fizhi_init_vars.F 119 FIZHI_INIT_VARS
     .                        drF(L-1)*hfacC(i,j,L-1,bi,bj)
./pkg/fizhi/fizhi_init_vars.F 130 FIZHI_INIT_VARS
         do L = 2,Nrphys+1
./pkg/fizhi/fizhi_init_vars.F 131 FIZHI_INIT_VARS
          pephy(i,j,L,bi,bj)=pephy(i,j,L-1,bi,bj)-dpphys0(i,j,L-1,bi,bj)
./pkg/fizhi/fizhi_init_vars.F 142 FIZHI_INIT_VARS
        do L = 1,Nrphys
./pkg/fizhi/fizhi_init_vars.F 145 FIZHI_INIT_VARS
         windphy(i,j,L,bi,bj) = 0.025 *
./pkg/fizhi/fizhi_init_vars.F 146 FIZHI_INIT_VARS
     .             log((pephy(i,j,1,bi,bj)-pephy(i,j,L+1,bi,bj))*10.)
./pkg/fizhi/fizhi_init_vars.F 166 FIZHI_INIT_VARS
        do L = 1,Nrphys
./pkg/fizhi/fizhi_init_vars.F 169 FIZHI_INIT_VARS
         uphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_init_vars.F 175 FIZHI_INIT_VARS
        do L = 1,Nrphys
./pkg/fizhi/fizhi_init_vars.F 178 FIZHI_INIT_VARS
         vphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_init_vars.F 184 FIZHI_INIT_VARS
        do L = 1,Nrphys
./pkg/fizhi/fizhi_init_vars.F 187 FIZHI_INIT_VARS
         thphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_init_vars.F 193 FIZHI_INIT_VARS
        do L = 1,Nrphys
./pkg/fizhi/fizhi_init_vars.F 196 FIZHI_INIT_VARS
         sphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_init_vars.F 202 FIZHI_INIT_VARS
        do L = 1,Nrphys
./pkg/fizhi/fizhi_init_vars.F 205 FIZHI_INIT_VARS
         duphy(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 206 FIZHI_INIT_VARS
         dvphy(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 207 FIZHI_INIT_VARS
         dthphy(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 208 FIZHI_INIT_VARS
         dsphy(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 214 FIZHI_INIT_VARS
        do L = 1,Nr
./pkg/fizhi/fizhi_init_vars.F 217 FIZHI_INIT_VARS
         guphy(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 218 FIZHI_INIT_VARS
         gvphy(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 219 FIZHI_INIT_VARS
         gthphy(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 220 FIZHI_INIT_VARS
         gsphy(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 234 FIZHI_INIT_VARS
         do L = 1,Nrphys
./pkg/fizhi/fizhi_init_vars.F 236 FIZHI_INIT_VARS
          tke(i,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 237 FIZHI_INIT_VARS
          xlmt(i,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 238 FIZHI_INIT_VARS
          khmt(i,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 259 FIZHI_INIT_VARS
        do L = 1,Nrphys
./pkg/fizhi/fizhi_init_vars.F 262 FIZHI_INIT_VARS
         swlz(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 263 FIZHI_INIT_VARS
         lwlz(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 264 FIZHI_INIT_VARS
         qliqavesw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 265 FIZHI_INIT_VARS
         qliqavelw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 266 FIZHI_INIT_VARS
         fccavesw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 267 FIZHI_INIT_VARS
         fccavelw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 268 FIZHI_INIT_VARS
         cldtot_sw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 269 FIZHI_INIT_VARS
         cldras_sw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 270 FIZHI_INIT_VARS
         cldlsp_sw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 271 FIZHI_INIT_VARS
         cldtot_lw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 272 FIZHI_INIT_VARS
         cldras_lw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 273 FIZHI_INIT_VARS
         cldlsp_lw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 60 FIZHI_INIT_VARS
       integer i, j, L, bi, bj, Lbotij
./pkg/fizhi/fizhi_lwrad.F 100 LWRIO
           cldtot(i,j,L)=min(1.0 _d 0,max(cldlw(i,j,L),fccave(i,j,L)/
./pkg/fizhi/fizhi_lwrad.F 102 LWRIO
           cldmxo(i,j,L) =  min( 1.0 _d 0,      clwmo(i,j,L) )
./pkg/fizhi/fizhi_lwrad.F 103 LWRIO
           lwlz(i,j,L) =  lwlz(i,j,L) + qliqave(i,j,L)/imstturb
./pkg/fizhi/fizhi_lwrad.F 108 LWRIO
        do L =1,lm
./pkg/fizhi/fizhi_lwrad.F 111 LWRIO
         cldtot(i,j,L) =  min( 1.0 _d 0,cldlw(i,j,L) )
./pkg/fizhi/fizhi_lwrad.F 112 LWRIO
         cldmxo(i,j,L) =  min( 1.0 _d 0,clwmo(i,j,L) )
./pkg/fizhi/fizhi_lwrad.F 144 LWRIO
       DO L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 146 LWRIO
        ADELPL(I,L) = convrt / ( ple(I,L+1)-ple(I,L) )
./pkg/fizhi/fizhi_lwrad.F 153 LWRIO
       DO L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 155 LWRIO
        CLRO(I,L) = min( 1.0 _d 0,clro(i,L) )
./pkg/fizhi/fizhi_lwrad.F 156 LWRIO
        CLMO(I,L) = min( 1.0 _d 0,clmo(i,L) )
./pkg/fizhi/fizhi_lwrad.F 163 LWRIO
      DO L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 165 LWRIO
      TZL(I,L) = TZL(I,L)*pk(I,L)
./pkg/fizhi/fizhi_lwrad.F 178 LWRIO
      do L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 180 LWRIO
      tauc(i,L,n) = tauc(i,L,n)*0.75
./pkg/fizhi/fizhi_lwrad.F 189 LWRIO
      do L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 191 LWRIO
      taual(i,L,n) = 0.
./pkg/fizhi/fizhi_lwrad.F 192 LWRIO
      ssaal(i,L,n) = 0.
./pkg/fizhi/fizhi_lwrad.F 193 LWRIO
      asyal(i,L,n) = 0.
./pkg/fizhi/fizhi_lwrad.F 201 LWRIO
      do L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 203 LWRIO
       cwc(i,L,n) = 0.
./pkg/fizhi/fizhi_lwrad.F 204 LWRIO
      reff(i,L,n) = 0.
./pkg/fizhi/fizhi_lwrad.F 223 LWRIO
       do L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 225 LWRIO
         dtrad(i,L) = (   flx(i,L)-   flx(i,L+1))*adelpl(i,L)
./pkg/fizhi/fizhi_lwrad.F 226 LWRIO
         tmpstrip(i,L) = flx(i,L)
./pkg/fizhi/fizhi_lwrad.F 227 LWRIO
         dtdtg(i,L) = ( dfdts(i,L)- dfdts(i,L+1))*adelpl(i,L)
./pkg/fizhi/fizhi_lwrad.F 228 LWRIO
        dtradc(i,L) = (flxclr(i,L)-flxclr(i,L+1))*adelpl(i,L)
./pkg/fizhi/fizhi_lwrad.F 262 LWRIO
      DO L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 264 LWRIO
      DTRAD (I,L) = ple(i,lm+1) * DTRAD (I,L)/pk(I,L)
./pkg/fizhi/fizhi_lwrad.F 265 LWRIO
      DTRADC(I,L) = ple(i,lm+1) * DTRADC(I,L)/pk(I,L)
./pkg/fizhi/fizhi_lwrad.F 266 LWRIO
       dtdtg(I,L) = ple(i,lm+1) * dtdtg (I,L)/pk(I,L)
./pkg/fizhi/fizhi_lwrad.F 285 LWRIO
      do L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 290 LWRIO
         tmpdiag(i,j) = tz(i,j,L)*pkz(i,j,L)
./pkg/fizhi/fizhi_lwrad.F 293 LWRIO
        call DIAGNOSTICS_FILL(tmpdiag,'TLW     ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_lwrad.F 299 LWRIO
         tmpdiag(i,j) = qz(i,j,L)*1000.
./pkg/fizhi/fizhi_lwrad.F 302 LWRIO
        call DIAGNOSTICS_FILL(tmpdiag,'SHRAD   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_lwrad.F 308 LWRIO
         tmpdiag(i,j) = oz(i,j,L)
./pkg/fizhi/fizhi_lwrad.F 311 LWRIO
        call DIAGNOSTICS_FILL(tmpdiag,'OZLW    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_lwrad.F 333 LWRIO
      do L = 1,lm
./pkg/fizhi/fizhi_lwrad.F 336 LWRIO
         fccave(i,j,L) = 0.
./pkg/fizhi/fizhi_lwrad.F 337 LWRIO
        qliqave(i,j,L) = 0.
./pkg/fizhi/fizhi_lwrad.F 97 LWRIO
        do L =1,lm
./pkg/fizhi/fizhi_moist.F 1039 MOISTIO
      do L = 1,lm
./pkg/fizhi/fizhi_moist.F 1046 MOISTIO
       plev = pl(indx,L)
./pkg/fizhi/fizhi_moist.F 1050 MOISTIO
       watnow = cldwater(i,j,L)
./pkg/fizhi/fizhi_moist.F 1052 MOISTIO
          cldras = min( max( cldras_lw(i,j,L)*cldras_mem,cpen(i,j,L)),
./pkg/fizhi/fizhi_moist.F 1057 MOISTIO
       cldlsp = min( max( cldlsp_lw(i,j,L)*cldlsp_mem,cldls(i,j,L)),
./pkg/fizhi/fizhi_moist.F 1065 MOISTIO
       lwlz(i,j,L) = ( nlwlz*lwlz(i,j,L) + watnow)/(nlwlz +1)
./pkg/fizhi/fizhi_moist.F 1066 MOISTIO
       cldtot_lw(i,j,L) = (nlwcld*cldtot_lw(i,j,L) + cldnow)/(nlwcld+1)
./pkg/fizhi/fizhi_moist.F 1067 MOISTIO
       cldlsp_lw(i,j,L) = (nlwcld*cldlsp_lw(i,j,L) + cldlsp)/(nlwcld+1)
./pkg/fizhi/fizhi_moist.F 1068 MOISTIO
       cldras_lw(i,j,L) = (nlwcld*cldras_lw(i,j,L) + cldras)/(nlwcld+1)
./pkg/fizhi/fizhi_moist.F 1073 MOISTIO
       watnow = cldwater(i,j,L)
./pkg/fizhi/fizhi_moist.F 1075 MOISTIO
          cldras = min( max(cldras_sw(i,j,L)*cldras_mem, cpen(i,j,L)),
./pkg/fizhi/fizhi_moist.F 1080 MOISTIO
       cldlsp = min( max(cldlsp_sw(i,j,L)*cldlsp_mem,cldls(i,j,L)),
./pkg/fizhi/fizhi_moist.F 1088 MOISTIO
       swlz(i,j,L) = ( nswlz*swlz(i,j,L) + watnow)/(nswlz +1)
./pkg/fizhi/fizhi_moist.F 1089 MOISTIO
       cldtot_sw(i,j,L) = (nswcld*cldtot_sw(i,j,L) + cldnow)/(nswcld+1)
./pkg/fizhi/fizhi_moist.F 1090 MOISTIO
       cldlsp_sw(i,j,L) = (nswcld*cldlsp_sw(i,j,L) + cldlsp)/(nswcld+1)
./pkg/fizhi/fizhi_moist.F 1091 MOISTIO
       cldras_sw(i,j,L) = (nswcld*cldras_sw(i,j,L) + cldras)/(nswcld+1)
./pkg/fizhi/fizhi_moist.F 1097 MOISTIO
       if( L.lt.midlevel ) cldhi (i,j) = max( cldnow,cldhi (i,j) )
./pkg/fizhi/fizhi_moist.F 1098 MOISTIO
       if( L.ge.midlevel  .and.
./pkg/fizhi/fizhi_moist.F 1099 MOISTIO
     .      L.lt.lowlevel ) cldmid(i,j) = max( cldnow,cldmid(i,j) )
./pkg/fizhi/fizhi_moist.F 1100 MOISTIO
       if( L.ge.lowlevel ) cldlow(i,j) = max( cldnow,cldlow(i,j) )
./pkg/fizhi/fizhi_moist.F 1104 MOISTIO
       cldtmp(i,j) =  cldtmp(i,j) + cldnow*pkzgather(i,L)
./pkg/fizhi/fizhi_moist.F 1105 MOISTIO
     .             *  ( tz(i,j,L) + dtmoist(i,j,L)*tmstp/pz(i,j) )
./pkg/fizhi/fizhi_moist.F 1231 RAS
      integer i,L,nc,ib,nt
./pkg/fizhi/fizhi_moist.F 1235 RAS
      do L = 1,k
./pkg/fizhi/fizhi_moist.F 1323 RAS
      do L=ib,k
./pkg/fizhi/fizhi_moist.F 1325 RAS
       cmass(i,L) = rhfrac(i)*cmass(i,L) * dt
./pkg/fizhi/fizhi_moist.F 1329 RAS
      do L=ib,k
./pkg/fizhi/fizhi_moist.F 1331 RAS
       cldmas(i,L) = cldmas(i,L) + cmass(i,L)
./pkg/fizhi/fizhi_moist.F 1339 RAS
      DO L=IB,K
./pkg/fizhi/fizhi_moist.F 1341 RAS
        thbef(I,L) = POI(I,L)
./pkg/fizhi/fizhi_moist.F 1342 RAS
        POI(I,L) = POI(I,L) + TCU(I,L) * DT * rhfrac(i)
./pkg/fizhi/fizhi_moist.F 1343 RAS
        QOI(I,L) = QOI(I,L) + QCU(I,L) * DT * rhfrac(i)
./pkg/fizhi/fizhi_moist.F 1347 RAS
      DO L=IB,K
./pkg/fizhi/fizhi_moist.F 1349 RAS
        UOI(I,L+nltop-1,NT)=UOI(I,L+nltop-1,NT)+UCU(I,L,NT)*DT*rhfrac(i)
./pkg/fizhi/fizhi_moist.F 1359 RAS
       do L = 1,k
./pkg/fizhi/fizhi_moist.F 1360 RAS
        if(tcu(i,L).ne.0.)ifound = ifound + 1
./pkg/fizhi/fizhi_moist.F 1364 RAS
        do L = 1,k
./pkg/fizhi/fizhi_moist.F 1365 RAS
         temp = TCU(I,L) * DT * rhfrac(i)
./pkg/fizhi/fizhi_moist.F 1377 RAS
      do L=k-1,1,-1
./pkg/fizhi/fizhi_moist.F 1379 RAS
         tx1(i)   = 100*(prs(i,L+1)-prs(i,L))/grav
./pkg/fizhi/fizhi_moist.F 1380 RAS
         cln(i,L) = min(1600*rains(i,L)/tx1(i),rasmax )
./pkg/fizhi/fizhi_moist.F 1661 CLOUD
      integer nt,km1,ic1,i,L,len1,len2,isav,len11,ii
./pkg/fizhi/fizhi_moist.F 1689 CLOUD
         DO 2050 L=1,K
./pkg/fizhi/fizhi_moist.F 1691 CLOUD
          PRH(I,L) = (PRJ(I,L+1)*PRS(I,L+1) - PRJ(I,L)*PRS(I,L))
./pkg/fizhi/fizhi_moist.F 1692 CLOUD
     *                            / ((PRS(I,L+1)-PRS(I,L)) * RKAPP1)
./pkg/fizhi/fizhi_moist.F 1696 CLOUD
         DO 2070 L=1,K
./pkg/fizhi/fizhi_moist.F 1698 CLOUD
          TX5(I) = POI(I,L) * PRH(I,L)
./pkg/fizhi/fizhi_moist.F 1699 CLOUD
          TX1(I) = (PRS(I,L) + PRS(I,L+1)) * 0.5
./pkg/fizhi/fizhi_moist.F 1702 CLOUD
          ALF(I,L) = TX2(I) - TX4(I) * TX5(I)
./pkg/fizhi/fizhi_moist.F 1703 CLOUD
          BET(I,L) = TX4(I) * PRH(I,L)
./pkg/fizhi/fizhi_moist.F 1704 CLOUD
          GAM(I,L) = 1.0 / ((1.0 + TX4(I)*ALBCP) * PRH(I,L))
./pkg/fizhi/fizhi_moist.F 1705 CLOUD
          PRI(I,L) = (CP/CMB2PA) / (PRS(I,L+1) - PRS(I,L))
./pkg/fizhi/fizhi_moist.F 1712 CLOUD
      DO 10 L=1,K
./pkg/fizhi/fizhi_moist.F 1714 CLOUD
      TCU(I,L) = 0.0
./pkg/fizhi/fizhi_moist.F 1715 CLOUD
      QCU(I,L) = 0.0
./pkg/fizhi/fizhi_moist.F 1716 CLOUD
      CMASS(I,L) = 0.0
./pkg/fizhi/fizhi_moist.F 1720 CLOUD
      do L=1,K
./pkg/fizhi/fizhi_moist.F 1722 CLOUD
      ucu(I,L,nt) = 0.0
./pkg/fizhi/fizhi_moist.F 1738 CLOUD
         DO 3703 L=KM1,IC1,-1
./pkg/fizhi/fizhi_moist.F 1740 CLOUD
         QS1(I)   = ALF(I,L) + BET(I,L)*POI(I,L)
./pkg/fizhi/fizhi_moist.F 1741 CLOUD
         QOL(I,L) = MIN(QS1(I)*RHMAX,QOI(I,L))
./pkg/fizhi/fizhi_moist.F 1743 CLOUD
         TEM1 = TX2(I) + PRJ(I,L+1) * POI(I,L) * CP
./pkg/fizhi/fizhi_moist.F 1744 CLOUD
         HOL(I,L) = TEM1 + QOL(I,L )* ALHL
./pkg/fizhi/fizhi_moist.F 1745 CLOUD
         HST(I,L) = TEM1 + QS1(I)   * ALHL
./pkg/fizhi/fizhi_moist.F 1747 CLOUD
         TX1(I)   = (PRJ(I,L+1) - PRJ(I,L)) * POI(I,L)
./pkg/fizhi/fizhi_moist.F 1748 CLOUD
         ETA(I,L) = ETA(I,L+1) + TX1(I)*CPBG
./pkg/fizhi/fizhi_moist.F 1777 CLOUD
      DO 160 L=IC,KM1
./pkg/fizhi/fizhi_moist.F 1779 CLOUD
      TX1(I) = TX1(I) + (HST(I,IC) - HOL(I,L)) * (ETA(I,L) - ETA(I,L+1))
./pkg/fizhi/fizhi_moist.F 1811 CLOUD
         DO L = 1,K
./pkg/fizhi/fizhi_moist.F 1813 CLOUD
           HST(I,L) = 0.0
./pkg/fizhi/fizhi_moist.F 1814 CLOUD
           QOL(I,L) = 0.0
./pkg/fizhi/fizhi_moist.F 1847 CLOUD
      DO 260 L=KM1,IC,-1
./pkg/fizhi/fizhi_moist.F 1849 CLOUD
      TX1(I) = ETA(IA(I),L)
./pkg/fizhi/fizhi_moist.F 1852 CLOUD
      ETA(I,L) = 1.0 + ALM(I) * TX1(I)
./pkg/fizhi/fizhi_moist.F 1872 CLOUD
         DO 380 L=KM1,IC1,-1
./pkg/fizhi/fizhi_moist.F 1875 CLOUD
         TEM = TX1(I) + (ETA(I,L) - ETA(I,L+1)) * HOL(II,L)
./pkg/fizhi/fizhi_moist.F 1877 CLOUD
         PCU(I) = PRJ(II,L+1) - PRH(II,L)
./pkg/fizhi/fizhi_moist.F 1878 CLOUD
         TEM1   = ETA(I,L+1) * PCU(I)
./pkg/fizhi/fizhi_moist.F 1881 CLOUD
         PCU(I) = PRH(II,L) - PRJ(II,L)
./pkg/fizhi/fizhi_moist.F 1882 CLOUD
         TEM1   = (TEM1 + ETA(I,L) * PCU(I)) * HST(II,L)
./pkg/fizhi/fizhi_moist.F 1885 CLOUD
         WFN(I) = WFN(I) + (TX1(I) - TEM1) * GAM(II,L)
./pkg/fizhi/fizhi_moist.F 1924 CLOUD
         DO L = 1,K
./pkg/fizhi/fizhi_moist.F 1926 CLOUD
           HST(I,L) = 0.0
./pkg/fizhi/fizhi_moist.F 1927 CLOUD
           QOL(I,L) = 0.0
./pkg/fizhi/fizhi_moist.F 1942 CLOUD
      DO 520 L=IC,K
./pkg/fizhi/fizhi_moist.F 1944 CLOUD
      TX1(I) = ETA(I2(I),L)
./pkg/fizhi/fizhi_moist.F 1947 CLOUD
      ETA(I,L) = TX1(I)
./pkg/fizhi/fizhi_moist.F 1981 CLOUD
      DO 635 L=KM1,IC,-1
./pkg/fizhi/fizhi_moist.F 1984 CLOUD
      TEM    = ETA(I,L) - ETA(I,L+1)
./pkg/fizhi/fizhi_moist.F 1985 CLOUD
      WLQ(I) = WLQ(I) + TEM * QOL(II,L)
./pkg/fizhi/fizhi_moist.F 1989 CLOUD
      DO L=KM1,IC,-1
./pkg/fizhi/fizhi_moist.F 1992 CLOUD
      TEM    = ETA(I,L) - ETA(I,L+1)
./pkg/fizhi/fizhi_moist.F 1993 CLOUD
      UHT(I,NT) = UHT(I,NT) + TEM *  UOI(II,L+nltop-1,NT)
./pkg/fizhi/fizhi_moist.F 2011 CLOUD
         DO 750 L=KM1,IC1,-1
./pkg/fizhi/fizhi_moist.F 2014 CLOUD
         TEM      = (POI(II,L-1) - POI(II,L)) * ETA(I,L)
./pkg/fizhi/fizhi_moist.F 2015 CLOUD
     *                            / (PRH(II,L) - PRH(II,L-1))
./pkg/fizhi/fizhi_moist.F 2017 CLOUD
         HOL(I,L)   = TEM * (PRJ(II,L)-PRH(II,L-1)) * PRH(II,L)
./pkg/fizhi/fizhi_moist.F 2018 CLOUD
     *                    *  PRI(II,L)  + HOL(I,L)
./pkg/fizhi/fizhi_moist.F 2019 CLOUD
         HOL(I,L-1) = TEM * (PRH(II,L)-PRJ(II,L)) * PRH(II,L-1)
./pkg/fizhi/fizhi_moist.F 2020 CLOUD
     *                                              * PRI(II,L-1)
./pkg/fizhi/fizhi_moist.F 2022 CLOUD
         AKM(I)   = AKM(I) - HOL(I,L)
./pkg/fizhi/fizhi_moist.F 2023 CLOUD
     *            * (ETA(I,L)   * (PRH(II,L)-PRJ(II,L)) +
./pkg/fizhi/fizhi_moist.F 2024 CLOUD
     *               ETA(I,L+1) * (PRJ(II,L+1)-PRH(II,L))) / PRH(II,L)
./pkg/fizhi/fizhi_moist.F 2067 CLOUD
         DO 840 L=KM1,IC1,-1
./pkg/fizhi/fizhi_moist.F 2071 CLOUD
         TX3(I) = (QOL(II,L-1) - QOL(II,L)) * ETA(I,L)
./pkg/fizhi/fizhi_moist.F 2074 CLOUD
         GMH(I,L) = HOL(I,L) + TX2(I)   * PRI(II,L) * (ALBCP*HALF)
./pkg/fizhi/fizhi_moist.F 2094 CLOUD
         DO 3725 L=KM1,IC1,-1
./pkg/fizhi/fizhi_moist.F 2097 CLOUD
         TX1(I) = TX1(I) + (ETA(I,L) - ETA(I,L+1)) * GMH(I,L)
./pkg/fizhi/fizhi_moist.F 2098 CLOUD
         TX2(I) = GAM(II,L-1) * (PRJ(II,L) - PRH(II,L-1))
./pkg/fizhi/fizhi_moist.F 2101 CLOUD
         IF (L .EQ. IC1) THEN
./pkg/fizhi/fizhi_moist.F 2109 CLOUD
     *          (TX2(I) + GAM(II,L)*(PRH(II,L)-PRJ(II,L)))
./pkg/fizhi/fizhi_moist.F 2171 CLOUD
      DO L = IC,K
./pkg/fizhi/fizhi_moist.F 2174 CLOUD
       if(L.lt.K)then
./pkg/fizhi/fizhi_moist.F 2175 CLOUD
       CMASS(II,L) =  ETA(I,L+1) * WFN(I) * ONEBG
./pkg/fizhi/fizhi_moist.F 2177 CLOUD
       CMASS(II,L) =  WFN(I) * ONEBG
./pkg/fizhi/fizhi_moist.F 2192 CLOUD
      DO 990 L=IC,K
./pkg/fizhi/fizhi_moist.F 2195 CLOUD
      TEM      = (GMH(I,L) - HOL(I,L)) * TX4(I)
./pkg/fizhi/fizhi_moist.F 2196 CLOUD
      TEM1     =  HOL(I,L) * TX5(I)
./pkg/fizhi/fizhi_moist.F 2198 CLOUD
      TCU(II,L) = TEM1 / PRH(II,L)
./pkg/fizhi/fizhi_moist.F 2199 CLOUD
      QCU(II,L) = TEM
./pkg/fizhi/fizhi_moist.F 2225 CLOUD
      DO 1020 L=KM1,IC1,-1
./pkg/fizhi/fizhi_moist.F 2228 CLOUD
      TEM = half*TX5(I) * PRI(II,L)
./pkg/fizhi/fizhi_moist.F 2230 CLOUD
      TX1(I) = (UOI(II,L-1+nltop-1,nt)-UOI(II,L+nltop-1,nt)) * ETA(I,L)
./pkg/fizhi/fizhi_moist.F 2235 CLOUD
      ucu(II,L,nt) = TX3(I)
./pkg/fizhi/fizhi_moist.F 223 MOISTIO
      DO L = lm-1,1,-1
./pkg/fizhi/fizhi_moist.F 232 MOISTIO
     1                   .and.(levpbl(i,j).eq.lm+1))levpbl(i,j)=L+1
./pkg/fizhi/fizhi_moist.F 2330 ACRITN
      integer i,L
./pkg/fizhi/fizhi_moist.F 2344 ACRITN
      do L=1,lma-1
./pkg/fizhi/fizhi_moist.F 2346 ACRITN
         if( pl(i).ge.p(L)   .and.
./pkg/fizhi/fizhi_moist.F 2347 ACRITN
     .       pl(i).le.p(L+1)) then
./pkg/fizhi/fizhi_moist.F 2348 ACRITN
         temp = ( pl(i)-p(L) )/( p(L+1)-p(L) )
./pkg/fizhi/fizhi_moist.F 2349 ACRITN
         acr(i) = a(L+1)*temp + a(L)*(1-temp)
./pkg/fizhi/fizhi_moist.F 2407 RNEVP
      integer i,L,n,nlaym1,irnlay,irnlm1
./pkg/fizhi/fizhi_moist.F 2438 RNEVP
       do L = 1,nlay
./pkg/fizhi/fizhi_moist.F 2440 RNEVP
          EVP9(i,L) = 0.
./pkg/fizhi/fizhi_moist.F 2441 RNEVP
         TEMP1(i,L) = 0.
./pkg/fizhi/fizhi_moist.F 2442 RNEVP
         TEMP2(i,L) = 0.
./pkg/fizhi/fizhi_moist.F 2443 RNEVP
         TEMP3(i,L) = 0.
./pkg/fizhi/fizhi_moist.F 2444 RNEVP
        CLSBTH(i,L) = 0.
./pkg/fizhi/fizhi_moist.F 2445 RNEVP
         cldlz(i,L) = 0.
./pkg/fizhi/fizhi_moist.F 2451 RNEVP
      DO L = 1,NLAY
./pkg/fizhi/fizhi_moist.F 2453 RNEVP
        TEMP2(I,L) = PL(I,L)*ZP001
./pkg/fizhi/fizhi_moist.F 2454 RNEVP
        TEMP2(I,L) = SQRT(TEMP2(I,L))
./pkg/fizhi/fizhi_moist.F 2460 RNEVP
      DO L = 1,NLAY
./pkg/fizhi/fizhi_moist.F 2462 RNEVP
      TEMP3(I,L) = GRAVITY*ZP01 / DP(I,L)
./pkg/fizhi/fizhi_moist.F 2468 RNEVP
      DO 100 L=1,NLAY
./pkg/fizhi/fizhi_moist.F 2471 RNEVP
        TEMP1(I,3) = TL(I,L)
./pkg/fizhi/fizhi_moist.F 2472 RNEVP
        TEMP1(I,4) = QL(I,L)
./pkg/fizhi/fizhi_moist.F 2480 RNEVP
       call QSAT ( temp1(i,3),pl(i,L),temp1(i,2),temp1(i,6),.true. )
./pkg/fizhi/fizhi_moist.F 248 MOISTIO
      do L = nsubmin,nltop,-1
./pkg/fizhi/fizhi_moist.F 2490 RNEVP
          EVP9(I,L) = (TEMP1(I,4) -  QL(I,L))/TEMP3(I,L)
./pkg/fizhi/fizhi_moist.F 2492 RNEVP
         cldlz(i,L) = rain(i,L)*temp3(i,L)
./pkg/fizhi/fizhi_moist.F 2493 RNEVP
       if(  tl(i,L).gt.tice-20.) then
./pkg/fizhi/fizhi_moist.F 2494 RNEVP
           water(i) =   water(i) + rain(i,L)
./pkg/fizhi/fizhi_moist.F 2496 RNEVP
         crystal(i) = crystal(i) + rain(i,L)
./pkg/fizhi/fizhi_moist.F 2510 RNEVP
       if( (evp9(i,L).gt.0.) .and. (crystal(i).gt.0.) ) then
./pkg/fizhi/fizhi_moist.F 2511 RNEVP
       iceevap(I) = EVP9(I,L)*fracice
./pkg/fizhi/fizhi_moist.F 2513 RNEVP
       EVP9(I,L)=EVP9(I,L)-iceevap(I)
./pkg/fizhi/fizhi_moist.F 2518 RNEVP
       if( (evp9(i,L).gt.0.) .and. (water(i).gt.0.) ) then
./pkg/fizhi/fizhi_moist.F 2519 RNEVP
       exparg    = ZM1P04*tmscl*((water(i)*RPHF*TEMP2(I,L))**ZP578)
./pkg/fizhi/fizhi_moist.F 251 MOISTIO
        if(levpbl(i,j).eq.L) then
./pkg/fizhi/fizhi_moist.F 2521 RNEVP
       watevap(I) = EVP9(I,L)*AREARAT*fracwat
./pkg/fizhi/fizhi_moist.F 2523 RNEVP
       EVP9(I,L)=EVP9(I,L)-watevap(I)
./pkg/fizhi/fizhi_moist.F 2527 RNEVP
       QL(I,L) = QL(I,L)+(iceevap(i)+watevap(i))*TEMP3(I,L)
./pkg/fizhi/fizhi_moist.F 2528 RNEVP
       TL(I,L) = TL(I,L)-(iceevap(i)+watevap(i))*TEMP3(I,L)*ELOCP
./pkg/fizhi/fizhi_moist.F 2542 RNEVP
      DO 200 L=1,NLAY
./pkg/fizhi/fizhi_moist.F 2544 RNEVP
        rainmax(i) = rhcrit(i,L)*evp9(i,L) +
./pkg/fizhi/fizhi_moist.F 2545 RNEVP
     .               ql(i,L)*(rhcrit(i,L)-1.)/temp3(i,L)
./pkg/fizhi/fizhi_moist.F 2548 RNEVP
           call QSAT( tl(i,L),pl(i,L),rh,dum,.false.)
./pkg/fizhi/fizhi_moist.F 2549 RNEVP
           rh = ql(i,L)/rh
./pkg/fizhi/fizhi_moist.F 2551 RNEVP
           if( rhcrit(i,L).eq.1.0 ) then
./pkg/fizhi/fizhi_moist.F 2554 RNEVP
           fact = min( 1.0 _d 0, alpha + (1.0-alpha)*( rh-rhcrit(i,L)) /
./pkg/fizhi/fizhi_moist.F 2555 RNEVP
     1                                          (1.0-rhcrit(i,L)) )
./pkg/fizhi/fizhi_moist.F 2559 RNEVP
           if( pl(i,L).ge.10.0 ) CLSBTH(I,L) = fact
./pkg/fizhi/fizhi_moist.F 2561 RNEVP
           QL(I,L) = QL(I,L)+rainmax(I)*TEMP3(I,L)
./pkg/fizhi/fizhi_moist.F 2562 RNEVP
           TL(I,L) = TL(I,L)-rainmax(I)*TEMP3(I,L)*ELOCP
./pkg/fizhi/fizhi_moist.F 2564 RNEVP
        cldlz(i,L) = cldlz(i,L) - rainmax(i)*temp3(i,L)
./pkg/fizhi/fizhi_moist.F 2570 RNEVP
          RPOW=(RLAR(I)*RPHF*TEMP2(I,L))**ZP578
./pkg/fizhi/fizhi_moist.F 2576 RNEVP
          QL(I,L) =    QL(I,L)+TEMP1(I,7)*TEMP3(I,L)
./pkg/fizhi/fizhi_moist.F 2577 RNEVP
          TL(I,L) =    TL(I,L)-TEMP1(I,7)*TEMP3(I,L)*ELOCP
./pkg/fizhi/fizhi_moist.F 2627 SRCLOUDS
      integer i,L
./pkg/fizhi/fizhi_moist.F 2645 SRCLOUDS
      do L = 1,irise
./pkg/fizhi/fizhi_moist.F 2647 SRCLOUDS
                  temp = th(i,L)*plk(i,L)
./pkg/fizhi/fizhi_moist.F 2648 SRCLOUDS
      call QSAT ( temp,pl(i,L),qs(i,L),dum,.false. )
./pkg/fizhi/fizhi_moist.F 2652 SRCLOUDS
      do L  = 2,irise
./pkg/fizhi/fizhi_moist.F 2654 SRCLOUDS
      rh    = q(i,L)/qs(i,L)
./pkg/fizhi/fizhi_moist.F 2656 SRCLOUDS
      rhcrit = rhc(i,L) - offset
./pkg/fizhi/fizhi_moist.F 2659 SRCLOUDS
      if(cloud(i,L).eq.  0.0 .and. ratio.gt.0.0 ) then
./pkg/fizhi/fizhi_moist.F 2660 SRCLOUDS
         cloud(i,L) = min( ratio,1.0 _d 0)
./pkg/fizhi/fizhi_moist.F 2678 CTEI
      integer i,L
./pkg/fizhi/fizhi_moist.F 2687 CTEI
      do L=lm,2,-1
./pkg/fizhi/fizhi_moist.F 2689 CTEI
          dq = q(i,L)+cldwat(i,L)-q(i,L-1)-cldwat(i,L-1)
./pkg/fizhi/fizhi_moist.F 2691 CTEI
      k = 1.0 + cpoel*plke(i,L)*( th(i,L)-th(i,L-1) ) / dq
./pkg/fizhi/fizhi_moist.F 2693 CTEI
      t = th(i,L)*plk(i,L)
./pkg/fizhi/fizhi_moist.F 2694 CTEI
      p = pl(i,L)
./pkg/fizhi/fizhi_moist.F 2708 CTEI
      cldfrc(i,L) = cldfrc(i,L)*f
./pkg/fizhi/fizhi_moist.F 2709 CTEI
      cldwat(i,L) = cldwat(i,L)*f
./pkg/fizhi/fizhi_moist.F 272 MOISTIO
      do L = 1,lm
./pkg/fizhi/fizhi_moist.F 283 MOISTIO
         thgather(indx,L) = tz(i,j,L)
./pkg/fizhi/fizhi_moist.F 284 MOISTIO
         shgather(indx,L) = qz(i,j,L,1)
./pkg/fizhi/fizhi_moist.F 285 MOISTIO
        pkegather(indx,L) = pkht(i,j,L)
./pkg/fizhi/fizhi_moist.F 286 MOISTIO
        pkzgather(indx,L) = pkl(i,j,L)
./pkg/fizhi/fizhi_moist.F 287 MOISTIO
        plegather(indx,L) = plze(i,j,L)
./pkg/fizhi/fizhi_moist.F 288 MOISTIO
        plzgather(indx,L) = plz(i,j,L)
./pkg/fizhi/fizhi_moist.F 289 MOISTIO
         dpgather(indx,L) = dpres(i,j,L)
./pkg/fizhi/fizhi_moist.F 295 MOISTIO
      do L = 1,lm
./pkg/fizhi/fizhi_moist.F 300 MOISTIO
        ugather(indx,L,nt) = qz(i,j,L,nt+ptracer)
./pkg/fizhi/fizhi_moist.F 308 MOISTIO
      do L = 1,lm
./pkg/fizhi/fizhi_moist.F 314 MOISTIO
        ugather(indx,L,ntracerin-ptracer+1) = uz(i,j,L)
./pkg/fizhi/fizhi_moist.F 315 MOISTIO
        ugather(indx,L,ntracerin-ptracer+2) = vz(i,j,L)
./pkg/fizhi/fizhi_moist.F 371 MOISTIO
      DO L = 1,lm
./pkg/fizhi/fizhi_moist.F 373 MOISTIO
             TH(I,L) = TH(I,L) * P0KAPPA
./pkg/fizhi/fizhi_moist.F 374 MOISTIO
         CLMAXO(I,L) = 0.
./pkg/fizhi/fizhi_moist.F 375 MOISTIO
         CLBOTH(I,L) = 0.
./pkg/fizhi/fizhi_moist.F 376 MOISTIO
         cldmas(I,L) = 0.
./pkg/fizhi/fizhi_moist.F 377 MOISTIO
        detrain(I,L) = 0.
./pkg/fizhi/fizhi_moist.F 381 MOISTIO
      do L = 1,lm
./pkg/fizhi/fizhi_moist.F 382 MOISTIO
      depths(L) = 0
./pkg/fizhi/fizhi_moist.F 386 MOISTIO
      do L = nsubmin,nltop,-1
./pkg/fizhi/fizhi_moist.F 387 MOISTIO
                       nindeces(L) = 0
./pkg/fizhi/fizhi_moist.F 389 MOISTIO
       if(pbl(i).eq.L) nindeces(L) = nindeces(L) + 1
./pkg/fizhi/fizhi_moist.F 391 MOISTIO
       if(nindeces(L).gt.0) then
./pkg/fizhi/fizhi_moist.F 393 MOISTIO
       depths(numdeps) = L
./pkg/fizhi/fizhi_moist.F 415 MOISTIO
       DO 601 L=NSUBCL,lm
./pkg/fizhi/fizhi_moist.F 417 MOISTIO
        TMP1(I,2) = TMP1(I,2) + (PLE(I,L+1)-PLE(I,L))*TH (I,L)/sp(i)
./pkg/fizhi/fizhi_moist.F 418 MOISTIO
        TMP1(I,3) = TMP1(I,3) + (PLE(I,L+1)-PLE(I,L))*SHL(I,L)/sp(i)
./pkg/fizhi/fizhi_moist.F 430 MOISTIO
          do  L = 1,lm
./pkg/fizhi/fizhi_moist.F 432 MOISTIO
          saveu(i,L,nt) = ul(i,L,nt)
./pkg/fizhi/fizhi_moist.F 438 MOISTIO
          DO L=NSUBCL,lm
./pkg/fizhi/fizhi_moist.F 440 MOISTIO
           TMP1(I,2) = TMP1(I,2)+(PLE(I,L+1)-PLE(I,L))*UL(I,L,NT)/sp(i)
./pkg/fizhi/fizhi_moist.F 451 MOISTIO
       DO 111 L=1,lm
./pkg/fizhi/fizhi_moist.F 453 MOISTIO
        TMP4(I,L) = PLE(I,L)
./pkg/fizhi/fizhi_moist.F 459 MOISTIO
       DO L=2,lm
./pkg/fizhi/fizhi_moist.F 461 MOISTIO
        TMP5(I,L) = PLKE(I,L)*P0KINV
./pkg/fizhi/fizhi_moist.F 503 MOISTIO
       do L = nsubcl, lm
./pkg/fizhi/fizhi_moist.F 504 MOISTIO
         rhcrit(i,L) = 1.
./pkg/fizhi/fizhi_moist.F 506 MOISTIO
       do L = 1, nsubcl-1
./pkg/fizhi/fizhi_moist.F 507 MOISTIO
        pcheck = pl(i,L)
./pkg/fizhi/fizhi_moist.F 509 MOISTIO
         rhcrit(i,L) = rhmin
./pkg/fizhi/fizhi_moist.F 512 MOISTIO
         rhcrit(i,L) = rhmin + (1.-rhmin)/(19.) *
./pkg/fizhi/fizhi_moist.F 522 MOISTIO
      do L = 1,lm
./pkg/fizhi/fizhi_moist.F 524 MOISTIO
        saveth(i,L) = th (i,L)
./pkg/fizhi/fizhi_moist.F 525 MOISTIO
        saveq (i,L) = shl(i,L)
./pkg/fizhi/fizhi_moist.F 526 MOISTIO
        PCPEN (i,L) = 0.
./pkg/fizhi/fizhi_moist.F 527 MOISTIO
        CLFRAC(i,L) = 0.
./pkg/fizhi/fizhi_moist.F 540 MOISTIO
       do L=nsubcl,lm
./pkg/fizhi/fizhi_moist.F 542 MOISTIO
        dum = dp(i,L)/(ple(i,lm+1)-ple(i,nsubcl))
./pkg/fizhi/fizhi_moist.F 543 MOISTIO
        cldmas(i,L) = cldmas(i,L-1) - dum*cldmas(i,nsubcl-1)
./pkg/fizhi/fizhi_moist.F 549 MOISTIO
       DO L=1,nsubcl
./pkg/fizhi/fizhi_moist.F 551 MOISTIO
        CVTH(I,L) =  (TH (I,L) - saveth(i,l))
./pkg/fizhi/fizhi_moist.F 552 MOISTIO
        CVQ (I,L) =  (SHL(I,L) - saveq (i,l))
./pkg/fizhi/fizhi_moist.F 555 MOISTIO
       DO L=nsubcl+1,lm
./pkg/fizhi/fizhi_moist.F 557 MOISTIO
        CVTH(I,L) = cvth(i,nsubcl)
./pkg/fizhi/fizhi_moist.F 558 MOISTIO
        CVQ (I,L) = cvq (i,nsubcl)
./pkg/fizhi/fizhi_moist.F 562 MOISTIO
       DO L=nsubcl+1,lm
./pkg/fizhi/fizhi_moist.F 564 MOISTIO
        TH (I,L) = saveth(i,l) + cvth(i,l)
./pkg/fizhi/fizhi_moist.F 565 MOISTIO
        SHL(I,L) = saveq (i,l) + cvq (i,l)
./pkg/fizhi/fizhi_moist.F 568 MOISTIO
       DO L=1,lm
./pkg/fizhi/fizhi_moist.F 570 MOISTIO
        CVTH(I,L) =  CVTH(I,L) *P0KINV*SP(I)*tminv
./pkg/fizhi/fizhi_moist.F 571 MOISTIO
        CVQ (I,L) =  CVQ (I,L) *SP(I)*tminv
./pkg/fizhi/fizhi_moist.F 578 MOISTIO
        DO L=1,nsubcl-1
./pkg/fizhi/fizhi_moist.F 580 MOISTIO
         CVU(I,L,nt) = ( UL(I,L,nt)-saveu(i,l,nt) )*sp(i)*tminv
./pkg/fizhi/fizhi_moist.F 583 MOISTIO
        DO L=nsubcl,lm
./pkg/fizhi/fizhi_moist.F 586 MOISTIO
          cvu(i,L,nt) = ( ul(i,nsubcl,nt)-usubcl(i,nt) ) *
./pkg/fizhi/fizhi_moist.F 587 MOISTIO
     .                     ( saveu(i,L,nt)/usubcl(i,nt) )*sp(i)*tminv
./pkg/fizhi/fizhi_moist.F 589 MOISTIO
          cvu(i,L,nt) = 0.0
./pkg/fizhi/fizhi_moist.F 599 MOISTIO
       do L=nltop,nsubcl
./pkg/fizhi/fizhi_moist.F 600 MOISTIO
       if( cvq(i,L).ne.0.0 ) lras = .true.
./pkg/fizhi/fizhi_moist.F 644 MOISTIO
      DO L=1,lm
./pkg/fizhi/fizhi_moist.F 646 MOISTIO
           TH(I,L) =  TH(I,L) + CVTH(I,L)*tmstp/SP(I)
./pkg/fizhi/fizhi_moist.F 647 MOISTIO
          SHL(I,L) = SHL(I,L) +  CVQ(I,L)*tmstp/SP(I)
./pkg/fizhi/fizhi_moist.F 648 MOISTIO
           TL(I,L) =  TH(I,L)*PLK(I,L)
./pkg/fizhi/fizhi_moist.F 649 MOISTIO
       saveth(I,L) =  th(I,L)
./pkg/fizhi/fizhi_moist.F 650 MOISTIO
       saveq (I,L) = SHL(I,L)
./pkg/fizhi/fizhi_moist.F 662 MOISTIO
      DO L=1,lm
./pkg/fizhi/fizhi_moist.F 665 MOISTIO
       TMP1(I,L) = sp(i) * (SHL(I,L)-saveq(I,L)) * tminv
./pkg/fizhi/fizhi_moist.F 667 MOISTIO
       CALL PSTBMP(TMP1(1,L),delqgather(1,L),ISTRIP,im*jm,1,NN)
./pkg/fizhi/fizhi_moist.F 670 MOISTIO
       TMP1(I,L) = sp(i) * ((TL(I,L)/PLK(I,L))-saveth(i,l)) * tminv
./pkg/fizhi/fizhi_moist.F 672 MOISTIO
       CALL PSTBMP(TMP1(1,L),deltgather(1,L),ISTRIP,im*jm,1,NN)
./pkg/fizhi/fizhi_moist.F 677 MOISTIO
        TMP1(I,L) = ((TL(I,L)/PLK(I,L))-saveth(i,l))*plk(i,l)*sday*tminv
./pkg/fizhi/fizhi_moist.F 679 MOISTIO
       call PASTE(tmp1(1,L),deltrnev(1,L),istrip,im*jm,1,NN)
./pkg/fizhi/fizhi_moist.F 682 MOISTIO
        TMP1(I,L) = (SHL(I,L)-saveq(I,L)) * 1000. * sday * tminv
./pkg/fizhi/fizhi_moist.F 684 MOISTIO
       call PASTE(tmp1(1,L),delqrnev(1,L),istrip,im*jm,1,NN)
./pkg/fizhi/fizhi_moist.F 694 MOISTIO
      do L=1,lm
./pkg/fizhi/fizhi_moist.F 696 MOISTIO
      srcld(i,L) = -clsbth(i,L)
./pkg/fizhi/fizhi_moist.F 703 MOISTIO
      do L=1,lm
./pkg/fizhi/fizhi_moist.F 705 MOISTIO
      srcld(i,L) = srcld(i,L)+clsbth(i,L)
./pkg/fizhi/fizhi_moist.F 778 MOISTIO
      do L = 1,lm
./pkg/fizhi/fizhi_moist.F 781 MOISTIO
       call BACK2GRD (deltgather(1,L),pblindex, dtmoist(1,1,L)  ,im*jm)
./pkg/fizhi/fizhi_moist.F 782 MOISTIO
       call BACK2GRD (delqgather(1,L),pblindex, dqmoist(1,1,L,1),im*jm)
./pkg/fizhi/fizhi_moist.F 783 MOISTIO
       call BACK2GRD (    cpen(1,1,L),pblindex,    cpen(1,1,L)  ,im*jm)
./pkg/fizhi/fizhi_moist.F 784 MOISTIO
       call BACK2GRD (   cldls(1,1,L),pblindex,   cldls(1,1,L)  ,im*jm)
./pkg/fizhi/fizhi_moist.F 785 MOISTIO
       call BACK2GRD (cldwater(1,1,L),pblindex,cldwater(1,1,L)  ,im*jm)
./pkg/fizhi/fizhi_moist.F 786 MOISTIO
       call BACK2GRD ( pkzgather(1,L),pblindex, pkzgather(1,L)  ,im*jm)
./pkg/fizhi/fizhi_moist.F 790 MOISTIO
       call BACK2GRD(tmpgather(1,L),pblindex,
./pkg/fizhi/fizhi_moist.F 791 MOISTIO
     .                                            tmpgather(1,L),im*jm)
./pkg/fizhi/fizhi_moist.F 792 MOISTIO
       call BACK2GRD(pkegather(1,L),pblindex,
./pkg/fizhi/fizhi_moist.F 793 MOISTIO
     .                                            pkegather(1,L),im*jm)
./pkg/fizhi/fizhi_moist.F 794 MOISTIO
       call BACK2GRD(deltrnev(1,L),pblindex,
./pkg/fizhi/fizhi_moist.F 795 MOISTIO
     .                                             deltrnev(1,L),im*jm)
./pkg/fizhi/fizhi_moist.F 796 MOISTIO
       call BACK2GRD(delqrnev(1,L),pblindex,
./pkg/fizhi/fizhi_moist.F 797 MOISTIO
     .                                             delqrnev(1,L),im*jm)
./pkg/fizhi/fizhi_moist.F 798 MOISTIO
       call BACK2GRD(cldsr(1,1,L),pblindex,
./pkg/fizhi/fizhi_moist.F 799 MOISTIO
     .                                              cldsr(1,1,L),im*jm)
./pkg/fizhi/fizhi_moist.F 805 MOISTIO
       do L = 1,lm
./pkg/fizhi/fizhi_moist.F 806 MOISTIO
       call BACK2GRD (delugather(1,L,nt),pblindex,
./pkg/fizhi/fizhi_moist.F 807 MOISTIO
     .                                 dqmoist(1,1,L,ptracer+nt),im*jm)
./pkg/fizhi/fizhi_moist.F 815 MOISTIO
      do L = 1,lm
./pkg/fizhi/fizhi_moist.F 816 MOISTIO
       call BACK2GRD (delugather(1,L,ntracerin-ptracer+1),pblindex,
./pkg/fizhi/fizhi_moist.F 817 MOISTIO
     .                                 dumoist(1,1,L),im*jm)
./pkg/fizhi/fizhi_moist.F 818 MOISTIO
       call BACK2GRD (delugather(1,L,ntracerin-ptracer+2),pblindex,
./pkg/fizhi/fizhi_moist.F 819 MOISTIO
     .                                 dvmoist(1,1,L),im*jm)
./pkg/fizhi/fizhi_moist.F 828 MOISTIO
      do L = 1,lm
./pkg/fizhi/fizhi_moist.F 831 MOISTIO
       dumoist(i,j,L) = dumoist(i,j,L) * tmpimjm(i,j)
./pkg/fizhi/fizhi_moist.F 832 MOISTIO
       dvmoist(i,j,L) = dvmoist(i,j,L) * tmpimjm(i,j)
./pkg/fizhi/fizhi_moist.F 857 MOISTIO
       do L=1,lm
./pkg/fizhi/fizhi_moist.F 860 MOISTIO
        tmpdiag(i,j) = cldsr(i,j,L)
./pkg/fizhi/fizhi_moist.F 863 MOISTIO
       call DIAGNOSTICS_FILL(tmpdiag,'CLDNP   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_moist.F 870 MOISTIO
       do L=1,lm
./pkg/fizhi/fizhi_moist.F 874 MOISTIO
        tmpdiag(i,j)=(dtmoist(i,j,L)*sday*pkzgather(indgath,L)/pz(i,j))
./pkg/fizhi/fizhi_moist.F 877 MOISTIO
       call DIAGNOSTICS_FILL(tmpdiag,'MOISTT  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_moist.F 884 MOISTIO
       do L=1,lm
./pkg/fizhi/fizhi_moist.F 887 MOISTIO
        tmpdiag(i,j)=(dqmoist(i,j,L,1)*sday*1000./pz(i,j))
./pkg/fizhi/fizhi_moist.F 890 MOISTIO
       call DIAGNOSTICS_FILL(tmpdiag,'MOISTQ  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_moist.F 897 MOISTIO
       do L=1,lm
./pkg/fizhi/fizhi_moist.F 900 MOISTIO
        tmpdiag(i,j)=dumoist(i,j,L)*sday
./pkg/fizhi/fizhi_moist.F 903 MOISTIO
       call DIAGNOSTICS_FILL(tmpdiag,'MOISTU  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_moist.F 910 MOISTIO
       do L=1,lm
./pkg/fizhi/fizhi_moist.F 913 MOISTIO
        tmpdiag(i,j)=dvmoist(i,j,L)*sday
./pkg/fizhi/fizhi_moist.F 916 MOISTIO
       call DIAGNOSTICS_FILL(tmpdiag,'MOISTV  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_moist.F 923 MOISTIO
       do L=1,lm
./pkg/fizhi/fizhi_moist.F 927 MOISTIO
        tmpdiag(i,j)=tmpgather(indgath,L)
./pkg/fizhi/fizhi_moist.F 930 MOISTIO
       call DIAGNOSTICS_FILL(tmpdiag,'CLDMAS  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_moist.F 937 MOISTIO
       do L=1,lm
./pkg/fizhi/fizhi_moist.F 941 MOISTIO
        tmpdiag(i,j)=pkegather(indgath,L)
./pkg/fizhi/fizhi_moist.F 944 MOISTIO
       call DIAGNOSTICS_FILL(tmpdiag,'DTRAIN  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_moist.F 951 MOISTIO
       do L=1,lm
./pkg/fizhi/fizhi_moist.F 955 MOISTIO
        tmpdiag(i,j)=deltrnev(indgath,L)
./pkg/fizhi/fizhi_moist.F 958 MOISTIO
       call DIAGNOSTICS_FILL(tmpdiag,'DTLS    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_moist.F 965 MOISTIO
       do L=1,lm
./pkg/fizhi/fizhi_moist.F 969 MOISTIO
        tmpdiag(i,j)=delqrnev(indgath,L)
./pkg/fizhi/fizhi_moist.F 972 MOISTIO
       call DIAGNOSTICS_FILL(tmpdiag,'DQLS    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_rayleigh.F 100 RAYLEIGH
      do L=1,lm
./pkg/fizhi/fizhi_rayleigh.F 105 RAYLEIGH
        tmpdiag(i,j) = rfu(i,j,L)*86400
./pkg/fizhi/fizhi_rayleigh.F 114 RAYLEIGH
        tmpdiag(i,j) = rfv(i,j,L)*86400
./pkg/fizhi/fizhi_rayleigh.F 123 RAYLEIGH
        tmpdiag(i,j) = rft(i,j,L)*86400
./pkg/fizhi/fizhi_rayleigh.F 44 RAYLEIGH
      integer  i,j,L
./pkg/fizhi/fizhi_rayleigh.F 64 RAYLEIGH
      do L=1,lm
./pkg/fizhi/fizhi_rayleigh.F 67 RAYLEIGH
        dz(i,j,L) = cpog * (pekap(i,j,L+1)-pekap(i,j,L)) * t(i,j,L) *
./pkg/fizhi/fizhi_rayleigh.F 68 RAYLEIGH
     .                (1.+virtcon*s(i,j,L))
./pkg/fizhi/fizhi_rayleigh.F 79 RAYLEIGH
      do L=lm-1,1,-1
./pkg/fizhi/fizhi_rayleigh.F 82 RAYLEIGH
        z(i,j,L) = z(i,j,L+1) + 0.5 * (dz(i,j,L)+dz(i,j,L+1))
./pkg/fizhi/fizhi_rayleigh.F 87 RAYLEIGH
      do L=1,lm
./pkg/fizhi/fizhi_rayleigh.F 90 RAYLEIGH
        rf(i,j,L) = dampcoef*(1+tanh((z(i,j,L)-50000.)/5000.))/86400.
./pkg/fizhi/fizhi_rayleigh.F 91 RAYLEIGH
        rfu(i,j,L) = - rf(i,j,L) * u(i,j,L)
./pkg/fizhi/fizhi_rayleigh.F 92 RAYLEIGH
        rfv(i,j,L) = - rf(i,j,L) * v(i,j,L)
./pkg/fizhi/fizhi_rayleigh.F 93 RAYLEIGH
        rft(i,j,L) = -(u(i,j,L)*rfu(i,j,L) + v(i,j,L)*rfv(i,j,L) )*cpinv
./pkg/fizhi/fizhi_rayleigh.F 94 RAYLEIGH
     .                        /pkap(i,j,L)
./pkg/fizhi/fizhi_step_diag.F 166 FIZHI_STEP_DIAG
      do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 173 FIZHI_STEP_DIAG
        tmpdiag(i,j) = (moistu (i,j,L,bi,bj)+turbu(i,j,L,bi,bj) )*86400
./pkg/fizhi/fizhi_step_diag.F 176 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'DIABU   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 184 FIZHI_STEP_DIAG
        tmpdiag(i,j) = (moistv (i,j,L,bi,bj)+turbv(i,j,L,bi,bj) )*86400
./pkg/fizhi/fizhi_step_diag.F 187 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'DIABV   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 196 FIZHI_STEP_DIAG
     .   ( turbt(i,j,L,bi,bj) + moistt(i,j,L,bi,bj) +
./pkg/fizhi/fizhi_step_diag.F 197 FIZHI_STEP_DIAG
     .      lwdt(i,j,L,bi,bj) +
./pkg/fizhi/fizhi_step_diag.F 198 FIZHI_STEP_DIAG
     .      dlwdtg(i,j,L,bi,bj)*(tgz(i,j,bi,bj)-tg0(i,j,bi,bj)) +
./pkg/fizhi/fizhi_step_diag.F 199 FIZHI_STEP_DIAG
     .      swdt(i,j,L,bi,bj)*radswt(i,j,bi,bj) )
./pkg/fizhi/fizhi_step_diag.F 200 FIZHI_STEP_DIAG
     .      * pk(i,j,L,bi,bj)*pinv(i,j)*86400
./pkg/fizhi/fizhi_step_diag.F 203 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'DIABT   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 212 FIZHI_STEP_DIAG
     . ( turbq(i,j,L,1,bi,bj) + moistq(i,j,L,1,bi,bj) ) * 
./pkg/fizhi/fizhi_step_diag.F 216 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'DIABQ   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 226 FIZHI_STEP_DIAG
     .            dlwdtg (i,j,L,bi,bj)*(tgz(i,j,bi,bj)-tg0(i,j,bi,bj)) )
./pkg/fizhi/fizhi_step_diag.F 230 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'RADLW   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 240 FIZHI_STEP_DIAG
     .            dlwdtg (i,j,L,bi,bj)*(tgz(i,j,bi,bj)-tg0(i,j,bi,bj)) )
./pkg/fizhi/fizhi_step_diag.F 244 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'LWCLR   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 257 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'RADSW   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 270 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'SWCLR   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 278 FIZHI_STEP_DIAG
        tmpdiag(i,j) = uphy(i,j,L)
./pkg/fizhi/fizhi_step_diag.F 281 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'UWND    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 289 FIZHI_STEP_DIAG
        tmpdiag(i,j) = vphy(i,j,L)
./pkg/fizhi/fizhi_step_diag.F 292 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'VWND    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 300 FIZHI_STEP_DIAG
        tmpdiag(i,j) = thphy(i,j,L)*pk(i,j,L,bi,bj)
./pkg/fizhi/fizhi_step_diag.F 303 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'TMPU    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 311 FIZHI_STEP_DIAG
        tmpdiag(i,j) = qq(i,j,L,bi,bj)
./pkg/fizhi/fizhi_step_diag.F 314 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'TKE     ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 322 FIZHI_STEP_DIAG
        tmpdiag(i,j) = sphy(i,j,L) * 1000.
./pkg/fizhi/fizhi_step_diag.F 325 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'SPHU    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 340 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 344 FIZHI_STEP_DIAG
     .             moistt(i,j,L,bi,bj)*pk(i,j,l,bi,bj)*dp(i,j,L,bi,bj)
./pkg/fizhi/fizhi_step_diag.F 364 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 368 FIZHI_STEP_DIAG
     .             turbt(i,j,L,bi,bj)*pk(i,j,l,bi,bj)*dp(i,j,L,bi,bj)
./pkg/fizhi/fizhi_step_diag.F 388 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 391 FIZHI_STEP_DIAG
        qbar(i,j) = qbar(i,j) + ( lwdt(i,j,L,bi,bj) +
./pkg/fizhi/fizhi_step_diag.F 392 FIZHI_STEP_DIAG
     .  dlwdtg(i,j,L,bi,bj)*(tgz(i,j,bi,bj)-tg0(i,j,bi,bj)) )
./pkg/fizhi/fizhi_step_diag.F 393 FIZHI_STEP_DIAG
     .             *pk(i,j,l,bi,bj)*dp(i,j,L,bi,bj)
./pkg/fizhi/fizhi_step_diag.F 413 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 417 FIZHI_STEP_DIAG
     .             swdt(i,j,L,bi,bj)*pk(i,j,l,bi,bj)*dp(i,j,L,bi,bj)
./pkg/fizhi/fizhi_step_diag.F 41 FIZHI_STEP_DIAG
      integer  i,j,L
./pkg/fizhi/fizhi_step_diag.F 439 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 443 FIZHI_STEP_DIAG
     .             sphy(i,j,L)*dp(i,j,L,bi,bj)
./pkg/fizhi/fizhi_swrad.F 114 SWRIO
      do L=1,lm
./pkg/fizhi/fizhi_swrad.F 117 SWRIO
      tdry(i,j,L) = tz(i,j,L)*pkz(i,j,L)
./pkg/fizhi/fizhi_swrad.F 177 SWRIO
        do L =1,lm
./pkg/fizhi/fizhi_swrad.F 180 SWRIO
           cldtot(i,j,L)=min(1.0 _d 0,max(cldsw(i,j,L),fccave(i,j,L)/
./pkg/fizhi/fizhi_swrad.F 182 SWRIO
           cldmxo(i,j,L)=min(1.0 _d 0,cswmo(i,j,L))
./pkg/fizhi/fizhi_swrad.F 183 SWRIO
           swlz(i,j,L)=swlz(i,j,L)+qliqave(i,j,L)/imstturb
./pkg/fizhi/fizhi_swrad.F 188 SWRIO
        do L =1,lm
./pkg/fizhi/fizhi_swrad.F 191 SWRIO
         cldtot(i,j,L) =  min( 1.0 _d 0,cldsw(i,j,L) )
./pkg/fizhi/fizhi_swrad.F 192 SWRIO
         cldmxo(i,j,L) =  min( 1.0 _d 0,cswmo(i,j,L) )
./pkg/fizhi/fizhi_swrad.F 201 SWRIO
      do L = 1,lm
./pkg/fizhi/fizhi_swrad.F 206 SWRIO
      if(     L.lt.mid_level ) then
./pkg/fizhi/fizhi_swrad.F 207 SWRIO
             cldhi (i,j) = max( cldtot(i,j,L),cldhi (i,j) )
./pkg/fizhi/fizhi_swrad.F 208 SWRIO
      elseif( L.ge.mid_level .and.
./pkg/fizhi/fizhi_swrad.F 209 SWRIO
     .        L.lt.low_level ) then
./pkg/fizhi/fizhi_swrad.F 210 SWRIO
             cldmid(i,j) = max( cldtot(i,j,L),cldmid(i,j) )
./pkg/fizhi/fizhi_swrad.F 211 SWRIO
      elseif( L.ge.low_level ) then
./pkg/fizhi/fizhi_swrad.F 212 SWRIO
             cldlow(i,j) = max( cldtot(i,j,L),cldlow(i,j) )
./pkg/fizhi/fizhi_swrad.F 239 SWRIO
       do L=1,lm
./pkg/fizhi/fizhi_swrad.F 242 SWRIO
        tmpdiag(i,j) = cswmo(i,j,L)
./pkg/fizhi/fizhi_swrad.F 245 SWRIO
       call DIAGNOSTICS_FILL(tmpdiag,'CLDRAS  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_swrad.F 250 SWRIO
       do L=1,lm
./pkg/fizhi/fizhi_swrad.F 253 SWRIO
        tmpdiag(i,j) = cldtot(i,j,L)
./pkg/fizhi/fizhi_swrad.F 256 SWRIO
       call DIAGNOSTICS_FILL(tmpdiag,'CLDTOT  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_swrad.F 273 SWRIO
       do L=1,lm
./pkg/fizhi/fizhi_swrad.F 276 SWRIO
        tmpdiag(i,j) = swlz(i,j,L) * 1.0e6
./pkg/fizhi/fizhi_swrad.F 279 SWRIO
       call DIAGNOSTICS_FILL(tmpdiag,'LZRAD   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_swrad.F 308 SWRIO
      do L = 1,lm
./pkg/fizhi/fizhi_swrad.F 311 SWRIO
      tau(i,j,L) = tautype(i,j,L,1)+tautype(i,j,L,2)+tautype(i,j,L,3)
./pkg/fizhi/fizhi_swrad.F 318 SWRIO
       do L=1,lm
./pkg/fizhi/fizhi_swrad.F 321 SWRIO
        tmpdiag(i,j) = tau(i,j,L)*100/(plze(i,j,L+1)-plze(i,j,L))
./pkg/fizhi/fizhi_swrad.F 324 SWRIO
       call DIAGNOSTICS_FILL(tmpdiag,'TAUAVE  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_swrad.F 330 SWRIO
       do L=1,lm
./pkg/fizhi/fizhi_swrad.F 333 SWRIO
        if( cldtot(i,j,L).ne.0.0 ) then
./pkg/fizhi/fizhi_swrad.F 334 SWRIO
        tmpdiag(i,j) = tau(i,j,L)*100/(plze(i,j,L+1)-plze(i,j,L))
./pkg/fizhi/fizhi_swrad.F 342 SWRIO
       call DIAGNOSTICS_FILL(tmpdiag,'TAUCLD  ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_swrad.F 343 SWRIO
       call DIAGNOSTICS_FILL(tmpdiag,'TAUCLDC ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_swrad.F 355 SWRIO
         do L = low_level,lm
./pkg/fizhi/fizhi_swrad.F 356 SWRIO
          taulow(i,j) = taulow(i,j) + tau(i,j,L)
./pkg/fizhi/fizhi_swrad.F 374 SWRIO
         do L = mid_level,low_level+1
./pkg/fizhi/fizhi_swrad.F 375 SWRIO
          taumid(i,j) = taumid(i,j) + tau(i,j,L)
./pkg/fizhi/fizhi_swrad.F 393 SWRIO
         do L = 1,mid_level+1
./pkg/fizhi/fizhi_swrad.F 394 SWRIO
          tauhi(i,j) = tauhi(i,j) + tau(i,j,L)
./pkg/fizhi/fizhi_swrad.F 42 SWRIO
      integer   i,j,L,nn,nsecf
./pkg/fizhi/fizhi_swrad.F 436 SWRIO
      do L= 1,lm
./pkg/fizhi/fizhi_swrad.F 439 SWRIO
      taua(i,L)   = 0.
./pkg/fizhi/fizhi_swrad.F 442 SWRIO
      tauc(i,L,1) = taul(i,L)/(1.+alf/(1-alf) * (reffi/reffw*0.80) )
./pkg/fizhi/fizhi_swrad.F 443 SWRIO
      tauc(i,L,2) = taul(i,L)/(1.+(1-alf)/alf * (reffw/reffi*1.25) )
./pkg/fizhi/fizhi_swrad.F 447 SWRIO
      tauc(i,L,1) = taul(i,L)
./pkg/fizhi/fizhi_swrad.F 448 SWRIO
      tauc(i,L,2) = 0.0
./pkg/fizhi/fizhi_swrad.F 452 SWRIO
      tauc(i,L,1) = 0.0
./pkg/fizhi/fizhi_swrad.F 453 SWRIO
      tauc(i,L,2) = taul(i,L)
./pkg/fizhi/fizhi_swrad.F 456 SWRIO
      reff(i,L,1) = reffi
./pkg/fizhi/fizhi_swrad.F 457 SWRIO
      reff(i,L,2) = reffw
./pkg/fizhi/fizhi_swrad.F 473 SWRIO
      alf = grav*(ple(i,Lm+1)-ptop)/(cp*dpstrip(i,L)*100)
./pkg/fizhi/fizhi_swrad.F 474 SWRIO
      dtsw (i,L) = alf*( flux   (i,L)-flux   (i,L+1) )/pk(i,L)
./pkg/fizhi/fizhi_swrad.F 475 SWRIO
      dtswc(i,L) = alf*( fluxclr(i,L)-fluxclr(i,L+1) )/pk(i,L)
./pkg/fizhi/fizhi_swrad.F 536 SWRIO
      do L = 1,lm
./pkg/fizhi/fizhi_swrad.F 539 SWRIO
         fccave(i,j,L) = 0.
./pkg/fizhi/fizhi_swrad.F 540 SWRIO
        qliqave(i,j,L) = 0.
./pkg/fizhi/fizhi_swrad.F 579 OPTHK
      integer  im,jm,lm,i,j,L
./pkg/fizhi/fizhi_swrad.F 595 OPTHK
      do L=1,lm
./pkg/fizhi/fizhi_swrad.F 598 OPTHK
         alf   =  min( max(( tl(i,j,L)-233.15)/20.,0. _d 0), 1. _d 0)
./pkg/fizhi/fizhi_swrad.F 599 OPTHK
         dp   =  ple(i,j,L+1)-ple(i,j,L)
./pkg/fizhi/fizhi_swrad.F 600 OPTHK
         tau(i,j,L,1)  = 0.0
./pkg/fizhi/fizhi_swrad.F 601 OPTHK
         tau(i,j,L,2)  = 0.0
./pkg/fizhi/fizhi_swrad.F 602 OPTHK
         tau(i,j,L,3)  = 0.0
./pkg/fizhi/fizhi_swrad.F 603 OPTHK
      if( cf(i,j,L).ne.0.0 ) then
./pkg/fizhi/fizhi_swrad.F 607 OPTHK
         fracls = ( cf(i,j,L)-cfm(i,j,L) )/cf(i,j,L)
./pkg/fizhi/fizhi_swrad.F 618 OPTHK
         tauice = max( 0.0002 _d 0, 0.002*min(500*lz(i,j,L)*1000,
./pkg/fizhi/fizhi_swrad.F 620 OPTHK
         tau(i,j,L,1) = fracls*(1-alf)*tauice*dp
./pkg/fizhi/fizhi_swrad.F 626 OPTHK
            tauh2o = max( 0.0020 _d 0, 0.200*min(200*lz(i,j,L)*1000,
./pkg/fizhi/fizhi_swrad.F 628 OPTHK
          tau(i,j,L,3) = fracls*alf*tauh2o*dp
./pkg/fizhi/fizhi_swrad.F 631 OPTHK
          if( lz(i,j,L).eq.0.0 ) then
./pkg/fizhi/fizhi_swrad.F 633 OPTHK
           tau(i,j,L,2) = fracls*alf*tauh2o*dp
./pkg/fizhi/fizhi_swrad.F 636 OPTHK
             tauh2o = max( 0.0020 _d 0, 0.120*min( 20*lz(i,j,L)*1000,
./pkg/fizhi/fizhi_swrad.F 638 OPTHK
           tau(i,j,L,3) = fracls*alf*tauh2o*dp
./pkg/fizhi/fizhi_swrad.F 644 OPTHK
         if( tl(i,j,L).gt.210.0 ) then
./pkg/fizhi/fizhi_swrad.F 646 OPTHK
         tau(i,j,L,3) = tau(i,j,L,3) + fraccu*tauras*dp
./pkg/fizhi/fizhi_swrad.F 649 OPTHK
         tau(i,j,L,1) = tau(i,j,L,1) + fraccu*tauras*dp
./pkg/fizhi/fizhi_turb.F 1222 TURBIO
      do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 1224 TURBIO
           duturb(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1225 TURBIO
           dvturb(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1226 TURBIO
           dtturb(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1227 TURBIO
           qqgrid(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1228 TURBIO
          qliqtmp(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1229 TURBIO
           fcctmp(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1233 TURBIO
           dqturb(i,1,L,nt) = 0.
./pkg/fizhi/fizhi_turb.F 1241 TURBIO
      call MSC2GRD(igrd,chfr,dumsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1242 TURBIO
     .                                              duturb(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1243 TURBIO
      call MSC2GRD(igrd,chfr,dvmsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1244 TURBIO
     .                                              dvturb(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1245 TURBIO
      call MSC2GRD(igrd,chfr,dtmsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1246 TURBIO
     .                                              dtturb(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1248 TURBIO
      call MSC2GRD(igrd,chfr,dqmsc(1,L,nt),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1249 TURBIO
     .                                           dqturb(1,1,L,nt),im,jm)
./pkg/fizhi/fizhi_turb.F 1251 TURBIO
      call MSC2GRD(igrd,chfr,  tke(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1252 TURBIO
     .                                              qqgrid(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1254 TURBIO
      call MSC2GRD(igrd,chfr, fccmsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1255 TURBIO
     .                                              fcctmp(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1256 TURBIO
      call MSC2GRD(igrd,chfr,qliqmsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1257 TURBIO
     .                                             qliqtmp(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1276 TURBIO
         fccavesw(i,1,L) =  fccavesw(i,1,L) +  fcctmp(i,1,L)
./pkg/fizhi/fizhi_turb.F 1277 TURBIO
         fccavelw(i,1,L) =  fccavelw(i,1,L) +  fcctmp(i,1,L)
./pkg/fizhi/fizhi_turb.F 1278 TURBIO
        qliqavelw(i,1,L) = qliqavelw(i,1,L) + qliqtmp(i,1,L)
./pkg/fizhi/fizhi_turb.F 1279 TURBIO
        qliqavesw(i,1,L) = qliqavesw(i,1,L) + qliqtmp(i,1,L)
./pkg/fizhi/fizhi_turb.F 1385 TURBIO
        do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 1389 TURBIO
         call DIAGNOSTICS_FILL(tmpdiag,'TURBT   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_turb.F 1613 TRBFLX
      integer i,iter,init,n,LL,L,Lp,Lp1,lmin,lminq,lminq1,ibit
./pkg/fizhi/fizhi_turb.F 1665 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1667 TRBFLX
         wu(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1668 TRBFLX
         wv(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1669 TRBFLX
         eu(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1670 TRBFLX
         et(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1674 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1676 TRBFLX
         XLDIAG(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1680 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1682 TRBFLX
         FREQDG(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1733 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1735 TRBFLX
        XLSAVE(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1736 TRBFLX
        KHSAVE(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1754 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1756 TRBFLX
        ADZ1(I,L) = (CP/GRAV)*(PLKE(I,L+1)-PLKE(I,L))
./pkg/fizhi/fizhi_turb.F 1757 TRBFLX
        ADZ1(I,L) = THV(I,L) * ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1758 TRBFLX
        DZ1TMP(I,L) = ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1764 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1766 TRBFLX
        ADZ2(I,L) = 0.5 * (ADZ1(I,L)+ADZ1(I,L+1))
./pkg/fizhi/fizhi_turb.F 1782 TRBFLX
      DO L =2,NLEVM1
./pkg/fizhi/fizhi_turb.F 1784 TRBFLX
        DZ3(I,L) = ALPHA * ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1796 TRBFLX
      DO L =2,NLEV
./pkg/fizhi/fizhi_turb.F 1798 TRBFLX
        TEMP(I,L) = VK * ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1805 TRBFLX
      L = NLEV - LL
./pkg/fizhi/fizhi_turb.F 1806 TRBFLX
      LP1 = L + 1
./pkg/fizhi/fizhi_turb.F 1808 TRBFLX
      VKZE(I,L) = VKZE(I,LP1) + TEMP(I,LP1)
./pkg/fizhi/fizhi_turb.F 1814 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1816 TRBFLX
        VKZM(I,L) = VKZE(I,L) - 0.5 * TEMP(I,L+1)
./pkg/fizhi/fizhi_turb.F 1821 TRBFLX
      DO 200 L = 1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1822 TRBFLX
       LP1 = L + 1
./pkg/fizhi/fizhi_turb.F 1824 TRBFLX
       FAC1 = DPSTR(I,L) / ( DPSTR(I,L) + DPSTR(I,LP1) )
./pkg/fizhi/fizhi_turb.F 1826 TRBFLX
       RHODZ2(I,L) = FAC1 * THV(I,LP1)
./pkg/fizhi/fizhi_turb.F 1827 TRBFLX
       RHODZ2(I,L) = RHODZ2(I,L) + FAC2 * THV(I,L)
./pkg/fizhi/fizhi_turb.F 1841 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1843 TRBFLX
        RHODZ2(I,L) = (RGAS*0.01) * RHODZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1844 TRBFLX
        TEMP(I,L)   = PLKE(I,L+1) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1845 TRBFLX
        RHODZ2(I,L) = TEMP(I,L) * RHODZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1846 TRBFLX
        RHODZ2(I,L) = PLE(I,L+1) / RHODZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1847 TRBFLX
        RHOZPK(I,L) = RHODZ2(I,L) * PLKE(I,L+1)
./pkg/fizhi/fizhi_turb.F 1848 TRBFLX
        RHODZ1(I,L) = (RGAS*0.01) * THV(I,L+1)
./pkg/fizhi/fizhi_turb.F 1849 TRBFLX
        TEMP(I,L)   = PLK(I,L+1) * ADZ1(I,L+1)
./pkg/fizhi/fizhi_turb.F 1850 TRBFLX
        RHODZ1(I,L) = TEMP(I,L) * RHODZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1851 TRBFLX
        RHODZ1(I,L) = PL(I,L+1) / RHODZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1862 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1864 TRBFLX
        FLXFPK(I,L) = PLE(I,L+1) - PLE(I,L)
./pkg/fizhi/fizhi_turb.F 1865 TRBFLX
        FLXFPK(I,L) = FLXFPK(I,L) * PLK(I,L)
./pkg/fizhi/fizhi_turb.F 1866 TRBFLX
        FLXFPK(I,L) = (GRAV*DTAU*0.01) / FLXFPK(I,L)
./pkg/fizhi/fizhi_turb.F 1878 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1880 TRBFLX
        FLXFAC(I,L) = FLXFPK(I,L) * PLK(I,L)
./pkg/fizhi/fizhi_turb.F 1895 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1897 TRBFLX
        FLXFCE(I,L) = PL(I,L+1) - PL(I,L)
./pkg/fizhi/fizhi_turb.F 1898 TRBFLX
        FLXFCE(I,L) = (GRAV*DTAU*0.01) / FLXFCE(I,L)
./pkg/fizhi/fizhi_turb.F 1906 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1908 TRBFLX
        ADZ1(I,L) = 1. / ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1914 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1916 TRBFLX
        ADZ2(I,L) = 1. / ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1927 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1929 TRBFLX
        DPK(I,L) = ( PLK(I,L+1)-PLK(I,L) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1941 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1943 TRBFLX
        Q(I,L) = 2. * QQ(I,L)
./pkg/fizhi/fizhi_turb.F 1944 TRBFLX
        Q(I,L) = SQRT( Q(I,L) )
./pkg/fizhi/fizhi_turb.F 2049 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2051 TRBFLX
        DU(I,L) = (   U(I,L)-  U(I,L+1) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 2052 TRBFLX
        DV(I,L) = (   V(I,L)-  V(I,L+1) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 2063 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2065 TRBFLX
         XL(I,L) = XLSAVE(I,L)
./pkg/fizhi/fizhi_turb.F 2075 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2077 TRBFLX
        DTH(I,L) = ( TH(I,L)-TH(I,L+1) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 2078 TRBFLX
        DSH(I,L) = ( SH(I,L)-SH(I,L+1) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 2079 TRBFLX
        TL(I,L) = TH(I,L)*PLK(I,L)
./pkg/fizhi/fizhi_turb.F 2094 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2096 TRBFLX
        BB(I,L) = FACEPS*SHSAT(I,L)/(TL(I,L)*TL(I,L))
./pkg/fizhi/fizhi_turb.F 2097 TRBFLX
        AA(I,L) = 1. / (1. + CLH * BB(I,L) )
./pkg/fizhi/fizhi_turb.F 2099 TRBFLX
        BB(I,L) = BB(I,L) * AA(I,L)
./pkg/fizhi/fizhi_turb.F 2100 TRBFLX
        SBAR(I,L) = AA(I,L) * (SH(I,L) - SHSAT(I,L))
./pkg/fizhi/fizhi_turb.F 2125 TRBFLX
      DO L =2,NLEVM1
./pkg/fizhi/fizhi_turb.F 2128 TRBFLX
        SSDEV(I,L) = XL(I,L-1)*(AA(I,L)*DSH(I,L-1)-
./pkg/fizhi/fizhi_turb.F 2129 TRBFLX
     1               BB(I,L)*plke(I,L)*DTH(I,L-1))
./pkg/fizhi/fizhi_turb.F 2130 TRBFLX
        SSDEV(I,L) = B2 * KHSAVE(I,L-1) * SSDEV(I,L-1) * SSDEV(I,L-1)
./pkg/fizhi/fizhi_turb.F 2131 TRBFLX
        SVAR(I,L) = SQRT(SSDEV(I,L))
./pkg/fizhi/fizhi_turb.F 2133 TRBFLX
        SSDEV(I,L) = XL(I,L)*(AA(I,L)*DSH(I,L)-
./pkg/fizhi/fizhi_turb.F 2134 TRBFLX
     1               BB(I,L)*plke(I,L+1)*DTH(I,L))
./pkg/fizhi/fizhi_turb.F 2135 TRBFLX
        SSDEV(I,L) = B2 * KHSAVE(I,L) * SSDEV(I,L) * SSDEV(I,L)
./pkg/fizhi/fizhi_turb.F 2136 TRBFLX
        TEMP(I,L) = SQRT(SSDEV(I,L))
./pkg/fizhi/fizhi_turb.F 2137 TRBFLX
        SVAR(I,L) = (1./2.) * (SVAR(I,L) + TEMP(I,L))
./pkg/fizhi/fizhi_turb.F 2138 TRBFLX
        IF ( SVAR(I,L).LT.Z1PEM25) SVAR(I,L) = Z1PEM25
./pkg/fizhi/fizhi_turb.F 2158 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2160 TRBFLX
        Q1M(I,L) = SBAR(I,L) / SVAR(I,L)
./pkg/fizhi/fizhi_turb.F 2161 TRBFLX
        FCC(I,L) = (1./2.) * ( 1. + ERRF( P5SR*Q1M(I,L) ) )
./pkg/fizhi/fizhi_turb.F 2162 TRBFLX
        SHL(I,L) =  FCC(I,L) * SBAR(I,L)
./pkg/fizhi/fizhi_turb.F 2163 TRBFLX
        ARG(I,L)  = (1./2.)*Q1M(I,L)*Q1M(I,L)
./pkg/fizhi/fizhi_turb.F 2164 TRBFLX
        IF(ARG(I,L).LE.ARGMAX)
./pkg/fizhi/fizhi_turb.F 2165 TRBFLX
     1        SHL(I,L) =  SHL(I,L)+RSQ2PI*SVAR(I,L)*EXP(-ARG(I,L))
./pkg/fizhi/fizhi_turb.F 2166 TRBFLX
        BETAT(I,L) = 1. + VIRTCON * SH(I,L) - VRT1CON * SHL(I,L)
./pkg/fizhi/fizhi_turb.F 2167 TRBFLX
        BETAW(I,L) = VIRTCON *
./pkg/fizhi/fizhi_turb.F 2168 TRBFLX
     1                ( TH(I,L) + (CLH/plk(I,L)) * SHL(I,L) )
./pkg/fizhi/fizhi_turb.F 2169 TRBFLX
        BETAL(I,L) = ( 1. + VIRTCON*SH(I,L) - TWO*VRT1CON*SHL(I,L) )
./pkg/fizhi/fizhi_turb.F 2170 TRBFLX
     1     * (CLH/plke(I,L+1)) - VRT1CON * TH(I,L)
./pkg/fizhi/fizhi_turb.F 2172 TRBFLX
        BETAT1(I,L) = BETAT(I,L) -
./pkg/fizhi/fizhi_turb.F 2173 TRBFLX
     1               BB(I,L)*plk(i,L) * FCC(I,L) * BETAL(I,L)
./pkg/fizhi/fizhi_turb.F 2174 TRBFLX
        BETAW1(I,L) = BETAW(I,L) + AA(I,L) * FCC(I,L) * BETAL(I,L)
./pkg/fizhi/fizhi_turb.F 2181 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2183 TRBFLX
        DTHV(I,L) =  (1./2.)*((BETAT1(I,L)+BETAT1(I,L+1))*DTH(I,L)
./pkg/fizhi/fizhi_turb.F 2184 TRBFLX
     1               + (BETAW1(I,L)+BETAW1(I,L+1))*DSH(I,L))
./pkg/fizhi/fizhi_turb.F 2206 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2208 TRBFLX
        STRT(I,L) = CP * DTHV(I,L) * DPK(I,L)
./pkg/fizhi/fizhi_turb.F 2209 TRBFLX
        DW2(I,L) = DU(I,L) * DU(I,L) + DV(I,L) * DV(I,L)
./pkg/fizhi/fizhi_turb.F 2210 TRBFLX
        IF ( DW2(I,L) .LE. 1.e-4 ) DW2(I,L) = 1.e-4
./pkg/fizhi/fizhi_turb.F 2211 TRBFLX
        RI(I,L) = STRT(I,L) / DW2(I,L)
./pkg/fizhi/fizhi_turb.F 2220 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2222 TRBFLX
        SRI(I,L) = RI(I,L)
./pkg/fizhi/fizhi_turb.F 2256 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2258 TRBFLX
        KH(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2259 TRBFLX
        KM(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2260 TRBFLX
        QQE(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2261 TRBFLX
        QE(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2262 TRBFLX
        P3(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2263 TRBFLX
        IBITSTB(I,L) = 0
./pkg/fizhi/fizhi_turb.F 2264 TRBFLX
        IF ( QQ(I,L) .GT. 1.e-8 ) THEN
./pkg/fizhi/fizhi_turb.F 2265 TRBFLX
         INTQ(I,L) = 1
./pkg/fizhi/fizhi_turb.F 2267 TRBFLX
         INTQ(I,L) = 0
./pkg/fizhi/fizhi_turb.F 2269 TRBFLX
        IF ( QQ(I,L).LE.1.e-8 ) THEN
./pkg/fizhi/fizhi_turb.F 2270 TRBFLX
         QQ(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2271 TRBFLX
         Q(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2308 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2310 TRBFLX
        QQ(I,L) = QQE(I,L)
./pkg/fizhi/fizhi_turb.F 2311 TRBFLX
        Q(I,L) = QE(I,L)
./pkg/fizhi/fizhi_turb.F 2333 TRBFLX
       DO L =LMIN,NLEVM1
./pkg/fizhi/fizhi_turb.F 2335 TRBFLX
         IF (  (IBITSTB(I,L).EQ.1)  .AND.
./pkg/fizhi/fizhi_turb.F 2336 TRBFLX
     1             ( Q(I,L) .LE. QE(I,L) ) ) THEN
./pkg/fizhi/fizhi_turb.F 2337 TRBFLX
          IBITSTB(I,L) = 1
./pkg/fizhi/fizhi_turb.F 2339 TRBFLX
          IBITSTB(I,L) = 0
./pkg/fizhi/fizhi_turb.F 2358 TRBFLX
       DO L =LMIN,NLEVM1
./pkg/fizhi/fizhi_turb.F 2360 TRBFLX
         IF(IBITSTB(I,L).EQ.1 ) THEN
./pkg/fizhi/fizhi_turb.F 2361 TRBFLX
          TEMP(I,L) = Q(I,L) / QE(I,L)
./pkg/fizhi/fizhi_turb.F 2362 TRBFLX
          KH(I,L) = TEMP(I,L) * KH(I,L)
./pkg/fizhi/fizhi_turb.F 2363 TRBFLX
          KM(I,L) = TEMP(I,L) * KM(I,L)
./pkg/fizhi/fizhi_turb.F 2365 TRBFLX
         TEMP(I,L) = 0.01 * QQE(I,L)
./pkg/fizhi/fizhi_turb.F 2366 TRBFLX
         IF((IBITSTB(I,L).EQ.1) .AND.
./pkg/fizhi/fizhi_turb.F 2367 TRBFLX
     1         ( QQ(I,L) .LE. TEMP(I,L) )) THEN
./pkg/fizhi/fizhi_turb.F 2368 TRBFLX
          QQ(I,L) = TEMP(I,L)
./pkg/fizhi/fizhi_turb.F 2369 TRBFLX
          Q(I,L) = 0.1 * QE(I,L)
./pkg/fizhi/fizhi_turb.F 2371 TRBFLX
         IF(IBITSTB(I,L).EQ.1 ) P3(I,L) = (2.*B3) *
./pkg/fizhi/fizhi_turb.F 2372 TRBFLX
     1             ( QE(I,L) - Q(I,L) )
./pkg/fizhi/fizhi_turb.F 2394 TRBFLX
       DO L =LMIN,NLEVM1
./pkg/fizhi/fizhi_turb.F 2396 TRBFLX
         P3(I,L) = P3(I,L) * DTAU / XL(I,L)
./pkg/fizhi/fizhi_turb.F 2397 TRBFLX
         TEMP(I,L) = QQE(I,L) * P3(I,L)
./pkg/fizhi/fizhi_turb.F 2398 TRBFLX
         XQ(I,L) = QQE(I,L) - QQ(I,L)
./pkg/fizhi/fizhi_turb.F 2409 TRBFLX
       DO L =LMIN,NLEVM1
./pkg/fizhi/fizhi_turb.F 2411 TRBFLX
         IF( ( (IBITSTB(I,L).EQ.1) .AND.
./pkg/fizhi/fizhi_turb.F 2412 TRBFLX
     1             ( XQ(I,L) .LT. TEMP(I,L) ) )
./pkg/fizhi/fizhi_turb.F 2414 TRBFLX
     3       ( (IBITSTB(I,L).EQ.0) .AND.
./pkg/fizhi/fizhi_turb.F 2415 TRBFLX
     4             ( XQ(I,L) .GT. TEMP(I,L) ) ) )
./pkg/fizhi/fizhi_turb.F 2416 TRBFLX
     5   P3(I,L) = XQ(I,L) / QQE(I,L)
./pkg/fizhi/fizhi_turb.F 2437 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2439 TRBFLX
        RIINIT(I,L) = RI(I,L)
./pkg/fizhi/fizhi_turb.F 2440 TRBFLX
        QQINIT(I,L) = QQ(I,L)
./pkg/fizhi/fizhi_turb.F 2446 TRBFLX
      do L = 1,nlev-1
./pkg/fizhi/fizhi_turb.F 2451 TRBFLX
         if(pl(i,L).le.150.) then
./pkg/fizhi/fizhi_turb.F 2452 TRBFLX
          qxlm(i,L) = min(qxlm(i,L),5. _d 0)
./pkg/fizhi/fizhi_turb.F 2465 TRBFLX
      DO L =LMINQ1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2467 TRBFLX
        RHOKDZ(I,L) = RHODZ1(I,L) * QXLM(I,L)
./pkg/fizhi/fizhi_turb.F 2478 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2480 TRBFLX
        KHSAVE(I,L) = KH(I,L)
./pkg/fizhi/fizhi_turb.F 2493 TRBFLX
       DO L =1,LMINQ1-1
./pkg/fizhi/fizhi_turb.F 2495 TRBFLX
         KM(I,L) = KMBG
./pkg/fizhi/fizhi_turb.F 2496 TRBFLX
         KH(I,L) = KHBG
./pkg/fizhi/fizhi_turb.F 2510 TRBFLX
      DO L =LMINQ1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2512 TRBFLX
        Q(I,L) = 2. * QQ(I,L)
./pkg/fizhi/fizhi_turb.F 2513 TRBFLX
        Q(I,L) = SQRT(Q(I,L))
./pkg/fizhi/fizhi_turb.F 2514 TRBFLX
        XQ(I,L)= XL(I,L) * Q(I,L)
./pkg/fizhi/fizhi_turb.F 2515 TRBFLX
        KM(I,L)=XQ(I,L)*KM(I,L)+KMBG
./pkg/fizhi/fizhi_turb.F 2516 TRBFLX
        KH(I,L)=XQ(I,L)*KH(I,L)+KHBG
./pkg/fizhi/fizhi_turb.F 2521 TRBFLX
      do L = 1,nlev-1
./pkg/fizhi/fizhi_turb.F 2527 TRBFLX
         if(pl(i,L).le.150.) then
./pkg/fizhi/fizhi_turb.F 2528 TRBFLX
          kh(i,L) = min(kh(i,L),5. _d 0)
./pkg/fizhi/fizhi_turb.F 2529 TRBFLX
          km(i,L) = min(km(i,L),5. _d 0)
./pkg/fizhi/fizhi_turb.F 2542 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2544 TRBFLX
        TEMP(I,L) = RHOZPK(I,L) * KH(I,L)
./pkg/fizhi/fizhi_turb.F 2545 TRBFLX
        DELTH(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2556 TRBFLX
      do L = 1,nlev
./pkg/fizhi/fizhi_turb.F 2558 TRBFLX
        dthdthg(i,L) = delth(i,L)
./pkg/fizhi/fizhi_turb.F 2561 TRBFLX
      do L = 1,nlev
./pkg/fizhi/fizhi_turb.F 2563 TRBFLX
        transth(i,L) = temp(i,L)
./pkg/fizhi/fizhi_turb.F 2573 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2575 TRBFLX
        RHOKDZ(I,L) = RHODZ2(I,L) * KH(I,L)
./pkg/fizhi/fizhi_turb.F 2576 TRBFLX
        DELSH(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2589 TRBFLX
      do L = 1,nlev
./pkg/fizhi/fizhi_turb.F 2591 TRBFLX
        dshdshg(i,L) = delsh(i,L)
./pkg/fizhi/fizhi_turb.F 2594 TRBFLX
      do L = 1,nlev
./pkg/fizhi/fizhi_turb.F 2596 TRBFLX
        transsh(i,L) = rhokdz(i,L)
./pkg/fizhi/fizhi_turb.F 259 TURBIO
      integer istnp1,istnlay,itrtrb,i,L,nn,nt
./pkg/fizhi/fizhi_turb.F 2620 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2622 TRBFLX
        RHOKDZ(I,L) = RHODZ2(I,L) * KM(I,L)
./pkg/fizhi/fizhi_turb.F 2633 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2635 TRBFLX
        WU(I,L) = WU(I,L) + RHOKDZ(I,L) * ( U(I,L+1) - U(I,L) )
./pkg/fizhi/fizhi_turb.F 2636 TRBFLX
        WV(I,L) = WV(I,L) + RHOKDZ(I,L) * ( V(I,L+1) - V(I,L) )
./pkg/fizhi/fizhi_turb.F 2647 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2649 TRBFLX
        IF ( QQ(I,L) .GT. QQMIN ) THEN
./pkg/fizhi/fizhi_turb.F 2650 TRBFLX
         IBITSTB(I,L) = 1
./pkg/fizhi/fizhi_turb.F 2652 TRBFLX
         IBITSTB(I,L) = 0
./pkg/fizhi/fizhi_turb.F 2654 TRBFLX
        IF( IBITSTB(I,L).EQ.1 ) FREQDG(I,L) = FREQDG(I,L) + aitr
./pkg/fizhi/fizhi_turb.F 2662 TRBFLX
       DO L = nlev-1,1,-1
./pkg/fizhi/fizhi_turb.F 2696 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2698 TRBFLX
        EU(I,L) = EU(I,L) + AITR*KM(I,L)
./pkg/fizhi/fizhi_turb.F 2699 TRBFLX
        ET(I,L) = ET(I,L) + AITR*KH(I,L)
./pkg/fizhi/fizhi_turb.F 2712 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2714 TRBFLX
         XLDIAG(I,L) = XLDIAG(I,L) + AITR*XL(I,L)
./pkg/fizhi/fizhi_turb.F 2726 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2728 TRBFLX
         XLSAVE(I,L) = XL(I,L)
./pkg/fizhi/fizhi_turb.F 2745 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2747 TRBFLX
        turbfcc(I,L) = turbfcc(I,L) + fcc(I,L) * aitr
./pkg/fizhi/fizhi_turb.F 2748 TRBFLX
        qliq(I,L) = qliq(I,L) + shl(I,L) * aitr
./pkg/fizhi/fizhi_turb.F 2759 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2761 TRBFLX
        WU(I,L) =  WU(I,L) * AITR
./pkg/fizhi/fizhi_turb.F 2762 TRBFLX
        WV(I,L) =  WV(I,L) * AITR
./pkg/fizhi/fizhi_turb.F 326 TURBIO
      do L=1,nlay
./pkg/fizhi/fizhi_turb.F 328 TURBIO
      sum0 = sum0 + tke(n,L)
./pkg/fizhi/fizhi_turb.F 3861 TRBLEN
      integer i,L,LL
./pkg/fizhi/fizhi_turb.F 3877 TRBLEN
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 3879 TRBLEN
         STBFCN(I,L) = STRT(I,L) - RICR * DW2(I,L)
./pkg/fizhi/fizhi_turb.F 3880 TRBLEN
         INT1(I,L) = 0
./pkg/fizhi/fizhi_turb.F 3881 TRBLEN
         IF( STBFCN(I,L).LE.0. ) INT1(I,L) = 1
./pkg/fizhi/fizhi_turb.F 3889 TRBLEN
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3891 TRBLEN
        INT2(I,L) = 0
./pkg/fizhi/fizhi_turb.F 3892 TRBLEN
        IF( (INT1(I,L).EQ.1) .NEQV. (INT1(I,L+1).EQ.1) ) INT2(I,L) = 1
./pkg/fizhi/fizhi_turb.F 3910 TRBLEN
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3912 TRBLEN
        XL0(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 3936 TRBLEN
      DO L =LMIN1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3938 TRBLEN
        INT2(I,L) = 0
./pkg/fizhi/fizhi_turb.F 3939 TRBLEN
        IF ( INT1(I,L).EQ.1 .OR. INT1(I,L+1).EQ.1 ) INT2(I,L) = 1
./pkg/fizhi/fizhi_turb.F 3940 TRBLEN
        IF ( INT2(I,L).EQ.1 )
./pkg/fizhi/fizhi_turb.F 3941 TRBLEN
     .       XL0(I,L) = (0.5+DZITRP(I,L)) * DZ3(I,L+1)
./pkg/fizhi/fizhi_turb.F 3953 TRBLEN
      DO L =LMIN1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3955 TRBLEN
        IF ( INT2(I,L).EQ.1 ) THEN
./pkg/fizhi/fizhi_turb.F 3956 TRBLEN
         XL0(I,L+1) = XL0(I,L+1) + ( (0.5-DZITRP(I,L)) * DZ3(I,L+1) )
./pkg/fizhi/fizhi_turb.F 3968 TRBLEN
      DO 500 L = LMINP,NLEVM1
./pkg/fizhi/fizhi_turb.F 3969 TRBLEN
           LM1 = L-1
./pkg/fizhi/fizhi_turb.F 3971 TRBLEN
      IF( INT1(I,LM1).EQ.1 ) XL0(I,L) = XL0(I,L) + XL0(I,LM1)
./pkg/fizhi/fizhi_turb.F 3984 TRBLEN
              L   = NLEVM1 + LMIN - LL
./pkg/fizhi/fizhi_turb.F 3985 TRBLEN
              LP1 = L+1
./pkg/fizhi/fizhi_turb.F 3988 TRBLEN
      IF( INT1(I,L)  .EQ.1 ) THEN
./pkg/fizhi/fizhi_turb.F 3989 TRBLEN
      XL0(I,L) = XL0(I,LP1)
./pkg/fizhi/fizhi_turb.F 3991 TRBLEN
      XL0(I,L) = XL0(I,L) + XL0(I,LP1)
./pkg/fizhi/fizhi_turb.F 4001 TRBLEN
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 4003 TRBLEN
        IF ( XL0(I,L).LT.XL0CNV ) XL0(I,L) = XL0CNV
./pkg/fizhi/fizhi_turb.F 4018 TRBLEN
      DO L =1,LMINQ1
./pkg/fizhi/fizhi_turb.F 4020 TRBLEN
        INT2(I,L) = 1 - INT1(I,L)
./pkg/fizhi/fizhi_turb.F 4037 TRBLEN
       DO L =LMINQ,NLEVM1
./pkg/fizhi/fizhi_turb.F 4039 TRBLEN
         IF( INT1(I,L).EQ.0 ) THEN
./pkg/fizhi/fizhi_turb.F 4040 TRBLEN
           XL0(I,L) =         Q(I,L) / XL0(I,L)
./pkg/fizhi/fizhi_turb.F 4041 TRBLEN
           XL0(I,L) =       XL0(I,L) * XL0(I,L) + 1.0E-20
./pkg/fizhi/fizhi_turb.F 4042 TRBLEN
           XL0(I,L) =    STBFCN(I,L) + XL0(I,L)
./pkg/fizhi/fizhi_turb.F 4043 TRBLEN
           XL0(I,L) = SQRT( XL0(I,L) )
./pkg/fizhi/fizhi_turb.F 4044 TRBLEN
           XL0(I,L) =         Q(I,L) / XL0(I,L)
./pkg/fizhi/fizhi_turb.F 4046 TRBLEN
         INT2(I,L)  = 0
./pkg/fizhi/fizhi_turb.F 4047 TRBLEN
         IF( XL0(I,L).LT.XL0MIN ) INT2(I,L)  = 1
./pkg/fizhi/fizhi_turb.F 4058 TRBLEN
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 4060 TRBLEN
         INT2(I,L) = 1 - INT1(I,L)
./pkg/fizhi/fizhi_turb.F 4061 TRBLEN
         IF ( INT2(I,L).EQ.1 ) XL0(I,L) = XL0MIN
./pkg/fizhi/fizhi_turb.F 4068 TRBLEN
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 4070 TRBLEN
        IF ( INT2(I,L).EQ.1 ) XL0(I,L) = XL0MIN
./pkg/fizhi/fizhi_turb.F 4083 TRBLEN
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 4085 TRBLEN
        XL(I,L) = XL0(I,L) * VKZE(I,L) / ( XL0(I,L)+VKZE(I,L) )
./pkg/fizhi/fizhi_turb.F 4104 TRBLEN
      DO L =LMINQ1,NLEVM1
./pkg/fizhi/fizhi_turb.F 4106 TRBLEN
        Q1(I,L) = Q(I,L)
./pkg/fizhi/fizhi_turb.F 4107 TRBLEN
        INT1(I,L) = 0
./pkg/fizhi/fizhi_turb.F 4108 TRBLEN
        IF (    Q(I,L).LE.Q(I,L+1) ) INT1(I,L) = 1
./pkg/fizhi/fizhi_turb.F 4109 TRBLEN
        IF ( INT1(I,L).EQ.1 )  THEN
./pkg/fizhi/fizhi_turb.F 4110 TRBLEN
           XL0(I,L) = XL0(I,L+1)
./pkg/fizhi/fizhi_turb.F 4111 TRBLEN
            Q1(I,L) =   Q(I,L+1)
./pkg/fizhi/fizhi_turb.F 4121 TRBLEN
      DO L =LMINQ1,NLEVM1
./pkg/fizhi/fizhi_turb.F 4123 TRBLEN
        QXLM(I,L) =   XL0(I,L)*VKZM(I,L)
./pkg/fizhi/fizhi_turb.F 4124 TRBLEN
     .                 / ( XL0(I,L)+VKZM(I,L) )
./pkg/fizhi/fizhi_turb.F 4125 TRBLEN
        QXLM(I,L) = CLMT53 * Q1(I,L)*QXLM(I,L)
./pkg/fizhi/fizhi_turb.F 4180 TRBITP
      integer I,L
./pkg/fizhi/fizhi_turb.F 4191 TRBITP
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 4193 TRBITP
        DZITRP(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 4219 TRBITP
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 4221 TRBITP
        IF ( INTCHG(I,L).EQ.1 ) THEN
./pkg/fizhi/fizhi_turb.F 4222 TRBITP
         DDD(I,L) = (         CP  *(DTHV(I,L+1)*DPK(I,L)
./pkg/fizhi/fizhi_turb.F 4223 TRBITP
     .                            + DTHV(I,L)*DPK(I,L+1)) )
./pkg/fizhi/fizhi_turb.F 4224 TRBITP
     .            - ( (2.*RICR)  * ( DU(I,L+1)* DU(I,L)
./pkg/fizhi/fizhi_turb.F 4225 TRBITP
     .                            +   DV(I,L+1)* DV(I,L)) )
./pkg/fizhi/fizhi_turb.F 4226 TRBITP
         AAA(I,L) = STBFCN(I,L) + STBFCN(I,L+1)
./pkg/fizhi/fizhi_turb.F 4227 TRBITP
         BBB(I,L) = STBFCN(I,L) - STBFCN(I,L+1)
./pkg/fizhi/fizhi_turb.F 4228 TRBITP
         CCC(I,L) =            1. / BBB(I,L)
./pkg/fizhi/fizhi_turb.F 4229 TRBITP
         DZITRP(I,L) = AAA(I,L) * CCC(I,L)
./pkg/fizhi/fizhi_turb.F 4230 TRBITP
         AAA(I,L) = AAA(I,L) - DDD(I,L)
./pkg/fizhi/fizhi_turb.F 4231 TRBITP
         DDD(I,L) =        (   DDD(I,L) *    DDD(I,L) )
./pkg/fizhi/fizhi_turb.F 4232 TRBITP
     .            - 4. * (STBFCN(I,L+1) * STBFCN(I,L) )
./pkg/fizhi/fizhi_turb.F 4233 TRBITP
         DDD(I,L) = DDD(I,L)*CCC(I,L)*CCC(I,L)
./pkg/fizhi/fizhi_turb.F 4234 TRBITP
         DDD(I,L) =    SQRT( DDD(I,L) )
./pkg/fizhi/fizhi_turb.F 4237 TRBITP
        IF( INTCHG(I,L).EQ.1 .AND. AAA(I,L).NE.0. ) THEN
./pkg/fizhi/fizhi_turb.F 4238 TRBITP
         DZITRP(I,L) = ( BBB(I,L)*(1.-DDD(I,L)) ) / AAA(I,L)
./pkg/fizhi/fizhi_turb.F 4241 TRBITP
        DZITRP(I,L) = 0.5 * DZITRP(I,L)
./pkg/fizhi/fizhi_turb.F 4496 TRBDIF
      integer I,L
./pkg/fizhi/fizhi_turb.F 4514 TRBDIF
      DO L =1,NLEVX
./pkg/fizhi/fizhi_turb.F 4516 TRBDIF
        CC(I,L+1) = RHOKDZ(I,L) * FLXFAC(I,L+1)
./pkg/fizhi/fizhi_turb.F 4523 TRBDIF
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 4525 TRBDIF
        BB(I,L) =   RHOKDZ(I,L) * FLXFAC(I,L)
./pkg/fizhi/fizhi_turb.F 4526 TRBDIF
        AA(I,L) = 1. + CC(I,L) + BB(I,L)
./pkg/fizhi/fizhi_turb.F 4535 TRBDIF
       DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 4537 TRBDIF
         AA(I,L) = AA(I,L) - XX2(I,L)
./pkg/fizhi/fizhi_turb.F 4568 TRBDIF
       DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 4570 TRBDIF
         IF ( XX1(I,L).LT.EPSL ) XX1(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 4584 VTRI0
      integer i,L,Lm1
./pkg/fizhi/fizhi_turb.F 4590 VTRI0
      DO 100 L = 2,K
./pkg/fizhi/fizhi_turb.F 4591 VTRI0
       LM1 = L - 1
./pkg/fizhi/fizhi_turb.F 4593 VTRI0
        C(I,L) = C(I,L) * A(I,LM1)
./pkg/fizhi/fizhi_turb.F 4594 VTRI0
        A(I,L) = 1. / ( A(I,L) - B(I,LM1) * C(I,L) )
./pkg/fizhi/fizhi_turb.F 4595 VTRI0
        F(I,L) = F(I,L) + F(I,LM1) * C(I,L)
./pkg/fizhi/fizhi_turb.F 4599 VTRI0
      DO 200 L = K,1,-1
./pkg/fizhi/fizhi_turb.F 4601 VTRI0
       Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) * A(I,L)
./pkg/fizhi/fizhi_turb.F 4614 VTRI1
      integer i,L
./pkg/fizhi/fizhi_turb.F 4616 VTRI1
      DO 200 L = K,1,-1
./pkg/fizhi/fizhi_turb.F 4618 VTRI1
        Y(I,L) = B(I,L) * Y(I,L+1) * A(I,L)
./pkg/fizhi/fizhi_turb.F 4632 VTRI2
      integer i,L
./pkg/fizhi/fizhi_turb.F 4634 VTRI2
      DO 100 L = 2,K
./pkg/fizhi/fizhi_turb.F 4636 VTRI2
        F(I,L) = F(I,L) + F(I,L-1) * C(I,L)
./pkg/fizhi/fizhi_turb.F 4640 VTRI2
      DO 200 L = K,1,-1
./pkg/fizhi/fizhi_turb.F 4642 VTRI2
        Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) *  A(I,L)
./pkg/fizhi/fizhi_turb.F 603 TURBIO
      do L = 1,nlay+1
./pkg/fizhi/fizhi_turb.F 605 TURBIO
      thv(i,L) = 1. + virtcon * sh(i,L)
./pkg/fizhi/fizhi_turb.F 606 TURBIO
      thv(i,L) = th(i,L) * thv(i,L)
./pkg/fizhi/fizhi_turb.F 614 TURBIO
      do L =1,nlay
./pkg/fizhi/fizhi_turb.F 616 TURBIO
      qliq(i,L) = 0.
./pkg/fizhi/fizhi_turb.F 617 TURBIO
      turbfcc(i,L) = 0.
./pkg/fizhi/fizhi_turb.F 946 TURBIO
      do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 948 TURBIO
      th(i,L) = th(i,L) + dthdthg(i,L)*(tc(i)-savetc(i))/pke(i,nlay+1)
./pkg/fizhi/fizhi_turb.F 955 TURBIO
      do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 957 TURBIO
      sh(i,L) = sh(i,L) + dshdshg(i,L)*(qa(i)-saveqa(i))
./pkg/fizhi/fizhi_turb.F 961 TURBIO
      do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 963 TURBIO
       sttflux(i,L) = transth(i,L) * (th(i,L+1)-th(i,L))
./pkg/fizhi/fizhi_turb.F 964 TURBIO
       stqflux(i,L) = transsh(i,L) * (sh(i,L+1)-sh(i,L))
./pkg/fizhi/fizhi_utils.F 101 STRIPINT
        DO 100 J=1,L
./pkg/fizhi/fizhi_utils.F 106 STRIPINT
        DO 200 J=1,L
./pkg/fizhi/fizhi_utils.F 118 PASTE
      SUBROUTINE PASTE(B,A,IB,IA,L,K)
./pkg/fizhi/fizhi_utils.F 120 PASTE
      integer ia,ib,L,K
./pkg/fizhi/fizhi_utils.F 121 PASTE
      _RL A(IA,L), B(IB,L)
./pkg/fizhi/fizhi_utils.F 129 PASTE
      DO 100 J=1,L
./pkg/fizhi/fizhi_utils.F 136 PSTBMP
      SUBROUTINE PSTBMP(B,A,IB,IA,L,K)
./pkg/fizhi/fizhi_utils.F 138 PSTBMP
      integer ia,ib,L,K
./pkg/fizhi/fizhi_utils.F 139 PSTBMP
      _RL A(IA,L), B(IB,L)
./pkg/fizhi/fizhi_utils.F 147 PSTBMP
      DO 100 J=1,L
./pkg/fizhi/fizhi_utils.F 154 STRINT
      SUBROUTINE STRINT(A,B,IA,IB,L,K)
./pkg/fizhi/fizhi_utils.F 156 STRINT
      integer ia,ib,L,K
./pkg/fizhi/fizhi_utils.F 157 STRINT
      INTEGER A(IA,L), B(IB,L)
./pkg/fizhi/fizhi_utils.F 166 STRINT
        DO 100 J=1,L
./pkg/fizhi/fizhi_utils.F 171 STRINT
        DO 200 J=1,L
./pkg/fizhi/fizhi_utils.F 446 PASTIT
      subroutine PASTIT(b,a,ib,ia,irun,L,k)
./pkg/fizhi/fizhi_utils.F 448 PASTIT
      integer ib,ia,L,k,irun,Lena,offset
./pkg/fizhi/fizhi_utils.F 456 PASTIT
      do 100 j=1,L
./pkg/fizhi/fizhi_utils.F 465 PSTBITINT
      integer ib,ia,L,k,irun,Lena,offset
./pkg/fizhi/fizhi_utils.F 474 PSTBITINT
      do 100 j=1,L
./pkg/fizhi/fizhi_utils.F 484 PSTBMPIT
      integer ib,ia,L,k,irun,Lena,offset
./pkg/fizhi/fizhi_utils.F 492 PSTBMPIT
      do 100 j=1,L
./pkg/fizhi/fizhi_utils.F 577 PASTE2GRD_OLD
      integer i,L,offset,Lena
./pkg/fizhi/fizhi_utils.F 583 PASTE2GRD_OLD
      do L = 1,levs
./pkg/fizhi/fizhi_utils.F 585 PASTE2GRD_OLD
       a(indx(i+offset-1),L) = a(indx(i+offset-1),L) + b(i,L)*chfr(i)
./pkg/fizhi/fizhi_utils.F 60 STRIP
      SUBROUTINE STRIP(A,B,IA,IB,L,K)
./pkg/fizhi/fizhi_utils.F 624 PASTE2GRD
      integer i,L,offset,Lena
./pkg/fizhi/fizhi_utils.F 62 STRIP
      integer ia,ib,L,K
./pkg/fizhi/fizhi_utils.F 633 PASTE2GRD
       do L= 1,levs
./pkg/fizhi/fizhi_utils.F 635 PASTE2GRD
        if( a(indx(i+offset-1),L).eq.undef .or. b(i,L).eq.undef ) then
./pkg/fizhi/fizhi_utils.F 636 PASTE2GRD
         a(indx(i+offset-1),L) = undef
./pkg/fizhi/fizhi_utils.F 638 PASTE2GRD
         a(indx(i+offset-1),L)=a(indx(i+offset-1),L) + b(i,L)*chfr(i)
./pkg/fizhi/fizhi_utils.F 63 STRIP
      _RL A(IA,L), B(IB,L)
./pkg/fizhi/fizhi_utils.F 643 PASTE2GRD
       do L= 1,levs
./pkg/fizhi/fizhi_utils.F 645 PASTE2GRD
        a(indx(i+offset-1),L)=a(indx(i+offset-1),L) + b(i,L)*chfr(i)
./pkg/fizhi/fizhi_utils.F 72 STRIP
        DO 100 J=1,L
./pkg/fizhi/fizhi_utils.F 77 STRIP
        DO 200 J=1,L
./pkg/fizhi/fizhi_utils.F 89 STRIPINT
      SUBROUTINE STRIPINT(A,B,IA,IB,L,K)
./pkg/fizhi/fizhi_utils.F 91 STRIPINT
      integer ia,ib,L,K
./pkg/fizhi/fizhi_utils.F 925 PKAPPA
      integer i,j,L
./pkg/fizhi/fizhi_utils.F 929 PKAPPA
      do L = 1,lm
./pkg/fizhi/fizhi_utils.F 92 STRIPINT
      INTEGER A(IA,L), B(IB,L)
./pkg/fizhi/fizhi_utils.F 932 PKAPPA
       pkz(i,j,L) = ( ple (i,j,l+1)*pkle(i,j,l+1)
./pkg/fizhi/fizhi_wrapper.F 101 FIZHI_WRAPPER
        do L = 1,Nr
./pkg/fizhi/fizhi_wrapper.F 104 FIZHI_WRAPPER
         vdynbef(i,j,L,bi,bj) = vvel(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 115 FIZHI_WRAPPER
        do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 118 FIZHI_WRAPPER
         thbef(i,j,L,bi,bj) = thphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 122 FIZHI_WRAPPER
        do L = 1,Nr
./pkg/fizhi/fizhi_wrapper.F 125 FIZHI_WRAPPER
         thdynbef(i,j,L,bi,bj) = theta(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 136 FIZHI_WRAPPER
        do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 139 FIZHI_WRAPPER
         sbef(i,j,L,bi,bj) = sphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 143 FIZHI_WRAPPER
        do L = 1,Nr
./pkg/fizhi/fizhi_wrapper.F 146 FIZHI_WRAPPER
         sdynbef(i,j,L,bi,bj) = salt(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 165 FIZHI_WRAPPER
         do L = 2,Nrphys+1
./pkg/fizhi/fizhi_wrapper.F 166 FIZHI_WRAPPER
          pephy(i,j,L,bi,bj)=pephy(i,j,L-1,bi,bj)-dpphys(i,j,L-1,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 169 FIZHI_WRAPPER
         do L = 1,Nrphys+1
./pkg/fizhi/fizhi_wrapper.F 170 FIZHI_WRAPPER
          pephy4fiz(i,j,Nrphys+2-L,bi,bj)=pephy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 179 FIZHI_WRAPPER
         do L = 1,Nr
./pkg/fizhi/fizhi_wrapper.F 180 FIZHI_WRAPPER
          pedyn(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_wrapper.F 194 FIZHI_WRAPPER
         do L = Lbotij+1,Nr+1
./pkg/fizhi/fizhi_wrapper.F 195 FIZHI_WRAPPER
          pedyn(i,j,L,bi,bj) = pedyn(i,j,L-1,bi,bj) -
./pkg/fizhi/fizhi_wrapper.F 196 FIZHI_WRAPPER
     .           drF(L-1)* rStarExpC(i,j,bi,bj)*hfacC(i,j,L-1,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 213 FIZHI_WRAPPER
        do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 214 FIZHI_WRAPPER
         tempphy(i,j,L,bi,bj) = pephy4fiz(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 235 FIZHI_WRAPPER
        L = Nrphys+1
./pkg/fizhi/fizhi_wrapper.F 239 FIZHI_WRAPPER
         tempLdiag(i,j,L) = pephy4fiz(i,j,L,bi,bj)/100.
./pkg/fizhi/fizhi_wrapper.F 242 FIZHI_WRAPPER
        do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 245 FIZHI_WRAPPER
         tempLdiag(i,j,L) = pephy4fiz(i,j,L,bi,bj)/100.
./pkg/fizhi/fizhi_wrapper.F 246 FIZHI_WRAPPER
         tempLdiag2(i,j,L) = thphy(i,j,L,bi,bj) *
./pkg/fizhi/fizhi_wrapper.F 247 FIZHI_WRAPPER
     .        (1.+0.609*sphy(i,j,L,bi,bj))
./pkg/fizhi/fizhi_wrapper.F 288 FIZHI_WRAPPER
         do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 289 FIZHI_WRAPPER
          tempphy(i,j,Nrphys+1-L,bi,bj)=duphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 297 FIZHI_WRAPPER
         do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 298 FIZHI_WRAPPER
          tempphy(i,j,Nrphys+1-L,bi,bj)=dvphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 306 FIZHI_WRAPPER
         do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 307 FIZHI_WRAPPER
          tempphy(i,j,Nrphys+1-L,bi,bj)=dthphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 315 FIZHI_WRAPPER
         do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 316 FIZHI_WRAPPER
          tempphy(i,j,Nrphys+1-L,bi,bj)=dsphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 53 FIZHI_WRAPPER
       integer i, j, L, Lbotij, bi, bj
./pkg/fizhi/fizhi_wrapper.F 73 FIZHI_WRAPPER
        do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 76 FIZHI_WRAPPER
         ubef(i,j,L,bi,bj) = uphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 80 FIZHI_WRAPPER
        do L = 1,Nr
./pkg/fizhi/fizhi_wrapper.F 83 FIZHI_WRAPPER
         udynbef(i,j,L,bi,bj) = uvel(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 94 FIZHI_WRAPPER
        do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 97 FIZHI_WRAPPER
         vbef(i,j,L,bi,bj) = vphy(i,j,L,bi,bj)
./pkg/fizhi/getpwhere.F 105 GETPWHERE
       WRITE(ioUnit,'(10F10.4)') (pref(L),L=1,ntry)
./pkg/fizhi/getpwhere.F 112 GETPWHERE
       DO L = ntry+1,Nrphys
./pkg/fizhi/getpwhere.F 113 GETPWHERE
        pref(L) = (Nrphys-L+0.5)*dpTop
./pkg/fizhi/getpwhere.F 116 GETPWHERE
       DO L = 1,Nr
./pkg/fizhi/getpwhere.F 117 GETPWHERE
        rC_dyn(L) = rC(L)*1. _d -2
./pkg/fizhi/getpwhere.F 121 GETPWHERE
       DO L = 1,Nr
./pkg/fizhi/getpwhere.F 122 GETPWHERE
         IF ( rC_dyn(L).GE.dptry_pedge(ntry+1) ) dynlev = L
./pkg/fizhi/getpwhere.F 124 GETPWHERE
       DO L = ntry+1,MIN(Nrphys,ntry+Nr-dynlev)
./pkg/fizhi/getpwhere.F 125 GETPWHERE
         pref(L) = rC_dyn(dynlev+L-ntry)
./pkg/fizhi/getpwhere.F 130 GETPWHERE
       do L = 2,Nr+1
./pkg/fizhi/getpwhere.F 131 GETPWHERE
        rF_edge(L) = rF_edge(L-1) - (drF(L-1)/100.)
./pkg/fizhi/getpwhere.F 133 GETPWHERE
       do L = 1,Nr
./pkg/fizhi/getpwhere.F 134 GETPWHERE
        rF_pmid(L) = (rF_edge(L) + rF_edge(L+1))/2.
./pkg/fizhi/getpwhere.F 138 GETPWHERE
       do L = 1,Nr
./pkg/fizhi/getpwhere.F 139 GETPWHERE
        if(rF_pmid(L).ge.pref(ntry)) dynlev = L
./pkg/fizhi/getpwhere.F 157 GETPWHERE
        do L = ntry+1,Nrphys-1
./pkg/fizhi/getpwhere.F 158 GETPWHERE
         pref(L) = rF_pmid(dynlev+L-ntry)
./pkg/fizhi/getpwhere.F 162 GETPWHERE
        do L = ntry+1,Nrphys-1
./pkg/fizhi/getpwhere.F 163 GETPWHERE
         pref(L) = rF_pmid(dynlev+L-ntry+1)
./pkg/fizhi/getpwhere.F 172 GETPWHERE
        DO L = 1,Nrphys
./pkg/fizhi/getpwhere.F 173 GETPWHERE
         IF ( pref(L).GE.pressures(n) ) tmplev = L
./pkg/fizhi/getpwhere.F 186 GETPWHERE
       WRITE(ioUnit,'(10F10.4)') (pref(L),L=ntry,Nrphys)
./pkg/fizhi/getpwhere.F 49 GETPWHERE
       integer n,L,dynlev
./pkg/fizhi/getpwhere.F 73 GETPWHERE
        do L = 1,ntry
./pkg/fizhi/getpwhere.F 74 GETPWHERE
         dptry(L) = dptry10(L)
./pkg/fizhi/getpwhere.F 78 GETPWHERE
        do L = 1,ntry
./pkg/fizhi/getpwhere.F 79 GETPWHERE
         dptry(L) = dptry40(L)
./pkg/fizhi/getpwhere.F 92 GETPWHERE
       do L = 1,ntry
./pkg/fizhi/getpwhere.F 93 GETPWHERE
        dptry_pedge(L+1) = dptry_pedge(L) - dptry(L)
./pkg/fizhi/getpwhere.F 95 GETPWHERE
       do L = 1,ntry
./pkg/fizhi/getpwhere.F 96 GETPWHERE
        pref(L) = (dptry_pedge(L) + dptry_pedge(L+1))/2.
./pkg/fizhi/slprs.F 32 SLPRS
      integer i,L
./pkg/fizhi/slprs.F 48 SLPRS
      do L = lm,1,-1
./pkg/fizhi/slprs.F 50 SLPRS
        if ( ple(i,L+1).ge.(ple(i,lm+1)-100.) ) then
./pkg/fizhi/slprs.F 51 SLPRS
         Ltop(i) = L
./pkg/fizhi/slprs.F 52 SLPRS
         tm(i) = tm(i) + thz(i,L)*(ple(i,L+1)-ple(i,L))
./pkg/fizhi/step_fizhi_corr.F 109 STEP_FIZHI_CORR
         do L = 2,Nrphys+1
./pkg/fizhi/step_fizhi_corr.F 110 STEP_FIZHI_CORR
          pephy(i,j,L,bi,bj)=pephy(i,j,L-1,bi,bj)-dpphys(i,j,L-1,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 120 STEP_FIZHI_CORR
        do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 123 STEP_FIZHI_CORR
         windphy(i,j,L,bi,bj) = 
./pkg/fizhi/step_fizhi_corr.F 124 STEP_FIZHI_CORR
     .     sqrt(uphy(i,j,Nrphys+1-L,bi,bj)*uphy(i,j,Nrphys+1-L,bi,bj) 
./pkg/fizhi/step_fizhi_corr.F 125 STEP_FIZHI_CORR
     .        + vphy(i,j,Nrphys+1-L,bi,bj)*vphy(i,j,Nrphys+1-L,bi,bj))
./pkg/fizhi/step_fizhi_corr.F 139 STEP_FIZHI_CORR
        do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 142 STEP_FIZHI_CORR
         tempphy(i,j,Nrphys+1-L,bi,bj) = uphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 148 STEP_FIZHI_CORR
        do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 151 STEP_FIZHI_CORR
         tempphy(i,j,Nrphys+1-L,bi,bj) = vphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 157 STEP_FIZHI_CORR
        do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 160 STEP_FIZHI_CORR
         tempphy(i,j,Nrphys+1-L,bi,bj) = thphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 166 STEP_FIZHI_CORR
        do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 169 STEP_FIZHI_CORR
         tempphy(i,j,Nrphys+1-L,bi,bj) = sphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 191 STEP_FIZHI_CORR
        do L = 1,Nr
./pkg/fizhi/step_fizhi_corr.F 194 STEP_FIZHI_CORR
        udyntemp(i,j,L,bi,bj)=uvel(i,j,L,bi,bj)-udyntemp(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 195 STEP_FIZHI_CORR
        vdyntemp(i,j,L,bi,bj)=vvel(i,j,L,bi,bj)-vdyntemp(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 196 STEP_FIZHI_CORR
        thdyntemp(i,j,L,bi,bj)=theta(i,j,L,bi,bj)-thdyntemp(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 197 STEP_FIZHI_CORR
        sdyntemp(i,j,L,bi,bj)=salt(i,j,L,bi,bj)-sdyntemp(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 219 STEP_FIZHI_CORR
        do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 222 STEP_FIZHI_CORR
         uphytemp(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 228 STEP_FIZHI_CORR
        do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 231 STEP_FIZHI_CORR
         vphytemp(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 237 STEP_FIZHI_CORR
        do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 240 STEP_FIZHI_CORR
         thphytemp(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 246 STEP_FIZHI_CORR
        do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 249 STEP_FIZHI_CORR
         sphytemp(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 280 STEP_FIZHI_CORR
      do L=1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 287 STEP_FIZHI_CORR
         tempij(i,j) = (uphy (i,j,L,bi,bj)-ubef(i,j,L,bi,bj) )
./pkg/fizhi/step_fizhi_corr.F 291 STEP_FIZHI_CORR
        call DIAGNOSTICS_FILL(tempij,'TENDUFIZ',L,1,3,bi,bj,myThid)
./pkg/fizhi/step_fizhi_corr.F 299 STEP_FIZHI_CORR
         tempij(i,j) = (vphy (i,j,L,bi,bj)-vbef(i,j,L,bi,bj) )
./pkg/fizhi/step_fizhi_corr.F 303 STEP_FIZHI_CORR
        call DIAGNOSTICS_FILL(tempij,'TENDVFIZ',L,1,3,bi,bj,myThid)
./pkg/fizhi/step_fizhi_corr.F 311 STEP_FIZHI_CORR
         tempij(i,j) = (thphy (i,j,L,bi,bj)-thbef(i,j,L,bi,bj) )
./pkg/fizhi/step_fizhi_corr.F 315 STEP_FIZHI_CORR
        call DIAGNOSTICS_FILL(tempij,'TENDTFIZ',L,1,3,bi,bj,myThid)
./pkg/fizhi/step_fizhi_corr.F 323 STEP_FIZHI_CORR
         tempij(i,j) = (sphy (i,j,L,bi,bj)-sbef(i,j,L,bi,bj) )
./pkg/fizhi/step_fizhi_corr.F 327 STEP_FIZHI_CORR
        call DIAGNOSTICS_FILL(tempij,'TENDQFIZ',L,1,3,bi,bj,myThid)
./pkg/fizhi/step_fizhi_corr.F 343 STEP_FIZHI_CORR
       do L=1,Nr
./pkg/fizhi/step_fizhi_corr.F 346 STEP_FIZHI_CORR
        udyntemp(i,j,L,bi,bj) = 
./pkg/fizhi/step_fizhi_corr.F 347 STEP_FIZHI_CORR
     .    (uvel(i,j,L,bi,bj)-udynbef(i,j,L,bi,bj))*dtinv - 
./pkg/fizhi/step_fizhi_corr.F 348 STEP_FIZHI_CORR
     .                                        guphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 349 STEP_FIZHI_CORR
        vdyntemp(i,j,L,bi,bj) = 
./pkg/fizhi/step_fizhi_corr.F 350 STEP_FIZHI_CORR
     .    (vvel(i,j,L,bi,bj)-vdynbef(i,j,L,bi,bj))*dtinv - 
./pkg/fizhi/step_fizhi_corr.F 351 STEP_FIZHI_CORR
     .                                        gvphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 362 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 365 STEP_FIZHI_CORR
        windphy(i,j,L,bi,bj) = 
./pkg/fizhi/step_fizhi_corr.F 366 STEP_FIZHI_CORR
     .     sqrt(uphy(i,j,Nrphys+1-L,bi,bj)*uphy(i,j,Nrphys+1-L,bi,bj) 
./pkg/fizhi/step_fizhi_corr.F 367 STEP_FIZHI_CORR
     .        + vphy(i,j,Nrphys+1-L,bi,bj)*vphy(i,j,Nrphys+1-L,bi,bj))
./pkg/fizhi/step_fizhi_corr.F 376 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 379 STEP_FIZHI_CORR
         uphytemp(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 385 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 388 STEP_FIZHI_CORR
         vphytemp(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 393 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 396 STEP_FIZHI_CORR
        uphytemp(i,j,L,bi,bj)=
./pkg/fizhi/step_fizhi_corr.F 397 STEP_FIZHI_CORR
     .    ( (uphy(i,j,L,bi,bj)-ubef(i,j,L,bi,bj))*dtinv
./pkg/fizhi/step_fizhi_corr.F 398 STEP_FIZHI_CORR
     .     - duphy(i,j,L,bi,bj) - uphytemp(i,j,L,bi,bj) ) * 86400. _d 0
./pkg/fizhi/step_fizhi_corr.F 399 STEP_FIZHI_CORR
        vphytemp(i,j,L,bi,bj)=
./pkg/fizhi/step_fizhi_corr.F 400 STEP_FIZHI_CORR
     .    ( (vphy(i,j,L,bi,bj)-vbef(i,j,L,bi,bj))*dtinv
./pkg/fizhi/step_fizhi_corr.F 401 STEP_FIZHI_CORR
     .     - dvphy(i,j,L,bi,bj) - vphytemp(i,j,L,bi,bj) ) * 86400. _d 0
./pkg/fizhi/step_fizhi_corr.F 425 STEP_FIZHI_CORR
       do L=1,Nr
./pkg/fizhi/step_fizhi_corr.F 428 STEP_FIZHI_CORR
        thdyntemp(i,j,L,bi,bj) = 
./pkg/fizhi/step_fizhi_corr.F 429 STEP_FIZHI_CORR
     .    (theta(i,j,L,bi,bj)-thdynbef(i,j,L,bi,bj))*dtinv - 
./pkg/fizhi/step_fizhi_corr.F 430 STEP_FIZHI_CORR
     .                                        gthphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 440 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 443 STEP_FIZHI_CORR
         thphytemp(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 448 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 451 STEP_FIZHI_CORR
        thphytemp(i,j,L,bi,bj)=
./pkg/fizhi/step_fizhi_corr.F 452 STEP_FIZHI_CORR
     .    ( (thphy(i,j,L,bi,bj)-thbef(i,j,L,bi,bj))*dtinv
./pkg/fizhi/step_fizhi_corr.F 453 STEP_FIZHI_CORR
     .     - dthphy(i,j,L,bi,bj) - thphytemp(i,j,L,bi,bj) ) * 86400. _d 0
./pkg/fizhi/step_fizhi_corr.F 471 STEP_FIZHI_CORR
       do L=1,Nr
./pkg/fizhi/step_fizhi_corr.F 474 STEP_FIZHI_CORR
        sdyntemp(i,j,L,bi,bj) = 
./pkg/fizhi/step_fizhi_corr.F 475 STEP_FIZHI_CORR
     .    (salt(i,j,L,bi,bj)-sdynbef(i,j,L,bi,bj))*dtinv - 
./pkg/fizhi/step_fizhi_corr.F 476 STEP_FIZHI_CORR
     .                                        gsphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 486 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 489 STEP_FIZHI_CORR
         sphytemp(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 494 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 497 STEP_FIZHI_CORR
        sphytemp(i,j,L,bi,bj)=
./pkg/fizhi/step_fizhi_corr.F 498 STEP_FIZHI_CORR
     .    ( (sphy(i,j,L,bi,bj)-sbef(i,j,L,bi,bj))*dtinv
./pkg/fizhi/step_fizhi_corr.F 499 STEP_FIZHI_CORR
     .     - dsphy(i,j,L,bi,bj) - sphytemp(i,j,L,bi,bj) ) * 86400. _d 0
./pkg/fizhi/step_fizhi_corr.F 56 STEP_FIZHI_CORR
       integer i, j, L, Lbotij, bi, bj
./pkg/fizhi/step_fizhi_corr.F 81 STEP_FIZHI_CORR
         do L = 1,Nr
./pkg/fizhi/step_fizhi_corr.F 82 STEP_FIZHI_CORR
          pedyn(i,j,L,bi,bj) = 0.
./pkg/fizhi/step_fizhi_corr.F 96 STEP_FIZHI_CORR
         do L = Lbotij+1,Nr+1
./pkg/fizhi/step_fizhi_corr.F 97 STEP_FIZHI_CORR
          pedyn(i,j,L,bi,bj) = pedyn(i,j,L-1,bi,bj) -
./pkg/fizhi/step_fizhi_corr.F 98 STEP_FIZHI_CORR
     .               drF(L-1)* rStarExpC(i,j,bi,bj)*hfacC(i,j,L-1,bi,bj)
./pkg/fizhi/step_fizhi_fg.F 32 STEP_FIZHI_FG
       integer i,j,L
./pkg/fizhi/step_fizhi_fg.F 53 STEP_FIZHI_FG
      do L = 1,Nrphys
./pkg/fizhi/step_fizhi_fg.F 56 STEP_FIZHI_FG
        tempij(i,j) = uphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_fg.F 61 STEP_FIZHI_FG
      do L = 1,Nrphys
./pkg/fizhi/step_fizhi_fg.F 64 STEP_FIZHI_FG
        tempij(i,j) = vphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_fg.F 69 STEP_FIZHI_FG
      do L = 1,Nrphys
./pkg/fizhi/step_fizhi_fg.F 72 STEP_FIZHI_FG
        tempij(i,j) = thphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_fg.F 77 STEP_FIZHI_FG
      do L = 1,Nrphys
./pkg/fizhi/step_fizhi_fg.F 80 STEP_FIZHI_FG
        tempij(i,j) = sphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_fg.F 92 STEP_FIZHI_FG
      do L = 1,Nrphys
./pkg/fizhi/step_fizhi_fg.F 95 STEP_FIZHI_FG
        tempij(i,j) = sphy(i,j,L,bi,bj)
./pkg/fizhi/step_physics.F 25 STEP_PHYSICS
       integer i,j,L
./pkg/fizhi/step_physics.F 27 STEP_PHYSICS
       do L = 1,Nrphys
./pkg/fizhi/step_physics.F 30 STEP_PHYSICS
         uphy(i,j,L,bi,bj) = uphy(i,j,L,bi,bj) + duphy(i,j,L,bi,bj) * dt
./pkg/fizhi/step_physics.F 31 STEP_PHYSICS
         vphy(i,j,L,bi,bj) = vphy(i,j,L,bi,bj) + dvphy(i,j,L,bi,bj) * dt
./pkg/fizhi/step_physics.F 32 STEP_PHYSICS
         thphy(i,j,L,bi,bj)=thphy(i,j,L,bi,bj) + dthphy(i,j,L,bi,bj)* dt
./pkg/fizhi/step_physics.F 33 STEP_PHYSICS
         sphy(i,j,L,bi,bj) = sphy(i,j,L,bi,bj) + dsphy(i,j,L,bi,bj) * dt
./pkg/fizhi/update_chemistry_exports.F 35 UPDATE_CHEMISTRY_EXPORTS
      integer i, j, L, LL, bi, bj
./pkg/fizhi/update_chemistry_exports.F 57 UPDATE_CHEMISTRY_EXPORTS
         do L = 2,Nrphys+1
./pkg/fizhi/update_chemistry_exports.F 58 UPDATE_CHEMISTRY_EXPORTS
         LL = Nrphys+2-L
./pkg/fizhi/update_chemistry_exports.F 59 UPDATE_CHEMISTRY_EXPORTS
         pephy(i,j,LL,bi,bj)=pephy(i,j,LL+1,bi,bj)-dpphys(i,j,L-1,bi,bj)
./pkg/fizhi/update_chemistry_exports.F 65 UPDATE_CHEMISTRY_EXPORTS
         do L = 1,Nrphys
./pkg/fizhi/update_chemistry_exports.F 66 UPDATE_CHEMISTRY_EXPORTS
          pphy(i,j,L,bi,bj)=(pephy(i,j,L+1,bi,bj)+pephy(i,j,L,bi,bj))
./pkg/fizhi/update_chemistry_exports.F 75 UPDATE_CHEMISTRY_EXPORTS
         do L = 1,Nrphys
./pkg/fizhi/update_chemistry_exports.F 76 UPDATE_CHEMISTRY_EXPORTS
          waterin(i,j,L) = sphy(i,j,L,bi,bj)
./pkg/fizhi/update_chemistry_exports.F 84 UPDATE_CHEMISTRY_EXPORTS
        do L = 1,nlevsoz
./pkg/fizhi/update_chemistry_exports.F 86 UPDATE_CHEMISTRY_EXPORTS
         oz1(j,L) = ozone(j,L,imns)*facm + ozone(j,L,ipls)*facp
./pkg/fizhi/update_chemistry_exports.F 90 UPDATE_CHEMISTRY_EXPORTS
        do L = 1,nlevsq
./pkg/fizhi/update_chemistry_exports.F 92 UPDATE_CHEMISTRY_EXPORTS
         strq1(j,L) = stratq(j,L,imns)*facm + stratq(j,L,ipls)*facp
./pkg/generic_advdiff/salt_fill.F 62 SALT_FILL
      INTEGER bi,bj,i,j,L,LM1
./pkg/generic_advdiff/salt_fill.F 85 SALT_FILL
        do L=Nr,2,-1
./pkg/generic_advdiff/salt_fill.F 86 SALT_FILL
         LM1 = L-1
./pkg/generic_advdiff/salt_fill.F 87 SALT_FILL
         dpratio= rC(L)/rC(LM1)
./pkg/generic_advdiff/salt_fill.F 90 SALT_FILL
          if( salt(i,j,L,bi,bj).lt.0.0  _d 0) then
./pkg/generic_advdiff/salt_fill.F 92 SALT_FILL
     .          salt(i,j,L,bi,bj)*dpratio
./pkg/generic_advdiff/salt_fill.F 93 SALT_FILL
           salt(i,j,L,bi,bj) = 0.0 _d 0
./pkg/gridalt/dyn2phys.F 102 DYN2PHYS
          dpkedyn = pedyn(i,j,L,bi,bj)-pedyn(i,j,L+1,bi,bj)
./pkg/gridalt/dyn2phys.F 117 DYN2PHYS
           if( (flg.eq.2) .and. (pedyn(i,j,L,bi,bj).lt.10000.)) then
./pkg/gridalt/dyn2phys.F 119 DYN2PHYS
     .   (qd-5. + (10.*(pedyn(i,j,L,bi,bj)-pphy(i,j,Lphy))/dpkedyn))/qd
./pkg/gridalt/dyn2phys.F 120 DYN2PHYS
           elseif( (flg.eq.2) .and. (pedyn(i,j,L,bi,bj).ge.10000.)) then
./pkg/gridalt/dyn2phys.F 59 DYN2PHYS
      integer  i,j,L,Lout1,Lout2,Lphy
./pkg/gridalt/dyn2phys.F 86 DYN2PHYS
      do L = 1,lmdyn
./pkg/gridalt/dyn2phys.F 90 DYN2PHYS
         qd = qdyn(i,j,L,bi,bj)
./pkg/gridalt/dyn2phys.F 92 DYN2PHYS
         if(L.ge.Lbot(i,j,bi,bj))then
./pkg/gridalt/dyn2phys.F 93 DYN2PHYS
          if(L.eq.Lbot(i,j,bi,bj)) then
./pkg/gridalt/dyn2phys.F 96 DYN2PHYS
           Lout1 = nlperdyn(i,j,L-1,bi,bj)
./pkg/gridalt/dyn2phys.F 98 DYN2PHYS
          Lout2 = nlperdyn(i,j,L,bi,bj)
./pkg/gridalt/gridalt_update.F 31 GRIDALT_UPDATE
       integer i,j,L
./pkg/gridalt/gridalt_update.F 52 GRIDALT_UPDATE
         do L = 1,Nrphys
./pkg/gridalt/gridalt_update.F 55 GRIDALT_UPDATE
           dpphys(i,j,L,bi,bj) = dpphys0(i,j,L,bi,bj) *
./pkg/gridalt/gridalt_update.F 73 GRIDALT_UPDATE
         do L = 1,Nrphys
./pkg/gridalt/gridalt_update.F 76 GRIDALT_UPDATE
           temparr(i,j,L,bi,bj) = dpphys(i,j,L,bi,bj) * 1. _d -2
./pkg/gridalt/make_phys_grid.F 108 MAKE_PHYS_GRID
       do L = 1,Nr
./pkg/gridalt/make_phys_grid.F 109 MAKE_PHYS_GRID
        dpstar_accum = dpstar_accum + drF(L)
./pkg/gridalt/make_phys_grid.F 171 MAKE_PHYS_GRID
        do L = Lbotij+1,Nr
./pkg/gridalt/make_phys_grid.F 177 MAKE_PHYS_GRID
           if((hfacC(i,j,L,bi,bj)*drF(L)).ge.deltap+dptry(Lnew))then
./pkg/gridalt/make_phys_grid.F 186 MAKE_PHYS_GRID
            if( (dptry(Lnew-1)+(hfacC(i,j,L,bi,bj)*drF(L)-deltap))
./pkg/gridalt/make_phys_grid.F 191 MAKE_PHYS_GRID
     .                             (hfacC(i,j,L,bi,bj)*drF(L))-deltap
./pkg/gridalt/make_phys_grid.F 192 MAKE_PHYS_GRID
             deltap = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 196 MAKE_PHYS_GRID
     .                      (hfacC(i,j,L,bi,bj)*drF(L)-deltap)
./pkg/gridalt/make_phys_grid.F 197 MAKE_PHYS_GRID
             deltap = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 206 MAKE_PHYS_GRID
           if(abs(deltap-hfacC(i,j,L-1,bi,bj)*drF(L-1)).gt.0.001)then
./pkg/gridalt/make_phys_grid.F 208 MAKE_PHYS_GRID
            dpphys(i,j,nlphys,bi,bj) = hfacC(i,j,L-1,bi,bj)*drF(L-1)
./pkg/gridalt/make_phys_grid.F 216 MAKE_PHYS_GRID
          if(abs(deltap-hfacC(i,j,L-1,bi,bj)*drF(L-1)).gt.0.001)then
./pkg/gridalt/make_phys_grid.F 218 MAKE_PHYS_GRID
           dpphys(i,j,nlphys,bi,bj) = hfacC(i,j,L-1,bi,bj)*drF(L-1) 
./pkg/gridalt/make_phys_grid.F 221 MAKE_PHYS_GRID
           dpphys(i,j,nlphys,bi,bj) = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 224 MAKE_PHYS_GRID
           dpphys(i,j,nlphys,bi,bj) = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 230 MAKE_PHYS_GRID
          dpphys(i,j,nlphys,bi,bj) = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 232 MAKE_PHYS_GRID
         nlperdyn(i,j,L,bi,bj) = nlphys
./pkg/gridalt/make_phys_grid.F 55 MAKE_PHYS_GRID
       integer i,j,L,Lbotij,Lnew
./pkg/gridalt/make_phys_grid.F 77 MAKE_PHYS_GRID
        do L = 1,ntry
./pkg/gridalt/make_phys_grid.F 78 MAKE_PHYS_GRID
         dptry(L) = dptry10(L)
./pkg/gridalt/make_phys_grid.F 83 MAKE_PHYS_GRID
        do L = 1,ntry
./pkg/gridalt/make_phys_grid.F 84 MAKE_PHYS_GRID
         dptry(L) = dptry40(L)
./pkg/gridalt/make_phys_grid.F 94 MAKE_PHYS_GRID
       do L = 1,Nr
./pkg/gridalt/make_phys_grid.F 97 MAKE_PHYS_GRID
         nlperdyn(i,j,L,bi,bj) = 0
./pkg/gridalt/phys2dyn.F 51 PHYS2DYN
      integer  i,j,L,Lout1,Lout1p1,Lout2,Lphy
./pkg/gridalt/phys2dyn.F 70 PHYS2DYN
      do L = 1,lmdyn
./pkg/gridalt/phys2dyn.F 74 PHYS2DYN
         qdyn(i,j,L,bi,bj) = 0.
./pkg/gridalt/phys2dyn.F 76 PHYS2DYN
         if(L.ge.Lbot(i,j,bi,bj))then
./pkg/gridalt/phys2dyn.F 77 PHYS2DYN
          if(L.eq.Lbot(i,j,bi,bj)) then
./pkg/gridalt/phys2dyn.F 80 PHYS2DYN
           Lout1 = nlperdyn(i,j,L-1,bi,bj)
./pkg/gridalt/phys2dyn.F 82 PHYS2DYN
          Lout2 = nlperdyn(i,j,L,bi,bj)
./pkg/gridalt/phys2dyn.F 86 PHYS2DYN
          dpkedyn = pedyn(i,j,L,bi,bj)-pedyn(i,j,L+1,bi,bj)
./pkg/gridalt/phys2dyn.F 95 PHYS2DYN
          qdyn(i,j,L,bi,bj) = sum
./pkg/mdsio/mdsio_buffertorl.F 34 MDS_BUFFERTORL
      INTEGER L
./pkg/mdsio/mdsio_buffertorl.F 39 MDS_BUFFERTORL
        DO L=1,narr
./pkg/mdsio/mdsio_buffertorl.F 40 MDS_BUFFERTORL
         arr(L) = vecBuf(L)
./pkg/mdsio/mdsio_buffertorl.F 43 MDS_BUFFERTORL
        DO L=1,narr
./pkg/mdsio/mdsio_buffertorl.F 44 MDS_BUFFERTORL
          vecBuf(L) = arr(L)
./pkg/mdsio/mdsio_buffertors.F 34 MDS_BUFFERTORS
      INTEGER L
./pkg/mdsio/mdsio_buffertors.F 39 MDS_BUFFERTORS
        DO L=1,narr
./pkg/mdsio/mdsio_buffertors.F 40 MDS_BUFFERTORS
         arr(L) = vecBuf(L)
./pkg/mdsio/mdsio_buffertors.F 43 MDS_BUFFERTORS
        DO L=1,narr
./pkg/mdsio/mdsio_buffertors.F 44 MDS_BUFFERTORS
          vecBuf(L) = arr(L)
./pkg/mdsio/mdsio_readvector.F 294 MDSREADVECTOR
           DO L=1,narr*nPx*nPy
./pkg/mdsio/mdsio_readvector.F 295 MDSREADVECTOR
            global(L) = xy_buffer_r4(L)
./pkg/mdsio/mdsio_readvector.F 306 MDSREADVECTOR
           DO L=1,narr*nPx*nPy
./pkg/mdsio/mdsio_readvector.F 307 MDSREADVECTOR
            global(L) = xy_buffer_r8(L)
./pkg/mdsio/mdsio_writevector.F 180 MDSWRITEVECTOR
           DO L=1,narr*nPx*nPy
./pkg/mdsio/mdsio_writevector.F 181 MDSWRITEVECTOR
            xy_buffer_r4(L) = global(L)
./pkg/mdsio/mdsio_writevector.F 192 MDSWRITEVECTOR
           DO L=1,narr*nPx*nPy
./pkg/mdsio/mdsio_writevector.F 193 MDSWRITEVECTOR
            xy_buffer_r8(L) = global(L)
./pkg/mdsio/mdsio_writevector.F 81 MDSWRITEVECTOR
      INTEGER K,L
./pkg/sphere/sphere.F 266 SHCROTATE
      INTEGER K, L, M
./pkg/sphere/sphere.F 281 SHCROTATE
            DO L = M,LMAX
./pkg/sphere/sphere.F 282 SHCROTATE
               K = L*L+L+1
./pkg/sphere/sphere.F 313 SHC2ZONE
      INTEGER LMAX1, J, K, L, M
./pkg/sphere/sphere.F 335 SHC2ZONE
         DO L = M,LMAX-1,2
./pkg/sphere/sphere.F 337 SHC2ZONE
            J    = L*L+L+1
./pkg/sphere/sphere.F 338 SHC2ZONE
            Q    = FACT(L)*P(K)
./pkg/sphere/sphere.F 346 SHC2ZONE
            J    = J+L+L+2
./pkg/sphere/sphere.F 347 SHC2ZONE
            Q    = FACT(L+1)*P(K)
./pkg/sphere/sphere.F 355 SHC2ZONE
         IF(MOD(LMAX-L,2).EQ.0) THEN
./pkg/sphere/sphere.F 404 SHC4ZONE
      INTEGER I, J, K, L, M
./pkg/sphere/sphere.F 445 SHC4ZONE
         DO L = M,LMAX-1,2
./pkg/sphere/sphere.F 446 SHC4ZONE
            K = L*L+L+1
./pkg/sphere/sphere.F 449 SHC4ZONE
            K = K+L+L+2
./pkg/sphere/sphere.F 506 FRSBASE
      INTEGER N, L
./pkg/sphere/sphere.F 524 FRSBASE
      DO L = I+5,J
./pkg/sphere/sphere.F 525 FRSBASE
      H(L)   = T*H(L-2)-H(L-4)
./pkg/sphere/sphere.F 557 HELMHOLTZ
      INTEGER K, L, N, M
./pkg/sphere/sphere.F 561 HELMHOLTZ
         DO L = 1,LMAX
./pkg/sphere/sphere.F 562 HELMHOLTZ
            X(L) = SNGL(1.0D0/DSQRT(DBLE(FLOAT(4*L*L-1))))
./pkg/sphere/sphere.F 563 HELMHOLTZ
            Y(L) = SNGL(DSQRT(DBLE(FLOAT(L))))
./pkg/sphere/sphere.F 564 HELMHOLTZ
            Y(L+LMAX) = SNGL(DSQRT(DBLE(FLOAT(L+LMAX))))
./pkg/sphere/sphere.F 565 HELMHOLTZ
            ISECT(L) = L*LMAX-L*(L-3)/2+1
./pkg/sphere/sphere.F 582 HELMHOLTZ
            L = N+M
./pkg/sphere/sphere.F 583 HELMHOLTZ
            AK   = X(L)*Y(L+M)*Y(N)
./pkg/sphere/sphere.F 590 HELMHOLTZ
               L = N+M
./pkg/sphere/sphere.F 591 HELMHOLTZ
               AK   = X(L)*Y(L+M)*Y(N)
./pkg/sphere/sphere.F 611 RECUR_Z
      INTEGER L
./pkg/sphere/sphere.F 614 RECUR_Z
      DO L = 1,LMAX
./pkg/sphere/sphere.F 615 RECUR_Z
         ZZ   = ZZ * DBLE(FLOAT(L+L+1))/DBLE(FLOAT(L+L))
./pkg/sphere/sphere.F 616 RECUR_Z
         Z(L) = SNGL(DSQRT(ZZ))
./pkg/zonal_filt/fftpack.F 1206 R8ADBG
      DO 120 L=2,IPPH
./pkg/zonal_filt/fftpack.F 1207 R8ADBG
         LC = IPP2-L
./pkg/zonal_filt/fftpack.F 1212 R8ADBG
            C2(IK,L) = CH2(IK,1)+AR1*CH2(IK,2)
./pkg/zonal_filt/fftpack.F 1225 R8ADBG
               C2(IK,L) = C2(IK,L)+AR2*CH2(IK,J)
./pkg/zonal_filt/fftpack.F 1642 R8ADFG
      DO 127 L=2,IPPH
./pkg/zonal_filt/fftpack.F 1643 R8ADFG
         LC = IPP2-L
./pkg/zonal_filt/fftpack.F 1648 R8ADFG
            CH2(IK,L) = C2(IK,1)+AR1*C2(IK,2)
./pkg/zonal_filt/fftpack.F 1661 R8ADFG
               CH2(IK,L) = CH2(IK,L)+AR2*C2(IK,J)
./pkg/zonal_filt/fftpack.F 279 RADBG
      DO 120 L=2,IPPH
./pkg/zonal_filt/fftpack.F 280 RADBG
         LC = IPP2-L
./pkg/zonal_filt/fftpack.F 285 RADBG
            C2(IK,L) = CH2(IK,1)+AR1*CH2(IK,2)
./pkg/zonal_filt/fftpack.F 298 RADBG
               C2(IK,L) = C2(IK,L)+AR2*CH2(IK,J)
./pkg/zonal_filt/fftpack.F 729 RADFG
      DO 127 L=2,IPPH
./pkg/zonal_filt/fftpack.F 730 RADFG
         LC = IPP2-L
./pkg/zonal_filt/fftpack.F 735 RADFG
            CH2(IK,L) = C2(IK,1)+AR1*C2(IK,2)
./pkg/zonal_filt/fftpack.F 748 RADFG
               CH2(IK,L) = CH2(IK,L)+AR2*C2(IK,J)
./verification/OpenAD/code_oad_all/ad_s_ifnblnk.F 23 OAD_S_IFNBLNK
      INTEGER L, LS
./verification/OpenAD/code_oad_all/ad_s_ifnblnk.F 28 OAD_S_IFNBLNK
      L=1
./verification/OpenAD/code_oad_all/ad_s_ifnblnk.F 29 OAD_S_IFNBLNK
      DO while ((L.le.LS).and.(string(L:L) .eq. ' '))
./verification/OpenAD/code_oad_all/ad_s_ifnblnk.F 30 OAD_S_IFNBLNK
        L=L+1
./verification/OpenAD/code_oad_all/ad_s_ifnblnk.F 32 OAD_S_IFNBLNK
      index = L
./verification/OpenAD/code_oad_all/ad_s_ilnblnk.F 23 OAD_S_ILNBLNK
      INTEGER L, LS
./verification/OpenAD/code_oad_all/ad_s_ilnblnk.F 28 OAD_S_ILNBLNK
      L=LS
./verification/OpenAD/code_oad_all/ad_s_ilnblnk.F 29 OAD_S_ILNBLNK
      DO while ((L.gt.0).and.(string(L:L) .eq. ' '))
./verification/OpenAD/code_oad_all/ad_s_ilnblnk.F 30 OAD_S_ILNBLNK
        L=L-1
./verification/OpenAD/code_oad_all/ad_s_ilnblnk.F 32 OAD_S_ILNBLNK
      index = L
./verification/fizhi-cs-aqualev20/code/fizhi_init_fixed.F 107 FIZHI_INIT_FIXED
       do L = 2,Nrphys+1
./verification/fizhi-cs-aqualev20/code/fizhi_init_fixed.F 108 FIZHI_INIT_FIXED
        pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100.
./verification/fizhi-cs-aqualev20/code/fizhi_init_fixed.F 112 FIZHI_INIT_FIXED
       do L = 1,Nrphys
./verification/fizhi-cs-aqualev20/code/fizhi_init_fixed.F 113 FIZHI_INIT_FIXED
        pressure(L)=(pressure0(Nrphys+2-L)+pressure0(Nrphys+1-L))/2.
./verification/fizhi-cs-aqualev20/code/fizhi_init_fixed.F 40 FIZHI_INIT_FIXED
       integer i,j,L,bi,bj
./verification/fizhi-gridalt-hs/code/do_fizhi.F 130 DO_FIZHI
      integer i,j,L
./verification/fizhi-gridalt-hs/code/do_fizhi.F 169 DO_FIZHI
       do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/do_fizhi.F 174 DO_FIZHI
         termP=0.5 _d 0*((pephy(i,j,L,bi,bj)+pephy(i,j,L+1,bi,bj))
./verification/fizhi-gridalt-hs/code/do_fizhi.F 177 DO_FIZHI
         duphy(i,j,L,bi,bj)= -kV*uphy(i,j,L,bi,bj)
./verification/fizhi-gridalt-hs/code/do_fizhi.F 178 DO_FIZHI
         dvphy(i,j,L,bi,bj)= -kV*vphy(i,j,L,bi,bj)
./verification/fizhi-gridalt-hs/code/do_fizhi.F 183 DO_FIZHI
         termP=0.5 _d 0*( pephy(i,j,L,bi,bj) + pephy(i,j,L+1,bi,bj) )
./verification/fizhi-gridalt-hs/code/do_fizhi.F 194 DO_FIZHI
          dthphy(i,j,L,bi,bj)=- kT*( thphy(I,J,L,bi,bj)-thetaEq )
./verification/fizhi-gridalt-hs/code/do_fizhi.F 196 DO_FIZHI
          dthphy(i,j,L,bi,bj)=0.
./verification/fizhi-gridalt-hs/code/do_fizhi.F 201 DO_FIZHI
         dsphy(i,j,L,bi,bj)=0.
./verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F 101 FIZHI_INIT_FIXED
       do L = 2,Nrphys+1
./verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F 102 FIZHI_INIT_FIXED
        pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100.
./verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F 106 FIZHI_INIT_FIXED
       do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F 107 FIZHI_INIT_FIXED
        pressure(L)=(pressure0(Nrphys+2-L)+pressure0(Nrphys+1-L))/2.
./verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F 40 FIZHI_INIT_FIXED
       integer i,j,L,bi,bj
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 110 FIZHI_INIT_VARS
         do L = 2,Nrphys+1
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 111 FIZHI_INIT_VARS
          pephy(i,j,L,bi,bj)=pephy(i,j,L-1,bi,bj)-dpphys0(i,j,L-1,bi,bj)
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 122 FIZHI_INIT_VARS
        do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 125 FIZHI_INIT_VARS
         windphy(i,j,L,bi,bj) = 0.025 *
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 126 FIZHI_INIT_VARS
     .             log((pephy(i,j,1,bi,bj)-pephy(i,j,L+1,bi,bj))*10.)
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 146 FIZHI_INIT_VARS
        do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 149 FIZHI_INIT_VARS
         uphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 155 FIZHI_INIT_VARS
        do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 158 FIZHI_INIT_VARS
         vphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 164 FIZHI_INIT_VARS
        do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 167 FIZHI_INIT_VARS
         thphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 174 FIZHI_INIT_VARS
        do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 177 FIZHI_INIT_VARS
         sphy(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 183 FIZHI_INIT_VARS
        do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 186 FIZHI_INIT_VARS
         duphy(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 187 FIZHI_INIT_VARS
         dvphy(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 188 FIZHI_INIT_VARS
         dthphy(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 189 FIZHI_INIT_VARS
         dsphy(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 195 FIZHI_INIT_VARS
        do L = 1,Nr
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 198 FIZHI_INIT_VARS
         guphy(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 199 FIZHI_INIT_VARS
         gvphy(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 200 FIZHI_INIT_VARS
         gthphy(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 201 FIZHI_INIT_VARS
         gsphy(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 215 FIZHI_INIT_VARS
         do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 217 FIZHI_INIT_VARS
          tke(i,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 218 FIZHI_INIT_VARS
          xlmt(i,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 219 FIZHI_INIT_VARS
          khmt(i,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 251 FIZHI_INIT_VARS
        do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 254 FIZHI_INIT_VARS
         swlz(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 255 FIZHI_INIT_VARS
         lwlz(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 256 FIZHI_INIT_VARS
         qliqavesw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 257 FIZHI_INIT_VARS
         qliqavelw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 258 FIZHI_INIT_VARS
         fccavesw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 259 FIZHI_INIT_VARS
         fccavelw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 260 FIZHI_INIT_VARS
         cldtot_sw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 261 FIZHI_INIT_VARS
         cldras_sw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 262 FIZHI_INIT_VARS
         cldlsp_sw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 263 FIZHI_INIT_VARS
         cldtot_lw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 264 FIZHI_INIT_VARS
         cldras_lw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 265 FIZHI_INIT_VARS
         cldlsp_lw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 48 FIZHI_INIT_VARS
       integer i, j, L, bi, bj, Lbotij
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 82 FIZHI_INIT_VARS
         do L = 1,Nr
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 83 FIZHI_INIT_VARS
          pedyn(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 97 FIZHI_INIT_VARS
         do L = Lbotij+1,Nr+1
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 98 FIZHI_INIT_VARS
          pedyn(i,j,L,bi,bj) = pedyn(i,j,L-1,bi,bj) -
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 99 FIZHI_INIT_VARS
     .                        drF(L-1)*hfacC(i,j,L-1,bi,bj)
./verification/flt_example/extra/cvprofiles.F 658 NO PROCEDURE
      INTEGER L, LS
./verification/flt_example/extra/cvprofiles.F 662 NO PROCEDURE
      DO 10 L = LS, 1, -1
./verification/flt_example/extra/cvprofiles.F 663 NO PROCEDURE
        IF ( string(L:L) .EQ. ' ' ) GOTO 10
./verification/flt_example/extra/cvprofiles.F 664 NO PROCEDURE
         ILNBLNK = L
./eesupp/src/eedie.F 38 EEDIE
      INTEGER I, L
./eesupp/src/eedie.F 61 EEDIE
       L = ILNBLNK(msgBuf)
./eesupp/src/eedie.F 62 EEDIE
       WRITE(*,*) msgBuf(1:L)
./optim/optim_main.F 119 NO PROCEDURE
      INTEGER L, LS
./optim/optim_main.F 123 NO PROCEDURE
      DO 10 L = 1, LS
./optim/optim_main.F 124 NO PROCEDURE
       IF ( string(L:L) .EQ. ' ' ) GOTO 10
./optim/optim_main.F 125 NO PROCEDURE
        IFNBLNK = L
./optim/optim_main.F 142 NO PROCEDURE
      INTEGER L, LS
./optim/optim_main.F 146 NO PROCEDURE
      DO 10 L = LS, 1, -1
./optim/optim_main.F 147 NO PROCEDURE
        IF ( string(L:L) .EQ. ' ' ) GOTO 10
./optim/optim_main.F 148 NO PROCEDURE
         ILNBLNK = L
./optim/utils.F 105 DATE
      INTEGER L, LS
./optim/utils.F 109 DATE
      DO 10 L = LS, 1, -1
./optim/utils.F 110 DATE
        IF ( string(L:L) .EQ. ' ' ) GOTO 10
./optim/utils.F 111 DATE
         ILNBLNK = L
./optim/utils.F 154 LCASE
      INTEGER   I, L
./optim/utils.F 157 LCASE
        L = INDEX(UPPER,string(I:I))
./optim/utils.F 158 LCASE
        IF ( L .EQ. 0 ) GOTO 10
./optim/utils.F 159 LCASE
         string(I:I) = LOWER(L:L)
./optim/utils.F 247 UCASE
      INTEGER   I, L
./optim/utils.F 250 UCASE
        L = INDEX(LOWER,string(I:I))
./optim/utils.F 251 UCASE
        IF ( L .EQ. 0 ) GOTO 10
./optim/utils.F 252 UCASE
          string(I:I) = UPPER(L:L)
./optim/utils.F 82 DATE
      INTEGER L, LS
./optim/utils.F 86 DATE
      DO 10 L = 1, LS
./optim/utils.F 87 DATE
       IF ( string(L:L) .EQ. ' ' ) GOTO 10
./optim/utils.F 88 DATE
        IFNBLNK = L
./eesupp/src/utils.F 104 DATE
      INTEGER L, LS
./eesupp/src/utils.F 109 DATE
      DO 10 L = 1, LS
./eesupp/src/utils.F 110 DATE
       IF ( string(L:L) .EQ. ' ' ) GOTO 10
./eesupp/src/utils.F 111 DATE
        IFNBLNK = L
./eesupp/src/utils.F 138 DATE
      INTEGER L, LS
./eesupp/src/utils.F 144 DATE
      DO 10 L = LS, 1, -1
./eesupp/src/utils.F 145 DATE
        IF ( string(L:L) .EQ. ' ' ) GOTO 10
./eesupp/src/utils.F 146 DATE
         ILNBLNK = L
./eesupp/src/utils.F 208 LCASE
      INTEGER   I, L
./eesupp/src/utils.F 213 LCASE
        L = INDEX(UPPER,string(I:I))
./eesupp/src/utils.F 214 LCASE
        IF ( L .EQ. 0 ) GOTO 10
./eesupp/src/utils.F 215 LCASE
         string(I:I) = LOWER(L:L)
./eesupp/src/utils.F 326 UCASE
      INTEGER   I, L
./eesupp/src/utils.F 331 UCASE
        L = INDEX(LOWER,string(I:I))
./eesupp/src/utils.F 332 UCASE
        IF ( L .EQ. 0 ) GOTO 10
./eesupp/src/utils.F 333 UCASE
          string(I:I) = UPPER(L:L)
./pkg/diagnostics/diagnostics_interp_p2p.F 139 DIAGNOSTICS_INTERP_P2P
      DO L=1,lm-1
./pkg/diagnostics/diagnostics_interp_p2p.F 140 DIAGNOSTICS_INTERP_P2P
       pkmin = pkz(1,L+1)
./pkg/diagnostics/diagnostics_interp_p2p.F 141 DIAGNOSTICS_INTERP_P2P
       pkmax = pkz(1,L)
./pkg/diagnostics/diagnostics_interp_p2p.F 144 DIAGNOSTICS_INTERP_P2P
         IF ( pkz(i,L+1).LT.pkmin ) pkmin = pkz(i,L+1)
./pkg/diagnostics/diagnostics_interp_p2p.F 145 DIAGNOSTICS_INTERP_P2P
         IF ( pkz(i,L)  .GT.pkmax ) pkmax = pkz(i,L)
./pkg/diagnostics/diagnostics_interp_p2p.F 150 DIAGNOSTICS_INTERP_P2P
           IF ( pk.LE.pkz(i,L) .AND. pk.GE.pkz(i,L+1) ) THEN
./pkg/diagnostics/diagnostics_interp_p2p.F 151 DIAGNOSTICS_INTERP_P2P
             temp = ( pk-pkz(i,L) ) / ( pkz(i,L+1)-pkz(i,L) )
./pkg/diagnostics/diagnostics_interp_p2p.F 153 DIAGNOSTICS_INTERP_P2P
             IF     ( qinp(i,L)  .NE.undef  .AND.
./pkg/diagnostics/diagnostics_interp_p2p.F 154 DIAGNOSTICS_INTERP_P2P
     &                qinp(i,L+1).NE.undef ) THEN
./pkg/diagnostics/diagnostics_interp_p2p.F 155 DIAGNOSTICS_INTERP_P2P
               qprs(i) = qinp(i,L+1)*temp + qinp(i,L)*(1.-temp)
./pkg/diagnostics/diagnostics_interp_p2p.F 156 DIAGNOSTICS_INTERP_P2P
             ELSEIF ( qinp(i,L+1).NE.undef  .AND. temp.GE.0.5 ) THEN
./pkg/diagnostics/diagnostics_interp_p2p.F 157 DIAGNOSTICS_INTERP_P2P
               qprs(i) = qinp(i,L+1)
./pkg/diagnostics/diagnostics_interp_p2p.F 158 DIAGNOSTICS_INTERP_P2P
             ELSEIF ( qinp(i,L)  .NE.undef  .AND. temp.LE.0.5 ) THEN
./pkg/diagnostics/diagnostics_interp_p2p.F 159 DIAGNOSTICS_INTERP_P2P
               qprs(i) = qinp(i,L)
./pkg/diagnostics/diagnostics_interp_p2p.F 76 DIAGNOSTICS_INTERP_P2P
      DO L=1,lm-1
./pkg/diagnostics/diagnostics_interp_p2p.F 77 DIAGNOSTICS_INTERP_P2P
       pkmin = pkz(1,L)
./pkg/diagnostics/diagnostics_interp_p2p.F 78 DIAGNOSTICS_INTERP_P2P
       pkmax = pkz(1,L+1)
./pkg/diagnostics/diagnostics_interp_p2p.F 81 DIAGNOSTICS_INTERP_P2P
         IF ( pkz(i,L)  .LT.pkmin ) pkmin = pkz(i,L)
./pkg/diagnostics/diagnostics_interp_p2p.F 82 DIAGNOSTICS_INTERP_P2P
         IF ( pkz(i,L+1).GT.pkmax ) pkmax = pkz(i,L+1)
./pkg/diagnostics/diagnostics_interp_p2p.F 87 DIAGNOSTICS_INTERP_P2P
           IF ( pk.GE.pkz(i,L) .AND. pk.LE.pkz(i,L+1) ) THEN
./pkg/diagnostics/diagnostics_interp_p2p.F 88 DIAGNOSTICS_INTERP_P2P
             temp = ( pk-pkz(i,L) ) / ( pkz(i,L+1)-pkz(i,L) )
./pkg/diagnostics/diagnostics_interp_p2p.F 90 DIAGNOSTICS_INTERP_P2P
             IF     ( qinp(i,L)  .NE.undef  .AND.
./pkg/diagnostics/diagnostics_interp_p2p.F 91 DIAGNOSTICS_INTERP_P2P
     &                qinp(i,L+1).NE.undef ) THEN
./pkg/diagnostics/diagnostics_interp_p2p.F 92 DIAGNOSTICS_INTERP_P2P
               qprs(i) = qinp(i,L+1)*temp + qinp(i,L)*(1.-temp)
./pkg/diagnostics/diagnostics_interp_p2p.F 93 DIAGNOSTICS_INTERP_P2P
             ELSEIF ( qinp(i,L+1).NE.undef  .AND. temp.GE.0.5 ) THEN
./pkg/diagnostics/diagnostics_interp_p2p.F 94 DIAGNOSTICS_INTERP_P2P
               qprs(i) = qinp(i,L+1)
./pkg/diagnostics/diagnostics_interp_p2p.F 95 DIAGNOSTICS_INTERP_P2P
             ELSEIF ( qinp(i,L)  .NE.undef  .AND. temp.LE.0.5 ) THEN
./pkg/diagnostics/diagnostics_interp_p2p.F 96 DIAGNOSTICS_INTERP_P2P
               qprs(i) = qinp(i,L)
./pkg/fizhi/AtoC.F 40 ATOC
       integer i, j, L, bi, bj
./pkg/fizhi/AtoC.F 53 ATOC
        do L = 1,numlevs
./pkg/fizhi/AtoC.F 56 ATOC
           tmpfld1(i,j) = fieldin1(i,j,L,bi,bj)*mask(i,j,L,bi,bj)
./pkg/fizhi/AtoC.F 57 ATOC
           tmpfld2(i,j) = fieldin2(i,j,L,bi,bj)*mask(i,j,L,bi,bj)
./pkg/fizhi/AtoC.F 64 ATOC
          if( (mask(i-1,j,L,bi,bj).ne.0.) .or.
./pkg/fizhi/AtoC.F 65 ATOC
     .                                  (mask(i,j,L,bi,bj).ne.0.) ) then
./pkg/fizhi/AtoC.F 70 ATOC
           fieldout1(i,j,L,bi,bj) = 
./pkg/fizhi/AtoC.F 72 ATOC
     .           ( mask(i-1,j,L,bi,bj) + mask(i,j,L,bi,bj) )
./pkg/fizhi/AtoC.F 74 ATOC
           fieldout1(i,j,L,bi,bj) = 0.
./pkg/fizhi/AtoC.F 76 ATOC
          if( (mask(i,j-1,L,bi,bj).ne.0.) .or.
./pkg/fizhi/AtoC.F 77 ATOC
     .                                  (mask(i,j,L,bi,bj).ne.0.) ) then
./pkg/fizhi/AtoC.F 82 ATOC
           fieldout2(i,j,L,bi,bj) = 
./pkg/fizhi/AtoC.F 84 ATOC
     .           ( mask(i,j,L,bi,bj) + mask(i,j-1,L,bi,bj) )
./pkg/fizhi/AtoC.F 86 ATOC
           fieldout2(i,j,L,bi,bj) = 0.
./pkg/fizhi/CtoA.F 42 CTOA
       integer i, j, L, bi, bj
./pkg/fizhi/CtoA.F 55 CTOA
        do L = 1,numlevs
./pkg/fizhi/CtoA.F 58 CTOA
           tmpfld1(i,j) = fieldin1(i,j,L,bi,bj)*mask1(i,j,L,bi,bj)
./pkg/fizhi/CtoA.F 59 CTOA
           tmpfld2(i,j) = fieldin2(i,j,L,bi,bj)*mask2(i,j,L,bi,bj)
./pkg/fizhi/CtoA.F 64 CTOA
          if( (mask1(i,j,L,bi,bj).ne.0.) .or. 
./pkg/fizhi/CtoA.F 65 CTOA
     .                               (mask1(i+1,j,L,bi,bj).ne.0.) ) then
./pkg/fizhi/CtoA.F 70 CTOA
           fieldout1(i,j,L,bi,bj) = 
./pkg/fizhi/CtoA.F 72 CTOA
     .            ( mask1(i,j,L,bi,bj) + mask1(i+1,j,L,bi,bj) )
./pkg/fizhi/CtoA.F 74 CTOA
           fieldout1(i,j,L,bi,bj) = 0.
./pkg/fizhi/CtoA.F 76 CTOA
          if( (mask2(i,j,L,bi,bj).ne.0.) .or. 
./pkg/fizhi/CtoA.F 77 CTOA
     .                               (mask2(i,j+1,L,bi,bj).ne.0.) ) then
./pkg/fizhi/CtoA.F 82 CTOA
           fieldout2(i,j,L,bi,bj) = 
./pkg/fizhi/CtoA.F 84 CTOA
     .            ( mask2(i,j,L,bi,bj) + mask2(i,j+1,L,bi,bj) )
./pkg/fizhi/CtoA.F 86 CTOA
           fieldout2(i,j,L,bi,bj) = 0.