Ignore:
Timestamp:
Feb 5, 2019 5:02:38 PM (5 years ago)
Author:
eckhard
Message:

inifor: bugfix: removed dependency on soilmoisture input files; added netcdf preprocessor flag

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/UTIL/inifor/src/inifor_transform.f90

    r3680 r3716  
    2626! -----------------
    2727! $Id$
     28! Include out-of-bounds error message in log
     29!
     30!
     31! 3680 2019-01-18 14:54:12Z knoop
    2832! Check if set of averaging columns is empty
    2933!
     
    777781          latpos = (palm_clat(i,j) - lat0) * cosmo_dyi
    778782
    779           IF (lonpos < 0.0 .OR. latpos < 0.0)  THEN
    780              PRINT *, " Error while finding neighbours: lonpos or latpos out of bounds!"
    781              PRINT *, "     (i,j) = (", i, ",",j,")"
    782              PRINT *, "      lonpos ", lonpos*TO_DEGREES, ", latpos ", latpos*TO_DEGREES
    783              PRINT *, "        lon0 ", lon0  *TO_DEGREES,   ", lat0   ", lat0*TO_DEGREES
    784              PRINT *, "    PALM lon ", palm_clon(i,j)*TO_DEGREES,   ", PALM lat ",palm_clat(i,j)*TO_DEGREES
    785              STOP
     783          IF (lonpos < 0.0_dp .OR. latpos < 0.0_dp)  THEN
     784             message = "lonpos or latpos out of bounds " //                    &
     785                "while finding interpolation neighbours!" // NEW_LINE(' ') //  &
     786                "          (i,j) = (" //                                       &
     787                TRIM(str(i)) // ", " // TRIM(str(j)) // ")" // NEW_LINE(' ') //&
     788                "          lonpos " // TRIM(real_to_str(lonpos*TO_DEGREES)) // &
     789                ", latpos " // TRIM(real_to_str(latpos*TO_DEGREES)) // NEW_LINE(' ') // &
     790                "          lon0 " // TRIM(real_to_str(lon0  *TO_DEGREES)) //   &
     791                ", lat0   " // TRIM(real_to_str(lat0*TO_DEGREES)) // NEW_LINE(' ') // &
     792                "          PALM lon " // TRIM(real_to_str(palm_clon(i,j)*TO_DEGREES)) // &
     793                ", PALM lat " // TRIM(real_to_str(palm_clat(i,j)*TO_DEGREES))
     794             CALL inifor_abort('find_horizontal_neighbours', message)
    786795          END IF
    787796
Note: See TracChangeset for help on using the changeset viewer.