@@ -237,10 +237,6 @@ def _windfield_holland(track, centroids, model='H08'):
237237 intensity [0 , close_centr ] = np .maximum (
238238 intensity [0 , close_centr ].todense (), v_full )
239239
240- # keep maximum instantaneous wind
241- intensity [0 , close_centr ] = np .maximum (
242- intensity [0 , close_centr ].todense (), v_full )
243-
244240 return intensity
245241
246242def _extra_rad_max_wind (track , ureg ):
@@ -355,13 +351,13 @@ def _bs_value(v_trans, penv, pcen, prepcen, lat, hol_xx, tint):
355351 """ Halland's 2008 b value computation.
356352
357353 Parameters:
358- v_trans (float): translational wind in m/s
359- penv (float): environmental pressure
360- pcen (float): central pressure
361- prepcen (float): previous central pressure
362- lat (float): latitude
354+ v_trans (float): translational wind ( m/s)
355+ penv (float): environmental pressure (hPa)
356+ pcen (float): central pressure (hPa)
357+ prepcen (float): previous central pressure (hPa)
358+ lat (float): latitude (degrees)
363359 hol_xx (float): Holland's xx value
364- tint (float): time step
360+ tint (float): time step (h)
365361
366362 Returns:
367363 float
@@ -418,8 +414,8 @@ def _vang_holland(track, i_node, r_arr, v_trans, model='H08'):
418414 pre_pcen = track .central_pressure .values [i_node - 1 ]
419415 if pre_pcen < 850 :
420416 pre_pcen = track .central_pressure .values [i_node ]
421- hol_b = _bs_value (v_trans , penv , pcen , pre_pcen , \
422- ycoord , hol_xx , track .time_step .values [i_node ])
417+ hol_b = _bs_value (v_trans , penv , pcen , pre_pcen ,
418+ ycoord , hol_xx , track .time_step .values [i_node ])
423419 else :
424420 # TODO H80: b=b_value(v_trans,vmax,penv,pcen,rho);
425421 raise NotImplementedError
0 commit comments