integer max_lev2, max_lev3
max_lev3=nbt_in/(nchklev_1*nchklev_2)+1
if(ilev_3.le.max_lev3) then
integer max_lev2, max_lev3, key_in
integer max_lev3
max_lev3=nTimeSteps/(nchklev_1*nchklev_2)+1