Changeset 4144 for palm/trunk/SOURCE/biometeorology_mod.f90
 Timestamp:
 Aug 6, 2019 9:11:47 AM (23 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/biometeorology_mod.f90
r4127 r4144 27 27 !  28 28 ! $Id$ 29 ! relational operators .EQ., .NE., etc. replaced by ==, /=, etc. 30 ! 31 ! 4127 20190730 14:47:10Z suehring 29 32 ! Output for bio_mrt added (merge from branch resler) 30 33 ! … … 4319 4322 ! 4320 4323 ! calculation of solar azimuth angle 4321 IF (woz .LE.12.0_wp) alpha = pi  acos( ( sin(thetasr) * sin( lat * dtor )  &4324 IF (woz <= 12.0_wp) alpha = pi  acos( ( sin(thetasr) * sin( lat * dtor )  & 4322 4325 sin( declination * dtor ) ) / ( cos(thetasr) * cos( lat * dtor ) ) ) 4323 IF (woz .GT. 12.0_wp)alpha = pi + acos( ( sin(thetasr) * sin( lat * dtor )  &4326 IF (woz > 12.0_wp) alpha = pi + acos( ( sin(thetasr) * sin( lat * dtor )  & 4324 4327 sin( declination * dtor ) ) / ( cos(thetasr) * cos( lat * dtor ) ) ) 4325 4328 saa = alpha / dtor … … 4347 4350 CALL uvem_solar_position 4348 4351 4349 IF (sza .GE.90) THEN4352 IF (sza >= 90) THEN 4350 4353 vitd3_exposure(:,:) = 0.0_wp 4351 4354 ELSE … … 4446 4449 IF (consider_obstructions ) THEN 4447 4450 obstruction_temp1 = building_obstruction_f%var_3d(:,j,i) 4448 IF ( obstruction_temp1(0) .NE. 9) THEN4451 IF ( obstruction_temp1(0) /= 9 ) THEN 4449 4452 DO pobi = 0, 44 4450 4453 DO bi = 0, 7 … … 4468 4471 4469 4472 obstruction_direct_beam = obstruction( nint(startpos_saa_float), nint( sza / 10.0_wp ) ) 4470 IF (sza .GE.89.99_wp) THEN4473 IF (sza >= 89.99_wp) THEN 4471 4474 sza = 89.99999_wp 4472 4475 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.