Symbol Description Uses
L Temps for string locations 2610



File Line number Procedure Code
./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.
./pkg/fizhi/do_fizhi.F 131 DO_FIZHI
      INTEGER i,j,L
./pkg/fizhi/do_fizhi.F 203 DO_FIZHI
      do L = 1,Nrphys
./pkg/fizhi/do_fizhi.F 206 DO_FIZHI
       u(i,j,L) = uphy(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 207 DO_FIZHI
       v(i,j,L) = vphy(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 208 DO_FIZHI
       t(i,j,L) = thphy(i,j,L,bi,bj)/p0kappa
./pkg/fizhi/do_fizhi.F 209 DO_FIZHI
       q(i,j,L,1) = sphy(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 210 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 211 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 212 DO_FIZHI
       ple(i,j,L,bi,bj) = pephy(i,j,L,bi,bj)/100.
./pkg/fizhi/do_fizhi.F 213 DO_FIZHI
       if (ple(i,j,L,bi,bj).gt.0.) then
./pkg/fizhi/do_fizhi.F 214 DO_FIZHI
        pkle(i,j,L,bi,bj) = ple(i,j,L,bi,bj) **kappa
./pkg/fizhi/do_fizhi.F 216 DO_FIZHI
        pkle(i,j,L,bi,bj) = 0.
./pkg/fizhi/do_fizhi.F 218 DO_FIZHI
       cldtot_lw(i,j,L) = cldtotlwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 219 DO_FIZHI
       cldlsp_lw(i,j,L) = cldlsplwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 220 DO_FIZHI
       cldras_lw(i,j,L) = cldraslwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 221 DO_FIZHI
       lwlz(i,j,L) = lwlzin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 222 DO_FIZHI
       qliqavelw(i,j,L) = qliqavelwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 223 DO_FIZHI
       fccavelw(i,j,L) = fccavelwin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 224 DO_FIZHI
       cldtot_sw(i,j,L) = cldtotswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 225 DO_FIZHI
       cldlsp_sw(i,j,L) = cldlspswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 226 DO_FIZHI
       cldras_sw(i,j,L) = cldrasswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 227 DO_FIZHI
       swlz(i,j,L) = swlzin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 228 DO_FIZHI
       qliqavesw(i,j,L) = qliqaveswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 229 DO_FIZHI
       fccavesw(i,j,L) = fccaveswin(i,j,L,bi,bj)
./pkg/fizhi/do_fizhi.F 268 DO_FIZHI
      do L = 1,Nrphys
./pkg/fizhi/do_fizhi.F 271 DO_FIZHI
       cldtotlwin(i,j,L,bi,bj) = cldtot_lw(i,j,L)
./pkg/fizhi/do_fizhi.F 272 DO_FIZHI
       cldlsplwin(i,j,L,bi,bj) = cldlsp_lw(i,j,L)
./pkg/fizhi/do_fizhi.F 273 DO_FIZHI
       cldraslwin(i,j,L,bi,bj) = cldras_lw(i,j,L)
./pkg/fizhi/do_fizhi.F 274 DO_FIZHI
       lwlzin(i,j,L,bi,bj) = lwlz(i,j,L)
./pkg/fizhi/do_fizhi.F 275 DO_FIZHI
       qliqavelwin(i,j,L,bi,bj) = qliqavelw(i,j,L)
./pkg/fizhi/do_fizhi.F 276 DO_FIZHI
       fccavelwin(i,j,L,bi,bj) = fccavelw(i,j,L)
./pkg/fizhi/do_fizhi.F 277 DO_FIZHI
       cldtotswin(i,j,L,bi,bj) = cldtot_sw(i,j,L)
./pkg/fizhi/do_fizhi.F 278 DO_FIZHI
       cldlspswin(i,j,L,bi,bj) = cldlsp_sw(i,j,L)
./pkg/fizhi/do_fizhi.F 279 DO_FIZHI
       cldrasswin(i,j,L,bi,bj) = cldras_sw(i,j,L)
./pkg/fizhi/do_fizhi.F 280 DO_FIZHI
       swlzin(i,j,L,bi,bj) = swlz(i,j,L)
./pkg/fizhi/do_fizhi.F 281 DO_FIZHI
       qliqaveswin(i,j,L,bi,bj) = qliqavesw(i,j,L)
./pkg/fizhi/do_fizhi.F 282 DO_FIZHI
       fccaveswin(i,j,L,bi,bj) = fccavesw(i,j,L)
./pkg/fizhi/do_fizhi.F 297 DO_FIZHI
      do L = 1,Nrphys
./pkg/fizhi/do_fizhi.F 300 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 301 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 302 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 303 DO_FIZHI
     &   lwdt(i,j,L,bi,bj) +
./pkg/fizhi/do_fizhi.F 304 DO_FIZHI
     &   dlwdtg(i,j,L,bi,bj) * (tgz(i,j,bi,bj)-tg0(i,j,bi,bj)) +
./pkg/fizhi/do_fizhi.F 305 DO_FIZHI
     &   swdt(i,j,L,bi,bj)*radswt(i,j,bi,bj) )*p0kappa ) / p(i,j,bi,bj)
./pkg/fizhi/do_fizhi.F 306 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 260 FIZHI_INIT_VARS
        do L = 1,Nrphys
./pkg/fizhi/fizhi_init_vars.F 263 FIZHI_INIT_VARS
         swlz(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 264 FIZHI_INIT_VARS
         lwlz(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 265 FIZHI_INIT_VARS
         qliqavesw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 266 FIZHI_INIT_VARS
         qliqavelw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 267 FIZHI_INIT_VARS
         fccavesw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 268 FIZHI_INIT_VARS
         fccavelw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 269 FIZHI_INIT_VARS
         cldtot_sw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 270 FIZHI_INIT_VARS
         cldras_sw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 271 FIZHI_INIT_VARS
         cldlsp_sw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 272 FIZHI_INIT_VARS
         cldtot_lw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 273 FIZHI_INIT_VARS
         cldras_lw(i,j,L,bi,bj) = 0.
./pkg/fizhi/fizhi_init_vars.F 274 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 163 FIZHI_STEP_DIAG
      do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 170 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 173 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'DIABU   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 181 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 184 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'DIABV   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 193 FIZHI_STEP_DIAG
     &   ( turbt(i,j,L,bi,bj) + moistt(i,j,L,bi,bj) +
./pkg/fizhi/fizhi_step_diag.F 194 FIZHI_STEP_DIAG
     &      lwdt(i,j,L,bi,bj) +
./pkg/fizhi/fizhi_step_diag.F 195 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 196 FIZHI_STEP_DIAG
     &      swdt(i,j,L,bi,bj)*radswt(i,j,bi,bj) )
./pkg/fizhi/fizhi_step_diag.F 197 FIZHI_STEP_DIAG
     &      * pk(i,j,L,bi,bj)*pinv(i,j)*86400
./pkg/fizhi/fizhi_step_diag.F 200 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'DIABT   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 209 FIZHI_STEP_DIAG
     & ( turbq(i,j,L,1,bi,bj) + moistq(i,j,L,1,bi,bj) ) *
./pkg/fizhi/fizhi_step_diag.F 213 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'DIABQ   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 223 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 227 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'RADLW   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 237 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 241 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'LWCLR   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 254 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'RADSW   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 267 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'SWCLR   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 275 FIZHI_STEP_DIAG
        tmpdiag(i,j) = uphy(i,j,L)
./pkg/fizhi/fizhi_step_diag.F 278 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'UWND    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 286 FIZHI_STEP_DIAG
        tmpdiag(i,j) = vphy(i,j,L)
./pkg/fizhi/fizhi_step_diag.F 289 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'VWND    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 297 FIZHI_STEP_DIAG
        tmpdiag(i,j) = thphy(i,j,L)*pk(i,j,L,bi,bj)
./pkg/fizhi/fizhi_step_diag.F 300 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'TMPU    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 308 FIZHI_STEP_DIAG
        tmpdiag(i,j) = qq(i,j,L,bi,bj)
./pkg/fizhi/fizhi_step_diag.F 311 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'TKE     ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 319 FIZHI_STEP_DIAG
        tmpdiag(i,j) = sphy(i,j,L) * 1000.
./pkg/fizhi/fizhi_step_diag.F 322 FIZHI_STEP_DIAG
       call DIAGNOSTICS_FILL(tmpdiag,'SPHU    ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_step_diag.F 337 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 341 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 361 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 365 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 385 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 388 FIZHI_STEP_DIAG
        qbar(i,j) = qbar(i,j) + ( lwdt(i,j,L,bi,bj) +
./pkg/fizhi/fizhi_step_diag.F 389 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 390 FIZHI_STEP_DIAG
     &             *pk(i,j,l,bi,bj)*dp(i,j,L,bi,bj)
./pkg/fizhi/fizhi_step_diag.F 410 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 414 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 436 FIZHI_STEP_DIAG
       do L=1,Nrphys
./pkg/fizhi/fizhi_step_diag.F 440 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 1130 TURBIO
      do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 1132 TURBIO
           duturb(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1133 TURBIO
           dvturb(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1134 TURBIO
           dtturb(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1135 TURBIO
           qqgrid(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1136 TURBIO
          qliqtmp(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1137 TURBIO
           fcctmp(i,1,L) = 0.
./pkg/fizhi/fizhi_turb.F 1141 TURBIO
           dqturb(i,1,L,nt) = 0.
./pkg/fizhi/fizhi_turb.F 1149 TURBIO
      call MSC2GRD(igrd,chfr,dumsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1150 TURBIO
     &                                              duturb(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1151 TURBIO
      call MSC2GRD(igrd,chfr,dvmsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1152 TURBIO
     &                                              dvturb(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1153 TURBIO
      call MSC2GRD(igrd,chfr,dtmsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1154 TURBIO
     &                                              dtturb(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1156 TURBIO
      call MSC2GRD(igrd,chfr,dqmsc(1,L,nt),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1157 TURBIO
     &                                           dqturb(1,1,L,nt),im,jm)
./pkg/fizhi/fizhi_turb.F 1159 TURBIO
      call MSC2GRD(igrd,chfr,  tke(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1160 TURBIO
     &                                              qqgrid(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1162 TURBIO
      call MSC2GRD(igrd,chfr, fccmsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1163 TURBIO
     &                                              fcctmp(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1164 TURBIO
      call MSC2GRD(igrd,chfr,qliqmsc(1,L),nchp,nchptot,fracland,
./pkg/fizhi/fizhi_turb.F 1165 TURBIO
     &                                             qliqtmp(1,1,L),im,jm)
./pkg/fizhi/fizhi_turb.F 1176 TURBIO
         fccavesw(i,1,L) =  fccavesw(i,1,L) +  fcctmp(i,1,L)
./pkg/fizhi/fizhi_turb.F 1177 TURBIO
         fccavelw(i,1,L) =  fccavelw(i,1,L) +  fcctmp(i,1,L)
./pkg/fizhi/fizhi_turb.F 1178 TURBIO
        qliqavelw(i,1,L) = qliqavelw(i,1,L) + qliqtmp(i,1,L)
./pkg/fizhi/fizhi_turb.F 1179 TURBIO
        qliqavesw(i,1,L) = qliqavesw(i,1,L) + qliqtmp(i,1,L)
./pkg/fizhi/fizhi_turb.F 1224 TURBIO
       do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 1228 TURBIO
        call DIAGNOSTICS_FILL(tmpdiag,'TURBT   ',L,1,3,bi,bj,myid)
./pkg/fizhi/fizhi_turb.F 1443 TRBFLX
      integer i,iter,init,n,LL,L,Lp,Lp1,lmin,lminq,lminq1,ibit
./pkg/fizhi/fizhi_turb.F 1472 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1474 TRBFLX
         wu(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1475 TRBFLX
         wv(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1476 TRBFLX
         eu(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1477 TRBFLX
         et(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1481 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1483 TRBFLX
         XLDIAG(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1487 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1489 TRBFLX
         FREQDG(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1536 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1538 TRBFLX
        XLSAVE(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1539 TRBFLX
        KHSAVE(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1552 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1554 TRBFLX
        ADZ1(I,L) = (CP/GRAV)*(PLKE(I,L+1)-PLKE(I,L))
./pkg/fizhi/fizhi_turb.F 1555 TRBFLX
        ADZ1(I,L) = THV(I,L) * ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1556 TRBFLX
        DZ1TMP(I,L) = ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1559 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1561 TRBFLX
        ADZ2(I,L) = 0.5 * (ADZ1(I,L)+ADZ1(I,L+1))
./pkg/fizhi/fizhi_turb.F 1574 TRBFLX
      DO L =2,NLEVM1
./pkg/fizhi/fizhi_turb.F 1576 TRBFLX
        DZ3(I,L) = ALPHA * ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1585 TRBFLX
      DO L =2,NLEV
./pkg/fizhi/fizhi_turb.F 1587 TRBFLX
        TEMP(I,L) = VK * ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1594 TRBFLX
      L = NLEV - LL
./pkg/fizhi/fizhi_turb.F 1595 TRBFLX
      LP1 = L + 1
./pkg/fizhi/fizhi_turb.F 1597 TRBFLX
      VKZE(I,L) = VKZE(I,LP1) + TEMP(I,LP1)
./pkg/fizhi/fizhi_turb.F 1600 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1602 TRBFLX
        VKZM(I,L) = VKZE(I,L) - 0.5 * TEMP(I,L+1)
./pkg/fizhi/fizhi_turb.F 1607 TRBFLX
      DO 200 L = 1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1608 TRBFLX
       LP1 = L + 1
./pkg/fizhi/fizhi_turb.F 1610 TRBFLX
       FAC1 = DPSTR(I,L) / ( DPSTR(I,L) + DPSTR(I,LP1) )
./pkg/fizhi/fizhi_turb.F 1612 TRBFLX
       RHODZ2(I,L) = FAC1 * THV(I,LP1)
./pkg/fizhi/fizhi_turb.F 1613 TRBFLX
       RHODZ2(I,L) = RHODZ2(I,L) + FAC2 * THV(I,L)
./pkg/fizhi/fizhi_turb.F 1616 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1618 TRBFLX
        RHODZ2(I,L) = (RGAS*0.01) * RHODZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1619 TRBFLX
        TEMP(I,L)   = PLKE(I,L+1) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1620 TRBFLX
        RHODZ2(I,L) = TEMP(I,L) * RHODZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1621 TRBFLX
        RHODZ2(I,L) = PLE(I,L+1) / RHODZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1622 TRBFLX
        RHOZPK(I,L) = RHODZ2(I,L) * PLKE(I,L+1)
./pkg/fizhi/fizhi_turb.F 1623 TRBFLX
        RHODZ1(I,L) = (RGAS*0.01) * THV(I,L+1)
./pkg/fizhi/fizhi_turb.F 1624 TRBFLX
        TEMP(I,L)   = PLK(I,L+1) * ADZ1(I,L+1)
./pkg/fizhi/fizhi_turb.F 1625 TRBFLX
        RHODZ1(I,L) = TEMP(I,L) * RHODZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1626 TRBFLX
        RHODZ1(I,L) = PL(I,L+1) / RHODZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1632 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1634 TRBFLX
        FLXFPK(I,L) = PLE(I,L+1) - PLE(I,L)
./pkg/fizhi/fizhi_turb.F 1635 TRBFLX
        FLXFPK(I,L) = FLXFPK(I,L) * PLK(I,L)
./pkg/fizhi/fizhi_turb.F 1636 TRBFLX
        FLXFPK(I,L) = (GRAV*DTAU*0.01) / FLXFPK(I,L)
./pkg/fizhi/fizhi_turb.F 1645 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1647 TRBFLX
        FLXFAC(I,L) = FLXFPK(I,L) * PLK(I,L)
./pkg/fizhi/fizhi_turb.F 1656 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1658 TRBFLX
        FLXFCE(I,L) = PL(I,L+1) - PL(I,L)
./pkg/fizhi/fizhi_turb.F 1659 TRBFLX
        FLXFCE(I,L) = (GRAV*DTAU*0.01) / FLXFCE(I,L)
./pkg/fizhi/fizhi_turb.F 1664 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1666 TRBFLX
        ADZ1(I,L) = 1. / ADZ1(I,L)
./pkg/fizhi/fizhi_turb.F 1669 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1671 TRBFLX
        ADZ2(I,L) = 1. / ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1679 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1681 TRBFLX
        DPK(I,L) = ( PLK(I,L+1)-PLK(I,L) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1689 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1691 TRBFLX
        Q(I,L) = 2. * QQ(I,L)
./pkg/fizhi/fizhi_turb.F 1692 TRBFLX
        Q(I,L) = SQRT( Q(I,L) )
./pkg/fizhi/fizhi_turb.F 1793 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1795 TRBFLX
        DU(I,L) = (   U(I,L)-  U(I,L+1) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1796 TRBFLX
        DV(I,L) = (   V(I,L)-  V(I,L+1) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1804 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1806 TRBFLX
         XL(I,L) = XLSAVE(I,L)
./pkg/fizhi/fizhi_turb.F 1811 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1813 TRBFLX
        DTH(I,L) = ( TH(I,L)-TH(I,L+1) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1814 TRBFLX
        DSH(I,L) = ( SH(I,L)-SH(I,L+1) ) * ADZ2(I,L)
./pkg/fizhi/fizhi_turb.F 1815 TRBFLX
        TL(I,L) = TH(I,L)*PLK(I,L)
./pkg/fizhi/fizhi_turb.F 1823 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1825 TRBFLX
        BB(I,L) = FACEPS*SHSAT(I,L)/(TL(I,L)*TL(I,L))
./pkg/fizhi/fizhi_turb.F 1826 TRBFLX
        AA(I,L) = 1. / (1. + CLH * BB(I,L) )
./pkg/fizhi/fizhi_turb.F 1828 TRBFLX
        BB(I,L) = BB(I,L) * AA(I,L)
./pkg/fizhi/fizhi_turb.F 1829 TRBFLX
        SBAR(I,L) = AA(I,L) * (SH(I,L) - SHSAT(I,L))
./pkg/fizhi/fizhi_turb.F 1840 TRBFLX
      DO L =2,NLEVM1
./pkg/fizhi/fizhi_turb.F 1843 TRBFLX
        SSDEV(I,L) = XL(I,L-1)*(AA(I,L)*DSH(I,L-1)-
./pkg/fizhi/fizhi_turb.F 1844 TRBFLX
     1               BB(I,L)*plke(I,L)*DTH(I,L-1))
./pkg/fizhi/fizhi_turb.F 1845 TRBFLX
        SSDEV(I,L) = B2 * KHSAVE(I,L-1) * SSDEV(I,L-1) * SSDEV(I,L-1)
./pkg/fizhi/fizhi_turb.F 1846 TRBFLX
        SVAR(I,L) = SQRT(SSDEV(I,L))
./pkg/fizhi/fizhi_turb.F 1848 TRBFLX
        SSDEV(I,L) = XL(I,L)*(AA(I,L)*DSH(I,L)-
./pkg/fizhi/fizhi_turb.F 1849 TRBFLX
     1               BB(I,L)*plke(I,L+1)*DTH(I,L))
./pkg/fizhi/fizhi_turb.F 1850 TRBFLX
        SSDEV(I,L) = B2 * KHSAVE(I,L) * SSDEV(I,L) * SSDEV(I,L)
./pkg/fizhi/fizhi_turb.F 1851 TRBFLX
        TEMP(I,L) = SQRT(SSDEV(I,L))
./pkg/fizhi/fizhi_turb.F 1852 TRBFLX
        SVAR(I,L) = (1./2.) * (SVAR(I,L) + TEMP(I,L))
./pkg/fizhi/fizhi_turb.F 1853 TRBFLX
        IF ( SVAR(I,L).LT.Z1PEM25) SVAR(I,L) = Z1PEM25
./pkg/fizhi/fizhi_turb.F 1856 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1858 TRBFLX
        Q1M(I,L) = SBAR(I,L) / SVAR(I,L)
./pkg/fizhi/fizhi_turb.F 1859 TRBFLX
        FCC(I,L) = (1./2.) * ( 1. + ERRF( P5SR*Q1M(I,L) ) )
./pkg/fizhi/fizhi_turb.F 1860 TRBFLX
        SHL(I,L) =  FCC(I,L) * SBAR(I,L)
./pkg/fizhi/fizhi_turb.F 1861 TRBFLX
        ARG(I,L)  = (1./2.)*Q1M(I,L)*Q1M(I,L)
./pkg/fizhi/fizhi_turb.F 1862 TRBFLX
        IF(ARG(I,L).LE.ARGMAX)
./pkg/fizhi/fizhi_turb.F 1863 TRBFLX
     1        SHL(I,L) =  SHL(I,L)+RSQ2PI*SVAR(I,L)*EXP(-ARG(I,L))
./pkg/fizhi/fizhi_turb.F 1864 TRBFLX
        BETAT(I,L) = 1. + VIRTCON * SH(I,L) - VRT1CON * SHL(I,L)
./pkg/fizhi/fizhi_turb.F 1865 TRBFLX
        BETAW(I,L) = VIRTCON *
./pkg/fizhi/fizhi_turb.F 1866 TRBFLX
     1                ( TH(I,L) + (CLH/plk(I,L)) * SHL(I,L) )
./pkg/fizhi/fizhi_turb.F 1867 TRBFLX
        BETAL(I,L) = ( 1. + VIRTCON*SH(I,L) - TWO*VRT1CON*SHL(I,L) )
./pkg/fizhi/fizhi_turb.F 1868 TRBFLX
     1     * (CLH/plke(I,L+1)) - VRT1CON * TH(I,L)
./pkg/fizhi/fizhi_turb.F 1870 TRBFLX
        BETAT1(I,L) = BETAT(I,L) -
./pkg/fizhi/fizhi_turb.F 1871 TRBFLX
     1               BB(I,L)*plk(i,L) * FCC(I,L) * BETAL(I,L)
./pkg/fizhi/fizhi_turb.F 1872 TRBFLX
        BETAW1(I,L) = BETAW(I,L) + AA(I,L) * FCC(I,L) * BETAL(I,L)
./pkg/fizhi/fizhi_turb.F 1875 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1877 TRBFLX
        DTHV(I,L) =  (1./2.)*((BETAT1(I,L)+BETAT1(I,L+1))*DTH(I,L)
./pkg/fizhi/fizhi_turb.F 1878 TRBFLX
     1               + (BETAW1(I,L)+BETAW1(I,L+1))*DSH(I,L))
./pkg/fizhi/fizhi_turb.F 1894 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 1896 TRBFLX
        STRT(I,L) = CP * DTHV(I,L) * DPK(I,L)
./pkg/fizhi/fizhi_turb.F 1897 TRBFLX
        DW2(I,L) = DU(I,L) * DU(I,L) + DV(I,L) * DV(I,L)
./pkg/fizhi/fizhi_turb.F 1898 TRBFLX
        IF ( DW2(I,L) .LE. 1.e-4 ) DW2(I,L) = 1.e-4
./pkg/fizhi/fizhi_turb.F 1899 TRBFLX
        RI(I,L) = STRT(I,L) / DW2(I,L)
./pkg/fizhi/fizhi_turb.F 1905 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1907 TRBFLX
        SRI(I,L) = RI(I,L)
./pkg/fizhi/fizhi_turb.F 1918 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1920 TRBFLX
        KH(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1921 TRBFLX
        KM(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1922 TRBFLX
        QQE(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1923 TRBFLX
        QE(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1924 TRBFLX
        P3(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1925 TRBFLX
        IBITSTB(I,L) = 0
./pkg/fizhi/fizhi_turb.F 1926 TRBFLX
        IF ( QQ(I,L) .GT. 1.e-8 ) THEN
./pkg/fizhi/fizhi_turb.F 1927 TRBFLX
         INTQ(I,L) = 1
./pkg/fizhi/fizhi_turb.F 1929 TRBFLX
         INTQ(I,L) = 0
./pkg/fizhi/fizhi_turb.F 1931 TRBFLX
        IF ( QQ(I,L).LE.1.e-8 ) THEN
./pkg/fizhi/fizhi_turb.F 1932 TRBFLX
         QQ(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1933 TRBFLX
         Q(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 1966 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 1968 TRBFLX
        QQ(I,L) = QQE(I,L)
./pkg/fizhi/fizhi_turb.F 1969 TRBFLX
        Q(I,L) = QE(I,L)
./pkg/fizhi/fizhi_turb.F 1983 TRBFLX
       DO L =LMIN,NLEVM1
./pkg/fizhi/fizhi_turb.F 1985 TRBFLX
         IF (  (IBITSTB(I,L).EQ.1)  .AND.
./pkg/fizhi/fizhi_turb.F 1986 TRBFLX
     1             ( Q(I,L) .LE. QE(I,L) ) ) THEN
./pkg/fizhi/fizhi_turb.F 1987 TRBFLX
          IBITSTB(I,L) = 1
./pkg/fizhi/fizhi_turb.F 1989 TRBFLX
          IBITSTB(I,L) = 0
./pkg/fizhi/fizhi_turb.F 1993 TRBFLX
       DO L =LMIN,NLEVM1
./pkg/fizhi/fizhi_turb.F 1995 TRBFLX
         IF(IBITSTB(I,L).EQ.1 ) THEN
./pkg/fizhi/fizhi_turb.F 1996 TRBFLX
          TEMP(I,L) = Q(I,L) / QE(I,L)
./pkg/fizhi/fizhi_turb.F 1997 TRBFLX
          KH(I,L) = TEMP(I,L) * KH(I,L)
./pkg/fizhi/fizhi_turb.F 1998 TRBFLX
          KM(I,L) = TEMP(I,L) * KM(I,L)
./pkg/fizhi/fizhi_turb.F 2000 TRBFLX
         TEMP(I,L) = 0.01 * QQE(I,L)
./pkg/fizhi/fizhi_turb.F 2001 TRBFLX
         IF((IBITSTB(I,L).EQ.1) .AND.
./pkg/fizhi/fizhi_turb.F 2002 TRBFLX
     1         ( QQ(I,L) .LE. TEMP(I,L) )) THEN
./pkg/fizhi/fizhi_turb.F 2003 TRBFLX
          QQ(I,L) = TEMP(I,L)
./pkg/fizhi/fizhi_turb.F 2004 TRBFLX
          Q(I,L) = 0.1 * QE(I,L)
./pkg/fizhi/fizhi_turb.F 2006 TRBFLX
         IF(IBITSTB(I,L).EQ.1 ) P3(I,L) = (2.*B3) *
./pkg/fizhi/fizhi_turb.F 2007 TRBFLX
     1             ( QE(I,L) - Q(I,L) )
./pkg/fizhi/fizhi_turb.F 2024 TRBFLX
       DO L =LMIN,NLEVM1
./pkg/fizhi/fizhi_turb.F 2026 TRBFLX
         P3(I,L) = P3(I,L) * DTAU / XL(I,L)
./pkg/fizhi/fizhi_turb.F 2027 TRBFLX
         TEMP(I,L) = QQE(I,L) * P3(I,L)
./pkg/fizhi/fizhi_turb.F 2028 TRBFLX
         XQ(I,L) = QQE(I,L) - QQ(I,L)
./pkg/fizhi/fizhi_turb.F 2031 TRBFLX
       DO L =LMIN,NLEVM1
./pkg/fizhi/fizhi_turb.F 2033 TRBFLX
         IF( ( (IBITSTB(I,L).EQ.1) .AND.
./pkg/fizhi/fizhi_turb.F 2034 TRBFLX
     1             ( XQ(I,L) .LT. TEMP(I,L) ) )
./pkg/fizhi/fizhi_turb.F 2036 TRBFLX
     3       ( (IBITSTB(I,L).EQ.0) .AND.
./pkg/fizhi/fizhi_turb.F 2037 TRBFLX
     4             ( XQ(I,L) .GT. TEMP(I,L) ) ) )
./pkg/fizhi/fizhi_turb.F 2038 TRBFLX
     5   P3(I,L) = XQ(I,L) / QQE(I,L)
./pkg/fizhi/fizhi_turb.F 2055 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2057 TRBFLX
        RIINIT(I,L) = RI(I,L)
./pkg/fizhi/fizhi_turb.F 2058 TRBFLX
        QQINIT(I,L) = QQ(I,L)
./pkg/fizhi/fizhi_turb.F 2064 TRBFLX
      do L = 1,nlev-1
./pkg/fizhi/fizhi_turb.F 2069 TRBFLX
         if(pl(i,L).le.150.) then
./pkg/fizhi/fizhi_turb.F 2070 TRBFLX
          qxlm(i,L) = min(qxlm(i,L),5. _d 0)
./pkg/fizhi/fizhi_turb.F 2079 TRBFLX
      DO L =LMINQ1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2081 TRBFLX
        RHOKDZ(I,L) = RHODZ1(I,L) * QXLM(I,L)
./pkg/fizhi/fizhi_turb.F 2089 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2091 TRBFLX
        KHSAVE(I,L) = KH(I,L)
./pkg/fizhi/fizhi_turb.F 2100 TRBFLX
       DO L =1,LMINQ1-1
./pkg/fizhi/fizhi_turb.F 2102 TRBFLX
         KM(I,L) = KMBG
./pkg/fizhi/fizhi_turb.F 2103 TRBFLX
         KH(I,L) = KHBG
./pkg/fizhi/fizhi_turb.F 2110 TRBFLX
      DO L =LMINQ1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2112 TRBFLX
        Q(I,L) = 2. * QQ(I,L)
./pkg/fizhi/fizhi_turb.F 2113 TRBFLX
        Q(I,L) = SQRT(Q(I,L))
./pkg/fizhi/fizhi_turb.F 2114 TRBFLX
        XQ(I,L)= XL(I,L) * Q(I,L)
./pkg/fizhi/fizhi_turb.F 2115 TRBFLX
        KM(I,L)=XQ(I,L)*KM(I,L)+KMBG
./pkg/fizhi/fizhi_turb.F 2116 TRBFLX
        KH(I,L)=XQ(I,L)*KH(I,L)+KHBG
./pkg/fizhi/fizhi_turb.F 2121 TRBFLX
      do L = 1,nlev-1
./pkg/fizhi/fizhi_turb.F 2127 TRBFLX
         if(pl(i,L).le.150.) then
./pkg/fizhi/fizhi_turb.F 2128 TRBFLX
          kh(i,L) = min(kh(i,L),5. _d 0)
./pkg/fizhi/fizhi_turb.F 2129 TRBFLX
          km(i,L) = min(km(i,L),5. _d 0)
./pkg/fizhi/fizhi_turb.F 2136 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2138 TRBFLX
        TEMP(I,L) = RHOZPK(I,L) * KH(I,L)
./pkg/fizhi/fizhi_turb.F 2139 TRBFLX
        DELTH(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2150 TRBFLX
      do L = 1,nlev
./pkg/fizhi/fizhi_turb.F 2152 TRBFLX
        dthdthg(i,L) = delth(i,L)
./pkg/fizhi/fizhi_turb.F 2155 TRBFLX
      do L = 1,nlev
./pkg/fizhi/fizhi_turb.F 2157 TRBFLX
        transth(i,L) = temp(i,L)
./pkg/fizhi/fizhi_turb.F 2161 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2163 TRBFLX
        RHOKDZ(I,L) = RHODZ2(I,L) * KH(I,L)
./pkg/fizhi/fizhi_turb.F 2164 TRBFLX
        DELSH(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 2177 TRBFLX
      do L = 1,nlev
./pkg/fizhi/fizhi_turb.F 2179 TRBFLX
        dshdshg(i,L) = delsh(i,L)
./pkg/fizhi/fizhi_turb.F 2182 TRBFLX
      do L = 1,nlev
./pkg/fizhi/fizhi_turb.F 2184 TRBFLX
        transsh(i,L) = rhokdz(i,L)
./pkg/fizhi/fizhi_turb.F 2204 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2206 TRBFLX
        RHOKDZ(I,L) = RHODZ2(I,L) * KM(I,L)
./pkg/fizhi/fizhi_turb.F 2211 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2213 TRBFLX
        WU(I,L) = WU(I,L) + RHOKDZ(I,L) * ( U(I,L+1) - U(I,L) )
./pkg/fizhi/fizhi_turb.F 2214 TRBFLX
        WV(I,L) = WV(I,L) + RHOKDZ(I,L) * ( V(I,L+1) - V(I,L) )
./pkg/fizhi/fizhi_turb.F 2217 TRBFLX
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2219 TRBFLX
        IF ( QQ(I,L) .GT. QQMIN ) THEN
./pkg/fizhi/fizhi_turb.F 2220 TRBFLX
         IBITSTB(I,L) = 1
./pkg/fizhi/fizhi_turb.F 2222 TRBFLX
         IBITSTB(I,L) = 0
./pkg/fizhi/fizhi_turb.F 2224 TRBFLX
        IF( IBITSTB(I,L).EQ.1 ) FREQDG(I,L) = FREQDG(I,L) + aitr
./pkg/fizhi/fizhi_turb.F 2232 TRBFLX
       DO L = nlev-1,1,-1
./pkg/fizhi/fizhi_turb.F 2260 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2262 TRBFLX
        EU(I,L) = EU(I,L) + AITR*KM(I,L)
./pkg/fizhi/fizhi_turb.F 2263 TRBFLX
        ET(I,L) = ET(I,L) + AITR*KH(I,L)
./pkg/fizhi/fizhi_turb.F 2273 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2275 TRBFLX
         XLDIAG(I,L) = XLDIAG(I,L) + AITR*XL(I,L)
./pkg/fizhi/fizhi_turb.F 2284 TRBFLX
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 2286 TRBFLX
         XLSAVE(I,L) = XL(I,L)
./pkg/fizhi/fizhi_turb.F 2297 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2299 TRBFLX
        turbfcc(I,L) = turbfcc(I,L) + fcc(I,L) * aitr
./pkg/fizhi/fizhi_turb.F 2300 TRBFLX
        qliq(I,L) = qliq(I,L) + shl(I,L) * aitr
./pkg/fizhi/fizhi_turb.F 2307 TRBFLX
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 2309 TRBFLX
        WU(I,L) =  WU(I,L) * AITR
./pkg/fizhi/fizhi_turb.F 2310 TRBFLX
        WV(I,L) =  WV(I,L) * AITR
./pkg/fizhi/fizhi_turb.F 250 TURBIO
      integer istnp1,istnlay,itrtrb,i,L,nn,nt
./pkg/fizhi/fizhi_turb.F 3407 TRBLEN
      integer i,L,LL
./pkg/fizhi/fizhi_turb.F 3418 TRBLEN
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 3420 TRBLEN
         STBFCN(I,L) = STRT(I,L) - RICR * DW2(I,L)
./pkg/fizhi/fizhi_turb.F 3421 TRBLEN
         INT1(I,L) = 0
./pkg/fizhi/fizhi_turb.F 3422 TRBLEN
         IF( STBFCN(I,L).LE.0. ) INT1(I,L) = 1
./pkg/fizhi/fizhi_turb.F 3426 TRBLEN
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3428 TRBLEN
        INT2(I,L) = 0
./pkg/fizhi/fizhi_turb.F 3429 TRBLEN
        IF( (INT1(I,L).EQ.1) .NEQV. (INT1(I,L+1).EQ.1) ) INT2(I,L) = 1
./pkg/fizhi/fizhi_turb.F 3444 TRBLEN
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3446 TRBLEN
        XL0(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 3464 TRBLEN
      DO L =LMIN1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3466 TRBLEN
        INT2(I,L) = 0
./pkg/fizhi/fizhi_turb.F 3467 TRBLEN
        IF ( INT1(I,L).EQ.1 .OR. INT1(I,L+1).EQ.1 ) INT2(I,L) = 1
./pkg/fizhi/fizhi_turb.F 3468 TRBLEN
        IF ( INT2(I,L).EQ.1 )
./pkg/fizhi/fizhi_turb.F 3469 TRBLEN
     &       XL0(I,L) = (0.5+DZITRP(I,L)) * DZ3(I,L+1)
./pkg/fizhi/fizhi_turb.F 3476 TRBLEN
      DO L =LMIN1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3478 TRBLEN
        IF ( INT2(I,L).EQ.1 ) THEN
./pkg/fizhi/fizhi_turb.F 3479 TRBLEN
         XL0(I,L+1) = XL0(I,L+1) + ( (0.5-DZITRP(I,L)) * DZ3(I,L+1) )
./pkg/fizhi/fizhi_turb.F 3491 TRBLEN
      DO 500 L = LMINP,NLEVM1
./pkg/fizhi/fizhi_turb.F 3492 TRBLEN
           LM1 = L-1
./pkg/fizhi/fizhi_turb.F 3494 TRBLEN
      IF( INT1(I,LM1).EQ.1 ) XL0(I,L) = XL0(I,L) + XL0(I,LM1)
./pkg/fizhi/fizhi_turb.F 3507 TRBLEN
              L   = NLEVM1 + LMIN - LL
./pkg/fizhi/fizhi_turb.F 3508 TRBLEN
              LP1 = L+1
./pkg/fizhi/fizhi_turb.F 3511 TRBLEN
      IF( INT1(I,L)  .EQ.1 ) THEN
./pkg/fizhi/fizhi_turb.F 3512 TRBLEN
      XL0(I,L) = XL0(I,LP1)
./pkg/fizhi/fizhi_turb.F 3514 TRBLEN
      XL0(I,L) = XL0(I,L) + XL0(I,LP1)
./pkg/fizhi/fizhi_turb.F 3521 TRBLEN
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 3523 TRBLEN
        IF ( XL0(I,L).LT.XL0CNV ) XL0(I,L) = XL0CNV
./pkg/fizhi/fizhi_turb.F 3535 TRBLEN
      DO L =1,LMINQ1
./pkg/fizhi/fizhi_turb.F 3537 TRBLEN
        INT2(I,L) = 1 - INT1(I,L)
./pkg/fizhi/fizhi_turb.F 3543 TRBLEN
       DO L =LMINQ,NLEVM1
./pkg/fizhi/fizhi_turb.F 3545 TRBLEN
         IF( INT1(I,L).EQ.0 ) THEN
./pkg/fizhi/fizhi_turb.F 3546 TRBLEN
           XL0(I,L) =         Q(I,L) / XL0(I,L)
./pkg/fizhi/fizhi_turb.F 3547 TRBLEN
           XL0(I,L) =       XL0(I,L) * XL0(I,L) + 1.0E-20
./pkg/fizhi/fizhi_turb.F 3548 TRBLEN
           XL0(I,L) =    STBFCN(I,L) + XL0(I,L)
./pkg/fizhi/fizhi_turb.F 3549 TRBLEN
           XL0(I,L) = SQRT( XL0(I,L) )
./pkg/fizhi/fizhi_turb.F 3550 TRBLEN
           XL0(I,L) =         Q(I,L) / XL0(I,L)
./pkg/fizhi/fizhi_turb.F 3552 TRBLEN
         INT2(I,L)  = 0
./pkg/fizhi/fizhi_turb.F 3553 TRBLEN
         IF( XL0(I,L).LT.XL0MIN ) INT2(I,L)  = 1
./pkg/fizhi/fizhi_turb.F 3561 TRBLEN
       DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3563 TRBLEN
         INT2(I,L) = 1 - INT1(I,L)
./pkg/fizhi/fizhi_turb.F 3564 TRBLEN
         IF ( INT2(I,L).EQ.1 ) XL0(I,L) = XL0MIN
./pkg/fizhi/fizhi_turb.F 3568 TRBLEN
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3570 TRBLEN
        IF ( INT2(I,L).EQ.1 ) XL0(I,L) = XL0MIN
./pkg/fizhi/fizhi_turb.F 3580 TRBLEN
      DO L =1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3582 TRBLEN
        XL(I,L) = XL0(I,L) * VKZE(I,L) / ( XL0(I,L)+VKZE(I,L) )
./pkg/fizhi/fizhi_turb.F 3592 TRBLEN
      DO L =LMINQ1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3594 TRBLEN
        Q1(I,L) = Q(I,L)
./pkg/fizhi/fizhi_turb.F 3595 TRBLEN
        INT1(I,L) = 0
./pkg/fizhi/fizhi_turb.F 3596 TRBLEN
        IF (    Q(I,L).LE.Q(I,L+1) ) INT1(I,L) = 1
./pkg/fizhi/fizhi_turb.F 3597 TRBLEN
        IF ( INT1(I,L).EQ.1 )  THEN
./pkg/fizhi/fizhi_turb.F 3598 TRBLEN
           XL0(I,L) = XL0(I,L+1)
./pkg/fizhi/fizhi_turb.F 3599 TRBLEN
            Q1(I,L) =   Q(I,L+1)
./pkg/fizhi/fizhi_turb.F 3604 TRBLEN
      DO L =LMINQ1,NLEVM1
./pkg/fizhi/fizhi_turb.F 3606 TRBLEN
        QXLM(I,L) =   XL0(I,L)*VKZM(I,L)
./pkg/fizhi/fizhi_turb.F 3607 TRBLEN
     &                 / ( XL0(I,L)+VKZM(I,L) )
./pkg/fizhi/fizhi_turb.F 3608 TRBLEN
        QXLM(I,L) = CLMT53 * Q1(I,L)*QXLM(I,L)
./pkg/fizhi/fizhi_turb.F 3662 TRBITP
      integer I,L
./pkg/fizhi/fizhi_turb.F 3670 TRBITP
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 3672 TRBITP
        DZITRP(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 3675 TRBITP
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 3677 TRBITP
        IF ( INTCHG(I,L).EQ.1 ) THEN
./pkg/fizhi/fizhi_turb.F 3678 TRBITP
         DDD(I,L) = (         CP  *(DTHV(I,L+1)*DPK(I,L)
./pkg/fizhi/fizhi_turb.F 3679 TRBITP
     &                            + DTHV(I,L)*DPK(I,L+1)) )
./pkg/fizhi/fizhi_turb.F 3680 TRBITP
     &            - ( (2.*RICR)  * ( DU(I,L+1)* DU(I,L)
./pkg/fizhi/fizhi_turb.F 3681 TRBITP
     &                            +   DV(I,L+1)* DV(I,L)) )
./pkg/fizhi/fizhi_turb.F 3682 TRBITP
         AAA(I,L) = STBFCN(I,L) + STBFCN(I,L+1)
./pkg/fizhi/fizhi_turb.F 3683 TRBITP
         BBB(I,L) = STBFCN(I,L) - STBFCN(I,L+1)
./pkg/fizhi/fizhi_turb.F 3684 TRBITP
         CCC(I,L) =            1. / BBB(I,L)
./pkg/fizhi/fizhi_turb.F 3685 TRBITP
         DZITRP(I,L) = AAA(I,L) * CCC(I,L)
./pkg/fizhi/fizhi_turb.F 3686 TRBITP
         AAA(I,L) = AAA(I,L) - DDD(I,L)
./pkg/fizhi/fizhi_turb.F 3687 TRBITP
         DDD(I,L) =        (   DDD(I,L) *    DDD(I,L) )
./pkg/fizhi/fizhi_turb.F 3688 TRBITP
     &            - 4. * (STBFCN(I,L+1) * STBFCN(I,L) )
./pkg/fizhi/fizhi_turb.F 3689 TRBITP
         DDD(I,L) = DDD(I,L)*CCC(I,L)*CCC(I,L)
./pkg/fizhi/fizhi_turb.F 3690 TRBITP
         DDD(I,L) =    SQRT( DDD(I,L) )
./pkg/fizhi/fizhi_turb.F 3693 TRBITP
        IF( INTCHG(I,L).EQ.1 .AND. AAA(I,L).NE.0. ) THEN
./pkg/fizhi/fizhi_turb.F 3694 TRBITP
         DZITRP(I,L) = ( BBB(I,L)*(1.-DDD(I,L)) ) / AAA(I,L)
./pkg/fizhi/fizhi_turb.F 3697 TRBITP
        DZITRP(I,L) = 0.5 * DZITRP(I,L)
./pkg/fizhi/fizhi_turb.F 3938 TRBDIF
      integer I,L
./pkg/fizhi/fizhi_turb.F 3949 TRBDIF
      DO L =1,NLEVX
./pkg/fizhi/fizhi_turb.F 3951 TRBDIF
        CC(I,L+1) = RHOKDZ(I,L) * FLXFAC(I,L+1)
./pkg/fizhi/fizhi_turb.F 3954 TRBDIF
      DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 3956 TRBDIF
        BB(I,L) =   RHOKDZ(I,L) * FLXFAC(I,L)
./pkg/fizhi/fizhi_turb.F 3957 TRBDIF
        AA(I,L) = 1. + CC(I,L) + BB(I,L)
./pkg/fizhi/fizhi_turb.F 3963 TRBDIF
       DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 3965 TRBDIF
         AA(I,L) = AA(I,L) - XX2(I,L)
./pkg/fizhi/fizhi_turb.F 3993 TRBDIF
       DO L =1,NLEV
./pkg/fizhi/fizhi_turb.F 3995 TRBDIF
         IF ( XX1(I,L).LT.EPSL ) XX1(I,L) = 0.
./pkg/fizhi/fizhi_turb.F 4012 VTRI0
      integer i,L,Lm1
./pkg/fizhi/fizhi_turb.F 4018 VTRI0
      DO 100 L = 2,K
./pkg/fizhi/fizhi_turb.F 4019 VTRI0
       LM1 = L - 1
./pkg/fizhi/fizhi_turb.F 4021 VTRI0
        C(I,L) = C(I,L) * A(I,LM1)
./pkg/fizhi/fizhi_turb.F 4022 VTRI0
        A(I,L) = 1. / ( A(I,L) - B(I,LM1) * C(I,L) )
./pkg/fizhi/fizhi_turb.F 4023 VTRI0
        F(I,L) = F(I,L) + F(I,LM1) * C(I,L)
./pkg/fizhi/fizhi_turb.F 4027 VTRI0
      DO 200 L = K,1,-1
./pkg/fizhi/fizhi_turb.F 4029 VTRI0
       Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) * A(I,L)
./pkg/fizhi/fizhi_turb.F 4044 VTRI1
      integer i,L
./pkg/fizhi/fizhi_turb.F 4046 VTRI1
      DO 200 L = K,1,-1
./pkg/fizhi/fizhi_turb.F 4048 VTRI1
        Y(I,L) = B(I,L) * Y(I,L+1) * A(I,L)
./pkg/fizhi/fizhi_turb.F 4064 VTRI2
      integer i,L
./pkg/fizhi/fizhi_turb.F 4066 VTRI2
      DO 100 L = 2,K
./pkg/fizhi/fizhi_turb.F 4068 VTRI2
        F(I,L) = F(I,L) + F(I,L-1) * C(I,L)
./pkg/fizhi/fizhi_turb.F 4072 VTRI2
      DO 200 L = K,1,-1
./pkg/fizhi/fizhi_turb.F 4074 VTRI2
        Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) *  A(I,L)
./pkg/fizhi/fizhi_turb.F 552 TURBIO
      do L = 1,nlay+1
./pkg/fizhi/fizhi_turb.F 554 TURBIO
      thv(i,L) = 1. + virtcon * sh(i,L)
./pkg/fizhi/fizhi_turb.F 555 TURBIO
      thv(i,L) = th(i,L) * thv(i,L)
./pkg/fizhi/fizhi_turb.F 563 TURBIO
      do L =1,nlay
./pkg/fizhi/fizhi_turb.F 565 TURBIO
      qliq(i,L) = 0.
./pkg/fizhi/fizhi_turb.F 566 TURBIO
      turbfcc(i,L) = 0.
./pkg/fizhi/fizhi_turb.F 854 TURBIO
      do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 856 TURBIO
      th(i,L) = th(i,L) + dthdthg(i,L)*(tc(i)-savetc(i))/pke(i,nlay+1)
./pkg/fizhi/fizhi_turb.F 863 TURBIO
      do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 865 TURBIO
      sh(i,L) = sh(i,L) + dshdshg(i,L)*(qa(i)-saveqa(i))
./pkg/fizhi/fizhi_turb.F 869 TURBIO
      do L = 1,nlay
./pkg/fizhi/fizhi_turb.F 871 TURBIO
       sttflux(i,L) = transth(i,L) * (th(i,L+1)-th(i,L))
./pkg/fizhi/fizhi_turb.F 872 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 289 FIZHI_WRAPPER
         do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 290 FIZHI_WRAPPER
          tempphy(i,j,Nrphys+1-L,bi,bj)=duphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 298 FIZHI_WRAPPER
         do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 299 FIZHI_WRAPPER
          tempphy(i,j,Nrphys+1-L,bi,bj)=dvphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 307 FIZHI_WRAPPER
         do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 308 FIZHI_WRAPPER
          tempphy(i,j,Nrphys+1-L,bi,bj)=dthphy(i,j,L,bi,bj)
./pkg/fizhi/fizhi_wrapper.F 316 FIZHI_WRAPPER
         do L = 1,Nrphys
./pkg/fizhi/fizhi_wrapper.F 317 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)
./pkg/fizhi/step_fizhi_corr.F 472 STEP_FIZHI_CORR
       do L=1,Nr
./pkg/fizhi/step_fizhi_corr.F 475 STEP_FIZHI_CORR
        sdyntemp(i,j,L,bi,bj) =
./pkg/fizhi/step_fizhi_corr.F 476 STEP_FIZHI_CORR
     .    (salt(i,j,L,bi,bj)-sdynbef(i,j,L,bi,bj))*dtinv -
./pkg/fizhi/step_fizhi_corr.F 477 STEP_FIZHI_CORR
     .                                        gsphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 487 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 490 STEP_FIZHI_CORR
         sphytemp(i,j,Nrphys+1-L,bi,bj) = tempphy(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 495 STEP_FIZHI_CORR
       do L = 1,Nrphys
./pkg/fizhi/step_fizhi_corr.F 498 STEP_FIZHI_CORR
        sphytemp(i,j,L,bi,bj)=
./pkg/fizhi/step_fizhi_corr.F 499 STEP_FIZHI_CORR
     .    ( (sphy(i,j,L,bi,bj)-sbef(i,j,L,bi,bj))*dtinv
./pkg/fizhi/step_fizhi_corr.F 500 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 109 MAKE_PHYS_GRID
       do L = 1,Nr
./pkg/gridalt/make_phys_grid.F 110 MAKE_PHYS_GRID
        dpstar_accum = dpstar_accum + drF(L)
./pkg/gridalt/make_phys_grid.F 172 MAKE_PHYS_GRID
        do L = Lbotij+1,Nr
./pkg/gridalt/make_phys_grid.F 178 MAKE_PHYS_GRID
           if((hfacC(i,j,L,bi,bj)*drF(L)).ge.deltap+dptry(Lnew))then
./pkg/gridalt/make_phys_grid.F 187 MAKE_PHYS_GRID
            if( (dptry(Lnew-1)+(hfacC(i,j,L,bi,bj)*drF(L)-deltap))
./pkg/gridalt/make_phys_grid.F 192 MAKE_PHYS_GRID
     &                             (hfacC(i,j,L,bi,bj)*drF(L))-deltap
./pkg/gridalt/make_phys_grid.F 193 MAKE_PHYS_GRID
             deltap = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 197 MAKE_PHYS_GRID
     &                      (hfacC(i,j,L,bi,bj)*drF(L)-deltap)
./pkg/gridalt/make_phys_grid.F 198 MAKE_PHYS_GRID
             deltap = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 207 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 209 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 217 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 219 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 222 MAKE_PHYS_GRID
           dpphys(i,j,nlphys,bi,bj) = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 225 MAKE_PHYS_GRID
           dpphys(i,j,nlphys,bi,bj) = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 231 MAKE_PHYS_GRID
          dpphys(i,j,nlphys,bi,bj) = hfacC(i,j,L,bi,bj)*drF(L)
./pkg/gridalt/make_phys_grid.F 233 MAKE_PHYS_GRID
         nlperdyn(i,j,L,bi,bj) = nlphys
./pkg/gridalt/make_phys_grid.F 56 MAKE_PHYS_GRID
       integer i,j,L,Lbotij,Lnew
./pkg/gridalt/make_phys_grid.F 78 MAKE_PHYS_GRID
        do L = 1,ntry
./pkg/gridalt/make_phys_grid.F 79 MAKE_PHYS_GRID
         dptry(L) = dptry10(L)
./pkg/gridalt/make_phys_grid.F 84 MAKE_PHYS_GRID
        do L = 1,ntry
./pkg/gridalt/make_phys_grid.F 85 MAKE_PHYS_GRID
         dptry(L) = dptry40(L)
./pkg/gridalt/make_phys_grid.F 95 MAKE_PHYS_GRID
       do L = 1,Nr
./pkg/gridalt/make_phys_grid.F 98 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/openad/ad_s_ifnblnk.F 26 OAD_S_IFNBLNK
      INTEGER L, LS
./pkg/openad/ad_s_ifnblnk.F 31 OAD_S_IFNBLNK
      L=1
./pkg/openad/ad_s_ifnblnk.F 32 OAD_S_IFNBLNK
      DO while ((L.le.LS).and.(string(L:L) .eq. ' '))
./pkg/openad/ad_s_ifnblnk.F 33 OAD_S_IFNBLNK
        L=L+1
./pkg/openad/ad_s_ifnblnk.F 35 OAD_S_IFNBLNK
      index = L
./pkg/openad/ad_s_ilnblnk.F 26 OAD_S_ILNBLNK
      INTEGER L, LS
./pkg/openad/ad_s_ilnblnk.F 31 OAD_S_ILNBLNK
      L=LS
./pkg/openad/ad_s_ilnblnk.F 32 OAD_S_ILNBLNK
      DO while ((L.gt.0).and.(string(L:L) .eq. ' '))
./pkg/openad/ad_s_ilnblnk.F 33 OAD_S_ILNBLNK
        L=L-1
./pkg/openad/ad_s_ilnblnk.F 35 OAD_S_ILNBLNK
      index = 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/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 131 DO_FIZHI
      INTEGER i,j,L
./verification/fizhi-gridalt-hs/code/do_fizhi.F 170 DO_FIZHI
       do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/do_fizhi.F 175 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 178 DO_FIZHI
         duphy(i,j,L,bi,bj)= -kV*uphy(i,j,L,bi,bj)
./verification/fizhi-gridalt-hs/code/do_fizhi.F 179 DO_FIZHI
         dvphy(i,j,L,bi,bj)= -kV*vphy(i,j,L,bi,bj)
./verification/fizhi-gridalt-hs/code/do_fizhi.F 184 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 195 DO_FIZHI
          dthphy(i,j,L,bi,bj)=- kT*( thphy(I,J,L,bi,bj)-thetaEq )
./verification/fizhi-gridalt-hs/code/do_fizhi.F 197 DO_FIZHI
          dthphy(i,j,L,bi,bj)=0.
./verification/fizhi-gridalt-hs/code/do_fizhi.F 202 DO_FIZHI
         dsphy(i,j,L,bi,bj)=0.
./verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F 107 FIZHI_INIT_FIXED
       do L = 2,Nrphys+1
./verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F 108 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 112 FIZHI_INIT_FIXED
       do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_fixed.F 113 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 252 FIZHI_INIT_VARS
        do L = 1,Nrphys
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 255 FIZHI_INIT_VARS
         swlz(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 256 FIZHI_INIT_VARS
         lwlz(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 257 FIZHI_INIT_VARS
         qliqavesw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 258 FIZHI_INIT_VARS
         qliqavelw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 259 FIZHI_INIT_VARS
         fccavesw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 260 FIZHI_INIT_VARS
         fccavelw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 261 FIZHI_INIT_VARS
         cldtot_sw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 262 FIZHI_INIT_VARS
         cldras_sw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 263 FIZHI_INIT_VARS
         cldlsp_sw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 264 FIZHI_INIT_VARS
         cldtot_lw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 265 FIZHI_INIT_VARS
         cldras_lw(i,j,L,bi,bj) = 0.
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 266 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
./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
./pkg/cheapaml/cheapaml_coare3_flux.F 196 CHEAPAML_COARE3_FLUX
     &    ' L,zu,zL,zt    =', L, zu, zL, zt
./pkg/cheapaml/cheapaml_coare3_flux.F 211 CHEAPAML_COARE3_FLUX
       L=zu/zL
./pkg/cheapaml/cheapaml_coare3_flux.F 212 CHEAPAML_COARE3_FLUX
       psu=psiu(zu/L)
./pkg/cheapaml/cheapaml_coare3_flux.F 213 CHEAPAML_COARE3_FLUX
       pst=psit(zt/L)
./pkg/cheapaml/cheapaml_coare3_flux.F 214 CHEAPAML_COARE3_FLUX
       usr=Du*xkar/(LOG(zu/zo)-psiu(zu/L))
./pkg/cheapaml/cheapaml_coare3_flux.F 215 CHEAPAML_COARE3_FLUX
       tsr=-(Dt)*xkar/(LOG(zt/zot)-psit(zt/L))
./pkg/cheapaml/cheapaml_coare3_flux.F 216 CHEAPAML_COARE3_FLUX
       qsr=-(Dq)*xkar/(LOG(zq/zoq)-psit(zq/L))
./pkg/cheapaml/cheapaml_coare3_flux.F 234 CHEAPAML_COARE3_FLUX
      q100=qs+qsr*(LOG(100. _d 0/zoq)-psit(100. _d 0/L))
./pkg/cheapaml/cheapaml_coare3_flux.F 245 CHEAPAML_COARE3_FLUX
      dEvapdTs = rhoa*usr*( xkar/(LOG(zq/zoq)-psit(zq/L)) )
./pkg/cheapaml/cheapaml_coare3_flux.F 57 CHEAPAML_COARE3_FLUX
      _RL tau,L,psu,pst,Bf
./eesupp/src/utils.F 106 DATE
      INTEGER L, LS
./eesupp/src/utils.F 111 DATE
      DO 10 L = 1, LS
./eesupp/src/utils.F 112 DATE
       IF ( string(L:L) .EQ. ' ' ) GOTO 10
./eesupp/src/utils.F 113 DATE
        IFNBLNK = L
./eesupp/src/utils.F 141 DATE
      INTEGER L, LS
./eesupp/src/utils.F 147 DATE
      DO 10 L = LS, 1, -1
./eesupp/src/utils.F 148 DATE
        IF ( string(L:L) .EQ. ' ' ) GOTO 10
./eesupp/src/utils.F 149 DATE
         ILNBLNK = L
./eesupp/src/utils.F 212 LCASE
      INTEGER   I, L
./eesupp/src/utils.F 216 LCASE
        L = INDEX(UPPER,string(I:I))
./eesupp/src/utils.F 217 LCASE
        IF ( L .EQ. 0 ) GOTO 10
./eesupp/src/utils.F 218 LCASE
         string(I:I) = LOWER(L:L)
./eesupp/src/utils.F 330 UCASE
      INTEGER   I, L
./eesupp/src/utils.F 334 UCASE
        L = INDEX(LOWER,string(I:I))
./eesupp/src/utils.F 335 UCASE
        IF ( L .EQ. 0 ) GOTO 10
./eesupp/src/utils.F 336 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)