C C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< C C INPUT DECK FOR TIMESGCM DIFFERENCE FIELDS PROCESSOR (timesdif): C C 11/19/93: Modified to use a cubic spline interpolation for vertical C slices (see, ~foster/timesdif/htint.f) C C Notes: C 1) Memory requirements of this program will increase with the C number of values given for PERTTIME, CNTRTIME, LOC, FIELDS, C UTLOC and UTLINE. If you find it is needing more than 6 MW, C it is recommended that you cut back on some options and C make multiple runs to accomplish the needed work. E.g., C UTLOC and UTLINE should be used without the single-ut plot options. C 2) I/O wait time can be significant when making UTLOC and UTLINE C plots for long model runs. C C PERTVOL = Perturbed volume(s) (24-char string) C PERTTIME = model day, hour, minute (e.g. '01:12:00') of perturbed history C CNTRVOL = Control volume(s) (24-char string) C CNTRTIME = model day, hour, minute (e.g. '08:12:00') of control history C (PERTTIME and CNTRTIME may have different day numbers, but must C be same UT C IDEN = density flag: C 0 -> mass mixing ratios C 1 -> convert species to number densities (cm3) C 2 -> convert species to number density mixing ratios C IDIF = difference flag: C 0 -> make raw diffs of all fields (simple subtraction) C 1 -> make percent diffs of species, raw diffs of other fields C COLOR = 1 -> color fill (otherwise monochrome) C C--------- mapped projections at selected zp or ht and current ut --------- C C CENGLB = glat,glon pairs for center of cylindrical equidistant map C projection (rectangular with lon on x-axis, lat on y-axis) C (typically CENGLB = 0.,0.) C CENMOLW = lat,lon pairs for center of mollweide map projection C POLLAT = perimeter latitude(s) for stereographic polar contours C (control hem with sign) C CENSATV = lat,lon pairs for center of satellite view projection C ICONT = 0/1 -> do not / do plot continental outlines 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 mapped fields, where appropriate C (e.g., densities) (applies to CENGLB,CENMOLW,POLLAT,CENSATV) C ZP = selected pressure levels for mapped projections C HTS = selected height surfaces for interpolation on mapped projections C C------------------- vertical slices at current ut ----------------------- C C LON = selected longitudes for contours with lat on x-axis, zp or ht C on y-axis (use 999. to get zonal means) C SLT = selected solar local times for contours with lat on x-axis, zp or ht C on y-axis C LAT = selected latitudes for contours with lon on x-axis, zp or ht on y-axis C LINELOC = vertical profiles (line plots) at LOC locations and current ut C LOG10SLI = 1 -> plot log10 of differences in vertical slices C (note: control y-axis of vertical slices with ZP_RANGE and/or HT_SCALE) C C------------------ time-dependent (ut on x-axis) ----------------------- C C UTLOC = contour ut on x-axis, zp or ht on y-axis (depending on ZP_RANGE C and/or HT_SCALE) at LOC locations C UTLINE = line plots with ut on x-axis, field on y-axis at ZP pressures C and/or HTS heights, at LOC locations (if UTLINE=1 and HTS is set, then C HT_SCALE must also be set and all HTS heights must be within the height C scale defined by HT_SCALE) C UTLAT = contour ut on x-axis, latitude on y-axis at selected zp or ht C and longitudes specified in ZPHT_LON C Note: it is recommended that if UTLOC and/or UTLINE are turned on, that C other plots (those plotted at each ut) be turned off. This is for C conservation of memory C C------------------------------------------------------------------------ C C LOC = lat,lon pairs of locations for LINELOC, UTLOC, and/or UTLINE C (use lat,999. for zonal means or 999.,999. for global means) C ZP_RANGE = minimum,maximum range of zp (y-axis) for LON, SLT, LAT, C LINELOC, and UTLOC plots (full range = -17,5.) C HT_RANGE = linear height scale (bottom,top,delta in km) for LON, SLT, and C LAT plots C FIELDS = fields to plot (0 or 1 for each) C (unvn is neutral wind vectors (arrows) alone (mapped projections only)) C (uivi is ion drift vectors (arrows) alone (mapped projections only)) C WRASCII = C 0 -> do not make ascii data file(s) C 1 -> make an ascii data file for each frame C (file names will be timesdif.fr001.dat, timesdif.fr002.dat, ...) C 2 -> make single ascii file containing data from each frame C (file name will be timesdif.dat) C DIRASCII = 56-char unix directory to which ascii data file(s) may be C remote copied (rcp from shavano to DIRASCII will be used), e.g.: C DIRASCII = 'sting.hao.ucar.edu:/d/foster/timesdif ' C MSSASCII = 56-char mss dir to which ascii data file(s) are to be disposed. C (this should be a directory path, not a file name, however the dir C need not exist on mss prior to this job) C e.g., if MSSASCII=/FOSTER/timesdif and IWRASCII=2, then new mss C file will be disposed as /FOSTER/timesdif/timesdif.dat C DIRCGM = 56-char remote directory to which metacode file (gmeta) may be C remote copied (rcp of gmeta from shavano to DIRCGM). C (note DIRCGM may be dir/file name) 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 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 C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< C CPERTVOL = '/ROBLE/RGR94/TSEQDB4 ' CPERTTIME = 8,0,0 CPERTVOL = '/ROBLE/RGR94/TSSLDB4 ' CPERTTIME = 6,0,0 CPERTVOL = '/ROBLE/RGR94/TSESTB7 ' CPERTTIME = 30,0,0 CPERTVOL = '/ROBLE/RGR94/TSEDYM2 ' CPERTTIME = 2,0,0 CPERTVOL = '/ROBLE/RGR94/TSESTD6 ' CPERTTIME = 5,0,0 CPERTVOL = '/ROBLE/RGR94/TSEQBA4 ' CPERTTIME = 6,0,0 PERTVOL = '/ROBLE/RGR94/T15NEW9 ' PERTTIME = 13,0,0 C CCNTRVOL = '/ROBLE/RGR93/TSSOL12 ' CCNTRTIME = 22,0,0 CCNTRVOL = '/ROBLE/RGR93/TSEBZ10 ' CCNTRTIME = 18,0,0 CCNTRVOL = '/ROBLE/RGR94/TSESTB6 ' CCNTRTIME = 29,0,0 CCNTRVOL = '/ROBLE/RGR94/TSEDYL2 ' CCNTRTIME = 2,0,0 CCNTRVOL = '/ROBLE/RGR94/TSEDIV6 ' CCNTRTIME = 5,0,0 CCNTRVOL = '/ROBLE/RGR93/TSEQBZ4 ' CCNTRTIME = 6,0,0 CNTRVOL = '/ROBLE/RGR93/TSEQBZ4 ' CNTRTIME = 6,0,0 WRASCII = 0 DIRASCII = 'ouray.hao.ucar.edu:/d/times ' CMSSASCII = '/ROBLE/timesdif ' CDIRCGM = 'ouray.hao.ucar.edu:/d/times/solstice.co2dbl.dif.day6 ' CDIRCGM = 'ouray.hao.ucar.edu:/d/times/equinox.sol.diurn.dif.day6 ' CDIRCGM = 'ouray.hao.ucar.edu:/d/times/timesdif.new.old.histories ' DIRCGM = 'vishnu.hao.ucar.edu:/d/foster/timesdif/timesdif.cgm ' CDIRCGM = 'ouray.hao.ucar.edu:/d/times/equinox.smax-smin.day1 ' C C IDEN=0 -> mass mixing ratios; IDEN=1 -> cm3; IDEN=2 -> gm/cm3 C IDEN = 1 C C IDIF = 0 -> make raw diffs of all fields (simple subtraction) C IDIF = 1 -> make percent diffs of species, raw diffs of other fields C IDIF = 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. CPOLLAT = -27.5 CCENSATV = 60.,999., -60.,999. ICONT = 0 TUV = 1 PUV = 1 LOG10 = 0 ZP = -15.,-13.,-11.,-9.,-7.,-4.,-1.,+3. HTS = 80.,100.,120. 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 = 999.,0.,-180. LON = 999. CLON = -147.4 CSLT = 0. CLAT = 2.5,17.5,42.5,67.5 CLAT = 2.5,17.5 CLAT = 65.1 CLOG10SLI = 0 C Poker Flat for spire CLOC = 65.,-147.4 CLOC = 18.3,-66.8, 22.3,114.1, 3.9,-38.0 C C LINELOC = vertical profiles (line plots) at LOC locations and current ut: LINELOC = 0 C C UTLOC = contour ut on x-axis, zp on y-axis at LOC locations UTLOC = 0 C C UTLINE = line plots with ut on x-axis, field on y-axis at ZP (not HTS) C at LOC locations UTLINE = 0 C C UTLAT = 0/1 for latitudes at which to contour ut on x-axis, field on y-axis C ZPHT_LON = zp/ht,lon pairs for UTLAT C UTLAT = 1 ZPHT_LON = -4.,-70., -4.,999., 2.,-70., 2.,999. C C Pressure and/or height for y-axes C ZP_RANGE = -17.,+5. CHT_SCALE = 30.,300.,5. CHT_SCALE = 30.,120.,2. C C tn un vn o2 o n4s noz co co2 h2o h2 hox o+ ch4 o21d FIELDS = 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, C no2 no o3 o1 oh ho2 h n2d ti te ne o2+ w z epot n2 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, C ui vi wi fof2 hmf2 rho unvn uivi 0, 0, 0, 0, 0, 0, 0, 0 C tn un vn o2 o n4s noz co co2 h2o h2 hox o+ ch4 o21d CFIELDS = 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C no2 no o3 o1 oh ho2 h n2d ti te ne o2+ w z epot n2 C 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, C ui vi wi fof2 hmf2 rho unvn uivi C 0, 0, 0, 0, 0, 0, 0, 0 CWRASCII = 0 CDIRASCII = 'ouray.hao.ucar.edu:/d/times ' CMSSASCII = '/ROBLE/timesdif ' CDIRCGM = 'ouray.hao.ucar.edu:/d/times/equinox.co2dbl.dif.day18a ' C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o n4s noz co co2 h2o h2 hox o+ ch4 o21d CCMIN = 1.,1.,1.,1., 1.,1.,1., 1., 1.,1., 1., 1., 1.,1.,1., C no2 no o3 o1 oh ho2 h n2d ti te ne o2+ w z epot n2 C 1., 1.,1.,1., 1.,1., 1.,1.,1.,1.,1.,1., 1.,1.,1., 1., C ui vi wi fof2 hmf2 rho unvn uivi C 1., 1., 1., 1., 1., 1., 1., 1. C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o n4s noz co co2 h2o h2 hox o+ ch4 o21d CCMAX = 0.,0.,0.,0., 0.,0.,0., 0., 0.,0., 0., 0., 0.,0.,0., C no2 no o3 o1 oh ho2 h n2d ti te ne o2+ w z epot n2 C 0., 0.,0.,0., 0.,0., 0.,0.,0.,0.,0.,0., 0.,0.,0., 0., C ui vi wi fof2 hmf2 rho unvn uivi C 0., 0., 0., 0., 0., 0., 0., 0. C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o n4s noz co co2 h2o h2 hox o+ ch4 o21d CCINT = 20.,10.,10.,0., 0.,0.,0., 0., 0.,0., 0., 0., 0.,0.,0., C no2 no o3 o1 oh ho2 h n2d ti te ne o2+ w z epot n2 C 0., 0.,0.,0., 0.,0., 0.,0.,0.,0.,0.,0., 0.,0.,0., 0., C ui vi wi fof2 hmf2 rho unvn uivi C 0., 0., 0., 0., 0., 0., 0., 0. ENDOFREAD