C Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc C C Input file for timesproc processor (zp = -17 to +5, kmx=45): C C HISTVOL = timegcm history volume(s) (24-char mss path per volume) C DDHHMM = model day,hour,min of each desired history (3 integers per history) C IDEN = density flag C 0 -> leave species as on history (most are mass mixing ratio) C 1 -> convert species to number densities (cm3) C 2 -> convert species to number density mixing ratios C MODELHTS = 1 -> use heights from model history C MODELHTS = 0 -> (default) use tn and mean mass with g0 to calculate heights C COLOR = 1 -> make color fill contours, monochrome otherwise C LABELS = 1 or 2 for minimal or normal labelling (use LABELS=1 when C metacode must be translated at <= 9600 baud)(when LABELS=1, only C monochrome plots will be made, and cmin,cmax,cint are not used) C PLTGLB = 1 -> cylindrical equidistant plots at ZP/HTS pressures/heights C SLTGLB = 0. -> cyl equid plots at ZP/HTS with midnight at center of x-axis C 12. -> cyl equid plots at ZP/HTS with noon at center of x-axis C (comment out for no cyl equid plots with slt x-axis) C ICONT = 1 -> add continental outlines to mapped projections, otherwise no C TUV = 1 -> add wind vectors to tn plots on mapped projections, otherwise no C ZUV = 1 -> add wind vectors to z plots on all mapped projections C PUV = 1 -> add ui+vi vectors to potential plots on mapped projections C (defaults for TUV and PUV = 1) C LOG10 = 1 -> plot log10 of fields, where appropriate (e.g., densities), C otherwise no (does not apply to vertical slices) C (default LOG10 = 0) C POLLAT = perimeter latitude(s) for stereographic polar contours at ZP/HTS C (comment out for no polar plots; control hem with sign) C CENSATV = lat,lon pairs of center of sat view projection(s) C (comment out for no sat view plots) C ZP = selected pressure levels for PLTGLB, SLTGLB, POLLAT, and CENSATV C HTS = selected heights for PLTGLB, SLTGLB, POLLAT, and CENSATV C (interpolation to constant height surfaces) C (ZP or HTS may be commented out if such plots are not desired) C PLTLON = 1 -> longitude slice contours (lat on x-axis, zp/ht on y-axis) C (y-axis determined by ZP_RANGE and HT_SCALE) C (note if 999. is used with PLTLON to get zonal means and VN is a C field to be plotted, then the stream function will also be contoured) C OVPOLLAT: if PLTLON=1, OVPOLLAT=latitudes at which to plot over the C pole on x-axis (e.g., OVPOLLAT=45 -> x-axis will be 45->90->45 C LON = selected longitudes for PLTLON (use 999. to get zonal means) C SLT = selected solar local time(s) for PLTLON C PLTLAT = 1 -> latitude slice contours (lon on x-axis, zp/ht on y-axis) C (y-axis determined by ZP_RANGE and HT_SCALE) C LAT = selected latitudes for PLTLAT C PLTLOC = 1 -> line plots with zp or height on y-axis, at LOC locations C (y-axis determined by ZP_RANGE and HT_SCALE) C LOC = lat,lon pairs of locations for PLTLOC (use lon=999. for zonal means) C PLGLBM = 1 -> line plots of global means (zp/ht on y-axis) C (y-axis determined by ZP_RANGE and HT_SCALE) C AMPHASE = n (n=1,2,3, or 4) -> contour amplitude and phase for first n C wave numbers. (e.g., AMPHASE=4 will make 8 frames per field), C lat on x-axis, zp and/or ht on y-axis, depending on ZP_RANGE and HT_SCALE. C WARNING: note that if you set AMPHASE=4, set all fields to be plotted, C and set ZP_RANGE and HT_SCALE, you will get 30x8x8=1920 frames! C ZP_RANGE = minimum,maximum range (-17 to +5) of zp for contours (y-axis) C by PLTLAT, PLTLON, and line plots by PLGLBM and PLTLOC C HT_SCALE = bot,top,del = bottom, top, and delta of desired linear height C scale (km) for interpolation by contours PLTLAT, PLTLON, C and line plots by PLTLOC and PLGLBM. (e.g., 60.,450.,10.) C LOGVERT = 1 -> plot log10 of field for vertical slices (PLTLON, PLTLAT, C PLTLOC, or PLTGLBM) C YAXRIGHT = flag for right-hand y-axes when left y-axis is ZP or HT (slices) C (10/31/94: currently implemented only for lon slices (PLTLON)) C 0 -> Make no extra right-hand axes C 1 -> If left axis is ZP, add right axis in height (km), C if left axis is HT, add right axis in pressure (mb) C 2 -> If left axis is ZP, add right axes in height (km) and pressure (mb) C if left axis is HT, add right axis in pressure (mb) (same as 1) C FIELDS = fields to plot (0 or 1) C (note epot (electric potential), and ion drifts are available only on C dynamo histories -- this is checked in the code: if potential or C drifts are requested on a non-dynamo history, a warning message C is printed and plots of these fields are not made) C E5577 is greenline volume emissions (5577A) C E6300 is redline volume emissions (6300A) C EO200 is O2 (0-0) band volume emissions C Notes on FIELDS: C 1. rho is o+o2+n2; any means or ht interp is done before the sum C 2. htot ("total hydrogen") is plotted only if iden=1 C (species obtained in cm3), but is plotted in volume mixing ratios C htot = (2*h2o+h+oh+ho2+4*ch4+2*h2)/(o2+o+n2) where species are cm3 C IONVEL = flag for calculation of ion drift vel (for ui,vi,wi only) C = 1 -> calculate ExB drifts C = 2 -> calculate magnetic ExB + effects of neutral atmos C = 3 -> calculate E0xB+unvn (same as ionvel=2, but epot=0.) C WRGRID = 0 -> do not make ascii file(s) at global grid C WRGRID = 1 -> make ascii file(s) of each field at each time over the C entire global grid C WRASCII = 0 -> do not make ascii data file(s) C WRASCII = 1 -> make an ascii data file for each frame C (file names will be timesproc.fr001.dat, timesproc.fr002.dat, ...) C WRASCII = 2 -> make single ascii file containing data from each frame C (file name will be timesproc.dat) C DIRASCII = 56-char remote directory to which ascii data file(s) may be C remote copied, e.g.: C e.g., DIRASCII = 'sting.hao.ucar.edu:/d/foster/timesproc ' C CMIN, CMAX, CINT = min, max values to contour, and contour interval C (if = 1.,0.,0., these values are chosen in software) Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc C C DR1=ut0-3, DR2=ut4-7, DR3=ut8-11, DR4=ut12-15, DR5=ut16-19, DR6=ut20-23 CHISTVOL = '/ROBLE/RGR95/TSEDYD7 ' CHISTVOL = '/ROBLE/RGR94/TSFCSS7 ' CDDHHMM = 30,0,0 CHISTVOL = '/ROBLE/RGR95/TSEYDR1 ','/ROBLE/RGR95/TSEYDR2 ', C '/ROBLE/RGR95/TSEYDR3 ','/ROBLE/RGR95/TSEYDR4 ', C '/ROBLE/RGR95/TSEYDR5 ','/ROBLE/RGR95/TSEYDR6 ', C '/ROBLE/RGR95/TSEYDR7 ' CDDHHMM = 0,2,0 C 0, 0, 0, 0, 1, 0, 0, 2, 0, 0, 3, 0, 0, 4, 0, 0, 5, 0, C 0, 6, 0, 0, 7, 0, 0, 8, 0, 0, 9, 0, 0,10, 0, 0,11, 0, C 0,12, 0, 0,13, 0, 0,14, 0, 0,15, 0, 0,16, 0, 0,17, 0, C 0,18, 0, 0,19, 0, 0,20, 0, 0,21, 0, 0,22, 0, 0,23, 0, 1,0,0 HISTVOL = '/ROBLE/RGR95/TSALP05 ' DDHHMM = 282,0,0 C C IDEN = 0 -> leave species as on history (most are mass mixing ratio) C IDEN = 1 -> convert species to number densities (cm3) C IDEN = 2 -> convert species to number density mixing ratios C IDEN = 1 C C MODELHTS = 1 -> use heights from model history C MODELHTS = 0 -> (default) tn and mean mass with g0 to calculate heights MODELHTS = 1 COLOR = 0 C LABELS = 1,2 for simple,normal contouring and labelling LABELS = 1 PLTGLB = 1 C C SLTGLB = 0. or 12. for cyl equid plots with midnight or noon C at center of x-axis (may be commented out) C CSLTGLB = 12. ICONT = 0 C C TUV = 1 -> add uv vectors to tn plots on all mapped projections C ZUV = 1 -> add uv vectors to z plots on all mapped projections C PUV = 1 -> add uivi vectors to epot plots on mapped projections TUV = 0 ZUV = 0 PUV = 1 UVMAX = 0. LOG10 = 0 CPOLLAT = 27.5,-27.5 CPOLLAT = 42.5,-42.5 C C CENSATV = lat,lon of center of sat view proj (use 999 for noon longitude) CCENSATV = 65.,999., -65.,999. C C ZP and HTS are used by all plots with a mapped projection C (i.e., pltglb (cyl equi), pollat (polar), and censatv (sat view)) C (ZP or HTS can be commented out) C CZP = -10.5 HTS = 35.,40.,45.,50.,55.,60.,65.,70.,75.,80.,85.,90.,95.,100., 105.,110.,115.,120. C C Longitude slices (ht or zp on y-axis, lat on x-axis at LON longitudes: C (Can use LON=999. to get zonal means) C PLTLON = 0 LON = -180.,180. CSLT = 0. C OVPOLLAT = if PLTLON=1, OVPOLLAT=latitudes at which to plot over the C pole on x-axis (e.g., OVPOLLAT=45 -> x-axis will be 45->90->45 COVPOLLAT = 30. C C Latitude slices (ht or zp on y-axis, lon on x-axis at LAT latitudes: C PLTLAT = 0 LAT = -87.5,87.5 C C Location line plots (vertical profiles): C can use 999. for lon part of LOC to get zonal means (e.g., LOC=60.,999.) C PLTLOC = 0 LOC = 2.5,-100., 2.5,-105., 2.5,-110., 2.5,-115., 2.5,-120., 2.5,-125. PLGLBM = 0 AMPHASE = 0 C C ZP_RANGE and HT_SCALE are used by all plots with zp/ht on y-axis C (i.e., pltlat, pltlon, pltloc, and pltglbm): C (ZP_RANGE or HT_SCALE can be commented out) C ZP_RANGE = -17.,-5. CHT_SCALE = 75.,105.,1. YAXRIGHT = 0 LOGVERT = 0 C C tn un vn o2 ox n4s noz co co2 h2o h2 hox o+ ch4 o21d FIELDS = 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C no2 no o3 o1 oh ho2 h n2d ti te ne o2+ w z epot n2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, C ui vi wi o/o2 o/n2 o2/n2 rho htot e5577 e6300 eo200 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0 C IONVEL = 3 WRGRID = 0 WRASCII = 2 DIRCGM = 'vishnu.hao.ucar.edu:/d/foster/timesproc/timesproc.cgm ' CDIRASCII = 'vishnu.hao.ucar.edu:/d/foster/timesproc/timesproc.dat ' DIRASCII = 'vishnu.hao:/d/foster/timesproc/connor/tseydr.ut2.dat ' C C tn un vn o2 o n4s noz co co2 h2o h2 hox o+ ch4 o21d CCMIN = 1.,1.,1.,1.,1., 1., 1.,1., 1.,1.25e-6,1., 1.,1., 1.,1., C no2 no o3 o1 oh ho2 h n2d ti te ne o2+ w z epot n2 C 1., 1.,1., 1.,1.,1.,1.,1.,1., 1.,1.,1., 1.,1., 1., 1., C ui vi wi o/o2 o/n2 o2/n2 rho htot C 1., 1., 1., 1., 1., 1., 1., 1., C e5577 e6300 eo200 C 1., 1., 1. C C tn un vn o2 o n4s noz co co2 h2o h2 hox o+ ch4 o21d CCMAX = 0.,0.,0.,0., 0., 0., 0.,0., 0.,2.3e-6,0., 0.,0., 0.,0., C no2 no o3 o1 oh ho2 h n2d ti te ne o2+ w z epot n2 C 0., 0., 0., 0.,0.,0.,0.,0., 0.,0.,0., 0., 0.,0.,0., 0., C ui vi wi o/o2 o/n2 o2/n2 rho htot C 0., 0., 0., 0., 0., 0., 0., 0., C e5577 e6300 eo200 C 0., 0., 0. C tn un vn o2 o n4s noz co co2 h2o h2 hox o+ ch4 o21d CCINT = 0.,0., 0.,0., 0.,0., 0.,0., 0.,5.e-8,0., .2,0., 0.,0., C no2 no o3 o1 oh ho2 h n2d ti te ne o2+ w z epot n2 C 0., 0.,0.,.3, 0., 0.,0.,0.,0.,0.,0.,0., 0.,0., 0., 0., C ui vi wi o/o2 o/n2 o2/n2 rho htot C 0., 0., 0., 0., 0., 0., 0., 0., C e5577 e6300 eo200 C 0., 0., 0. C ENDOFREAD