Changeset 3528 for palm/trunk
- Timestamp:
- Nov 15, 2018 7:07:11 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/radiation_model_mod.f90
r3524 r3528 28 28 ! ----------------- 29 29 ! $Id$ 30 ! Add an epsilon value to compare values in if statement to fix possible 31 ! precsion related errors in raytrace routines. 32 ! 33 ! 3524 2018-11-14 13:36:44Z raasch 30 34 ! missing cpp-directives added 31 35 ! … … 7084 7088 INTEGER(iwp) :: ip !< number of processor where gridbox reside 7085 7089 INTEGER(iwp) :: ig !< 1D index of gridbox in global 2D array 7086 REAL(wp) :: lad_s_target !< recieved lad_s of particular grid box 7090 7091 REAL(wp) :: eps = 1E-10_wp !< epsilon for value comparison 7092 REAL(wp) :: lad_s_target !< recieved lad_s of particular grid box 7087 7093 7088 7094 ! … … 7167 7173 ENDIF 7168 7174 7169 IF ( nextdist >= distance )EXIT7175 IF ( ABS(distance - nextdist) < eps ) EXIT 7170 7176 lastdist = nextdist 7171 7177 dimnext(seldim) = dimnext(seldim) + dimdelta(seldim) … … 7286 7292 INTEGER(iwp) :: nly !< maximum plant canopy height 7287 7293 INTEGER(iwp) :: ntrack 7288 REAL(wp) :: zbottom, ztop !< urban surface boundary in real numbers 7289 REAL(wp) :: zorig !< z coordinate of ray column entry 7290 REAL(wp) :: zexit !< z coordinate of ray column exit 7291 REAL(wp) :: qdist !< ratio of real distance to z coord difference 7292 REAL(wp) :: dxxyy !< square of real horizontal distance 7293 REAL(wp) :: curtrans !< transparency of current PC box crossing 7294 7294 7295 INTEGER(iwp) :: zb0 7295 7296 INTEGER(iwp) :: zb1 … … 7304 7305 INTEGER(MPI_ADDRESS_KIND) :: wdisp !< RMA window displacement 7305 7306 #endif 7307 7308 REAL(wp) :: eps = 1E-10_wp !< epsilon for value comparison 7309 REAL(wp) :: zbottom, ztop !< urban surface boundary in real numbers 7310 REAL(wp) :: zorig !< z coordinate of ray column entry 7311 REAL(wp) :: zexit !< z coordinate of ray column exit 7312 REAL(wp) :: qdist !< ratio of real distance to z coord difference 7313 REAL(wp) :: dxxyy !< square of real horizontal distance 7314 REAL(wp) :: curtrans !< transparency of current PC box crossing 7315 7316 7306 7317 7307 7318 yxorigin(:) = origin(2:3) … … 7422 7433 ENDIF 7423 7434 7424 IF ( nextdist >= distance) EXIT7435 IF ( ABS(distance - nextdist) < eps ) EXIT 7425 7436 7426 7437 IF ( rad_angular_discretization .AND. calc_svf ) THEN
Note: See TracChangeset
for help on using the changeset viewer.