30-day-ave.sh - Produce a thirty-day time average. --------------------------------------------------- The next job produces a thirty day time average for Case A, using do-loop notation to specify the days list. Note, that if only one case is processed, Case A is used for case dependent keywords. The HPROJ ICP generates 28 horizontal projection plots (4 fields 7 pressure levels), while the ZONAVG ICP produces 4 zonal mean cross-section plots. C C Example 4. - 30-day-ave.sh - Time average of 30 days. (.03295 GAUS) C TITLEA = 'Example 4. - 30-day-ave.sh - 30 DAY AVG (Days 20-49) ' TAPESA = '/CSM/ccm2/414/hist/h0003' C C NINTAPA automatically increments the input tape number when needed. C DAYSA list is in do-loop notation. C NINTAPA = 3,4,1 DAYSA = 20.0,49.0,1.0 C FIELDA1 = 'T','U','V','Q' PRESSLE = 1000.,850.,700.,500.,300.,200.,100. C C Request time averaging. C TIMAVGA = 'YES' C C Produce plots by using the horizontal projection option HPROJ. C HPROJ = 'RECT' C C Set the scale factor for meridional cross-section plots of Q. C MXSCAL = 'Q',1.E3 C C Compute, plot and print the zonal averages, send to laser printer. C ZONAVG = 'YES' ZAVGPRN = 'YES' DPLTMF = 'MP' C SAVHSTA = '/USERNAME/ccm2/414/hist/example4' PWDHSTA = 'WRPASS' C C BPHSTA = 'YES' allows blocked points to be written out. C SAVMHST specifies that the time averaged data should be written out. C BPHSTA = 'YES' SAVMHST = 'TAV' C ENDOFDATA -------------------------------------------------------------- Example 4 is the format for time averaging a series of model days. Here, the first input Mass Store file is /CSM/ccm2/414/hist/h0003 . NINTAPA automatically increments the input tape number when the next input tape is needed. Here, NINTAPA = 3,4,1 requests a total of 3 tapes whose names are formed by numerically incrementing the last 4 characters of the file specified by the ICP TAPESA by 1 when needed. The same effect could be obtained by specifying the three tape names with the TAPESA ICP. Do-loop notation is used for the DAYSA specification. Time averaging of all the fields is specified via the TIMAVGA = 'YES' ICP. The only difference between processing of individual days and generating a time average of the individual days is the value of TIMAVGA. DPLTMF = 'MP' sends the plots of the time averaged data to the SCD laser printer. The data is saved on a Mass Store file for future use by invoking the SAVHSTA keyword. PWDHSTA assigns a write password to the MSS file. The ICPs TIMAVGA = 'YES' and SAVMHST = 'TAV' are both required to ensure that only the time averaged data is to be saved to the tape. The default for SAVMHST is 'TSR' , which writes each individual time sample to the output tape, regardless of the value of TIMAVGA .