parameter ( nidow = 2*nchklev_1*nchklev_2*nchklev_3 )
do ilev_3 = 1,nchklev_3
integer nchklev_3
parameter( nchklev_3 = 150 )
parameter( nchklev_3 = 10 )
parameter( nchklev_3 = 60 )
parameter( nchklev_3 = 5 )
parameter( nchklev_3 = 3 )
parameter( nchklev_3 = 36 )
parameter( nchklev_3 = 48 )
parameter( nchklev_3 = 24 )
parameter( nchklev_3 = 50 )
parameter( nchklev_3 = 40 )
parameter( nchklev_3 = 80 )
parameter( nchklev_3 = 2 )
max_lev4=nTimeSteps/(nchklev_1*nchklev_2*nchklev_3)+1
& + (ilev_4 - 1)*nchklev_3*nchklev_2*nchklev_1
parameter( nchklev_3 = 20 )
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