Ignore:
Timestamp:
Aug 13, 2008 5:09:26 PM (16 years ago)
Author:
letzel
Message:
  • changed: allow 100 spectra levels instead of 10 for consistency with define_netcdf_header
File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/data_output_spectra.f90

    r144 r189  
    44! Actual revisions:
    55! -----------------
     6! allow 100 spectra levels instead of 10 for consistency with
     7! define_netcdf_header
    68! user-defined spectra
    79!
     
    270272    INTEGER                  :: cranz, labforx = 3, labfory = 3, legpos = 3, &
    271273                                timodex = 1
    272     INTEGER, DIMENSION(1:10) :: cucol = 1, klist = 999999, lstyle = 0
     274    INTEGER, DIMENSION(1:100):: cucol = 1, klist = 999999, lstyle = 0
    273275    LOGICAL                  :: datleg = .TRUE., grid = .TRUE., &
    274276                                lclose = .TRUE., rand = .TRUE., &
     
    278280    REAL                     :: gwid = 0.1, rlegfak = 0.7, uxmin, uxmax, &
    279281                                uymin, uymax
    280     REAL, DIMENSION(1:10)    :: lwid = 0.6
    281     REAL, DIMENSION(10)      :: uyma, uymi
     282    REAL, DIMENSION(1:100)   :: lwid = 0.6
     283    REAL, DIMENSION(100)     :: uyma, uymi
    282284
    283285    NAMELIST /RAHMEN/  cranz, datleg, rtext, swap
     
    301303!-- different quantities u, v, w, pt and q)
    302304    DO  k = 1, n_sp_x
    303        IF ( k < 10 )  THEN
     305       IF ( k < 100 )  THEN
    304306          IF ( pr == 3 )  THEN
    305307             WRITE ( 82, 100 )  '#', k, header_char( pr ),        &
     
    351353!-- Calculate klist-values from the available comp_spectra_level values
    352354    i = 1; k = 1
    353     DO WHILE ( i <= 10  .AND.  plot_spectra_level(i) /= 999999 )
     355    DO WHILE ( i <= 100  .AND.  plot_spectra_level(i) /= 999999 )
    354356       DO WHILE ( k <= n_sp_x  .AND. &
    355357                  plot_spectra_level(i) >= comp_spectra_level(k) )
     
    398400100 FORMAT (A,I1,1X,A,1X,I4,'m ',A)
    399401101 FORMAT (A,I2,1X,A,1X,I4,'m ',A)
    400 102 FORMAT (E15.7,10(1X,E15.7))
     402102 FORMAT (E15.7,100(1X,E15.7))
    401403103 FORMAT ('NEXT')
    402404104 FORMAT ('time averaged over',F7.1,' s')
     
    429431    INTEGER                  :: cranz, labforx = 3, labfory = 3, legpos = 3, &
    430432                                timodex = 1
    431     INTEGER, DIMENSION(1:10) :: cucol = 1, klist = 999999, lstyle = 0
     433    INTEGER, DIMENSION(1:100):: cucol = 1, klist = 999999, lstyle = 0
    432434    LOGICAL                  :: datleg = .TRUE., grid = .TRUE., &
    433435                                lclose = .TRUE., rand = .TRUE., &
     
    437439    REAL                     :: gwid = 0.1, rlegfak = 0.7, uxmin, uxmax, &
    438440                                uymin, uymax
    439     REAL, DIMENSION(1:10)    :: lwid = 0.6
    440     REAL, DIMENSION(10)      :: uyma, uymi
     441    REAL, DIMENSION(1:100)   :: lwid = 0.6
     442    REAL, DIMENSION(100)     :: uyma, uymi
    441443
    442444    NAMELIST /RAHMEN/  cranz, datleg, rtext, swap
     
    460462!-- different quantities u, v, w, pt and q)
    461463    DO  k = 1, n_sp_y
    462        IF ( k < 10 )  THEN
     464       IF ( k < 100 )  THEN
    463465          IF ( pr == 3 ) THEN
    464466             WRITE ( 84, 100 )  '#', k, header_char( pr ),        &
     
    510512!-- Calculate klist-values from the available comp_spectra_level values
    511513    j = 1; k = 1
    512     DO WHILE ( j <= 10  .AND.  plot_spectra_level(j) /= 999999 )
     514    DO WHILE ( j <= 100  .AND.  plot_spectra_level(j) /= 999999 )
    513515       DO WHILE ( k <= n_sp_y  .AND. &
    514516                  plot_spectra_level(j) >= comp_spectra_level(k) )
     
    557559100 FORMAT (A,I1,1X,A,1X,I4,'m ',A)
    558560101 FORMAT (A,I2,1X,A,1X,I4,'m ',A)
    559 102 FORMAT (E15.7,10(1X,E15.7))
     561102 FORMAT (E15.7,100(1X,E15.7))
    560562103 FORMAT ('NEXT')
    561563104 FORMAT ('time averaged over',F7.1,' s')
Note: See TracChangeset for help on using the changeset viewer.