Warning, stokes.F constants are defined for pure CO2.
 Cannot find required file "run.def"
   (which should contain some input parameters
    along with the following line:
    INCLUDEDEF=callphys.def
    )
  ... might as well stop here ...