Time units in simulation

I am hoping that someone could give me some insight about my problem. Currently I want to simulate the evolution of a gaussian wavepacket in a time-dependent potential using the Crank-Nicolson scheme. I have already the potential and wavefunction in dimensionless units. However I have the doubt of which time and frequency units I need to use. I want to run my simulation for a total time of 200 ps, and the frequency that I use is 3 THz. However my question is, can I just put the number 3 (instead of 3^12) as my frequency in THz, and my simulation time just as 200 (which will be in ps)? I am posting my python code with comments only on the relevant part of the question.

import scipy as sp
import numpy as np
from scipy import integrate, sparse, linalg,special
import scipy.sparse.linalg
import pylab as pl
ax=-10. #initial point
bx=10. #final point
tf=200 #final time in picoseconds
nx= 500#numero de puntos en el x-grid
Lx=bx-ax #total length
delx= Lx/nx #x-step size=Lx/nx
dt= 1.e-3#T/mt # time step delta t
omega=3# frequency is 3THz
mt=int((tf/omega)/dt)#
print (' MT = ', mt)
nux= 1j*dt/(4.0*(delx**2))

gridx = sp.zeros(nx) #space grid
igridx = sp.array(range(nx)) # time grid
psi = sp.zeros(nx)
pot = sp.zeros(nx)
gridx = delx*(igridx - nx/2)

###############These are just parameters#############
sigma=0.1
alpha= 2.614
D=1
R=2.95
r_e=.9699
R_e=R-2*r_e
r_0=(1./alpha)*np.arccosh(0.5*np.exp(alpha*R_e/2.0))
x=alpha*gridx
X_e=alpha*R_e
E_barrier=D*( 1-2*np.exp(-alpha*R_e/2.) )**2

psi=(2*sp.pi*sigma**2)**(-0.25) *np.exp( -(x+alpha*r_0)**2 /(4*sigma**2) )

print(nux)
####### A+######################
Asup=sp.empty(nx,dtype=complex)
Asub=sp.empty(nx,dtype=complex)
Asup.fill(-nux)
Asub.fill(-nux)
####### matrix A-########
asup=sp.empty(nx,dtype=complex)
asub=sp.empty(nx,dtype=complex)
asup.fill(nux)
asub.fill(nux)