Changeset 3001 for palm/trunk

Timestamp:

Apr 20, 2018 12:27:13 PM (7 years ago)

Author:

suehring

Message:

Bugfix, replace MERGE function by an if-condition in the anterpolation in order to avoid floating-point exceptions.

File:

• palm/trunk/SOURCE/pmc_interface_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix, replace MERGE function by an if-condition in the anterpolation (in 
+! order to avoid floating-point exceptions).
+! 
+! 2997 2018-04-19 13:35:17Z suehring
 ! Mask topography grid points in anterpolation
 ! 
@@ -5793 +5797 @@
 !
 !--             In case all child grid points are inside topography, i.e. 
-!--             cellsum is zero, also parent solution would have zero values at
-!--             that point, which may cause problems in particular for the 
-!--             temperature. Therefore, in case cellsum is zero, keep the 
-!--             parent solution at this point.  
-                fc(kk,jj,ii) = MERGE( ( 1.0_wp - fra ) * fc(kk,jj,ii) +        &
-                                      fra * cellsum  /                         &
-                                      REAL( ijkfc(kk,jj,ii), KIND = wp ),      &
-                                      fc(kk,jj,ii),                            &
-                                      cellsum /= 0.0_wp )       
+!--             ijkfc and cellsum are zero, also parent solution would have 
+!--             zero values at that grid point, which may cause problems in 
+!--             particular for the temperature. Therefore, in case cellsum is 
+!--             zero, keep the parent solution at this point. 
+                IF ( ijkfc(kk,jj,ii) /= 0 )  THEN
+                   fc(kk,jj,ii) = ( 1.0_wp - fra ) * fc(kk,jj,ii) +            &
+                                    fra * cellsum  /                           &
+                                    REAL( ijkfc(kk,jj,ii), KIND=wp )
+                ENDIF  
 
              ENDDO