# # QSUB -s /bin/csh # QSUB -eo # QSUB -lT 300 # QSUB -q reg # QSUB -nr # QSUB -lM 3MW ja cd $TMPDIR set job = tigcmproc batchname $job.$$ if (! -e $job.ymp) then echo 'Doing msread of executable...' msread $job.ymp /FOSTER/bin/$job endif chmod u+x $job.ymp $job.ymp << 'EOF' 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 HISTVOL = '/ROBLE/RGR95/TEESN01 ' DDHHMM = 5,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 CENGLB = 0.,0. CCENGLB = 0.,999. CCENMOLW = 60.,999., -60.,999. CPOLLAT = -37.5 CCENSATV = 60.,999.,-60.,999. ICONT = 1 TUV = 1 PUV = 1 UVMAX = 0. CUVMAX = 0. LOG10 = 0 ZP = 1.,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 CLON = -180.,-75.,0. CSLT = 12. CLAT = 57.5,0. CLOC = 40.1,-105.5, -60.,999., 999.,999. CZP_RANGE = -7.,5. HT_SCALE = 100.,600.,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, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, C o/(o2+n2) 0 IONVEL = 1 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+ CMIN = 480., 1., 1., 1., 1., 1., 1., 1., C n2d ti te ne o2+ w z pot 1., 1., 1., 1., 1., 1., 1., -18000., C n2 ui vi wi fof2 hmf2 rho unvn 1., -50., -50., -65., 1., 100., 1., 1., C uivi o/o2 o/n2 o2/n2 o/(o2+n2) 1., 1., 1., 1., 1. C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o1 n4s no o+ CMAX = 1640., 0., 0., 0., 0., 0., 0., 0., C n2d ti te ne o2+ w z pot 0., 0., 0., 0., 0., 0., 0., 14000., C n2 ui vi wi fof2 hmf2 rho unvn 0., 50., 50., 35., 20., 675., 0., 0., C uivi o/o2 o/n2 o2/n2 o/(o2+n2) 0., 0., 0., 0., 0. C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o1 n4s no o+ CINT = 40., 0., 0., 0., 0., 0., 0., 0., C n2d ti te ne o2+ w z pot 0., 0., 0., .2, 0., 0., 0., 2000., C n2 ui vi wi fof2 hmf2 rho unvn 0., 10., 10., 5., 1., 25., 0., 0., C uivi o/o2 o/n2 o2/n2 o/(o2+n2) 0., 0., 0., 0., 0. ENDOFREAD 'EOF' ja -st