diagUnits = trcUnits//flxUnits
CHARACTER*9 flxUnits, wUnits
flxUnits = '.Pa.m^2/s'
flxUnits = '.m^3/s '
WRITE(diagUnits,'(2A)') 'g/kg',flxUnits
WRITE(diagUnits,'(2A)') 'psu',flxUnits
CHARACTER*10 flxUnits
flxUnits = '.Pa.m^2/s '
WRITE(diagUnits,'(2A)') 'K',flxUnits
WRITE(diagUnits,'(2A)') 'degC',flxUnits