Ignore:
Timestamp:
Oct 30, 2007 12:12:24 PM (14 years ago)
Author:
letzel
Message:

prognostic_equations include the respective wall_*flux in the parameter list of
calls of diffusion_s. Same as before, only the values of wall_heatflux(0:4)
can be assigned. At present, wall_humidityflux, wall_qflux, wall_salinityflux,
wall_scalarflux are kept zero. diffusion_s uses the respective wall_*flux
instead of wall_heatflux. This update serves two purposes:

  • it avoids errors in calculations with humidity/scalar/salinity and prescribed

non-zero wall_heatflux,

  • it prepares PALM for a possible assignment of wall fluxes of

humidity/scalar/salinity in a future release.

Bugfix: assignment of fluxes at walls

Updates to documentation:
chapter_4.2.html#mode_dvrp
chapter_3.5.4.html#time_series

Default for mrun options -q and -n is "sla3" for lctit. Queues bes1 and bes2
removed. DOC/misc/Tsubame.html updated.

Modified default paths (/work/...) for lctit in .mrun.config.default

File:
1 edited

Legend:

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

    r110 r129  
    44! Actual revisions:
    55! -----------------
    6 !
     6! add wall_*flux to parameter list of calls of diffusion_s
    77!
    88! Former revisions:
     
    369369!--       Tendency terms
    370370          IF ( scalar_advec == 'bc-scheme' )  THEN
    371              CALL diffusion_s( i, j, ddzu, ddzw, kh, pt, shf, tswst, tend )
     371             CALL diffusion_s( i, j, ddzu, ddzw, kh, pt, shf, tswst, &
     372                               wall_heatflux, tend )
    372373          ELSE
    373374             IF ( tsc(2) == 2.0  .OR.  timestep_scheme(1:5) == 'runge' )  THEN
     
    383384             THEN
    384385                CALL diffusion_s( i, j, ddzu, ddzw, kh_m, pt_m, shf_m, &
    385                                   tswst_m, tend )
     386                                  tswst_m, wall_heatflux, tend )
    386387             ELSE
    387                 CALL diffusion_s( i, j, ddzu, ddzw, kh, pt, shf, tswst, tend )
     388                CALL diffusion_s( i, j, ddzu, ddzw, kh, pt, shf, tswst, &
     389                                  wall_heatflux, tend )
    388390             ENDIF
    389391          ENDIF
     
    471473             IF ( scalar_advec == 'bc-scheme' )  THEN
    472474                CALL diffusion_s( i, j, ddzu, ddzw, kh, sa, saswsb, saswst, &
    473                                   tend )
     475                                  wall_salinityflux, tend )
    474476             ELSE
    475477                IF ( tsc(2) == 2.0  .OR.  timestep_scheme(1:5) == 'runge' ) THEN
     
    483485                ENDIF
    484486                CALL diffusion_s( i, j, ddzu, ddzw, kh, sa, saswsb, saswst, &
    485                                   tend )
     487                                  wall_salinityflux, tend )
    486488             ENDIF
    487489       
     
    563565!--          Tendency-terms
    564566             IF ( scalar_advec == 'bc-scheme' )  THEN
    565                 CALL diffusion_s( i, j, ddzu, ddzw, kh, q, qsws, qswst, tend )
     567                CALL diffusion_s( i, j, ddzu, ddzw, kh, q, qsws, qswst, &
     568                                  wall_qflux, tend )
    566569             ELSE
    567570                IF ( tsc(2) == 2.0  .OR.  timestep_scheme(1:5) == 'runge' ) THEN
     
    577580                THEN
    578581                   CALL diffusion_s( i, j, ddzu, ddzw, kh_m, q_m, qsws_m, &
    579                                      qswst_m, tend )
     582                                     qswst_m, wall_qflux, tend )
    580583                ELSE
    581584                   CALL diffusion_s( i, j, ddzu, ddzw, kh, q, qsws, qswst, &
    582                                      tend )
     585                                     wall_qflux, tend )
    583586                ENDIF
    584587             ENDIF
     
    974977          THEN
    975978             CALL diffusion_s( i, j, ddzu, ddzw, kh_m, pt_m, shf_m, &
    976                                tswst_m, tend )
    977           ELSE
    978              CALL diffusion_s( i, j, ddzu, ddzw, kh, pt, shf, tswst, tend )
     979                               tswst_m, wall_heatflux, tend )
     980          ELSE
     981             CALL diffusion_s( i, j, ddzu, ddzw, kh, pt, shf, tswst, &
     982                               wall_heatflux, tend )
    979983          ENDIF
    980984
     
    10331037             ENDIF
    10341038             CALL diffusion_s( i, j, ddzu, ddzw, kh, sa, saswsb, saswst, &
    1035                                tend )
     1039                               wall_salinityflux, tend )
    10361040
    10371041             CALL user_actions( i, j, 'sa-tendency' )
     
    10871091             THEN
    10881092                CALL diffusion_s( i, j, ddzu, ddzw, kh_m, q_m, qsws_m, &
    1089                                   qswst_m, tend )
     1093                                  qswst_m, wall_qflux, tend )
    10901094             ELSE
    10911095                CALL diffusion_s( i, j, ddzu, ddzw, kh, q, qsws, qswst, &
    1092                                   tend )
     1096                                  wall_qflux, tend )
    10931097             ENDIF
    10941098       
     
    14871491!-- pt-tendency terms with no communication
    14881492    IF ( scalar_advec == 'bc-scheme' )  THEN
    1489        CALL diffusion_s( ddzu, ddzw, kh, pt, shf, tswst, tend )
     1493       CALL diffusion_s( ddzu, ddzw, kh, pt, shf, tswst, wall_heatflux, &
     1494                         tend )
    14901495    ELSE
    14911496       IF ( tsc(2) == 2.0  .OR.  timestep_scheme(1:5) == 'runge' )  THEN
     
    14991504       ENDIF
    15001505       IF ( tsc(2) == 2.0  .AND.  timestep_scheme(1:8) == 'leapfrog' )  THEN
    1501           CALL diffusion_s( ddzu, ddzw, kh_m, pt_m, shf_m, tswst_m, tend )
     1506          CALL diffusion_s( ddzu, ddzw, kh_m, pt_m, shf_m, tswst_m, &
     1507                            wall_heatflux, tend )
    15021508       ELSE
    1503           CALL diffusion_s( ddzu, ddzw, kh, pt, shf, tswst, tend )
     1509          CALL diffusion_s( ddzu, ddzw, kh, pt, shf, tswst, wall_heatflux, &
     1510                            tend )
    15041511       ENDIF
    15051512    ENDIF
     
    15891596
    15901597!
    1591 !--    Scalar/q-tendency terms with no communication
     1598!--    sa-tendency terms with no communication
    15921599       IF ( scalar_advec == 'bc-scheme' )  THEN
    1593           CALL diffusion_s( ddzu, ddzw, kh, sa, saswsb, saswst, tend )
     1600          CALL diffusion_s( ddzu, ddzw, kh, sa, saswsb, saswst, &
     1601                            wall_salinityflux, tend )
    15941602       ELSE
    15951603          IF ( tsc(2) == 2.0  .OR.  timestep_scheme(1:5) == 'runge' )  THEN
     
    16021610             ENDIF
    16031611          ENDIF
    1604           CALL diffusion_s( ddzu, ddzw, kh, sa, saswsb, saswst, tend )
     1612          CALL diffusion_s( ddzu, ddzw, kh, sa, saswsb, saswst, &
     1613                            wall_salinityflux, tend )
    16051614       ENDIF
    16061615       
     
    16891698!--    Scalar/q-tendency terms with no communication
    16901699       IF ( scalar_advec == 'bc-scheme' )  THEN
    1691           CALL diffusion_s( ddzu, ddzw, kh, q, qsws, qswst, tend )
     1700          CALL diffusion_s( ddzu, ddzw, kh, q, qsws, qswst, wall_qflux, tend )
    16921701       ELSE
    16931702          IF ( tsc(2) == 2.0  .OR.  timestep_scheme(1:5) == 'runge' )  THEN
     
    17011710          ENDIF
    17021711          IF ( tsc(2) == 2.0  .AND.  timestep_scheme(1:8) == 'leapfrog' )  THEN
    1703              CALL diffusion_s( ddzu, ddzw, kh_m, q_m, qsws_m, qswst_m, tend )
    1704           ELSE
    1705              CALL diffusion_s( ddzu, ddzw, kh, q, qsws, qswst, tend )
     1712             CALL diffusion_s( ddzu, ddzw, kh_m, q_m, qsws_m, qswst_m, &
     1713                               wall_qflux, tend )
     1714          ELSE
     1715             CALL diffusion_s( ddzu, ddzw, kh, q, qsws, qswst, &
     1716                               wall_qflux, tend )
    17061717          ENDIF
    17071718       ENDIF
Note: See TracChangeset for help on using the changeset viewer.