Changeset 2071
- Timestamp:
- Nov 17, 2016 11:22:14 AM (8 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/ls_forcing_mod.f90
r2038 r2071 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! 22 ! Small bugfix (Resler) 23 23 ! 24 24 ! Former revisions: … … 217 217 ENDDO 218 218 219 shf_surf = shf_surf * heatflux_input_conversion(nzb)220 qsws_surf = qsws_surf * waterflux_input_conversion(nzb)221 222 219 IF ( time_surf(1) > end_time ) THEN 223 220 WRITE ( message_string, * ) 'No time dependent surface variables in ',& … … 226 223 CALL message( 'ls_forcing', 'PA0371', 0, 0, 0, 6, 0 ) 227 224 lsf_surf = .FALSE. 225 ELSE 226 shf_surf = shf_surf * heatflux_input_conversion(nzb) 227 qsws_surf = qsws_surf * waterflux_input_conversion(nzb) 228 228 ENDIF 229 229 -
palm/trunk/SOURCE/urban_surface_mod.f90
r2032 r2071 21 21 ! Current revisions: 22 22 ! ------------------ 23 ! 23 ! Small bugfix (Resler) 24 24 ! 25 25 ! Former revisions: … … 422 422 INTEGER(iwp), PARAMETER :: icsurf = 5 !< Surface skin layer heat capacity (J mâ2 Kâ1 ) 423 423 INTEGER(iwp), PARAMETER :: ithick = 6 !< thickness of the surface (wall, roof, land) ( m ) 424 INTEGER(iwp), PARAMETER :: irhoC = 7 !< volumetric heat capacity rho _ocean*C of the material ( J mâ3 Kâ1 )424 INTEGER(iwp), PARAMETER :: irhoC = 7 !< volumetric heat capacity rho*C of the material ( J mâ3 Kâ1 ) 425 425 INTEGER(iwp), PARAMETER :: ilambdah = 8 !< thermal conductivity λH of the wall (W mâ1 Kâ1 ) 426 426 CHARACTER(12), DIMENSION(:), ALLOCATABLE :: surface_type_names !< names of wall types (used only for reports) … … 2889 2889 IF ( plant_canopy ) THEN 2890 2890 pchf_prep(:) = r_d * (hyp(nzub:nzut) / 100000.0_wp)**0.286_wp & 2891 / (cp * hyp(nzub:nzut) * dx*dy*dz) !< equals to 1 / (rho _ocean* c_p * Vbox * T)2891 / (cp * hyp(nzub:nzut) * dx*dy*dz) !< equals to 1 / (rho * c_p * Vbox * T) 2892 2892 ENDIF 2893 2893 … … 3035 3035 IF ( surf(id, isurfsrc) >= isky ) THEN 3036 3036 !-- Diffuse rad from boundary surfaces. See comments for svf above. 3037 pcbinsw(ipcgb) = pcbinsw(ipcgb) + svf(1,isvf) * svf(2,isvf) * rad_sw_in_diff(j,i)3037 pcbinsw(ipcgb) = pcbinsw(ipcgb) + csf(1,icsf) * csf(2,icsf) * rad_sw_in_diff(j,i) 3038 3038 !-- canopy shading is applied only to shortwave, therefore no absorbtion for lw 3039 3039 !-- pcbinlw(ipcgb) = pcbinlw(ipcgb) + svf(1,isvf) * rad_lw_in_diff(j,i) … … 3047 3047 pc_abs_frac = 1._wp - exp(pc_abs_eff * lad_s(k,j,i)) 3048 3048 pcbinsw(ipcgb) = pcbinsw(ipcgb) & 3049 + rad_sw_in_dir(j, i) * pc_box_area * svf(2,isvf) * pc_abs_frac3049 + rad_sw_in_dir(j, i) * pc_box_area * csf(2,icsf) * pc_abs_frac 3050 3050 ENDIF 3051 3051 ENDDO … … 3496 3496 ENDIF 3497 3497 #if defined( __parallel ) 3498 CALL MPI_B ARRIER( comm2d, ierr )3498 CALL MPI_Barrier( comm2d, ierr ) 3499 3499 #endif 3500 3500 ENDDO … … 3866 3866 ENDIF 3867 3867 3868 !-- volumetric heat capacity rho _ocean*C of the wall ( J mâ3 Kâ1 )3868 !-- volumetric heat capacity rho*C of the wall ( J mâ3 Kâ1 ) 3869 3869 rho_c_wall(:,l) = surface_params(irhoC, ip) 3870 3870 … … 3933 3933 pt1 = pt(k,j,i) 3934 3934 3935 !-- calculate rho _ocean* cp coefficient at surface layer3935 !-- calculate rho * cp coefficient at surface layer 3936 3936 rho_cp = cp * hyp(k) / ( r_d * pt1 * exn(k) ) 3937 3937 … … 4246 4246 ENDIF 4247 4247 #if defined( __parallel ) 4248 CALL MPI_B ARRIER( comm2d, ierr )4248 CALL MPI_Barrier( comm2d, ierr ) 4249 4249 #endif 4250 4250 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.