c subroutine setflab(istiegcm,iden,flab,nftot,ionvel) logical istiegcm character*56 flab(nftot),flab_mmr(nftot),flab_cm3(nftot), + flab_ndmr(nftot),lab c c Define long field labels flab, with units according to iden: c c IDEN = 0 -> leave species in mass mixing ratio c IDEN = 1 -> convert species to number densities (cm3) c IDEN = 2 -> convert to number density mixing ratio c c data labshort c + /'TN ','UN ','VN ','O2 ','O1 ', c + 'N4S ','NO ','O+ ','N2D ','TI ', c + 'TE ','NE ','O2+ ','W ','Z ', c + 'POT ','N2 ','UI ','VI ','WI ', c + 'FOF2 ','HMF2 ','RHO ','UN+VN ','UI+VI ', c + 'O/O2 ','O/N2 ','O2/N2 ','O/O2+N2 '/ c c Labels for mass mixing ratios (IDEN=0): do i=1,nftot flab_mmr(i) = + ' ' enddo write(flab_mmr(1),"('NEUTRAL TEMPERATURE (DEG K)')") write(flab_mmr(2),"('NEUTRAL ZONAL WIND (M/S)')") write(flab_mmr(3),"('NEUTRAL MERIDIONAL WIND (M/S)')") write(flab_mmr(4),"('O2 MASS MIXING RATIO')") write(flab_mmr(5),"('O MASS MIXING RATIO')") write(flab_mmr(6),"('N4S MASS MIXING RATIO')") write(flab_mmr(7),"('NO MASS MIXING RATIO')") write(flab_mmr(8),"('O+ NUMBER DENSITY (CM-3)')") write(flab_mmr(9),"('N2D MASS MIXING RATIO')") write(flab_mmr(10),"('ION TEMPERATURE (DEG K)')") write(flab_mmr(11),"('ELECTRON TEMPERATURE (DEG K)')") write(flab_mmr(12),"('ELECTRON DENSITY (CM-3)')") write(flab_mmr(13),"('O2+ NUMBER DENSITY (CM-3)')") write(flab_mmr(14),"('VERTICAL WIND (PLUS UP) (M/S)')") write(flab_mmr(15),"('HEIGHT (KM)')") write(flab_mmr(16),"('ELECTRIC POTENTIAL (VOLTS)')") write(flab_mmr(17),"('N2 MASS MIXING RATIO')") write(flab_mmr(18),"('ZONAL ION DRIFT (M/S)')") write(flab_mmr(19),"('MERIDIONAL ION DRIFT (M/S)')") write(flab_mmr(20),"('VERTICAL ION DRIFT (M/S)')") write(flab_mmr(21),"('FOF2 (MHz)')") write(flab_mmr(22),"('HMF2 (KM)')") write(flab_mmr(23),"('TOTAL DENSITY MASS MIXING RATIO', + ' (O+O2+N2)')") write(flab_mmr(24),"('UN+VN VECTORS (M/S)')") write(flab_mmr(25),"('UI+VI VECTORS (M/S)')") write(flab_mmr(26),"('MASS MIXING RATIO O/O2')") write(flab_mmr(27),"('MASS MIXING RATIO O/N2')") write(flab_mmr(28),"('MASS MIXING RATIO O2/N2')") write(flab_mmr(29),"('MASS MIXING RATIO O/(O2+N2)')") c C Labels for number densities (IDEN=1): do i=1,nftot flab_cm3(i) = + ' ' enddo write(flab_cm3(1),"('NEUTRAL TEMPERATURE (DEG K)')") write(flab_cm3(2),"('NEUTRAL ZONAL WIND (M/S)')") write(flab_cm3(3),"('NEUTRAL MERIDIONAL WIND (M/S)')") write(flab_cm3(4),"('O2 NUMBER DENSITY (CM-3)')") write(flab_cm3(5),"('O NUMBER DENSITY (CM-3)')") write(flab_cm3(6),"('N4S NUMBER DENSITY (CM-3)')") write(flab_cm3(7),"('NO NUMBER DENSITY (CM-3)')") write(flab_cm3(8),"('O+ NUMBER DENSITY (CM-3)')") write(flab_cm3(9),"('N2D NUMBER DENSITY (CM-3)')") write(flab_cm3(10),"('ION TEMPERATURE (DEG K)')") write(flab_cm3(11),"('ELECTRON TEMPERATURE (DEG K)')") write(flab_cm3(12),"('ELECTRON DENSITY (CM-3)')") write(flab_cm3(13),"('O2+ NUMBER DENSITY (CM-3)')") write(flab_cm3(14),"('VERTICAL WIND (PLUS UP) (M/S)')") write(flab_cm3(15),"('HEIGHT (KM)')") write(flab_cm3(16),"('ELECTRIC POTENTIAL (VOLTS)')") write(flab_cm3(17),"('N2 NUMBER DENSITY (CM-3)')") write(flab_cm3(18),"('ZONAL ION DRIFT (M/S)')") write(flab_cm3(19),"('MERIDIONAL ION DRIFT (M/S)')") write(flab_cm3(20),"('VERTICAL ION DRIFT (M/S)')") write(flab_cm3(21),"('FOF2 (MHz)')") write(flab_cm3(22),"('HMF2 (KM)')") write(flab_cm3(23),"('TOTAL DENSITY (O+O2+N2) (CM-3)')") write(flab_cm3(24),"('UN+VN VECTORS (M/S)')") write(flab_cm3(25),"('UI+VI VECTORS (M/S)')") write(flab_cm3(26),"('RATIO O/O2 (CM-3)')") write(flab_cm3(27),"('RATIO O/N2 (CM-3)')") write(flab_cm3(28),"('RATIO O2/N2 (CM-3)')") write(flab_cm3(29),"('RATIO O/(O2+N2) (CM-3)')") c c Define ion drift field labels according to ionvel flag: if (ionvel.le.1) then write(flab_mmr(18),"('ZONAL ION DRIFT (M/S) (ExB)')") write(flab_mmr(19),"('MERIDIONAL ION DRIFT (M/S) (ExB)')") write(flab_mmr(20),"('VERTICAL ION DRIFT (M/S) (ExB)')") c write(flab_cm3(18),"('ZONAL ION DRIFT (M/S) (ExB)')") write(flab_cm3(19),"('MERIDIONAL ION DRIFT (M/S) (ExB)')") write(flab_cm3(20),"('VERTICAL ION DRIFT (M/S) (ExB)')") c write(flab_ndmr(18),"('ZONAL ION DRIFT (M/S) (ExB)')") write(flab_ndmr(19),"('MERIDIONAL ION DRIFT (M/S) (ExB)')") write(flab_ndmr(20),"('VERTICAL ION DRIFT (M/S) (ExB)')") else write(flab_mmr(18),"('MAGNETIC ZONAL ION DRIFT (M/S) ', + '(mag ExB+unvn)')") write(flab_mmr(19),"('MAGNETIC MERIDIONAL ION DRIFT (M/S) ', + '(mag ExB+unvn)')") write(flab_mmr(20),"('MAGNETIC VERTICAL ION DRIFT (M/S) ', + '(mag ExB+unvn)')") c write(flab_cm3(18),"('MAGNETIC ZONAL ION DRIFT (M/S) ', + '(mag ExB+unvn)')") write(flab_cm3(19),"('MAGNETIC MERIDIONAL ION DRIFT (M/S) ', + '(mag ExB+unvn)')") write(flab_cm3(20),"('MAGNETIC VERTICAL ION DRIFT (M/S) ', + '(mag ExB+unvn)')") c write(flab_ndmr(18),"('MAGNETIC ZONAL ION DRIFT (M/S) ', + '(mag ExB+unvn)')") write(flab_ndmr(19),"('MAGNETIC MERIDIONAL ION DRIFT (M/S) ', + '(mag ExB+unvn)')") write(flab_ndmr(20),"('MAGNETIC VERTICAL ION DRIFT (M/S) ', + '(mag ExB+unvn)')") endif c c Labels for mass density (gm/cm3) (IDEN=2): do i=1,nftot flab_ndmr(i) = + ' ' enddo write(flab_ndmr(1),"('NEUTRAL TEMPERATURE (DEG K)')") write(flab_ndmr(2),"('NEUTRAL ZONAL WIND (M/S)')") write(flab_ndmr(3),"('NEUTRAL MERIDIONAL WIND (M/S)')") write(flab_ndmr(4),"('O2 MASS DENSITY (GM/CM3)')") write(flab_ndmr(5),"('O MASS DENSITY (GM/CM3)')") write(flab_ndmr(6),"('N4S MASS DENSITY (GM/CM3)')") write(flab_ndmr(7),"('NO MASS DENSITY (GM/CM3)')") write(flab_ndmr(8),"('O+ NUMBER DENSITY (CM-3)')") write(flab_ndmr(9),"('N2D MASS DENSITY (GM/CM3)')") write(flab_ndmr(10),"('ION TEMPERATURE (DEG K)')") write(flab_ndmr(11),"('ELECTRON TEMPERATURE (DEG K)')") write(flab_ndmr(12),"('ELECTRON DENSITY (CM-3)')") write(flab_ndmr(13),"('O2+ NUMBER DENSITY (CM-3)')") write(flab_ndmr(14),"('VERTICAL WIND (PLUS UP) (M/S)')") write(flab_ndmr(15),"('HEIGHT (KM)')") write(flab_ndmr(16),"('ELECTRIC POTENTIAL (VOLTS)')") write(flab_ndmr(17),"('N2 MASS DENSITY (GM/CM3)')") write(flab_ndmr(18),"('ZONAL ION DRIFT (M/S)')") write(flab_ndmr(19),"('MERIDIONAL ION DRIFT (M/S)')") write(flab_ndmr(20),"('VERTICAL ION DRIFT (M/S)')") write(flab_ndmr(21),"('FOF2 (MHz)')") write(flab_ndmr(22),"('HMF2 (KM)')") write(flab_ndmr(23),"('TOTAL DENSITY (O+O2+N2) (GM/CM3)')") write(flab_ndmr(24),"('UN+VN VECTORS (M/S)')") write(flab_ndmr(25),"('UI+VI VECTORS (M/S)')") write(flab_ndmr(26),"('RATIO O/O2 (GM/CM3)')") write(flab_ndmr(27),"('RATIO O/N2 (GM/CM3)')") write(flab_ndmr(28),"('RATIO O2/N2 (GM/CM3)')") write(flab_ndmr(29),"('RATIO O/(O2+N2) (GM/CM3)')") c c Define flab to be used by plotting routines: do i=1,nftot flab(i) = + ' ' enddo if (iden.eq.0) then do i=1,nftot flab(i) = flab_mmr(i) enddo elseif (iden.eq.1) then do i=1,nftot flab(i) = flab_cm3(i) enddo elseif (iden.eq.2) then do i=1,nftot flab(i) = flab_ndmr(i) enddo else write(6,"('>>> setflab: bad iden=',i3)") iden endif c c Add TIGCM or TIEGCM at beginning of each label: do i=1,nftot len = lenstr(flab(i)) if (istiegcm) then if (len+7.le.56) then lab = flab(i) flab(i) = + ' ' write(flab(i)(1:7),"('TIEGCM ')") flab(i) = flab(i)(1:7)//lab(1:len) endif else if (len+6.le.56) then lab = flab(i) flab(i) = + ' ' write(flab(i)(1:6),"('TIGCM ')") flab(i) = flab(i)(1:6)//lab(1:len) endif endif enddo c write(6,"('setflab: iden=',i2,' istiegcm=',l1,' flab=')") c + iden,istiegcm c do i=1,nftot c write(6,"(a)") flab(i) c enddo c return end