Changeset 3716 for palm/trunk/UTIL/inifor/src/inifor_transform.f90
- Timestamp:
- Feb 5, 2019 5:02:38 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/UTIL/inifor/src/inifor_transform.f90
r3680 r3716 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Include out-of-bounds error message in log 29 ! 30 ! 31 ! 3680 2019-01-18 14:54:12Z knoop 28 32 ! Check if set of averaging columns is empty 29 33 ! … … 777 781 latpos = (palm_clat(i,j) - lat0) * cosmo_dyi 778 782 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) 786 795 END IF 787 796
Note: See TracChangeset
for help on using the changeset viewer.