Changeset 4365 for palm/trunk/SOURCE/init_3d_model.f90
- Timestamp:
- Jan 8, 2020 2:29:20 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_3d_model.f90
r4360 r4365 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Overwrite u_init, v_init, pt_init, q_init and s_init with hom for all 28 ! cyclic_fill-cases, not only for turbulent_inflow = .TRUE. 29 ! 30 ! 4360 2020-01-07 11:25:50Z suehring 27 31 ! Introduction of wall_flags_total_0, which currently sets bits based on static 28 32 ! topography information used in wall_flags_static_0 … … 1122 1126 #endif 1123 1127 ENDDO 1124 ! 1125 !-- In case of cyclic fill, call init_surfaces a second time, so that 1126 !-- surface properties such as heat fluxes are initialized as prescribed. 1127 IF ( TRIM( initializing_actions ) == 'cyclic_fill' ) & 1128 1129 1130 IF ( TRIM( initializing_actions ) == 'cyclic_fill' ) THEN 1131 1132 ! 1133 !-- In case of cyclic fill, call init_surfaces a second time, so that 1134 !-- surface properties such as heat fluxes are initialized as prescribed. 1128 1135 CALL init_surfaces 1129 1136 1137 ! 1138 !-- Overwrite u_init, v_init, pt_init, q_init and s_init with the 1139 !-- horizontally mean (hom) vertical profiles from the end 1140 !-- of the prerun, because these profiles shall be used as the reference 1141 !-- state for the rayleigh damping and the pt_damping. This is especially 1142 !-- important for the use of large_scale_subsidence, because the 1143 !-- reference temperature in the free atmosphere changes in time. 1144 u_init(:) = hom_sum(:,1,0) 1145 v_init(:) = hom_sum(:,2,0) 1146 pt_init(:) = hom_sum(:,4,0) 1147 IF ( humidity ) & 1148 q_init(:) = hom_sum(:,41,0) 1149 IF ( passive_scalar ) & 1150 s_init(:) = hom_sum(:,115,0) 1151 ENDIF 1130 1152 ! 1131 1153 !-- In case of complex terrain and cyclic fill method as initialization, … … 1161 1183 !-- profiles from the prerun. Alternatively, prescribed profiles 1162 1184 !-- for u,v-components can be used. 1163 !-- Overwrite u_init, v_init, pt_init, q_init and s_init with the1164 !-- (temporally) and horizontally averaged vertical profiles from the end1165 !-- of the prerun, because these profiles shall be used as the basic state1166 !-- for the rayleigh damping and the pt_damping.1167 1185 ALLOCATE( mean_inflow_profiles(nzb:nzt+1,1:num_mean_inflow_profiles) ) 1168 1186 … … 1172 1190 ELSE 1173 1191 mean_inflow_profiles(:,1) = hom_sum(:,1,0) ! u 1174 u_init(:) = hom_sum(:,1,0)1175 1192 mean_inflow_profiles(:,2) = hom_sum(:,2,0) ! v 1176 v_init(:) = hom_sum(:,2,0)1177 1193 ENDIF 1178 1194 mean_inflow_profiles(:,4) = hom_sum(:,4,0) ! pt 1179 pt_init(:) = hom_sum(:,4,0)1180 1195 IF ( humidity ) & 1181 1196 mean_inflow_profiles(:,6) = hom_sum(:,41,0) ! q 1182 q_init(:) = hom_sum(:,41,0)1183 1197 IF ( passive_scalar ) & 1184 1198 mean_inflow_profiles(:,7) = hom_sum(:,115,0) ! s 1185 s_init(:) = hom_sum(:,115,0) 1199 1186 1200 ! 1187 1201 !-- In case of complex terrain, determine vertical displacement at inflow
Note: See TracChangeset
for help on using the changeset viewer.