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 CHISTVOL = '/ROBLE/RGR94/TISP01 ','/ROBLE/RGR94/TISP02 ', C '/ROBLE/RGR94/TISP03 ' CHISTVOL = '/ROBLE/RGR94/TESP01 ','/ROBLE/RGR94/TESP02 ', C '/ROBLE/RGR94/TESP03 ' CHISTVOL = '/ROBLE/RGR94/TESP21 ','/ROBLE/RGR94/TESP22 ', C '/ROBLE/RGR94/TESP23 ' HISTVOL = '/ROBLE/RGR93/TESP01 ','/ROBLE/RGR93/TESP02 ', '/ROBLE/RGR93/TESP03 ' CDDHHMM = C 1, 0, 0, 1, 1, 0, 1, 2, 0, 1, 3, 0, 1, 4, 0, 1, 5, 0, 1,6,0 DDHHMM = 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 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 or ht on y-axis: PLUTVERT = 1 LOC = -87.5,0., -67.65,62.9, -75.5,-26.6, -43.98,170.42, 40.1,-105.5 LOCNAME = 'SOUTH POLE ','MAWSON ','HALLEY ', 'MT. JOHN ','FRITZ PEAK ' BILIN = 1 LOGLOC = 1 ZP_RANGE = -7.,5. CHT_SCALE = 100.,170.,1. 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 CZPHT_LON = 2.,125. CZPHT_LON = 2.,-55. LOGUTLAT = 1 DOPPLER = 1 C C tn un vn o2 o1 n4s no o+ n2d ti te ne o2+ w z FIELDS = 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, C pot n2 ui vi wi fof2 hmf2 rho o/(o2+n2) o/n2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 IONVEL = 2 WRASCII = 0 DIRASCII = 'vishnu.hao.ucar.edu:/d/foster/tigcmloc ' DIRCGM = 'vishnu.hao.ucar.edu:/d/foster/tigcmloc/tigcmloc.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