locUnitTr  = '(kg/kg) '

         locUnitTr  = '(kg/m^2)'

         ilnb = LEN(locUnitTr)

         locUnitTr = SItrUnit(iTr)(1:ilnb)

       diagUnits = locUnitTr//'        '

       diagUnits = locUnitTr//'        '

      CHARACTER*8   locUnitTr