parameter ( nidow = 2*nchklev_1*nchklev_2*nchklev_3 )
& + (ilev_3 - 1)*nchklev_2*nchklev_1
max_lev3=nbt_in/(nchklev_1*nchklev_2)+1
do ilev_2 = 1,nchklev_2
integer nchklev_2
parameter( nchklev_2 = 15 )
parameter( nchklev_2 = 120 )
parameter( nchklev_2 = 10 )
parameter( nchklev_2 = 22 )
parameter( nchklev_2 = 2 )
parameter( nchklev_2 = 50 )
parameter( nchklev_2 = 3 )
parameter( nchklev_2 = 30 )
parameter( nchklev_2 = 8 )
parameter( nchklev_2 = 100 )
parameter( nchklev_2 = 5 )
parameter( nchklev_2 = 6 )
max_lev4=nTimeSteps/(nchklev_1*nchklev_2*nchklev_3)+1
max_lev3=nTimeSteps/(nchklev_1*nchklev_2)+1
& + (ilev_4 - 1)*nchklev_3*nchklev_2*nchklev_1
if (nchklev_1*nchklev_2 .lt. nTimeSteps) then
& nchklev_1*nchklev_2
if (nchklev_1*nchklev_2*nchklev_3*nchklev_4 .lt. nTimeSteps) then
& nchklev_1*nchklev_2*nchklev_3*nchklev_4
if (nchklev_1*nchklev_2*nchklev_3 .lt. nTimeSteps) then
& nchklev_1*nchklev_2*nchklev_3
iWh=tapeMaxCounter*MAX(nchklev_2,nchklev_3)/nWh+1