_RL MapFact(FineGridMax)
& layers_G, MapFact, dZZf
TatU = MapFact(kk) *
& (1-MapFact(kk)) *
TatV = MapFact(kk) *
MapFact(kk) = 1.0 - (( ZZc(kk) - Zc(k) ) / drC(k+1))
& ', MapFact(',kk,')=',MapFact(kk),
MapFact(kk) = 1.0 _d 0
MapFact(kk) = 0.0 _d 0