[FDTD] codes for vlf propagation on earth-ionosphere waveguide

hi there

i am currently working on fdtd programming to model vlf propagation on earth-ionosphere waveguide during daytime to detect ionospheric anomalies before major aerthquake. i have go through some paper from berenger, 2002 and cummer,2000 but still did not have complete parameters to use in the codes. does anyone here have experience wrote the code? can i use it as reference? thank you in advance
Thanks Report