C C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< C C INPUT DECK FOR TI(E)GCM DIFFERENCE FIELDS PROCESSOR (tigcmdif): 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. E.g., reading hourly histories C over a 6-day period for both perturbed and control cases consumed C over 1.5 hours of locked i/o wait time (0.5 gau in prem queue)! C These jobs should be run in the econ queue. C C PERTVOL = Perturbed volume(s) C PERTTIME = model day, hour, minute (e.g. '01:12:00') of perturbed history C CNTRVOL = Control volume(s) 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 mass densities (gm/cm3) 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 UVMAX = magnitude of maximum vector (m/s) 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 (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 ZPHT_LON = pairs of (zp or ht),lon for use with UTLAT. If first number C of a pair is <= 5., it is assumed to be a pressure, otherwise is C taken as a height (km). Second number of a pair may be 999. for C zonal means. C Note: it is recommended that if UTLOC, UTLINE, or UTLAT are turned on, C that C other plots (those plotted at each ut) be turned off. C This is for 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 LOCNAME = 16-char location name. Optional, but if given, there must be C one name per lat,lon pair in LOC (will appear on plot labels) C (blank or dummy names may be used for zonal or global mean locs C BILIN = 1 -> do bilinear interpolation to LOC (also to lon in lon slices, C lat in lat slices), otherwise use nearest grid points C ZP_RANGE = minimum,maximum range of zp (y-axis) for LON, SLT, LAT, C LINELOC, and UTLOC plots 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 (vpar = neutral wind parallel to magnetic field) C (wpar = vertical component of vpar) 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 tigcmdif.fr001.dat, tigcmdif.fr002.dat, ...) C 2 -> make single ascii file containing data from each frame C (file name will be tigcmdif.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/tigcmdif ' 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/tigcmdif and IWRASCII=2, then new mss C file will be disposed as /FOSTER/tigcmdif/tigcmdif.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 PERTVOL = '/EMERY/TI95/M92A28 ' PERTTIME = 28,00,00, 28,01,00, 28,02,00, 28,03,00, 28,04,00, 28,05,00, 28,06,00 C CNTRVOL = '/EMERY/TI95/M92C28 ' CNTRTIME = 28,00,00, 28,01,00, 28,02,00, 28,03,00, 28,04,00, 28,05,00, 28,06,00 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, 27.5 CCENSATV = 60.,999., -60.,999. ICONT = 1 TUV = 1 PUV = 1 UVMAX = 0. LOG10 = 0 CZP = -4.,2. CHTS = 120.,180.,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 = 0.,-70. CSLT = 12. CLAT = 0.,60. C C 18.3 -66.8 -8.0 ARECIBO C 67.0 -51.0 -39.0 SONDRESTROM C 69.6 19.2 1.0 EISCAT C 40.1 -105.5 13.0 FRITZ PEAK C 42.6 -71.5 -14.0 MILLSTONE C CLOC = 18.3,-66.8, 67.,-51., 69.6,19.2, 40.1,-105.5, 42.6,-71.5 C CLOCNAME = 'ARECIBO ','SONDRESTROM ','EISCAT ', C 'FRITZ PEAK ','MILLSTONE ' C LOC = 74.7, -94.9, 64.6, 18.8, 45.4, -75.9, 54.7, 20.6, 56.7, 61.1 LOCNAME = 'Resolute Bay ','Lycksele ','Ottawa ', 'Kaliningrad ','Sverdlovsk ' C BILIN = 1 C LINELOC = vertical profiles (line plots) at LOC locations and each ut: LINELOC = 0 C C UTLOC = contour ut on x-axis, zp or ht on y-axis at LOC locations UTLOC = 1 C C UTLINE = line plots with ut on x-axis, field on y-axis at ZP or HTS C at LOC locations (if HTS is set, HT_SCALE must also be set, C and all HTS must within range of HT_SCALE) 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 = 0 ZPHT_LON = -4.,-70., 2.,-70., -4.,999., 2.,999. CZPHT_LON = 120.,-70., 120.,999., 180.,-70., 180.,999., C 300.,-70., 300.,999. C C Pressure and/or height for y-axes C 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 = 0, 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 vpar wpar o/n2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, C n2/o 1 WRASCII = 0 DIRASCII = 'vishnu.hao.ucar.edu:/d/foster/tigcmdif ' CMSSASCII = '/FOSTER/tigcmdif ' DIRCGM = 'vishnu.hao.ucar.edu:/d/foster/tigcmdif/tigcmdif.cgm ' CDIRCGM = 'vishnu.hao.ucar.edu:/d/foster/tigcmdif/tesmxmb.utlat.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 C 1., 1., 1., 1., 1., 1., 1., C unvn uivi C 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 unvn uivi C 0., 0. C C23456789012345678901234567890123456789012345678901234567890123456789012 C tn un vn o2 o1 n4s no o+ CCINT = 10., 10., 10., 0., 0., 0., 0., 0., C n2d ti te ne o2+ w z pot C 0., 0., 0., .2, 0., 0., 0., 500., C n2 ui vi wi fof2 hmf2 rho C 0., 0., 0., 1., 0., 0., 0., C unvn uivi vpar wpar C 0., 0., .6, .6 ENDOFREAD