added a test script for forcing function estimation