.py

# VENTO U / V

# criando um 'vento_dirrad' com a direcao do vento em radiano
dado['vento_dirrad'] = dado['vento_dir'] + 180
dado['vento_dirrad'] = np.deg2rad(dado['vento_dirrad'])

# decompondo a vento_velocidade em u e v
#--------------------------------------------------------------
f = dado.vento_dirrad == 0
dado.loc[f,'vento_dir_u'] = np.cos(dado.vento_dirrad[f])
dado.loc[f,'vento_dir_v'] = np.sin(dado.vento_dirrad[f])

# angulo entre 0 e 90 1Q
f = ((0< dado.vento_dirrad)&(dado.vento_dirrad <= pi/2))
dado.loc[f,'vento_dir_u'] = np.cos(dado.vento_dirrad[f])
dado.loc[f,'vento_dir_v'] = np.sin(dado.vento_dirrad[f])

# 90 <angulo<=180 2Q
f = (pi/2 < dado.vento_dirrad) & (dado.vento_dirrad<= pi)
dado.loc[f,'vento_dir_u'] = np.cos(dado.vento_dirrad[f]-(pi/2))
dado.loc[f,'vento_dir_v'] = -np.sin(dado.vento_dirrad[f]-(pi/2))

# 180 <angulo<= 270 3Q
f = ((pi < dado.vento_dirrad) & (dado.vento_dirrad <= (3*pi)/2))
dado.loc[f,'vento_dir_u'] = -np.sin(dado.vento_dirrad[f]-(pi))
dado.loc[f,'vento_dir_v'] = -np.cos(dado.vento_dirrad[f]-(pi))

# 270 < angulo <= 359
f = (((3*pi)/2< dado.vento_dirrad) & (dado.vento_dirrad <= 2*pi))
dado.loc[f,'vento_dir_u'] = -np.cos(dado.vento_dirrad[f]-(3*pi/2))
dado.loc[f,'vento_dir_v'] = np.sin(dado.vento_dirrad[f]-(3*pi/2))
#---------------------------------------------------------------

del f