c subroutine setflab(iden,flab,nftot) character*56 lab,flab(nftot) c write(flab(1),"('NEUTRAL TEMPERATURE (DEG K)')") write(flab(2),"('NEUTRAL ZONAL WIND (M/S)')") write(flab(3),"('NEUTRAL MERIDIONAL WIND (M/S)')") write(flab(4),"('O2')") write(flab(5),"('OX')") write(flab(6),"('N4S')") write(flab(7),"('NOZ')") write(flab(8),"('CO')") write(flab(9),"('CO2')") write(flab(10),"('H2O')") write(flab(11),"('H2')") write(flab(12),"('HOX')") write(flab(13),"('O+')") write(flab(14),"('CH4')") write(flab(15),"('O21D')") write(flab(16),"('NO2')") write(flab(17),"('NO')") write(flab(18),"('O3')") write(flab(19),"('O')") write(flab(20),"('OH')") write(flab(21),"('HO2')") write(flab(22),"('H')") write(flab(23),"('N2D')") write(flab(24),"('ION TEMPERATURE (DEG K)')") write(flab(25),"('ELECTRON TEMPERATURE (DEG K)')") write(flab(26),"('ELECTRON DENSITY (CM-3)')") write(flab(27),"('O2+ NUMBER DENSITY (CM-3)')") write(flab(28),"('VERTICAL WIND (PLUS UP) (M/S)')") write(flab(29),"('HEIGHT (KM)')") write(flab(30),"('ELECTRIC POTENTIAL (VOLTS)')") write(flab(31),"('ZONAL ION DRIFT (M/S)')") write(flab(32),"('MERIDIONAL ION DRIFT (M/S)')") write(flab(33),"('VERTICAL ION DRIFT (M/S)')") write(flab(34),"('N2')") write(flab(35),"('E5577 GREENLINE EMISSION (photons/cm3/sec)')") write(flab(36),"('E6300 REDLINE EMISSION (photons/cm3/sec)')") write(flab(37),"('O2 ATMOS (0-0) BAND EMISSION (photons/cm3/sec)' + )") c c Add TIGCM or TIEGCM at beginning of each label: do i=1,nftot len = lenstr(flab(i)) if (len+9.le.56) then lab = flab(i) flab(i) = + ' ' write(flab(i)(1:9),"('TIMESGCM ')") flab(i) = flab(i)(1:9)//lab(1:len) endif enddo c return end