7c7
<       use params_module,only: dz,nlonp4
---
>       use params_module,only: dz,nlonp4,spval
74a75,77
> ! Lei
>      |  qout,      ! electron heating
>      |  lossout,   ! electron heating
115a119
>       real :: bb
136a141,148
> ! 8/28/06 btf: mods by Jiuhou Lei:
>        if (abs(rlatm(i,lat))-pi/4.5 >= 0.) then
>           bb = 5.0
>         else
>           bb=abs(rlatm(i,lat))*180/pi
>           bb = 5.0 -4.3e-3*bb*bb + 6.85
>         endif 
> !
138,139c150,153
<         fed(i) = ( -5.0e+7*f107te*a(i)-4.0e+7*f107te)*1.2
<         fen(i) = fed(i)/2.
---
> ! 8/28/06 btf: mods by Jiuhou:
>         fed(i) = -4.0e+7*f107te*a(i)-2.0e+7*f107te   
>         fen(i) = fed(i)/bb
> !
393a408,410
> ! by Lei
>           qout(k,i) = qe(k,i)*qtot(k,i)*evergs
> 
503a521,523
> ! by Lei
>           lossout(k,i) =  loss_en(k,i)*(te(k,i)-tn(k,i))
>      |		        + loss_ei(k,i)*(te(k,i)-ti(k,i))
506a527,531
> !
> ! 8/28/06 btf: Jiuhou added EHEAT and ELOSS:
> !     write(*,*)'come test tei',lossout(k,i), qout(k,i)
> !     call addfsech('EHEAT'  ,' ',' ',qout    ,lon0,lon1,nk,nkm1,lat)
> !     call addfsech('ELOSS'  ,' ',' ',lossout ,lon0,lon1,nk,nkm1,lat)    
588c613,615
< ! Te is not defined at lev1 (only up to lev1-1)
---
> ! 10/4/05: put spval in top level:
>       te_out(lev1,:) = spval
> !
604a632,635
> !
> ! 10/4/05: put spval in top level:
>       ti_out(lev1,:) = spval
> !