C Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc C C INPUT FILE FOR TIGCMPROC CRAY PROCESSOR (zp = -7 to +5, kmx=25): C (binary executable is on mss /FOSTER/bin/tigcmproc) C (source is on sting:~foster/tigcmproc) C C This is a "snapshot" processor, i.e., all plots are at a selected C ut (from a selected history). Multiple histories may be specified C in a single run. C C This processor will accept tigcm or tiegcm histories, however: C If electric potential or ion drift plots are requested (pot,ui,vi,wi), C and the history is tigcm (not tiegcm), a notice will be printed C to stdout (unit 6), and potential and/or ion drifts will not be plotted. C C All values below DDHHMM except FIELDS may be commented out to get the C default value, or to suppress the plotting option. Always comment out C in column one with an upper case C. C All uncommented keyword = value lexical read lines must start in C column two and be <= 72 columns per line. (you can use the ruler at C CMIN,CMAX,CINT for reference). New lines may be used if a comma C appears at the end of the previous line. C Also note that tabs cannot not be used for white space in lexical reads C C HISTVOL = tigcm or tiegcm history volume(s) (24-char mss path per volume) C e.g., '/ROBLE/RGR92/LTCS01 ' C DDHHMM = model day,hour,min of each desired history C (3 comma-separated integers per history, e.g., 3,0,0, 3,1,0) C (can use 99 as one of the integers to determine histories on volume; C 99 will not be found, and times found will be printed to stdout) C (HISTVOL and DDHHMM are required inputs) C IDEN = 0 -> leave species as on history (most are mass mixing ratios) C IDEN = 1 -> convert species to number densities (cm3) C IDEN = 2 -> convert species to mass density (gm/cm3) C (except ne, o+, and o2+ which remain number density) C (default IDEN = 1) C Note o+, ne, and o2+ are number densities on the history, and will C be left as such (cm-3) even when IDEN = 0 or 2 C COLOR = 1 -> make color fill contours, monochrome otherwise C (default COLOR = 0) C C --------------- CONTOURS ON MAPPED PROJECTIONS ---------------- C C Note all lat,lon values and ZP and HTS are floating point values, C whereas ICONT, TUV, PUV, and LOG10 are integers C CENGLB = lat,lon of center of cylindrical equidistant projection(s) C at ZP/HTS. Can use CENGLB=0.,999. to get local noon longitude at C center of x-axis. For "normal" global plots use CENGLB = 0.,0. C X-axis labels are adjusted for lon.ne.0. For lat.ne.0., the map itself C is labelled (NP,SP,GM,ID) and gridded, and map axes are not labelled. C Note extra slt x-axis rounds local time values so may be inaccurate C for certain center longitudes C (Comment out for no cylindrical equidistant plots) C CENMOLW = lat,lon of center of mollweide projection(s) at ZP/HTS C (vectors (TUV and PUV) not recommended with this projection) C (use 999. as lon to get local noon) (comment out for no molw) C CENSATV = lat,lon of center of sat view projection(s) at ZP/HTS C (use 999. as lon to get local noon) (comment out for no satv) C POLLAT = perimeter latitudes(s) for stereographic (polar) projections C (comment out for no polar plots) (this is a local time dial) C Note: nearest grid points to above lat,lon values will be used C (i.e., no interpolation between grid points) C ICONT = 1 -> add continental outlines to mapped projections, otherwise no C TUV = 1 -> add un+vn vectors to tn plots on mapped projections C PUV = 1 -> add ui+vi vectors to potential plots on mapped projections C (defaults for TUV and PUV = 1) C UVMAX = magnitude of vector scale arrow (m/s) (for arrow vector plots) C LOG10 = 1 -> plot log10 of fields, where appropriate (e.g., densities), C otherwise no (does not apply to vertical slices) (default LOG10 = 0) C ZP = selected pressure levels for mapped projections (-7. -> +5.) C (may be commented out for no plots at selected zp) C (nearest grid point is used, i.e., no interpolation in pressure) C HTS = selected heights for interpolation on mapped projections C (may be commented out for no plots at height surfaces) C C ----------------- VERTICAL SLICE CONTOURS ---------------- C C (note all vertical slice input values are floating point) C LON = longitude(s) for lon slices (use 999. to get zonal means) C (zp and/or ht on y-axis, latitude on x-axis) (may be commented out) C SLT = solar local time(s) for lon slices C (zp and/or ht on y-axis, latitude on x-axis) (may be commented out) C LAT = latitude(s) for lat slices C (zp and/or ht on y-axis, longitude on x-axis)(may be commented out) C LOC = lat,lon pairs of locations at which line plots are made with C zp and/or height on y-axis, depending on ZP_RANGE and HT_SCALE C (can use LOC=999.,999. to get global means) (may be commented out) C ZP_RANGE = bot,top of zp pressure range for y-axis of vertical slices C (within -7. to +5.) (may be commented out) C HT_SCALE = bot,top,del = bottom, top, and delta of height scale (km) for C y-axis of vertical slices (e.g., 100.,500.,25.) (may be commented out) C C Note: log10 of field is always plotted on vertical slices where appropriate C (e.g., densities, not temp or winds) C C --------------- FIELDS, FILES, AND CONTOUR VALUES ---------------- C C FIELDS = fields to plot (integers 0 or 1) C (this is required, and must provide a value for each field) C (rho is total density o2 + o1 + n2) (unvn and uivi are vector C plots or neutral and ion velocities without contours) 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 WRGRID = 1 or 2 -> make ascii data files of fields at entire 3-d 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 tigcmproc.fr001.dat, tigcmproc.fr002.dat, ...) C WRASCII = 2 -> make single ascii file containing data from each frame C (file name will be tigcmproc.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/tigcmproc ' C MSSASCII = 56-char mss path to which global ascii files (WRGRID) C may be disposed. C DIRCGM = 56-char remote directory to which metacode file (gmeta) may be C remote copied. C CMIN = contour minimums (if > CMAX, will use min of plotted array) C CMAX = contour maximums (if < CMIN, will use max of plotted array) C CINT = contour intervals (if = 0., will choose interval in software) C (CMIN,CMAX,CINT are floating point values) C (defaults for CMIN,CMAX,CINT = 1.,0.,0. for all fields) C C Note CMIN must be < CMAX, and CINT > 0. for CMIN,CMAX to take effect C however, CINT > 0. will take effect regardless of CMIN,CMAX C CMIN,CMAX,CINT apply to all plots. Remember that these will change C radically on mapped projections depending on LOG10. Note also that a C CINT that is appropriate for mapped projections may not be appropriate C for vertical slices. C C Please direct comments, bugs, and suggestions for additional options to: C Ben Foster (hao) (foster@ncar.ucar.edu) 303-497-1595 C Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc C CHISTVOL = '/ROBLE/RGR93/TESP03 ' CDDHHMM = 2,0,0 CHISTVOL = '/ROBLE/RGR94/TESP01 ','/ROBLE/RGR94/TESP02 ', C '/ROBLE/RGR94/TESP03 ' CDDHHMM = C 1, 0, 0, 1, 1, 0, 1, 2, 0, 1, 3, 0, 1, 4, 0, 1, 5, 0, C 1, 6, 0, 1, 7, 0, 1, 8, 0, 1, 9, 0, 1,10, 0, 1,11, 0, C 1,12, 0, 1,13, 0, 1,14, 0, 1,15, 0, 1,16, 0, 1,17, 0, C 1,18, 0, 1,19, 0, 1,20, 0, 1,21, 0, 1,22, 0, 1,23, 0, 2,0,0 C CHISTVOL = '/EMERY/HIST95/J93NSC ' CDDHHMM = 25,13,15 CHISTVOL = '/EMERY/HIST93/M92BG2 ' CDDHHMM = 27,1,0, 27,2,0 HISTVOL = '/ROBLE/RGR94/TETWJ3 ','/ROBLE/RGR94/TETWJ4 ' DDHHMM = 3,0,0, 4,0,0 C C IDEN=0 -> mass mixing ratios; IDEN=1 -> cm3; IDEN=2 -> gm/cm3 C IDEN = 1 COLOR = 0 C C Mapped projections: C (use CENxxx = lat,999. for local noon longitude) C CCENGLB = 0.,0. CCENGLB = 0.,999. CCENMOLW = 60.,999., -60.,999. POLLAT = -42.5, 42.5 CCENSATV = 60.,999.,-60.,999. ICONT = 1 TUV = 1 PUV = 1 UVMAX = 0. CUVMAX = 0. LOG10 = 0 ZP = -4.,2. CHTS = 300. C C Vertical Slices and location line plots: C (use LON = 999. for zonal means) C (use LOC = 999.,999. for global means) C LON = 0.,-70. CSLT = 12. CLAT = 57.5,0. LOC = 40.1,-105.5, -60.,999., 999.,999. ZP_RANGE = -7.,5. CHT_SCALE = 100.,500.,10. C C tn un vn o2 o1 n4s no o+ n2d ti te ne o2+ w z FIELDS = 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C pot n2 ui vi wi fof2 hmf2 rho unvn uivi o/o2 o/n2 o2/n2 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, C o/(o2+n2) 0 IONVEL = 2 WRGRID = 0 WRASCII = 0 DIRASCII = 'vishnu.hao.ucar.edu:/d/foster/tigcmproc ' C C MSSASCII currently for WRGRID only (7/8/93) CMSSASCII = '/FOSTER/tigcmproc ' DIRCGM = 'vishnu.hao.ucar.edu:/d/foster/tigcmproc/tigcmproc.cgm ' C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o1 n4s no o+ CCMIN = 1., 1., 1., 1., 1., 1., 1., 1., C n2d ti te ne o2+ w z pot C 1., 1., 1., 1., 1., 1., 1., -2000., C n2 ui vi wi fof2 hmf2 rho unvn C 1., -50., -50., -50., 1., 1., 1., 1., C uivi o/o2 o/n2 o2/n2 o/(o2+n2) C 1., 1., 1., 1., 1. C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o1 n4s no o+ CCMAX = 0., 0., 0., 0., 0., 0., 0., 0., C n2d ti te ne o2+ w z pot C 0., 0., 0., 0., 0., 0., 0., 2000., C n2 ui vi wi fof2 hmf2 rho unvn C 0., 50., 50., 50., 0., 0., 0., 0., C uivi o/o2 o/n2 o2/n2 o/(o2+n2) C 0., 0., 0., 0., 0. C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o1 n4s no o+ CCINT = 25., 0., 0., 0., 0., 0., 0., 0., C n2d ti te ne o2+ w z pot C 0., 0., 0., .2, 0., 0., 0., 400., C n2 ui vi wi fof2 hmf2 rho unvn C 0., 10., 10., 10., 0., 0., 0., 0., C uivi o/o2 o/n2 o2/n2 o/(o2+n2) C 0., 0., 0., 0., 0. ENDOFREAD