C Ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc C C INPUT FILE FOR TIGCMLOC CRAY PROCESSOR (zp = -7 to +5, kmx=25): C (binary executable is on mss /FOSTER/bin/tigcmloc) C (source is on sting:~foster/tigcmloc and mss /FOSTER/src/tigcmloc.tar) C C This is a multi-ut processor, i.e., most plots are made with ut on x-axis 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 Do not use tabs as 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 ','/ROBLE/RGR92/LTCS02 ' 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 C COLOR = 1 -> make color fill contours, monochrome otherwise C (default COLOR = 0) C PLUTVERT = 1 -> contour with ut on x-axis, zp/ht on y-axis at LOC C locations. ZP_RANGE and/or HT_SCALE are used to determine y-axis C (set PLUTVERT=0 for no such contours) C LOC = lat,lon pairs of selected locations. May use lat,999. to get C zonal means at selected lat, or 999.,999. to get global means C BILIN = 1 -> do bilinear interpolation to desired location; if BILIN=0, C use nearest grid point. (No interp is done for global means. C If doing zonal means and BILIN=1, interp is done to desired lat) C LOGLOC = 1 -> plot log10 of density fields at locations, otherwise no C ZP_RANGE = bottom,top range of pressures for y-axis (-7 -> +5) C (e.g.: ZP_RANGE= -4.,2.) C (used for PLUTVERT contours and UTLINE line plots) C HT_SCALE = bottom,top,delta heights for y-axis, e.g., 100.,500.,10. C (used for PLUTVERT contours and UTLINE line plots) C ZPHTLINE = selected zp(s) and/or ht(s) at which to make line plots at C LOC locations with ut on x-axis, field on y-axis C (is assumed a height if > 5., otherwise is assumed zp) C UTLINE = selected times (integer triplets) at which to make line plots C at LOC locations with field on x-axis, zp/ht on y-axis C (must be subset of DDHHMM) (ZP_RANGE and/or HT_SCALE will be C used for y-axis) C ZPHT_LON = pairs of zp/ht,lon for contours with ut on x-axis, lat on y-axis C (first of each pair may be pressure or height, 2nd of each pair may be C a longitude, or 999. for zonal means) C DOPPLER = 1 -> will calculate doppler t,u,v, and brightness, and make C line plots at LOC locations with ut on x=axis (if ZPHTLINE > 0). Also C calculate volume emission rate, and contour with ut on x-axis, zp/ht C on y-axis if (PLUTVERT > 0) (see doptuv.f and solred.f in libutil.a) C C --------------- FIELDS 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) C (fof2 and hmf2 will not be plotted where zp/ht is on y-axis) C IONVEL = flag for calculation of ion drift vel (for ui,vi,wi only) C = 1 -> calculate ExB drifts C = 2 -> caluclate ExB + effects of neutral atmos 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/tigcmloc ' C DIRCGM = 56-char remote directory to which metacode file (gmeta) may be C remote copied. C 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 HISTVOL = '/ROBLE/RGR94/TESMX1 ','/ROBLE/RGR94/TESMX2 ', '/ROBLE/RGR94/TESMX3 ','/ROBLE/RGR94/TESMX4 ', '/ROBLE/RGR94/TESMX5 ','/ROBLE/RGR94/TESMX6 ', '/ROBLE/RGR94/TESMX7 ','/ROBLE/RGR94/TESMX8 ', '/ROBLE/RGR94/TESMX9 ','/ROBLE/RGR94/TESMX10 ' DDHHMM = 0, 0,0, 0, 1,0, 0, 2,0, 0, 3,0, 0, 4,0, 0, 5,0, 0, 6,0, 0, 7,0, 0, 8,0, 0, 9,0, 0,10,0, 0,11,0, 0,12,0, 0,13,0, 0,14,0, 0,15,0, 0,16,0, 0,17,0, 0,18,0, 0,19,0, 0,20,0, 0,21,0, 0,22,0, 0,23,0, 1, 0,0, 1, 1,0, 1, 2,0, 1, 3,0, 1, 4,0, 1, 5,0, 1, 6,0, 1, 7,0, 1, 8,0, 1, 9,0, 1,10,0, 1,11,0, 1,12,0, 1,13,0, 1,14,0, 1,15,0, 1,16,0, 1,17,0, 1,18,0, 1,19,0, 1,20,0, 1,21,0, 1,22,0, 1,23,0, 2, 0,0, 2, 1,0, 2, 2,0, 2, 3,0, 2, 4,0, 2, 5,0, 2, 6,0, 2, 7,0, 2, 8,0, 2, 9,0, 2,10,0, 2,11,0, 2,12,0, 2,13,0, 2,14,0, 2,15,0, 2,16,0, 2,17,0, 2,18,0, 2,19,0, 2,20,0, 2,21,0, 2,22,0, 2,23,0, 3, 0,0, 3, 1,0, 3, 2,0, 3, 3,0, 3, 4,0, 3, 5,0, 3, 6,0, 3, 7,0, 3, 8,0, 3, 9,0, 3,10,0, 3,11,0, 3,12,0, 3,13,0, 3,14,0, 3,15,0, 3,16,0, 3,17,0, 3,18,0, 3,19,0, 3,20,0, 3,21,0, 3,22,0, 3,23,0, 4, 0,0, 4, 1,0, 4, 2,0, 4, 3,0, 4, 4,0, 4, 5,0, 4, 6,0, 4, 7,0, 4, 8,0, 4, 9,0, 4,10,0, 4,11,0, 4,12,0, 4,13,0, 4,14,0, 4,15,0, 4,16,0, 4,17,0, 4,18,0, 4,19,0, 4,20,0, 4,21,0, 4,22,0, 4,23,0, 5,0,0 C C IDEN=0 -> mass mixing ratios; IDEN=1 -> cm3; IDEN=2 -> gm/cm3 C IDEN = 1 COLOR = 0 C C PLUTVERT = 1 -> : contour ut on x-axis, zp and/or ht on y-axis: PLUTVERT = 0 LOC = 999.,999. CLOCNAME = ' ' BILIN = 1 LOGLOC = 1 ZP_RANGE = -7.,5. HT_SCALE = 100.,600.,10. CZPHTLINE = -6.,-4.,-2.,0. CUTLINE = 3,0,0, 3,6,0, 3,12,0, 3,18,0, 4,0,0 C C ZPHT_LON = pairs of zp/ht,lon for contours with ut on x-axis, lat on y-axis C (first of each pair may be pressure or height, 2nd of each pair may be C a longitude or 999. for zonal means) C ZPHT_LON = -4.,999., 2.,999. CZPHT_LON = 2.,-55. LOGUTLAT = 1 DOPPLER = 0 C C tn un vn o2 o1 n4s no o+ n2d ti te ne o2+ w z FIELDS = 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, C pot n2 ui vi wi fof2 hmf2 rho o/(o2+n2) o/n2 0, 1, 0, 0, 0, 0, 0, 1, 0, 0 IONVEL = 2 WRASCII = 0 DIRASCII = 'vishnu.hao.ucar.edu:/d/foster/tigcmloc ' CDIRCGM = 'vishnu.hao.ucar.edu:/d/foster/tigcmloc/tigcmloc.cgm ' DIRCGM = 'vishnu.hao.ucar.edu:/d/foster/tigcmloc/TESMX.utlat.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., 1., C n2 ui vi wi fof2 hmf2 rho o/(o2+n2) C 1., 1., 1., 1., 1., 1., 1., 1., C o/n2 C 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., 0., C n2 ui vi wi fof2 hmf2 rho o/(o2+n2) C 0., 0., 0., 0., 0., 0., 0., 0., C o/n2 C 0. C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o1 n4s no o+ CCINT = 0., 0., 0., 0., 0., 0., 0., 0., C n2d ti te ne o2+ w z pot C 0., 0., 0., .2, 0., 0., 0., 0., C n2 ui vi wi fof2 hmf2 rho o/(o2+n2) C 0., 0., 0., 0., 0., 0., 0., 0., C o/n2 C 0. ENDOFREAD