Stan: Below is the new tigcmproc script with wrascii, as we discussed. You will probably want to set WRASCII=2, comment out DIRASCII and DIRCGM, then put ftp's in your irje script to transfer the tigcmproc.dat and gmeta files to your home machine. (You could, of course, submit via MIGS from sting and let it do the rcp's, but you probably prefer to use irje from your own machine). (You might also consider changing your home login to stans because then you could use rcp/rsh, etc.) Note that memory required for this job varies according to the number of fields you are plotting. Keep an eye on HWM mem, and adjust QSUB -lM accordingly. As I said, the wrascii bit is brand new, so please let me know if you have problems or suggestions... --Ben # # QSUB -s /bin/csh # QSUB -eo # QSUB -lT 600 # QSUB -q prem # QSUB -nr # QSUB -lM 3MW ja cd $TMPDIR set job = tigcmproc batchname $job.$$ msread $job /FOSTER/bin/$job.dev chmod u+x $job pshell $job << '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 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 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 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 DIRCGM = 56-char remote directory to which metacode file (gmeta) may be C remote copied. C (if DIRASCII or DIRCGM are commented out, no rcp will be done) 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/RGR92/DYEM03 ' HISTVOL = '/ROBLE/RGR91/AFSY05 ' CHISTVOL = '/ROBLE/RGR92/DYEM03 ','/ROBLE/RGR91/AFSY05 ' CDDHHMM = 3,0,0, 3,1,0 DDHHMM = 3,4,0, 3,5,0 CDDHHMM = 3,0,0, 3,4,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. CCENMOLW = 60.,999., -60.,999. POLLAT = 30., -30. CENSATV = 60.,999. ICONT = 1 TUV = 1 PUV = 1 LOG10 = 0 ZP = -4. HTS = 200. 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.,999. CSLT = 0.,12. CLAT = 62.5 LOC = 62.5,-70., 999.,999. ZP_RANGE = -5.,2. 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, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C pot n2 ui vi wi fof2 hmf2 rho unvn uivi 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 WRASCII = 2 DIRASCII = 'sting.hao.ucar.edu:/d/foster/tigcmproc ' DIRCGM = 'sting.hao.ucar.edu:/d/foster/tigcmproc ' 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 unvn C 1., 1., 1., 1., 1., 1., 1., 1., C uivi 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 unvn C 0., 0., 0., 0., 0., 0., 0., 0., C uivi 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 unvn C 0., 0., 0., 0., 0., 0., 0., 0., C uivi C 0. ENDOFREAD 'EOF' ja -st