diff --git a/science/Makefile b/science/Makefile
index 00f49e72bc6e..a670fb648dff 100644
--- a/science/Makefile
+++ b/science/Makefile
@@ -172,6 +172,7 @@
     SUBDIR += nwchem-data
     SUBDIR += octopus
     SUBDIR += openbabel
+    SUBDIR += openems
     SUBDIR += openkim
     SUBDIR += openkim-models
     SUBDIR += openmx
diff --git a/science/openems/Makefile b/science/openems/Makefile
new file mode 100644
index 000000000000..083f233dcd3d
--- /dev/null
+++ b/science/openems/Makefile
@@ -0,0 +1,35 @@
+PORTNAME=	openems
+DISTVERSIONPREFIX=	v
+DISTVERSION=	0.0.35-71
+DISTVERSIONSUFFIX=	-g4c24b6e
+CATEGORIES=	science
+
+MAINTAINER=	yuri@FreeBSD.org
+COMMENT=	Electromagnetic field solver using the EC-FDTD method
+
+LICENSE=	GPLv3
+LICENSE_FILE=	${WRKSRC}/COPYING
+
+LIB_DEPENDS=	libboost_system.so:devel/boost-libs \
+		libCSXCAD.so:cad/csxcad \
+		libfparser-4.5.so:math/fparser \
+		libhdf5.so:science/hdf5 \
+		libsz.so:science/szip \
+		libtinyxml.so:textproc/tinyxml \
+		libvtkCommonCore-9.0.so:math/vtk9
+
+USES=		cmake pkgconfig
+USE_LDCONFIG=	yes
+
+USE_GITHUB=	yes
+GH_ACCOUNT=	thliebig
+GH_PROJECT=	openEMS
+
+CXXFLAGS+=	-I${LOCALBASE}/include/fparser # fparser.pc probably returns a wrong value
+
+DATADIR=	${PREFIX}/share/openEMS
+
+post-install:
+	@${REINPLACE_CMD} -i '' -e 's|#!/bin/bash|#!/bin/sh|' ${STAGEDIR}${PREFIX}/bin/openEMS.sh
+
+.include <bsd.port.mk>
diff --git a/science/openems/distinfo b/science/openems/distinfo
new file mode 100644
index 000000000000..64eb8ebc0c96
--- /dev/null
+++ b/science/openems/distinfo
@@ -0,0 +1,3 @@
+TIMESTAMP = 1629334571
+SHA256 (thliebig-openEMS-v0.0.35-71-g4c24b6e_GH0.tar.gz) = 33690c506c52d8139e8e92734a4478f323156538aafa1e9260c3a4a68e6b563e
+SIZE (thliebig-openEMS-v0.0.35-71-g4c24b6e_GH0.tar.gz) = 1726662
diff --git a/science/openems/pkg-descr b/science/openems/pkg-descr
new file mode 100644
index 000000000000..68a45a169430
--- /dev/null
+++ b/science/openems/pkg-descr
@@ -0,0 +1,8 @@
+openEMS is a free and open electromagnetic field solver using the FDTD method.
+Matlab or Octave are used as an easy and flexible scripting interface.
+
+It features:
+* fully 3D Cartesian and cylindrical coordinates graded mesh
+* multi-threading, SIMD (SSE) and MPI support for high speed FDTD
+
+WWW: https://openems.de/start/
diff --git a/science/openems/pkg-plist b/science/openems/pkg-plist
new file mode 100644
index 000000000000..a886fef35d2f
--- /dev/null
+++ b/science/openems/pkg-plist
@@ -0,0 +1,132 @@
+bin/nf2ff
+bin/openEMS
+bin/openEMS.sh
+include/openEMS/nf2ff.h
+include/openEMS/openems.h
+include/openEMS/openems_global.h
+lib/libnf2ff.so
+lib/libnf2ff.so.0
+lib/libnf2ff.so.0.1.0
+lib/libopenEMS.so
+lib/libopenEMS.so.0
+lib/libopenEMS.so.0.0.35
+%%DATADIR%%/matlab/AR_estimate.m
+%%DATADIR%%/matlab/Add2Queue.m
+%%DATADIR%%/matlab/AddCPWPort.m
+%%DATADIR%%/matlab/AddCircWaveGuidePort.m
+%%DATADIR%%/matlab/AddCoaxialPort.m
+%%DATADIR%%/matlab/AddCurvePort.m
+%%DATADIR%%/matlab/AddLumpedPort.m
+%%DATADIR%%/matlab/AddMRStub.m
+%%DATADIR%%/matlab/AddMSLPort.m
+%%DATADIR%%/matlab/AddPML.m
+%%DATADIR%%/matlab/AddRectWaveGuidePort.m
+%%DATADIR%%/matlab/AddStripLinePort.m
+%%DATADIR%%/matlab/AddWaveGuidePort.m
+%%DATADIR%%/matlab/AnalyzeNF2FF.m
+%%DATADIR%%/matlab/CalcNF2FF.m
+%%DATADIR%%/matlab/CheckQueue.m
+%%DATADIR%%/matlab/ConvertHDF5_VTK.m
+%%DATADIR%%/matlab/CreateNF2FFBox.m
+%%DATADIR%%/matlab/DFT_time2freq.m
+%%DATADIR%%/matlab/DelayFidelity.m
+%%DATADIR%%/matlab/Dump2VTK.m
+%%DATADIR%%/matlab/DumpFF2VTK.m
+%%DATADIR%%/matlab/FFT_time2freq.m
+%%DATADIR%%/matlab/FindFreeSSH.m
+%%DATADIR%%/matlab/FinishQueue.m
+%%DATADIR%%/matlab/GetField_Interpolation.m
+%%DATADIR%%/matlab/GetField_Range.m
+%%DATADIR%%/matlab/GetField_SubSampling.m
+%%DATADIR%%/matlab/GetField_TD2FD.m
+%%DATADIR%%/matlab/InitCylindricalFDTD.m
+%%DATADIR%%/matlab/InitFDTD.m
+%%DATADIR%%/matlab/InitQueue.m
+%%DATADIR%%/matlab/PlotHDF5FieldData.m
+%%DATADIR%%/matlab/ReadHDF5Attribute.m
+%%DATADIR%%/matlab/ReadHDF5Dump.m
+%%DATADIR%%/matlab/ReadHDF5FieldData.m
+%%DATADIR%%/matlab/ReadHDF5Mesh.m
+%%DATADIR%%/matlab/ReadUI.m
+%%DATADIR%%/matlab/ResultsQueue.m
+%%DATADIR%%/matlab/RunOpenEMS.m
+%%DATADIR%%/matlab/RunOpenEMS_MPI.m
+%%DATADIR%%/matlab/RunOpenEMS_Parallel.m
+%%DATADIR%%/matlab/SetBoundaryCond.m
+%%DATADIR%%/matlab/SetCustomExcite.m
+%%DATADIR%%/matlab/SetDiracExcite.m
+%%DATADIR%%/matlab/SetGaussExcite.m
+%%DATADIR%%/matlab/SetSinusExcite.m
+%%DATADIR%%/matlab/SetStepExcite.m
+%%DATADIR%%/matlab/SetupMPI.m
+%%DATADIR%%/matlab/Tutorials/Bent_Patch_Antenna.m
+%%DATADIR%%/matlab/Tutorials/CRLH_Extraction.m
+%%DATADIR%%/matlab/Tutorials/CRLH_LeakyWaveAnt.m
+%%DATADIR%%/matlab/Tutorials/Circ_Waveguide.m
+%%DATADIR%%/matlab/Tutorials/Conical_Horn_Antenna.m
+%%DATADIR%%/matlab/Tutorials/CreateCRLH.m
+%%DATADIR%%/matlab/Tutorials/CylindricalWave_CC.m
+%%DATADIR%%/matlab/Tutorials/Dipole_SAR.m
+%%DATADIR%%/matlab/Tutorials/Helical_Antenna.m
+%%DATADIR%%/matlab/Tutorials/Horn_Antenna.m
+%%DATADIR%%/matlab/Tutorials/MRI_LP_Birdcage.m
+%%DATADIR%%/matlab/Tutorials/MRI_Loop_Coil.m
+%%DATADIR%%/matlab/Tutorials/MSL_NotchFilter.m
+%%DATADIR%%/matlab/Tutorials/Parallel_Plate_Waveguide.m
+%%DATADIR%%/matlab/Tutorials/Patch_Antenna_Array.m
+%%DATADIR%%/matlab/Tutorials/Patch_Antenna_Phased_Array.m
+%%DATADIR%%/matlab/Tutorials/RCS_Sphere.m
+%%DATADIR%%/matlab/Tutorials/RadarUWBTutorial.m
+%%DATADIR%%/matlab/Tutorials/Rect_Waveguide.m
+%%DATADIR%%/matlab/Tutorials/Simple_Patch_Antenna.m
+%%DATADIR%%/matlab/Tutorials/StripLine2MSL.m
+%%DATADIR%%/matlab/Tutorials/readme
+%%DATADIR%%/matlab/WriteHDF5.m
+%%DATADIR%%/matlab/WriteOpenEMS.m
+%%DATADIR%%/matlab/calcLumpedPort.m
+%%DATADIR%%/matlab/calcPort.m
+%%DATADIR%%/matlab/calcTLPort.m
+%%DATADIR%%/matlab/calcWGPort.m
+%%DATADIR%%/matlab/calc_ypar.m
+%%DATADIR%%/matlab/examples/__deprecated__/MSL2.m
+%%DATADIR%%/matlab/examples/antennas/Bi_Quad_Antenna.m
+%%DATADIR%%/matlab/examples/antennas/Patch_Antenna.m
+%%DATADIR%%/matlab/examples/antennas/Patch_Antenna_Array.m
+%%DATADIR%%/matlab/examples/antennas/infDipol.m
+%%DATADIR%%/matlab/examples/antennas/inverted_f.m
+%%DATADIR%%/matlab/examples/optimizer/optimizer_asco.m
+%%DATADIR%%/matlab/examples/optimizer/optimizer_simfun.m
+%%DATADIR%%/matlab/examples/other/Helix.m
+%%DATADIR%%/matlab/examples/other/LumpedElement.m
+%%DATADIR%%/matlab/examples/other/Metamaterial_PlaneWave_Drude.m
+%%DATADIR%%/matlab/examples/other/PML_reflection_analysis.m
+%%DATADIR%%/matlab/examples/other/PlaneWave.m
+%%DATADIR%%/matlab/examples/other/gauss_excitation_test.m
+%%DATADIR%%/matlab/examples/other/resistance_sheet.m
+%%DATADIR%%/matlab/examples/transmission_lines/CPW_Line.m
+%%DATADIR%%/matlab/examples/transmission_lines/Finite_Stripline.m
+%%DATADIR%%/matlab/examples/transmission_lines/MSL.m
+%%DATADIR%%/matlab/examples/transmission_lines/MSL_Losses.m
+%%DATADIR%%/matlab/examples/transmission_lines/Stripline.m
+%%DATADIR%%/matlab/examples/transmission_lines/directional_coupler.m
+%%DATADIR%%/matlab/examples/waveguide/Circ_Waveguide.m
+%%DATADIR%%/matlab/examples/waveguide/Circ_Waveguide_CylinderCoords.m
+%%DATADIR%%/matlab/examples/waveguide/Coax.m
+%%DATADIR%%/matlab/examples/waveguide/Coax_CylinderCoords.m
+%%DATADIR%%/matlab/examples/waveguide/Coax_Cylindrical_MG.m
+%%DATADIR%%/matlab/examples/waveguide/Rect_Waveguide.m
+%%DATADIR%%/matlab/h5readatt_octave.cc
+%%DATADIR%%/matlab/harminv.m
+%%DATADIR%%/matlab/optimize.m
+%%DATADIR%%/matlab/optimizer_asco_sim.m
+%%DATADIR%%/matlab/physical_constants.m
+%%DATADIR%%/matlab/plotFF3D.m
+%%DATADIR%%/matlab/plotFFdB.m
+%%DATADIR%%/matlab/plotRefl.m
+%%DATADIR%%/matlab/polarFF.m
+%%DATADIR%%/matlab/private/ReadNF2FF.m
+%%DATADIR%%/matlab/private/invoke_openEMS.m
+%%DATADIR%%/matlab/queue_addProcess.m
+%%DATADIR%%/matlab/queue_checkProcess.m
+%%DATADIR%%/matlab/queue_delProcess.m
+%%DATADIR%%/matlab/setup.m