parameter ( nidow = 2*nchklev_1*nchklev_2*nchklev_3 )
do ilev_1 = 1,nchklev_1
iloop = (ilev_2 - 1)*nchklev_1 + ilev_1
& + (ilev_3 - 1)*nchklev_2*nchklev_1
max_lev3=nbt_in/(nchklev_1*nchklev_2)+1
max_lev2=nbt_in/nchklev_1+1
integer nchklev_1
parameter( nchklev_1 = 3 )
parameter( nchklev_1 = 50 )
parameter( nchklev_1 = 20 )
parameter( nchklev_1 = 5 )
parameter( nchklev_1 = 96 )
parameter( nchklev_1 = 48 )
parameter( nchklev_1 = 24 )
parameter( nchklev_1 = 1 )
parameter( nchklev_1 = 6 )
parameter( nchklev_1 = 36 )
parameter( nchklev_1 = 12 )
max_lev3=nTimeSteps/(nchklev_1*nchklev_2)+1
max_lev2=nTimeSteps/nchklev_1+1
max_lev4=nTimeSteps/(nchklev_1*nchklev_2*nchklev_3)+1
& + (ilev_4 - 1)*nchklev_3*nchklev_2*nchklev_1
parameter( nchklev_1 = 30 )
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
nchklev_1/2/3