CALL SEQUENCE...
       spldiss  IN OUT INPUT OUTPUT [ -n VAL1 -a VAL2 -b VAL3 ]
 
 
      DESCRIPTION...
       spldiss - COMPUTES SPECTRAL LEITH DISSIPATION FUNCTION
 
 
      AUTHOR   - T.Shepherd, June 01/79.
      LAST REVISION: $Header: Autres-corrections-d-usages-de-MAXW-NWDS-modifs-dans-GGSIGN-SIGNAG @ 2018-06-05 17:00:06 -0400 (dugas) $
 
 
      PURPOSE - COMPUTES THE LEITH DISSIPATION FUNCTION L(M,N) FOR EACH FIELD
                IN A GLOBAL SPECTRAL FILE. THIS FUNCTION IS DEFINED AS:
 
                      L(M,N) = ETA3*F(N) FOR N>NZ
 
                WHERE,
 
                      F(N)   = 4*(N-NZ)/NSTAR
                      NZ     = 0.55*NSTAR,
                      ETA3   = THE CUBIC ROOT OF THE ENSTROPHY FLUX RATE,
                AND   NSTAR  = THE TRUNCATION WAVENUMBER.
 
      NOTE    - OUTPUTED RESULTS ARE ON AS MANY LEVELS AND FOR AS MANY
                STEPS AS THE INPUT FILE HAS.
 
 
      INPUT FILE...
       IN  = GLOBAL SPECTRAL FIELDS
 
      OUTPUT FILE...
       OUT = CORRESPONDING FILE OF LEITH DISSIPATION FOR EACH WAVENUMBER
 
 
      INPUT PARAMETERS...
 
       NSTAR (n) = THE TRUNCATION WAVENUMBER.
       ETA3  (a) = THE CUBIC ROOT OF THE ENSTROPHY FLUX RATE,
       TEST  (b) = 0, THE NEGATIVE VALUES IN THE DISSIPATION FUNCTION ARE SET
                      TO ZERO.
                 <>0, (OTHERWISE) THE NEGATIVE VALUES ARE RETAINED.
 
       READ( 5,5010 ) NSTAR,ETA3,TEST
 5010  FORMAT(10X,I5,E10.0,I5)
 
 
      EXAMPLE OF INPUT CARD...
 
   SPLDISS   15       1.0    0
 
         0123456789012345678901234567890123456789012345678901234567890
         1         2         3         4         5         6         7
 
      EXIT CONDITIONS...
        0 NO PROBLEMS
        1 ERROR READING FILE IN
        2 ERROR READING INPUT DIRECTIVES