Hi, Perhaps you've caught this already, but currently in hydro_1d.F90, there is this statement: ipres = dBaseKeyNumber('ipres') I think it should be: ipres = dBaseKeyNumber('pres') However, "ipres" is not used in the rest of the routine, so this error is inconsequential to its function. David