Fixed bug in heat transport diagnostic method in EBM code.

climlab · Apr 9, 2015 · 12857ec · 12857ec
1 parent f45c284
commit 12857ec
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/climlab/model/ebm.py b/climlab/model/ebm.py
@@ -78,8 +78,8 @@ def global_mean_temperature(self):
 
     def inferred_heat_transport( self ):
         '''Returns the inferred heat transport (in PW) by integrating the TOA energy imbalance from pole to pole.'''
-        phi = self.phi
-        energy_in = self.diagnostics['net_radiation']
+        phi = np.deg2rad(self.lat)
+        energy_in = np.squeeze(self.diagnostics['net_radiation'])
         return ( 1E-15*2* np.math.pi*const.a**2 * 
             integrate.cumtrapz(np.cos(phi)*energy_in,
             x=phi, initial=0. ) )