CALL SEQUENCE... spcsum IN OUT INPUT OUTPUT [-kind VAL1 [VAL2] \ -t1 VAL3 -name VAL4 -lv1 VAL5 \ -i VAL6 -j VAL7 -k VAL8 \ -label VAL9] DESCRIPTION... spcsum - PRINTS SPECTRAL SUMS OVER M OR N AUTHORS - S.Lambert AND T.Shepherd, MAY 13/83. LAST REVISION: $Header: /home/cvs/RMNLIB/utils/rdiag/lspgm/spcsum.ptn,v 4.20 2014/09/25 20:20:53 dugas Exp $ PURPOSE - SUMS A COMPLEX-VALUED SPECTRAL FILE SC(N,M) EITHER OVER THE PARAMETER M, PRINTING OUT THE RESULT SUM(N), OR OVER N, WITH THE RESULT SUM(M). IT ALSO PRODUCES A CROSS-SECTION FILE OF SPECTRAL DENSITY AS A FUNCTION OF WAVENUMBER AND PRESSURE. NOTE - IT CAN BE USED AS A DISPLAY PROGRAM FOR SPCDNS. INPUT FILE... IN = GLOBAL SPECTRAL COEFF FILE OUTPUT FILE... OUT = CROSS-SECTION OF HEIGHT VS. WAVENUMBER. INPUT PARAMETERS... SUMTYP (kind-1)= 'M', CAUSES SUM OVER N PRODUCING A FUNCTION OF M. = 'N', CAUSES SUM OVER M PRODUCING A FUNCTION OF N. VKIND (kind-2)= VERTICAL LEVEL TYPE. ONLY USE THIS WITH IN RPN STD FILES AND WHEN LEVX CONTAINS DECODED VALUES. FURTHERMORE, THIS PARAMETRE CAN ONLY BE SPECIFIED IN COMMAND-LINE MODE. THE CURRENTLY SUPPORTED VALUES ARE: 'HAUTEUR' , heights (m) with respect to sea level 'SIGMA' , classical terrain following (0.0 -> 1.0) 'ETA' , first version of the GEM ETA coordinate 'PRESSION' , pressure (mb) 'ARBITRAIRE', arbitrary level (no units apply) 'HAUTEURS' , heights (M) with respect to the surface 'HYBRIDE' , hybrid GEM coordinate 'THETA' , potential temperature coordinate 'GALCHEN' , GalChen generalized coordinate LABEL (label) = 256 CHARACTER LABEL PRINTED UNDER DISPLAY. NSTEP (t1) \ NAME (name) = RECORD LABEL WORDS 2,3,4 IDENTIFYING THE FIELD TO SUM. LEVEL (lv1) / IX (i) = 1, CAUSES THE CROSS-SECTION TO CONTAIN LOG10 OF THE SPECTRAL DENSITIES INSTEAD OF THE SPECTRAL DENSITIES THEMSELVES (WAVE NUMBER ZERO IS OMITTED). NWAV (j) = 1, CAUSES THE NUMBER OF WAVES ON FILE OUT TO BE TRUNCATED TO AN EVEN NUMBER. (THE DEFAULT IS 0; NO TRUNCATION). IPNT (k) = 1, PRINTS THE VALUES OF THE SPECTRUM, 0, NO PRINTING. NOTES: 1) IF PARAMETRES ARE TO BE PASSED VIA THE COMMAND LINE, ONLY ONE SET OF SUCH PARAMETRES CAN BE PROCESSED BY SPCSUM. 2) THE DEFAULT VALUES FOR BOTH NSTEP AND LEVEL IMPLY THAT THEY ARE NOT REALLY USED FOR THE DATA SELECTION. 3) A) IF NAME=ALL (THE DEFAULT), ALL RECORDS THAT SATISFY THE CRITERIA ON NSTEP AND LEVEL WILL BE TREATED. B) IF NAME=NEXT, ONLY THE FIRST RECORD IS TREATED AND THE VALUES OF NSTEP,LEVEL ARE IGNORED. C) IF NAME=SKIP, THE NEXT RECORD IS SKIPPED. THIS IS INVALID IN COMMAND-LINE MODE. READ( 5,5005 ) SUMTYP *ONLY ONCE* READ( 5,5015 ) LABEL *ONLY ONCE* READ( 5,5010 ) NSTEP,NAME,LEVEL,IX,NWAV,IPNT 5005 FORMAT(10X,4X,A1) 5010 FORMAT(10X,I10,1X,A4,4I5) 5015 FORMAT(64A4) EXAMPLES OF INPUT CARDS... SPCSUM M DISPLAY LABEL GOES HERE 36 VORT 500 1 0 1 NOTE THAT BLANK INPUT IS EQUIVALENCED TO THE FOLLOWING 3 CARDS SPCSUM N -1 ALL -1 0 1 0 0123456789012345678901234567890123456789012345678901234567890 1 2 3 4 5 6 7 EXIT CONDITIONS... 0 NO PROBLEMS 1 ERROR READING FILE IN 2 ERROR READING SUMTYP 3 ERROR READING LABEL 4 ERROR DECODING NSTEP PARAMETERS 5 NO MORE RECORDS TO SKIP (NAM WAS 'SKIP')