(un)trustworthy tools -- a problem...
  Home FAQ Contact Sign in
comp.lang.fortran only
 
Advanced search
POPULAR GROUPS

more...

 Up
(un)trustworthy tools -- a problem...         

Group: comp.lang.fortran · Group Profile
Author: Carlie J. Coats
Date: Nov 21, 2007 08:39

On a server I'm using for my work, I'm using"

#52 which idb
/opt/intel/idbe/9.1.051/bin/idb
#53 which ifort
/opt/intel/fce/9.1.051/bin/ifort
#54 ifort -V
Intel(R) Fortran Compiler for Intel(R) EM64T-based applications,
Version 9.1 Build 20070510 Package ID: l_fc_c_9.1.051
Copyright (C) 1985-2007 Intel Corporation. All rights reserved.

with debug-and-check-everything compiler flags:

FFLAGS = -C -g -O0 -traceback -debug-parameters all -ftrapuv \
-check bounds -check uninit -check noarg_temp_created \
-ip -tpp7 -ftz

with relevant code segment below, where I'm going into the
"CALL ADVHS" line, so that at that point CON(:,:) and CWORK(:,C,:)
should match, for second subscript in the range [1,NROWS].

But they don't, as shown by the debugger-output extracts following...
This is very disconcerting ;-(

Any ideas?

-- Carlie J. Coats, Jr.
Baron Advanced Meteorological Systems, LLC.

--------------- code follows ---------------------------------------

...
REAL CON( NDYNSD,0:MXDIM )
REAL CWORK( NDYNSD, NCOLS, NROWS)
...
DO C = 1, NCOLS ! Y-direction

DO R = 1, NROWS+1
VEL(R) = VWIND( C,R,L ) + VWIND( C+1,R,L )
END DO

R = 0
MFAC = SBDY( NDYNSD,C,L )*METRIC( C,R,L )
CON(NDYNSD,R) = MFAC
DO V = 1, TR_DYN_HI
DFAC = CFAC * MFAC * DYNWT( V )
CON(V,R) = DFAC * SBDY( V,C,L )
END DO
DFAC = METRIC( C,R,L )
DO V = AE_DYN_LO, AE_DYN_HI
CON(V,C) = DFAC * SBDY( V,C,L )
END DO

DO R = 1, NROWS
DO V = 1, NDYNSD
CON(V,R) = CWORK( V,C,R )
END DO
END DO

R = NROWS+1
MFAC = NBDY( NDYNSD,C,L )*METRIC( C,R,L )
CON(NDYNSD,R) = MFAC
DO V = 1, TR_DYN_HI
DFAC = CFAC * MFAC * DYNWT( V )
CON(V,R) = DFAC * NBDY( V,C,L )
END DO
DFAC = METRIC( C,R,L )
DO V = AE_DYN_LO, AE_DYN_HI
CON(V,C) = DFAC * NBDY( V,C,L )
END DO

CALL ADVHS( NROWS, CON, VEL, DTS, DX2 )

DO R = 1, NROWS
DO V = 1, NDYNSD
CDYNS( V,C,R,L ) = CON(V,R)
END DO
END DO

END DO ! end loop on columns C

----------------------- debugger output ----------------------------

I'm having a failure at C=152, for which the debugger is telling me
that these values _don't_ match -- badly so!

(idb) print c
152
(idb) print ncols,nrows,ndynsd
163 160 79

(idb) print con(:,152)
(1,152) 1048.13696
(2,152) 11313.2256
(3,152) 745341.875
(4,152) 998.225708
(5,152) 73.2032166
(6,152) 2329.1936
(7,152) 3419.1167
(8,152) 8.98403156e-24
(9,152) 46.5838661
(10,152) 30.6122551
(11,152) 7187.2251
(12,152) 2.79503232e-23
(13,152) 7.05412817
(14,152) 5.98935461
(15,152) 226.264496
(16,152) 2.32919339e-23
(17,152) 2.12954832e-23
(18,152) 3.26087078e-23
(19,152) 4.52528985e-23
(20,152) 4.52528985e-23
(21,152) 1.53061277e-23
(22,152) 1.99645155e-23
(23,152) 56113.8203
(24,152) 306.122559
(25,152) 199.645126
(26,152) 739632.938
(27,152) 2.06299985e-23
(28,152) 3.59361262e-23
(29,152) 1.56388708e-23
(30,152) 2.62866108e-23
(31,152) 805.235352
(32,152) 5.65661232e-24
(33,152) 562.024536
(34,152) 1.26617152e-12
(35,152) 1.53522822e-12
(36,152) 4.99112863e-23
(37,152) 4.99112863e-23
(38,152) 4.99112863e-23
(39,152) 4.99112863e-23
(40,152) 11.2259274
(41,152) 5.88953202e-23
(42,152) 5.88953202e-23
(43,152) 80.9332809
(44,152) 6185.65771
(45,152) 177.746674
(46,152) 434.116211
(47,152) 13.8310785
(48,152) 9.77517084e-06
(49,152) 1.75146738e-07
(50,152) -809.561584
(51,152) -0.225888312
(52,152) -3008.32007
(53,152) -2.84895444
(54,152) -230555.188
(55,152) -64.2384109
(56,152) -1139.53516
(57,152) -1.02539742
(58,152) -4335.09521
(59,152) 561.864136
(60,152) 5053.53955
(61,152) -174522.188
(62,152) -2568.03491
(63,152) 11.4412622
(64,152) -0.331968963
(65,152) 1889.46985
(66,152) 2640.65796
(67,152) 369.057983
(68,152) -9385.83105
(69,152) -3.0564332
(70,152) -137380
(71,152) -36.5769768
(72,152) -847.055054
(73,152) -0.239837781
(74,152) 2.71012882e+12
(75,152) -1.35481e+12
(76,152) 16244672
(77,152) 0.0112076588
(78,152) -1.45581281
(79,152) 9637.13672
(idb) print cwork(:,c,152)
(1,152,152) 1048.13696
(2,152,152) 11313.2256
(3,152,152) 745341.875
(4,152,152) 998.225708
(5,152,152) 73.2032166
(6,152,152) 2329.1936
(7,152,152) 3419.1167
(8,152,152) 8.98403156e-24
(9,152,152) 46.5838661
(10,152,152) 30.6122551
(11,152,152) 7187.2251
(12,152,152) 2.79503232e-23
(13,152,152) 7.05412817
(14,152,152) 5.98935461
(15,152,152) 226.264496
(16,152,152) 2.32919339e-23
(17,152,152) 2.12954832e-23
(18,152,152) 3.26087078e-23
(19,152,152) 4.52528985e-23
(20,152,152) 4.52528985e-23
(21,152,152) 1.53061277e-23
(22,152,152) 1.99645155e-23
(23,152,152) 56113.8203
(24,152,152) 306.122559
(25,152,152) 199.645126
(26,152,152) 739632.938
(27,152,152) 2.06299985e-23
(28,152,152) 3.59361262e-23
(29,152,152) 1.56388708e-23
(30,152,152) 2.62866108e-23
(31,152,152) 805.235352
(32,152,152) 5.65661232e-24
(33,152,152) 562.024536
(34,152,152) 1.26617152e-12
(35,152,152) 1.53522822e-12
(36,152,152) 4.99112863e-23
(37,152,152) 4.99112863e-23
(38,152,152) 4.99112863e-23
(39,152,152) 4.99112863e-23
(40,152,152) 11.2259274
(41,152,152) 5.88953202e-23
(42,152,152) 5.88953202e-23
(43,152,152) 80.9332809
(44,152,152) 6030.29004
(45,152,152) 180.078339
(46,152,152) 9.15134348e-27
(47,152,152) 9.15134348e-27
(48,152,152) 0
(49,152,152) 5.12089871e-10
(50,152,152) 9.00042761e-27
(51,152,152) 9.00042761e-27
(52,152,152) 9.00042761e-27
(53,152,152) 9.00042761e-27
(54,152,152) 9.00042761e-27
(55,152,152) 9.00042761e-27
(56,152,152) 9.00042761e-27
(57,152,152) 9.00042761e-27
(58,152,152) 9.00042761e-27
(59,152,152) 285.963348
(60,152,152) 2572.31738
(61,152,152) 9.00042761e-27
(62,152,152) 9.00042761e-27
(63,152,152) 92.8395844
(64,152,152) 3.08911074e-09
(65,152,152) 5234.54883
(66,152,152) 7315.61572
(67,152,152) 1023.21887
(68,152,152) 0.383342326
(69,152,152) 9.00042761e-27
(70,152,152) 9.00042761e-27
(71,152,152) 9.00042761e-27
(72,152,152) 9.00042761e-27
(73,152,152) 9.00042761e-27
(74,152,152) 2.1532423e+12
(75,152,152) 9.00042761e-27
(76,152,152) 432170080
(77,152,152) 0.00680456636
(78,152,152) 9.00042761e-27
(79,152,152) 9637.13672
4 Comments
diggit! del.icio.us! reddit!