From: Geoff Crowley <[email-masked]>
Date: Thu, 23 Feb 1995 18:44:21 -0700

95.02.23 18:13:37 IDENTIFIER 87279.shavano NAME HC2949 JOBID 16967
shavano: UNICOS 8.0.2.1, IOS 8.0.2
=======================================================================
THIS SYSTEM IS FOR THE USE OF AUTHORIZED USERS ONLY.  INDIVIDUALS USING
THIS COMPUTER SYSTEM WITHOUT AUTHORITY, OR IN EXCESS OF THEIR AUTHORITY,
ARE SUBJECT TO HAVING ALL THEIR ACTIVITIES ON THIS SYSTEM MONITORED AND
RECORDED BY SYSTEM PERSONNEL.  IN THE COURSE OF MONITORING INDIVIDUALS
IMPROPERLY USING THIS SYSTEM, OR IN THE COURSE OF SYSTEM MAINTENANCE,
THE ACTIVITIES OF AUTHORIZED USERS MAY ALSO BE MONITORED.
=======================================================================
System testing 2/28/95 05:30-09:00 MST
=======================================================================
news: CF9010.news CrayTools13.news nag 8.0.notes 7.0.5.notes SCC40.news
 CrayLibs10.news emulator10.news packing gpfa imsl lapack UNICOS70.news
 SCC30.news CF7750.news
Principal Project GAU use:
                 ---Project GAUs----  --Division GAUs--  ---Program GAUs--
YYMMDD  Project  Allocated      Used     Mon  3_Mon_avg     Mon  3_Mon_avg
======  ======== =========  ========  ======  =========  ======  =========
950222  37241000     99.00     55.33     N/A        N/A     N/A        N/A
Thu Feb 23 18:13:38 
Thu Feb 23 18:13:38 + exit 0
Thu Feb 23 18:13:38 + /usr/spool/nqs/scripts/++KIj+++++0+++
Batch output name set to 'timesdif.84327'
           C                                                                               
           C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<               
           C                                                                               
           C  INPUT DECK FOR TIMESGCM DIFFERENCE FIELDS PROCESSOR (timesdif):              
           C  (executable on ncar mss /FOSTER/bin/timesdif)                                
           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 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 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 = 'vishnu.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) ([email-masked]) 303-497-1595                          
           C                                                                               
           C<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<               
           C                                                                               
            PERTVOL = '/ROBLE/RGR94/TS10DS6   '                                            
            PERTTIME = 23,0,0                                                              
           C                                                                               
            CNTRVOL = '/ROBLE/RGR94/TS10D12   '                                            
            CNTRTIME = 23,0,0                                                              
           C                                                                               
           C IDEN=0 -> mass mixing ratios; IDEN=1 -> cm3; IDEN=2 -> gm/cm3                 
           C                                                                               
            IDEN = 0                                                                       
           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 = 0                                                                       
            COLOR = 0                                                                      
           C                                                                               
           C Mapped projections:                                                           
           C (use CENxxx = lat,999. for local noon longitude)                              
           C                                                                               
            CENGLB = 0.,0.                                                                 
           CCENMOLW = 60.,999., -60.,999.                                                  
           CPOLLAT = -27.5                                                                 
           CCENSATV = 60.,999., -60.,999.                                                  
            ICONT = 1                                                                      
            TUV = 0                                                                        
            PUV = 1                                                                        
            LOG10 = 0                                                                      
            ZP = -6.,-13.                                                                  
           CHTS = 60.,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                                                                               
            LON = 999.,0.,-180.                                                            
           CSLT = 0.                                                                       
           CLAT = 2.5,42.5,67.5                                                            
            LOG10SLI = 0                                                                   
           C Arecibo, Hong Kong, Fortaleza                                                 
           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 Pressure and/or height for y-axes                                             
           C                                                                               
            ZP_RANGE = -17.,5.                                                             
           CHT_SCALE = 50.,120.,5.                                                         
           C                                                                               
           C          tn un vn o2  o n4s noz co co2 h2o h2 hox o+ ch4 o21d                 
            FIELDS =   1, 0, 0, 1, 1,  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               
                       0,   0, 1, 1, 0,  0, 0, 0, 0,  0, 0, 0,  0, 0, 0,  1,               
           C          ui  vi  wi fof2 hmf2 rho unvn uivi                                   
                       0,  0,  0,  0,   0,  1,   0,   0                                    
            WRASCII = 0                                                                    
            DIRASCII = 'cedar.hao.ucar.edu:/d/crowley/timesdif                  '          
           CMSSASCII = '/FOSTER/timesdif                                        '          
            DIRCGM =   'crowley.hao.ucar.edu:/d/crowley/timesdif/timesdif.cgm   '          
           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 =    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.                                    
            ENDOFREAD                                                                      
getinp: zprange=  -17.00    5.00 izprange=  1 45
getvol: acquiring history volume /ROBLE/RGR94/TS10DS6 from mss
        to disk file RGR94.TS10DS6 ...
getvol: successful msread of mss file /ROBLE/RGR94/TS10DS6 to 
  unicos file RGR94.TS10DS6
assgn: successful assign of unit  8 to file RGR94.TS10DS6
 
Gettgcm: Searching for  23: 0: 0 Found  22:15: 0 vol /ROBLE/RGR94/TS10DS6 (it= 1)
Gettgcm: Searching for  23: 0: 0 Found  22:16: 0 vol /ROBLE/RGR94/TS10DS6 (it= 2)
Gettgcm: Searching for  23: 0: 0 Found  22:17: 0 vol /ROBLE/RGR94/TS10DS6 (it= 3)
Gettgcm: Searching for  23: 0: 0 Found  22:18: 0 vol /ROBLE/RGR94/TS10DS6 (it= 4)
Gettgcm: Searching for  23: 0: 0 Found  22:19: 0 vol /ROBLE/RGR94/TS10DS6 (it= 5)
Gettgcm: Searching for  23: 0: 0 Found  22:20: 0 vol /ROBLE/RGR94/TS10DS6 (it= 6)
Gettgcm: Searching for  23: 0: 0 Found  22:21: 0 vol /ROBLE/RGR94/TS10DS6 (it= 7)
Gettgcm: Searching for  23: 0: 0 Found  22:22: 0 vol /ROBLE/RGR94/TS10DS6 (it= 8)
Gettgcm: Searching for  23: 0: 0 Found  22:23: 0 vol /ROBLE/RGR94/TS10DS6 (it= 9)
Gettgcm: Found history  23: 0: 0 on vol /ROBLE/RGR94/TS10DS6 (iv= 1, it=10)
 
This file contains secondary histories from a timesgcm run (without dynamo)
(Secondary histories contain only selected fields. They may be read by
  post-model processors but may not be used to restart the model)
 
The secondary histories in this file contain 13 fields,
  each with index > 0 in the following list:
 
   tn   un   vn   o2   ox  n4s  noz   co  co2  h2o   h2  hox   o+  ch4 o21d
    1    2    3    4    0    0    0    0    0    0    0    0    0    0    0
  no2   no   o3   o1   oh  ho2    h  n2d   ti   te   ne  o2+    w    z epot
    0    5    6    7    8    0    0    0    9   10   11    0   12   13    0
 
>>> Warning: field OX    not available on this timesgcm secondary history
gettgcm: allocated space for fields array: nlon=73 nzp=45 nlat=36 nflds= 6
The following fields will be defined:
TN      O2      O3      O1      Z       N2      
 
------------------------- TGCM HEADER SUMMARY --------------------------
 
HISTORY VOLUME = /ROBLE/RGR94/TS10DS6
MODEL TIME =  23: 0: 0
 
 iter nday nhr nmin =  8280    23     0     0
 label(1-5) = TGCM17 NMC BNDRYZ START SEPT.20,1985    
 date =      1993       21
 output(1) = ROBLE   RGR94   TS10D11 
 output(2) = ROBLE   RGR94   TS10D12 
 start =    21    0    0
 stp =    23    0    0    0    0    0
 hist =     0   12    0    0    0    0
 sav =     0   12    0
 step =   240
 mag =    -74.500   127.000    79.000   -70.000
 difhor,iuivi =   1  1
 sdtide =   0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
            0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
 dtide =   0.1600E+04  0.1150E+02
 ipower,aurora,dispos =   0  1  1
 data = FOSTER  TIMEDATVSOLM01Z 
 source =                         
 sourct =    21    0    0
 rdate = 11/01/94
 naurp =    3
 hp =    21.94
 cp =    58.15
 byimf =     0.00
 colfac =     1.50
 f107d =   101.78
 f107a =   130.49
 
------------------------------------------------------------------------
 
link command=ln RGR94.TS10DS6 /usr/tmp/TIGCM
timesdif: found perturbed history  23: 0: 0 on vol /ROBLE/RGR94/TS10DS6 it=  1 istimes=F isdyn=T
getvol: acquiring history volume /ROBLE/RGR94/TS10D12 from mss
        to disk file RGR94.TS10D12 ...
getvol: successful msread of mss file /ROBLE/RGR94/TS10D12 to 
  unicos file RGR94.TS10D12
assgn: successful assign of unit  9 to file RGR94.TS10D12
 
Gettgcm: Found history  23: 0: 0 on vol /ROBLE/RGR94/TS10D12 (iv= 1, it= 1)
         (this history is from a timesgcm model run (without dynamo))
gettgcm: allocated space for fields array: nlon=73 nzp=45 nlat=36 nflds= 7
The following fields will be defined:
TN      O2      OX      O3      O1      Z       N2      
 
------------------------- TGCM HEADER SUMMARY --------------------------
 
HISTORY VOLUME = /ROBLE/RGR94/TS10D12
MODEL TIME =  23: 0: 0
 
 iter nday nhr nmin =  8280    23     0     0
 label(1-5) = TGCM17 NMC BNDRYZ START SEPT.20,1985    
 date =      1993       21
 output(1) = ROBLE   RGR94   TS10D11 
 output(2) = ROBLE   RGR94   TS10D12 
 start =    21    0    0
 stp =    23    0    0    0    0    0
 hist =     0   12    0    0    0    0
 sav =     0   12    0
 step =   240
 mag =    -74.500   127.000    79.000   -70.000
 difhor,iuivi =   1  1
 sdtide =   0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
            0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
 dtide =   0.1600E+04  0.1150E+02
 ipower,aurora,dispos =   0  1  1
 data = FOSTER  TIMEDATVSOLM01Z 
 source =                         
 sourct =    21    0    0
 rdate = 11/01/94
 naurp =    3
 hp =    21.94
 cp =    58.15
 byimf =     0.00
 colfac =     1.50
 f107d =   101.78
 f107a =   130.49
 
------------------------------------------------------------------------
 
link command=ln RGR94.TS10D12 /usr/tmp/TIGCM
timesdif: found control history  23: 0: 0 on vol /ROBLE/RGR94/TS10D12 it=  1 istimes=F isdyn=T
 
>>> Perturbed and Control histories must have the same fields
    Fields on perturbed histories = 
  1  0  0  1  0  0  0  0  0  0  0  0  0  0  0  0  0  1  1  0  0  0  0  0  0  0  0  0  1  0  1
    Fields on control histories = 
  1  0  0  1  1  0  0  0  0  0  0  0  0  0  0  0  0  1  1  0  0  0  0  0  0  0  0  0  1  0  1
 STOP fields
 STOP executed at line 124 in Fortran routine 'TIMESDIF'
 CP: 1.196s,  Wallclock: 392.363s
 HWM mem: 3358034, HWM stack: 27095, Stack overflows: 0

Job Accounting - Summary Report
===============================

Job Accounting File Name         : /usr/tmp/nqs.+++++0KIj/.jacct16967
Operating System                 : sn1036 shavano 8.0.2 cbr.15 CRAY Y-MP
User Name (ID)                   : crowley (6891)
Group Name (ID)                  : univ (1000)
Account Name (ID)                : 37241000 (14404)
    GAUs Allocated               :      99.0000
    GAUs Used, as of 02/22/95    :      55.3300
    Division Monthly Use         :   N/A
    Program Monthly Use          :   N/A
Job Name (ID)                    : HC2949 (16967)
Report Starts                    : 02/23/95 18:13:38
Report Ends                      : 02/23/95 18:20:19
Elapsed Time                     :          401      Seconds
User CPU Time                    :            1.2313 Seconds
System CPU Time                  :            0.5578 Seconds
I/O Wait Time (Locked)           :           29.2657 Seconds
I/O Wait Time (Unlocked)         :            0.3294 Seconds
CPU Time Memory Integral         :            3.5689 Mword-seconds
SDS Time Memory Integral         :            0.0000 Mword-seconds
I/O Wait Time Memory Integral    :           71.6911 Mword-seconds
Data Transferred                 :           20.0442 MWords
Maximum job memory used          :            3.5405 MWords
Maximum process memory used      :            3.2061 MWords
Maximum SDS used                 :            0.0000 MWords
Logical I/O Requests             :          923
Physical I/O Requests            :          930
Number of Commands               :            9

                                                 GAU Components
                                            Job Charge        :     0.00100
CPU Hours            :     0.00050          CPU Charge        :     0.00047
Avg Memory (mwd)     :     1.99485          Memory Charge     :     0.00001
Disk Activity (mwd)  :    20.04417          Disk Act Charge   :     0.00335

Charge before Queue Factor       :                              0.00483 GAUs
      (Excluding MSS/NTWK/TAGS) 
Multiplier for prem Queue        :            1.50 
Charged against Allocation       :                              0.00724 GAUs
logout