Ignore:
Timestamp:
Sep 18, 2017 1:37:24 PM (4 years ago)
Author:
suehring
Message:

Bugfixes concerning top fluxes and TKE production

File:
1 edited

Legend:

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

    r2378 r2478  
    2525! -----------------
    2626! $Id$
     27! Bugfixes in initializing model top
     28!
     29! 2378 2017-08-31 13:57:27Z suehring
    2730! Bugfix in write restart data
    2831!
     
    14921495!--          Initialize top heat flux
    14931496             IF ( constant_top_heatflux )                                      &
    1494                 surf%shf = top_heatflux * heatflux_input_conversion(nzt+1)
     1497                surf%shf(num_h) = top_heatflux * heatflux_input_conversion(nzt+1)
    14951498!
    14961499!--          Initialization in case of a coupled model run
    14971500             IF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
    1498                 surf%shf = 0.0_wp
     1501                surf%shf(num_h) = 0.0_wp
    14991502             ENDIF
    15001503!
    15011504!--          Prescribe latent heat flux at the top     
    15021505             IF ( humidity )  THEN
    1503              surf%qsws = 0.0_wp
     1506             surf%qsws(num_h) = 0.0_wp
    15041507                IF ( cloud_physics  .AND.  microphysics_morrison ) THEN
    1505                    surf%ncsws = 0.0_wp
    1506                    surf%qcsws = 0.0_wp
     1508                   surf%ncsws(num_h) = 0.0_wp
     1509                   surf%qcsws(num_h) = 0.0_wp
    15071510                ENDIF
    15081511                IF ( cloud_physics  .AND.  microphysics_seifert ) THEN
    1509                    surf%nrsws = 0.0_wp
    1510                    surf%qrsws = 0.0_wp
     1512                   surf%nrsws(num_h) = 0.0_wp
     1513                   surf%qrsws(num_h) = 0.0_wp
    15111514                ENDIF
    15121515             ENDIF
     
    15141517!--          Prescribe top scalar flux
    15151518             IF ( passive_scalar .AND. constant_top_scalarflux )               &
    1516                 surf%ssws = top_scalarflux
     1519                surf%ssws(num_h) = top_scalarflux
    15171520!
    15181521!--          Prescribe top salinity flux
    1519              IF ( ocean .AND. constant_top_salinityflux)                          &
    1520                 surf%sasws = top_salinityflux
     1522             IF ( ocean .AND. constant_top_salinityflux)                       &
     1523                surf%sasws(num_h) = top_salinityflux
    15211524!
    15221525!--          Initialization in case of a coupled model run
    15231526             IF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
    1524                 surf%shf = 0.0_wp
     1527                surf%shf(num_h) = 0.0_wp
    15251528             ENDIF
    15261529!
    15271530!--          Top momentum fluxes
    1528              surf%usws = top_momentumflux_u * momentumflux_input_conversion(nzt+1)
    1529              surf%vsws = top_momentumflux_v * momentumflux_input_conversion(nzt+1)
     1531             surf%usws(num_h) = top_momentumflux_u *                           &
     1532                                            momentumflux_input_conversion(nzt+1)
     1533             surf%vsws(num_h) = top_momentumflux_v *                           &
     1534                                            momentumflux_input_conversion(nzt+1)
    15301535!
    15311536!--          Increment surface indices
Note: See TracChangeset for help on using the changeset viewer.