Ignore:
Timestamp:
Nov 28, 2011 9:47:19 AM (10 years ago)
Author:
raasch
Message:

scalar quantities can be excluded from Rayleigh damping; bugfix for long lines in configuration file with more than 300 characters

File:
1 edited

Legend:

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

    r737 r785  
    44! Current revisions:
    55! -----------------
     6! new factor rdf_sc allows separate Rayleigh damping of scalars
    67!
    78! Former revisions:
     
    511512                                     sbt * tend(k,j,i) + tsc(3) * tpt_m(k,j,i) &
    512513                                   ) -                                         &
    513                            tsc(5) * rdf(k) * ( pt(k,j,i) - pt_init(k) )
     514                           tsc(5) * rdf_sc(k) * ( pt(k,j,i) - pt_init(k) )
    514515          ENDDO
    515516
     
    602603                                     sbt * tend(k,j,i) + tsc(3) * tsa_m(k,j,i) &
    603604                                      ) -                                      &
    604                               tsc(5) * rdf(k) * ( sa(k,j,i) - sa_init(k) )
     605                              tsc(5) * rdf_sc(k) * ( sa(k,j,i) - sa_init(k) )
    605606                IF ( sa_p(k,j,i) < 0.0 )  sa_p(k,j,i) = 0.1 * sa(k,j,i)
    606607             ENDDO
     
    723724                                      sbt * tend(k,j,i) + tsc(3) * tq_m(k,j,i) &
    724725                                     ) -                                       &
    725                              tsc(5) * rdf(k) * ( q(k,j,i) - q_init(k) )
     726                             tsc(5) * rdf_sc(k) * ( q(k,j,i) - q_init(k) )
    726727                IF ( q_p(k,j,i) < 0.0 )  q_p(k,j,i) = 0.1 * q(k,j,i)
    727728             ENDDO
     
    12161217                               tsc(2) * tend(k,j,i) + tsc(3) * tpt_m(k,j,i) &
    12171218                                   ) -                                      &
    1218                            tsc(5) * rdf(k) * ( pt(k,j,i) - pt_init(k) )
     1219                           tsc(5) * rdf_sc(k) * ( pt(k,j,i) - pt_init(k) )
    12191220          ENDDO
    12201221
     
    12661267                               tsc(2) * tend(k,j,i) + tsc(3) * tsa_m(k,j,i) &
    12671268                                      ) -                                   &
    1268                              tsc(5) * rdf(k) * ( sa(k,j,i) - sa_init(k) )
     1269                             tsc(5) * rdf_sc(k) * ( sa(k,j,i) - sa_init(k) )
    12691270                IF ( sa_p(k,j,i) < 0.0 )  sa_p(k,j,i) = 0.1 * sa(k,j,i)
    12701271             ENDDO
     
    13461347                                tsc(2) * tend(k,j,i) + tsc(3) * tq_m(k,j,i) &
    13471348                                     ) -                                    &
    1348                              tsc(5) * rdf(k) * ( q(k,j,i) - q_init(k) )
     1349                             tsc(5) * rdf_sc(k) * ( q(k,j,i) - q_init(k) )
    13491350                IF ( q_p(k,j,i) < 0.0 )  q_p(k,j,i) = 0.1 * q(k,j,i)
    13501351             ENDDO
     
    18461847                                     sbt * tend(k,j,i) + tsc(3) * tpt_m(k,j,i) &
    18471848                                   ) -                                         &
    1848                            tsc(5) * rdf(k) * ( pt(k,j,i) - pt_init(k) )
     1849                           tsc(5) * rdf_sc(k) * ( pt(k,j,i) - pt_init(k) )
    18491850          ENDDO
    18501851       ENDDO
     
    19401941                                     sbt * tend(k,j,i) + tsc(3) * tsa_m(k,j,i) &
    19411942                                      ) -                                      &
    1942                               tsc(5) * rdf(k) * ( sa(k,j,i) - sa_init(k) )
     1943                              tsc(5) * rdf_sc(k) * ( sa(k,j,i) - sa_init(k) )
    19431944                IF ( sa_p(k,j,i) < 0.0 )  sa_p(k,j,i) = 0.1 * sa(k,j,i)
    19441945             ENDDO
     
    20652066                                      sbt * tend(k,j,i) + tsc(3) * tq_m(k,j,i) &
    20662067                                     ) -                                       &
    2067                              tsc(5) * rdf(k) * ( q(k,j,i) - q_init(k) )
     2068                             tsc(5) * rdf_sc(k) * ( q(k,j,i) - q_init(k) )
    20682069                IF ( q_p(k,j,i) < 0.0 )  q_p(k,j,i) = 0.1 * q(k,j,i)
    20692070             ENDDO
Note: See TracChangeset for help on using the changeset viewer.