* 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.