Changeset 2256 for palm


Ignore:
Timestamp:
Jun 7, 2017 1:58:08 PM (7 years ago)
Author:
suehring
Message:

change default value of zeta_min; enable scalar/heat/water-fluxes at downward-facing surfaces, bugfix Makefile

Location:
palm/trunk/SOURCE
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/Makefile

    r2238 r2256  
    2525# -----------------
    2626# $Id$
     27# Remove ring dependency in init_pegrid
     28#
     29# 2238 2017-05-31 16:49:16Z suehring
    2730# Bugfix, further missing dependency on surface_mod
    2831#
     
    493496init_masks.o: modules.o mod_kinds.o netcdf_interface_mod.o
    494497init_ocean.o: modules.o eqn_state_seawater.o mod_kinds.o
    495 init_pegrid.o: modules.o mod_kinds.o spectra_mod.o pmc_interface_mod.o
     498init_pegrid.o: modules.o mod_kinds.o spectra_mod.o
    496499init_pt_anomaly.o: modules.o mod_kinds.o
    497500init_rankine.o: modules.o mod_kinds.o
  • palm/trunk/SOURCE/modules.f90

    r2233 r2256  
    2525! -----------------
    2626! $Id$
     27! Change default value of zeta_min to -20
     28! Increase dimension for wall_heatflux, etc.
     29!
     30! 2233 2017-05-30 18:08:54Z suehring
    2731!
    2832! 2232 2017-05-30 17:47:52Z suehring
     
    13311335    REAL(wp) ::  wall_adjustment_factor = 1.8_wp               !<
    13321336    REAL(wp) ::  zeta_max = 20.0_wp                            !< Maximum value of zeta = z/L
    1333     REAL(wp) ::  zeta_min = -9990.0_wp                         !< Minimum value of zeta = z/L
     1337    REAL(wp) ::  zeta_min = -20.0_wp                           !< Minimum value of zeta = z/L
    13341338    REAL(wp) ::  z_max_do2d = -1.0_wp                          !<
    13351339    REAL(wp) ::  z0h_factor = 1.0_wp                           !<
     
    13641368    REAL(wp) ::  volume_flow_area(1:3) = 0.0_wp                    !<
    13651369    REAL(wp) ::  volume_flow_initial(1:3) = 0.0_wp                 !<
    1366     REAL(wp) ::  wall_heatflux(0:4) = 0.0_wp                       !<
    1367     REAL(wp) ::  wall_humidityflux(0:4) = 0.0_wp                   !<
    1368     REAL(wp) ::  wall_salinityflux(0:4) = 0.0_wp                   !<
    1369     REAL(wp) ::  wall_scalarflux(0:4) = 0.0_wp                     !<
     1370    REAL(wp) ::  wall_heatflux(0:5) = 0.0_wp                       !<
     1371    REAL(wp) ::  wall_humidityflux(0:5) = 0.0_wp                   !<
     1372    REAL(wp) ::  wall_salinityflux(0:5) = 0.0_wp                   !<
     1373    REAL(wp) ::  wall_scalarflux(0:5) = 0.0_wp                     !<
    13701374    REAL(wp) ::  subs_vertical_gradient(10) = 0.0_wp               !<
    13711375    REAL(wp) ::  subs_vertical_gradient_level(10) = -9999999.9_wp  !<
  • palm/trunk/SOURCE/surface_mod.f90

    r2233 r2256  
    2424! Former revisions:
    2525! -----------------
    26 
    27 ! Revision 1.1  2016/12/20 10:06:06  suehring
     26! $Id$
     27! Enable heating at downward-facing surfaces
     28!
     29! 2233 2017-05-30 18:08:54Z suehring
    2830! Initial revision
    2931!
     
    13511353                   ENDIF
    13521354!
    1353 !--             Set heat-flux at downward-facing surfaces to zero
     1355!--             Set heat-flux at downward-facing surfaces
    13541356                ELSE
    1355                    surf%shf(num_h) = 0.0_wp
     1357                   surf%shf(num_h) = wall_heatflux(5) *                        &
     1358                                             heatflux_input_conversion(k)
    13561359                ENDIF
    13571360
     
    13691372                      ENDIF
    13701373                   ELSE
    1371                       surf%qsws(num_h) = 0.0_wp
     1374                      surf%qsws(num_h) = wall_humidityflux(5) *                &
     1375                                             heatflux_input_conversion(k)
    13721376                   ENDIF
    13731377                ENDIF
     
    13791383
    13801384                         IF ( k-1 /= 0 )                                       &
    1381                             surf%ssws(num_h) = wall_scalarflux(0) *            &
    1382                                                   waterflux_input_conversion(k-1)
     1385                            surf%ssws(num_h) = wall_scalarflux(0)
    13831386
    13841387                      ELSE
     
    13861389                      ENDIF
    13871390                   ELSE
    1388                       surf%ssws(num_h) = 0.0_wp
     1391                      surf%ssws(num_h) = wall_scalarflux(5)
    13891392                   ENDIF
    13901393                ENDIF
Note: See TracChangeset for help on using the changeset viewer.