30,31c30
< integer, save :: ncid = -1
< integer norbs_id, nspin_id, nnzs_id, scf_step_id
---
> integer ncid, norbs_id, nspin_id, nnzs_id, scf_step_id
103d101
<       if ( ncid > -1 ) call check( nf90_close(ncid))