216,218d215
< 
< ! 05/15 BA Emery: added to get mean energy from 2012 TR 
<        real :: alfk,alf6,alf21,rrot6,rrot21,rt6,rt21,rd6,rd21,rh6,rh21
328,330c325
< !      h2=2.5+0.025*max(power,55.)+0.01*min(0.,power-55.)
< ! 05/14 BA Emery: h2 based on Emery 2012 TR w limit of 4.25 deg from Fig 19b
<       h2 = min(4.25,2.87 + 0.15*plevel)
---
>       h2=2.5+0.025*max(power,55.)+0.01*min(0.,power-55.)
353,357d347
<       if (irdpot==1) then
< !  Revise arad=crad+2deg with theta0 of calccloc (offc=1.1deg and dskofc=0deg)
<         arad(isouth) = theta0(isouth)/dtr + 2.
<         arad(inorth) = theta0(inorth)/dtr + 2.
<       endif
366,380d355
< ! 05/14 BA Emery:  set alfa0=0 for Emery 2012 TR 2*alfa eq A7 w limit 4.5keV for 6,21MLT
<       alfa0 = 0.
< ! Remember 2*alfa = mean energy, so set alfk,alf6,alf21 to half the values for mean energy
<       alfk = 0.50/2.
<       alf6 = min(4.5,0.36 + 0.48 * plevel)/2.
<       alf21 = min(4.5,max(1.00,-1.75+0.69*plevel))/2.
< !  Convert MLT from hours to degrees to radians, and degrees to radians
<       rrot6 = 6.00 * h2deg * dtr
<       rrot21 = 21.00 * h2deg * dtr
<       rt6 = 7.00 * h2deg * dtr
<       rt21 = 4.00 * h2deg * dtr
<       rd6 = -4.0 * dtr
<       rd21 = 4.0 * dtr
<       rh6 = 7.0 * dtr
<       rh21 = 10.0 * dtr
389,393c364,365
< !      roth = 0.81 - 0.06 * plevel
< !      rote = 0.17 - 0.04 * plevel
< ! 05/14 BA Emery: rote,h based on Emery 2012 TR w much reduced ranges from Figs 17c,19c
<       rote = max(-2.0,0.5 - 0.15 * plevel)
<       roth = max(4.0,8.0 - 0.30 * plevel)
---
>       roth = 0.81 - 0.06 * plevel
>       rote = 0.17 - 0.04 * plevel
778,785c750,756
< ! 1/02 bf: this alfa not currently in use: 05/14 BA Emery use in tiegcm_superdarn
<          alfa(i+2) = alfk+
< !!     |      alf6*exp(-((dtheta(i)-rd6-rd6v*cos(alon(i)) )/rh6)**2)* ! no rd6v in 2012 TR
<      |      alf6*exp(-((dtheta(i)-rd6)/rh6)**2)*
<      |      exp(-(atan2(sin(alon(i)-rrot6),cos(alon(i)-rrot6))/rt6)**2)+
<      |      alf21*exp(-((dtheta(i)-rd21)/rh21)**2)*
<      |      exp(-(atan2(sin(alon(i)-rrot21),cos(alon(i)-rrot21))/
<      |      rt21)**2)
---
> ! 1/02 bf: this alfa not currently in use:
> !         alfa(i+2) = alfk+
> !    |      alf6*exp(-((dtheta(i)-rd6-rd6v*cos(alon(i)) )/rh6)**2)*
> !    |      exp(-(atan2(sin(alon(i)-rrot6),cos(alon(i)-rrot6))/rt6)**2)+
> !    |      alf21*exp(-((dtheta(i)-rd21)/rh21)**2)*
> !    |      exp(-(atan2(sin(alon(i)-rrot21),cos(alon(i)-rrot21))/
> !    |      rt21)**2)