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 C TIEGCM runs for spring 1994 AGU (hourly histories for model days 0->5): C C /ROBLE/RGR94/TESMB1-10: Equinox solar min (base case) C /ROBLE/RGR94/TESMX1-10: Equinox solar max (perturbed case) C C History Volume Mday:ut Jday Date C C /RGR94/TESMx1 00:00->00:12 94080 3/21/94 C /RGR94/TESMx2 00:13->01:01 94080-94081 3/21-22/94 C /RGR94/TESMx3 01:02->01:14 94081 3/22/94 C /RGR94/TESMx4 01:15->02:03 94081-94082 3/22-23/94 C /RGR94/TESMx5 02:04->02:16 94082 3/23/94 C /RGR94/TESMx6 02:17->03:05 94082-94083 3/23-24/94 C /RGR94/TESMx7 03:06->03:18 94083 3/24/94 C /RGR94/TESMx8 03:19->04:07 94083-94084 3/24-25/94 C /RGR94/TESMx9 04:08->04:20 94084 3/25/94 C /RGR94/TESMx10 04:21->05:00 94084-94085 3/25-26/94 C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< C PERTVOL = '/ROBLE/RGR94/TESMX1 ' CPERTVOL = '/ROBLE/RGR94/TESMX1 ','/ROBLE/RGR94/TESMX2 ', C '/ROBLE/RGR94/TESMX3 ','/ROBLE/RGR94/TESMX4 ', C '/ROBLE/RGR94/TESMX5 ','/ROBLE/RGR94/TESMX6 ', C '/ROBLE/RGR94/TESMX7 ','/ROBLE/RGR94/TESMX8 ', C '/ROBLE/RGR94/TESMX9 ','/ROBLE/RGR94/TESMX10 ' PERTTIME = C 0, 0,0, 0, 1,0 0, 0,0, 0, 1,0, 0, 2,0, 0, 3,0, 0, 4,0, 0, 5,0, 0, 6,0 CPERTTIME = C 0, 0,0, 0, 1,0, 0, 2,0, 0, 3,0, 0, 4,0, 0, 5,0, 0, 6,0, C 0, 7,0, 0, 8,0, 0, 9,0, 0,10,0, 0,11,0, 0,12,0, 0,13,0, C 0,14,0, 0,15,0, 0,16,0, 0,17,0, 0,18,0, 0,19,0, 0,20,0, C 0,21,0, 0,22,0, 0,23,0, C 1, 0,0, 1, 1,0, 1, 2,0, 1, 3,0, 1, 4,0, 1, 5,0, 1, 6,0, C 1, 7,0, 1, 8,0, 1, 9,0, 1,10,0, 1,11,0, 1,12,0, 1,13,0, C 1,14,0, 1,15,0, 1,16,0, 1,17,0, 1,18,0, 1,19,0, 1,20,0, C 1,21,0, 1,22,0, 1,23,0, C 2, 0,0, 2, 1,0, 2, 2,0, 2, 3,0, 2, 4,0, 2, 5,0, 2, 6,0, C 2, 7,0, 2, 8,0, 2, 9,0, 2,10,0, 2,11,0, 2,12,0, 2,13,0, C 2,14,0, 2,15,0, 2,16,0, 2,17,0, 2,18,0, 2,19,0, 2,20,0, C 2,21,0, 2,22,0, 2,23,0, C 3, 0,0, 3, 1,0, 3, 2,0, 3, 3,0, 3, 4,0, 3, 5,0, 3, 6,0, C 3, 7,0, 3, 8,0, 3, 9,0, 3,10,0, 3,11,0, 3,12,0, 3,13,0, C 3,14,0, 3,15,0, 3,16,0, 3,17,0, 3,18,0, 3,19,0, 3,20,0, C 3,21,0, 3,22,0, 3,23,0, C 4, 0,0, 4, 1,0, 4, 2,0, 4, 3,0, 4, 4,0, 4, 5,0, 4, 6,0, C 4, 7,0, 4, 8,0, 4, 9,0, 4,10,0, 4,11,0, 4,12,0, 4,13,0, C 4,14,0, 4,15,0, 4,16,0, 4,17,0, 4,18,0, 4,19,0, 4,20,0, C 4,21,0, 4,22,0, 4,23,0, 5,0,0 C CNTRVOL = '/ROBLE/RGR94/TESMB1 ' CCNTRVOL = '/ROBLE/RGR94/TESMB1 ','/ROBLE/RGR94/TESMB2 ', C '/ROBLE/RGR94/TESMB3 ','/ROBLE/RGR94/TESMB4 ', C '/ROBLE/RGR94/TESMB5 ','/ROBLE/RGR94/TESMB6 ', C '/ROBLE/RGR94/TESMB7 ','/ROBLE/RGR94/TESMB8 ', C '/ROBLE/RGR94/TESMB9 ','/ROBLE/RGR94/TESMB10 ' CNTRTIME = C 0, 0,0, 0, 1,0 0, 0,0, 0, 1,0, 0, 2,0, 0, 3,0, 0, 4,0, 0, 5,0, 0, 6,0 CCNTRTIME = C 0, 0,0, 0, 1,0, 0, 2,0, 0, 3,0, 0, 4,0, 0, 5,0, 0, 6,0, C 0, 7,0, 0, 8,0, 0, 9,0, 0,10,0, 0,11,0, 0,12,0, 0,13,0, C 0,14,0, 0,15,0, 0,16,0, 0,17,0, 0,18,0, 0,19,0, 0,20,0, C 0,21,0, 0,22,0, 0,23,0, C 1, 0,0, 1, 1,0, 1, 2,0, 1, 3,0, 1, 4,0, 1, 5,0, 1, 6,0, C 1, 7,0, 1, 8,0, 1, 9,0, 1,10,0, 1,11,0, 1,12,0, 1,13,0, C 1,14,0, 1,15,0, 1,16,0, 1,17,0, 1,18,0, 1,19,0, 1,20,0, C 1,21,0, 1,22,0, 1,23,0, C 2, 0,0, 2, 1,0, 2, 2,0, 2, 3,0, 2, 4,0, 2, 5,0, 2, 6,0, C 2, 7,0, 2, 8,0, 2, 9,0, 2,10,0, 2,11,0, 2,12,0, 2,13,0, C 2,14,0, 2,15,0, 2,16,0, 2,17,0, 2,18,0, 2,19,0, 2,20,0, C 2,21,0, 2,22,0, 2,23,0, C 3, 0,0, 3, 1,0, 3, 2,0, 3, 3,0, 3, 4,0, 3, 5,0, 3, 6,0, C 3, 7,0, 3, 8,0, 3, 9,0, 3,10,0, 3,11,0, 3,12,0, 3,13,0, C 3,14,0, 3,15,0, 3,16,0, 3,17,0, 3,18,0, 3,19,0, 3,20,0, C 3,21,0, 3,22,0, 3,23,0, C 4, 0,0, 4, 1,0, 4, 2,0, 4, 3,0, 4, 4,0, 4, 5,0, 4, 6,0, C 4, 7,0, 4, 8,0, 4, 9,0, 4,10,0, 4,11,0, 4,12,0, 4,13,0, C 4,14,0, 4,15,0, 4,16,0, 4,17,0, 4,18,0, 4,19,0, 4,20,0, C 4,21,0, 4,22,0, 4,23,0, 5,0,0 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 ZP = -4.,2. CHTS = 100.,200.,300.,400. 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.,-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 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 = 0 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 = 1 ZPHT_LON = -4.,-70., 2.,-70., -4.,999., 2.,999. C C Pressure and/or height for y-axes C CZP_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, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, C pot n2 ui vi wi fof2 hmf2 rho unvn uivi vpar wpar 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0 WRASCII = 0 DIRASCII = 'vishnu.hao.ucar.edu:/d/foster/tigcmdif ' CMSSASCII = '/FOSTER/tigcmdif ' CDIRCGM = 'vishnu.hao.ucar.edu:/d/foster/tigcmdif/tigcmdif.cgm ' DIRCGM = '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