_RL dZZf(FineGridMax)
INTEGER MapIndex(FineGridMax), CellIndex(FineGridMax)
_RL MapFact(FineGridMax)
INTEGER Nlayers, FineGridFact, FineGridMax, layers_maxNum
PARAMETER( FineGridMax = Nr * FineGridFact )
_RL ZZf(FineGridMax+1)
_RL ZZc(FineGridMax)
NZZ = FineGridMax