Ignore:
Timestamp:
Apr 13, 2018 11:52:24 AM (3 years ago)
Author:
suehring
Message:

Bugfixes in initalization of land-surface model as well as timeseries data output in case of elevated model surfaces.

File:
1 edited

Legend:

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

    r2963 r2968  
    2525! -----------------
    2626! $Id$
     27! Bugfix in initialization in case of elevated model surface
     28!
     29! 2963 2018-04-12 14:47:44Z suehring
    2730! - In initialization of surface types, consider the case that surface_fractions
    2831!   is not given in static input file.
     
    23742377
    23752378       USE indices,                                                            &
    2376            ONLY:  nx, ny
     2379           ONLY:  nx, ny, topo_min_level
    23772380
    23782381       USE pmc_interface,                                                      &
     
    24082411       IF (  .NOT.  cloud_physics )  THEN
    24092412          CALL calc_mean_profile( pt, 4 )
    2410           rho_surface = surface_pressure * 100.0_wp / ( r_d * hom(nzb+1,1,4,0) * exn )
     2413          rho_surface = surface_pressure * 100.0_wp / ( r_d * hom(topo_min_level+1,1,4,0) * exn )
    24112414       ENDIF
    24122415
     
    24242427       tt_soil_h_m%var_2d    = 0.0_wp
    24252428       tm_soil_h_m%var_2d    = 0.0_wp
    2426        tm_liq_h_m%var_1d  = 0.0_wp
    2427        surf_lsm_h%c_liq = 0.0_wp
     2429       tm_liq_h_m%var_1d     = 0.0_wp
     2430       surf_lsm_h%c_liq      = 0.0_wp
    24282431
    24292432       surf_lsm_h%ghf = 0.0_wp
     
    24432446          tt_soil_v_m(l)%var_2d    = 0.0_wp
    24442447          tm_soil_v_m(l)%var_2d    = 0.0_wp
    2445           tm_liq_v_m(l)%var_1d  = 0.0_wp
    2446           surf_lsm_v(l)%c_liq = 0.0_wp
     2448          tm_liq_v_m(l)%var_1d     = 0.0_wp
     2449          surf_lsm_v(l)%c_liq      = 0.0_wp
    24472450
    24482451          surf_lsm_v(l)%ghf = 0.0_wp
Note: See TracChangeset for help on using the changeset viewer.