Dear ti(e)gcm data user: Below is a script for an updated ti(e)gcm "snapshot" cray processor. This program makes plots at single ut's from ti(e)gcm histories existing on the ncar mass store system. It may be submitted to shavano from hao machines via MIGS with: "nrnet shjob scriptname", or "shjob scriptname", if shjob is aliased as "alias shjob nrnet shjob \!*". Those of you not on MIGS, or submitting from non-hao machines must change the script accordingly, e.g., to use irje. (The script may also be submitted interactively from shavano). Contact me if you need help. Before submitting the job: 1. Change any QSUB values you want (e.g., -q and -lT) (-lM can be 3MW) 2. Fill in the metapath script variable with the machine and directory where you want the metacode returned, for example: set metapath = sting:/d/foster/tigcmproc or set metapath = ouray.hao.ucar.edu:/d/roble/tigcmproc (The name of the file sent to this directory will be tigcmproc.cgm) Note that with certain input values, several hundred plots may be generated, producing large metacode files (e.g., 348 monochrome plots generated a 13.8Mb file). Estimate the number of plots you should get before running the job, and be sure the metapath directory has enough space. 3. READ ALL COMMENTS FOR THE LEXICAL READS, and change the input keyword = value pairs for the options you want. Output listings will be returned to your default MIGS directory (usually ~/ntwk) with file name tigcmproc.xxxxx, where xxxxx is the PID from the cray job. Please report bugs and comments to me. P.S.: if you do not wish to receive notices of tgcm model processors, such as this one, please let me know and I will remove you from the mailing list. --Ben Foster (ncar/hao) foster@ncar.ucar.edu (303)-497-1595 # # QSUB -s /bin/csh # QSUB -eo # QSUB -lT 600 # QSUB -q prem # QSUB -nr # QSUB -lM 3MW ja cd $TMPDIR set job = tigcmproc set metapath = msread $job /FOSTER/bin/$job 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 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) 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 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 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 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/RGR92/DYEM03 ' CHISTVOL = '/ROBLE/RGR91/AFSY05 ' CHISTVOL = '/ROBLE/RGR92/DYEM03 ','/ROBLE/RGR91/AFSY05 ' DDHHMM = 3,0,0 CDDHHMM = 3,4,0 CDDHHMM = 3,0,0, 3,4,0 C C IDEN=0 -> mass mixing ratios; IDEN=1 -> cm3; IDEN=2 -> cm3 mixing ratios C IDEN = 1 COLOR = 1 C C Mapped projections: C (use CENxxx = lat,999. for local noon longitude) C CENGLB = 0.,0. CCENMOLW = 60.,999., -60.,999. CPOLLAT = 30., -30. CCENSATV = 60.,999., -60.,999. ICONT = 1 TUV = 1 PUV = 1 LOG10 = 1 ZP = -4.,2. CHTS = 200.,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 = 999. CSLT = 0.,12. CLAT = 60. LOC = 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, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C pot n2 ui vi wi fof2 hmf2 rho 1, 0, 0, 0, 0, 0, 0, 0 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 C 1., 1., 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., 0., C n2 ui vi wi fof2 hmf2 rho C 0., 0., 0., 0., 0., 0., 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 C 0., 0., 0., 0., 0., 0., 0. ENDOFREAD 'EOF' batchname $job.$$ rcp gmeta $metapath ja -st