 Timestamp:
 Apr 8, 2019 11:52:36 AM (2 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/land_surface_model_mod.f90
r3856 r3868 25 25 !  26 26 ! $Id$ 27 ! More strict limitation of roughness length when it is in the order of the 28 ! vertical grid spacing 29 ! 30 ! 3856 20190403 11:06:59Z suehring 27 31 ! Bugfix in lsm_init in case no surfacefractions are provided 28 32 ! … … 4879 4883 ! necessary. 4880 4884 DO m = 1, surf_lsm_h%ns 4881 IF ( surf_lsm_h%z0(m) > =surf_lsm_h%z_mo(m) ) THEN4885 IF ( surf_lsm_h%z0(m) > 0.5_wp * surf_lsm_h%z_mo(m) ) THEN 4882 4886 4883 surf_lsm_h%z0(m) = 0. 9_wp * surf_lsm_h%z_mo(m)4887 surf_lsm_h%z0(m) = 0.5_wp * surf_lsm_h%z_mo(m) 4884 4888 4885 4889 WRITE( message_string, * ) 'z0 exceeds surfacelayer height ' // & … … 4890 4894 0, 0, 0, 6, 0 ) 4891 4895 ENDIF 4892 IF ( surf_lsm_h%z0h(m) > =surf_lsm_h%z_mo(m) ) THEN4896 IF ( surf_lsm_h%z0h(m) > 0.5_wp * surf_lsm_h%z_mo(m) ) THEN 4893 4897 4894 surf_lsm_h%z0h(m) = 0. 9_wp * surf_lsm_h%z_mo(m)4895 surf_lsm_h%z0q(m) = 0. 9_wp * surf_lsm_h%z_mo(m)4898 surf_lsm_h%z0h(m) = 0.5_wp * surf_lsm_h%z_mo(m) 4899 surf_lsm_h%z0q(m) = 0.5_wp * surf_lsm_h%z_mo(m) 4896 4900 4897 4901 WRITE( message_string, * ) 'z0h exceeds surfacelayer height ' // & … … 4906 4910 DO l = 0, 3 4907 4911 DO m = 1, surf_lsm_v(l)%ns 4908 IF ( surf_lsm_v(l)%z0(m) > =surf_lsm_v(l)%z_mo(m) ) THEN4912 IF ( surf_lsm_v(l)%z0(m) > 0.5_wp * surf_lsm_v(l)%z_mo(m) ) THEN 4909 4913 4910 surf_lsm_v(l)%z0(m) = 0. 9_wp * surf_lsm_v(l)%z_mo(m)4914 surf_lsm_v(l)%z0(m) = 0.5_wp * surf_lsm_v(l)%z_mo(m) 4911 4915 4912 4916 WRITE( message_string, * ) 'z0 exceeds surfacelayer height '//& … … 4918 4922 0, 0, 0, 6, 0 ) 4919 4923 ENDIF 4920 IF ( surf_lsm_v(l)%z0h(m) > =surf_lsm_v(l)%z_mo(m) ) THEN4924 IF ( surf_lsm_v(l)%z0h(m) > 0.5_wp * surf_lsm_v(l)%z_mo(m) ) THEN 4921 4925 4922 surf_lsm_v(l)%z0h(m) = 0. 9_wp * surf_lsm_v(l)%z_mo(m)4923 surf_lsm_v(l)%z0q(m) = 0. 9_wp * surf_lsm_v(l)%z_mo(m)4926 surf_lsm_v(l)%z0h(m) = 0.5_wp * surf_lsm_v(l)%z_mo(m) 4927 surf_lsm_v(l)%z0q(m) = 0.5_wp * surf_lsm_v(l)%z_mo(m) 4924 4928 4925 4929 WRITE( message_string, * ) 'z0h exceeds surfacelayer height '//&
Note: See TracChangeset
for help on using the changeset viewer.