_RL layers_G(nLayers+1)
& layers_G, MapFact, dZZf, layers_bounds
& layers_G, layers_taveFreq, layers_diagFreq,
layers_G(k) = UNSET_RL
layers_bounds(k,1) = layers_G(k)