Changeset 2478 for palm/trunk/SOURCE


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

Bugfixes concerning top fluxes and TKE production

Location:
palm/trunk/SOURCE
Files:
3 edited

Legend:

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

    r2365 r2478  
    2525! -----------------
    2626! $Id$
     27! Bugfix, correct flag for use_top
     28!
     29! 2365 2017-08-21 14:59:59Z kanani
    2730! Vertical nesting implemented (SadiqHuq)
    2831!
     
    15541557          wall_flags_0(:,j,i) = IBSET( wall_flags_0(:,j,i), 9 )
    15551558          IF ( use_top_fluxes )                                                &
    1556              wall_flags_0(nzt:nzt+1,j,i) = IBCLR( wall_flags_0(nzt:nzt+1,j,i), 9 )
     1559             wall_flags_0(nzt+1,j,i) = IBCLR( wall_flags_0(nzt+1,j,i), 9 )
    15571560
    15581561          DO k = nzb+1, nzt
  • palm/trunk/SOURCE/production_e.f90

    r2329 r2478  
    2525! -----------------
    2626! $Id$
     27! Bugfix, consider case where no constant-flux layer and no surfaces fluxes
     28! are used
     29!
     30! 2329 2017-08-03 14:24:56Z knoop
    2731! Bugfix: added division by density as kinematic fluxes are needed
    2832!
     
    432436             ENDDO
    433437
    434           ELSEIF ( use_surface_fluxes )  THEN
     438          ELSE
    435439
    436440             DO  j = nys, nyn
     
    11821186          ENDDO
    11831187
    1184        ELSEIF ( use_surface_fluxes )  THEN
     1188       ELSE
    11851189!
    11861190!--       Calculate TKE production by shear. Here, no additional
     
    12351239!--                the bottom and top surface layer
    12361240                   DO  k = nzb+1, nzt
     1241 
    12371242                      tend(k,j,i) = tend(k,j,i) +                              &
    12381243                                    kh(k,j,i) * g / rho_reference *            &
     
    12411246                                      MERGE( 1.0_wp, 0.0_wp,                   &
    12421247                                             BTEST( wall_flags_0(k,j,i), 0 )   &
    1243                                            )                                  
     1248                                           )                   
    12441249                   ENDDO
    12451250
     
    13231328                                      MERGE( 1.0_wp, 0.0_wp,                   &
    13241329                                             BTEST( wall_flags_0(k,j,i), 0 )   &
    1325                                            )
     1330                                           )           
    13261331                   ENDDO
    13271332
  • 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.