* p0 in TIME-GCM output: 0.0005 mb seems to be incorrect: it should be 5.e-7 mb

conversion 

      pInt(1,ilev) = p0/1000._r8 * exp(-zilev(ilev)) /100._r8 

should be *100._r8, not /100._r8


*When determining OPBot and ionBot, there seems to be an inconsistency between units again (pMid is in Pa, while the threshold values
are in hPa).

! ionBot should correspond to about 60 km or 0.5 hPa. 
! oPBot should correspond to about 130 km or 0.00001 hPa.
!
    do ilev = 1, nlev

      if (pMid(1,ilev) >= 0.00001) oPBot = ilev - 1
      if (pMid(1,ilev) >= 0.5) ionBot = ilev - 1

    enddo

And the logic seems to be incorrect. (should be <=)


* All plasma densities (Op, O2p, NOp, Np, N2p and e) are number densities in TIME-GCM output.

It doesn't seem mmrN2 is defined anywhere (neither are some other mmr's)


* Does the code have geoMagLat etc? Seems to be 0 now

* Found a mistake in formulating the UBC: C31 *(ni_0-ni_2)/2/DelZ1 + C41 ni_1 = Flux, was mistaken as ...(ni_0-ni_1)...

* C41/C31 should be defined separately, because the sin(I) ^ 2 will cancel out. Now defined c43r, and then define C41 as c43r*C31.

* The top values of Te, Ti, and Op in TIME-GCM don't seem right.

* TIME-GCM output of B (mag field) is Gauss, not nano tesla. What are the E and B units in WACCM?

* DelTpDelZ should be on midpoints (as c1), but it's not the case now.
