PageRenderTime 70ms CodeModel.GetById 32ms RepoModel.GetById 0ms app.codeStats 1ms

/other/netcdf_write_matrix/src/nf_test/test_read.F

https://github.com/jbeezley/wrf-fire
FORTRAN Legacy | 1061 lines | 836 code | 107 blank | 118 comment | 224 complexity | 36e8bf933c8c38b4d31123f042cf73bb MD5 | raw file
Possible License(s): AGPL-1.0 
    

  1. C*********************************************************************
    2. 

  2. C   Copyright 1996, UCAR/Unidata
    3. 

  3. C   See netcdf/COPYRIGHT file for copying and redistribution conditions.
    4. 

  4. C   $Id: test_read.F,v 1.11 2005/11/07 15:13:32 ed Exp $
    5. 

  5. C*********************************************************************
    6. 

  6. 

  7. C Test nf_strerror.
    8. 

  8. C    Try on a bad error status.
    9. 

  9. C    Test for each defined error status.
    10. 

  10. C
    11. 

  11.         subroutine test_nf_strerror()
    12. 

  12.         implicit        none
    13. 

  13. #include "tests.inc"
    14. 

  14.         integer         number_of_messages
    15. 

  15.         parameter       (number_of_messages = 27)
    16. 

  16. 

  17.         integer         i
    18. 

  18.         integer         status(number_of_messages)
    19. 

  19.         character*80    message
    20. 

  20.         character*80    msg(number_of_messages)
    21. 

  21. 

  22.         data    status(1)  / NF_NOERR/
    23. 

  23.         data    status(2)  / NF_EBADID /
    24. 

  24.         data    status(3)  / NF_EEXIST /
    25. 

  25.         data    status(4)  / NF_EINVAL /
    26. 

  26.         data    status(5)  / NF_EPERM /
    27. 

  27.         data    status(6)  / NF_ENOTINDEFINE /
    28. 

  28.         data    status(7)  / NF_EINDEFINE /
    29. 

  29.         data    status(8)  / NF_EINVALCOORDS /
    30. 

  30.         data    status(9)  / NF_EMAXDIMS /
    31. 

  31.         data    status(10) / NF_ENAMEINUSE /
    32. 

  32.         data    status(11) / NF_ENOTATT /
    33. 

  33.         data    status(12) / NF_EMAXATTS /
    34. 

  34.         data    status(13) / NF_EBADTYPE /
    35. 

  35.         data    status(14) / NF_EBADDIM /
    36. 

  36.         data    status(15) / NF_EUNLIMPOS /
    37. 

  37.         data    status(16) / NF_EMAXVARS /
    38. 

  38.         data    status(17) / NF_ENOTVAR /
    39. 

  39.         data    status(18) / NF_EGLOBAL /
    40. 

  40.         data    status(19) / NF_ENOTNC /
    41. 

  41.         data    status(20) / NF_ESTS /
    42. 

  42.         data    status(21) / NF_EMAXNAME /
    43. 

  43.         data    status(22) / NF_EUNLIMIT /
    44. 

  44.         data    status(23) / NF_ENORECVARS /
    45. 

  45.         data    status(24) / NF_ECHAR /
    46. 

  46.         data    status(25) / NF_EEDGE /
    47. 

  47.         data    status(26) / NF_ESTRIDE /
    48. 

  48.         data    status(27) / NF_EBADNAME /
    49. 

  49. 

  50.         data    msg(1)  / 'No error' /
    51. 

  51.         data    msg(2)  / 'Not a netCDF id' /
    52. 

  52.         data    msg(3)  / 'netCDF file exists && NC_NOCLOBBER' /
    53. 

  53.         data    msg(4)  / 'Invalid argument' /
    54. 

  54.         data    msg(5)  / 'Write to read only' /
    55. 

  55.         data    msg(6)  / 'Operation not allowed in data mode' /
    56. 

  56.         data    msg(7)  / 'Operation not allowed in define mode' /
    57. 

  57.         data    msg(8)  / 'Index exceeds dimension bound' /
    58. 

  58.         data    msg(9)  / 'NC_MAX_DIMS exceeded' /
    59. 

  59.         data    msg(10) / 'String match to name in use' /
    60. 

  60.         data    msg(11) / 'Attribute not found' /
    61. 

  61.         data    msg(12) / 'NC_MAX_ATTRS exceeded' /
    62. 

  62.         data    msg(13)
    63. 

  63.      +      / 'Not a netCDF data type or _FillValue type mismatch' /
    64. 

  64.         data    msg(14) / 'Invalid dimension id or name' /
    65. 

  65.         data    msg(15) / 'NC_UNLIMITED in the wrong index' /
    66. 

  66.         data    msg(16) / 'NC_MAX_VARS exceeded' /
    67. 

  67.         data    msg(17) / 'Variable not found' /
    68. 

  68.         data    msg(18) / 'Action prohibited on NC_GLOBAL varid' /
    69. 

  69.         data    msg(19) / 'Not a netCDF file' /
    70. 

  70.         data    msg(20) / 'In Fortran, string too short' /
    71. 

  71.         data    msg(21) / 'NC_MAX_NAME exceeded' /
    72. 

  72.         data    msg(22) / 'NC_UNLIMITED size already in use' /
    73. 

  73.         data    msg(23) / 'nc_rec op when there are no record vars' /
    74. 

  74.         data    msg(24) / 'Attempt to convert between text & numbers' /
    75. 

  75.         data    msg(25) / 'Start+count exceeds dimension bound' /
    76. 

  76.         data    msg(26) / 'Illegal stride' /
    77. 

  77.         data    msg(27) 
    78. 

  78.      +      / 'Attribute or variable name contains illegal characters' /
    79. 

  79. 

  80. C       /* Try on a bad error status */
    81. 

  81.         message = nf_strerror(-666)!/* should fail */
    82. 

  82.         if (message .ne. 'Unknown Error')
    83. 

  83.      +      call errorc('nf_strerror on bad error status returned: ',
    84. 

  84.      +          message)
    85. 

  85. 

  86. C       /* Try on each legitimate error status */
    87. 

  87.         do 1, i=1, number_of_messages
    88. 

  88.             message = nf_strerror(status(i))
    89. 

  89.             if (message .ne. msg(i))
    90. 

  90.      +          call error('nf_strerror() should return "' // msg(i) //
    91. 

  91.      +                     '"')
    92. 

  92. 1       continue
    93. 

  93.         end
    94. 

  94. 

  95. 

  96. C Test nf_open.
    97. 

  97. C If in read-only section of tests,
    98. 

  98. C    Try to open a non-existent netCDF file, check error return.
    99. 

  99. C    Open a file that is not a netCDF file, check error return.
    100. 

  100. C    Open a netCDF file with a bad mode argument, check error return.
    101. 

  101. C    Open a netCDF file with NF_NOWRITE mode, try to write, check error.
    102. 

  102. C    Try to open a netcdf twice, check whether returned netcdf ids different.
    103. 

  103. C If in writable section of tests,
    104. 

  104. C    Open a netCDF file with NF_WRITE mode, write something, close it.
    105. 

  105. C On exit, any open netCDF files are closed.
    106. 

  106.         subroutine test_nf_open()
    107. 

  107.         implicit        none
    108. 

  108. #include "tests.inc"
    109. 

  109.         integer err
    110. 

  110.         integer ncid
    111. 

  111.         integer ncid2
    112. 

  112.         
    113. 

  113. C       /* Try to open a nonexistent file */
    114. 

  114.         err = nf_open('tooth-fairy.nc', NF_NOWRITE, ncid)!/* should fail */
    115. 

  115.         if (err .eq. NF_NOERR)
    116. 

  116.      +      call error('nf_open of nonexistent file should have failed')
    117. 

  117.         if (.not. NF_ISSYSERR(err))
    118. 

  118.      +      call error(
    119. 

  119.      +  'nf_open of nonexistent file should have returned system error')
    120. 

  120. 

  121. C       /* Open a file that is not a netCDF file. */
    122. 

  122.         err = nf_open('test_read.o', NF_NOWRITE, ncid)!/* should fail */
    123. 

  123.         if (err .ne. NF_ENOTNC)
    124. 

  124.      +      call errore('nf_open of non-netCDF file: ', err)
    125. 

  125. 

  126. C       /* Open a netCDF file in read-only mode, check that write fails */
    127. 

  127.         err = nf_open(testfile, NF_NOWRITE, ncid)
    128. 

  128.         if (err .ne. 0)
    129. 

  129.      +      call errore('nf_open: ', err)
    130. 

  130.         err = nf_redef(ncid)    !/* should fail */
    131. 

  131.         if (err .ne. NF_EPERM)
    132. 

  132.      +      call error('nf_redef of read-only file should fail')
    133. 

  133. C       /* Opened OK, see if can open again and get a different netCDF ID */
    134. 

  134.         err = nf_open(testfile, NF_NOWRITE, ncid2)
    135. 

  135.         if (err .ne. 0) then
    136. 

  136.             call errore('nf_open: ', err)
    137. 

  137.         else
    138. 

  138.             err = nf_close(ncid2)
    139. 

  139.         end if
    140. 

  140.         if (ncid2 .eq. ncid)
    141. 

  141.      +      call error(
    142. 

  142.      +  'netCDF IDs for first and second nf_open calls should differ')
    143. 

  143. 

  144.         if (.not. readonly) then        !/* tests using netCDF scratch file */
    145. 

  145.             err = nf_create(scratch, NF_NOCLOBBER, ncid2)
    146. 

  146.             if (err .ne. 0) then
    147. 

  147.                 call errore('nf_create: ', err)
    148. 

  148.             else 
    149. 

  149.                 err = nf_close(ncid2)
    150. 

  150.             end if
    151. 

  151.             err = nf_open(scratch, NF_WRITE, ncid2)
    152. 

  152.             if (err .ne. 0) then
    153. 

  153.                 call errore('nf_open: ', err)
    154. 

  154.             else 
    155. 

  155.                 err = nf_close(ncid2)
    156. 

  156.             end if
    157. 

  157.             err = nf_delete(scratch)
    158. 

  158.             if (err .ne. 0) 
    159. 

  159.      +          call errorc('delete of scratch file failed', scratch)
    160. 

  160.         end if
    161. 

  161. 

  162.         err = nf_close(ncid)
    163. 

  163.         if (err .ne. 0)
    164. 

  164.      +      call errore('nf_close: ', err)
    165. 

  165.         end
    166. 

  166. 

  167. 

  168. C 
    169. 

  169. C Test nf_close.
    170. 

  170. C    Try to close a netCDF file twice, check whether second close fails.
    171. 

  171. C    Try on bad handle, check error return.
    172. 

  172. C    Try in define mode and data mode.
    173. 

  173. C/
    174. 

  174.         subroutine test_nf_close()
    175. 

  175.         implicit        none
    176. 

  176. #include "tests.inc"
    177. 

  177.         integer ncid
    178. 

  178.         integer err
    179. 

  179. 

  180.         err = nf_open(testfile, NF_NOWRITE, ncid)
    181. 

  181.         if (err .ne. 0)
    182. 

  182.      +      call errore('nf_open: ', err)
    183. 

  183. 

  184. C       /* Close a netCDF file twice, second time should fail */
    185. 

  185.         err = nf_close(ncid)
    186. 

  186.         if (err .ne. 0)
    187. 

  187.      +      call errore('nf_close failed: ', err)
    188. 

  188.         err = nf_close(ncid)
    189. 

  189.         if (err .ne. NF_EBADID)
    190. 

  190.      +      call error('nf_close of closed file should have failed')
    191. 

  191.         
    192. 

  192. C       /* Try with a bad netCDF ID */
    193. 

  193.         err = nf_close(BAD_ID)!/* should fail */
    194. 

  194.         if (err .ne. NF_EBADID)
    195. 

  195.      +      call errore(
    196. 

  196.      +          'nf_close with bad netCDF ID returned wrong error: ', 
    197. 

  197.      +          err)
    198. 

  198. 

  199. C       /* Close in data mode */
    200. 

  200.         err = nf_open(testfile, NF_NOWRITE, ncid)
    201. 

  201.         if (err .ne. 0)
    202. 

  202.      +      call errore('nf_open: ', err)
    203. 

  203.         err = nf_close(ncid)
    204. 

  204.         if (err .ne. 0)
    205. 

  205.      +      call errore('nf_close in data mode failed: ', err)
    206. 

  206. 

  207.         if (.not. readonly) then        !/* tests using netCDF scratch file */
    208. 

  208.             err = nf_create(scratch, NF_NOCLOBBER, ncid)
    209. 

  209.             if (err .ne. 0) 
    210. 

  210.      +          call errore('nf_create: ', err)
    211. 

  211.             err = nf_close(ncid)
    212. 

  212.             if (err .ne. 0)
    213. 

  213.      +          call errore('nf_close in define mode: ', err)
    214. 

  214.             err = nf_delete(scratch)
    215. 

  215.             if (err .ne. 0)
    216. 

  216.      +          call errorc('delete of scratch file failed: ', 
    217. 

  217.      +              scratch)
    218. 

  218.         end if
    219. 

  219.         end
    220. 

  220. 

  221. 

  222. C Test nf_inq.
    223. 

  223. C    Try on bad handle, check error return.
    224. 

  224. C    Try in data mode, check returned values.
    225. 

  225. C    Try asking for subsets of info.
    226. 

  226. C If in writable section of tests,
    227. 

  227. C    Try in define mode, after adding an unlimited dimension, variable.
    228. 

  228. C On exit, any open netCDF files are closed.
    229. 

  229.         subroutine test_nf_inq()
    230. 

  230.         implicit        none
    231. 

  231. #include "tests.inc"
    232. 

  232.         integer ncid
    233. 

  233.         integer ncid2                   !/* for scratch netCDF dataset */
    234. 

  234.         integer ndims                   !/* number of dimensions */
    235. 

  235.         integer nvars                   !/* number of variables */
    236. 

  236.         integer ngatts                  !/* number of global attributes */
    237. 

  237.         integer recdim                  !/* id of unlimited dimension */
    238. 

  238.         integer err
    239. 

  239.         integer ndims0
    240. 

  240.         integer nvars0
    241. 

  241.         integer ngatts0
    242. 

  242.         integer recdim0
    243. 

  243.         integer did
    244. 

  244.         integer vid
    245. 

  245. 

  246.         err = nf_open(testfile, NF_NOWRITE, ncid)
    247. 

  247.         if (err .ne. 0)
    248. 

  248.      +      call errore('nf_open: ', err)
    249. 

  249.         
    250. 

  250. C       /* Try on bad handle */
    251. 

  251.         err = nf_inq(BAD_ID, ndims, nvars, ngatts, recdim)
    252. 

  252.         if (err .ne. NF_EBADID)
    253. 

  253.      +      call errore('bad ncid: ', err)
    254. 

  254.         
    255. 

  255.         err = nf_inq(ncid, ndims, nvars, ngatts, recdim)
    256. 

  256.         if (err .ne. 0) then
    257. 

  257.             call errore('nf_inq: ', err)
    258. 

  258.         else if (ndims .ne. NDIMS) then
    259. 

  259.             call errori('nf_inq: wrong number of dimensions returned: ',
    260. 

  260.      +                  ndims)
    261. 

  261.         else if (nvars .ne. NVARS) then
    262. 

  262.             call errori('nf_inq: wrong number of variables returned: ',
    263. 

  263.      +                  nvars)
    264. 

  264.         else if (ngatts .ne. NGATTS) then
    265. 

  265.             call errori(
    266. 

  266.      +          'nf_inq: wrong number of global atts returned: ',
    267. 

  267.      +          ngatts)
    268. 

  268.         else if (recdim .ne. RECDIM) then
    269. 

  269.             call errori('nf_inq: wrong record dimension ID returned: ',
    270. 

  270.      +                  recdim)
    271. 

  271.         end if
    272. 

  272. 

  273.         if (.not. readonly) then        !/* tests using netCDF scratch file */
    274. 

  274.             err = nf_create(scratch, NF_NOCLOBBER, ncid2)
    275. 

  275.             if (err .ne. 0) then
    276. 

  276.                 call errore('nf_create: ', err)
    277. 

  277.             else                !/* add dim, var, gatt, check inq */
    278. 

  278.                 err = nf_enddef(ncid2) !/* enter data mode */
    279. 

  279.                 err = nf_inq(ncid2, ndims0, nvars0, ngatts0, recdim0)
    280. 

  280.                 if (err .ne. 0)
    281. 

  281.      +              call errore('nf_inq: ', err)
    282. 

  282.                 err = nf_redef(ncid2) !/* enter define mode */
    283. 

  283. C               /* Check that inquire still works in define mode */
    284. 

  284.                 err = nf_inq(ncid2, ndims, nvars, ngatts, recdim)
    285. 

  285.                 if (err .ne. 0) then
    286. 

  286.                     call errore('nf_inq in define mode: ', err)
    287. 

  287.                 else if (ndims .ne. ndims0) then
    288. 

  288.                     call errori('nf_inq in define mode: ndims wrong, ',
    289. 

  289.      +                          ndims)
    290. 

  290.                 else if (nvars .ne. nvars0) then
    291. 

  291.                     call errori('nf_inq in define mode: nvars wrong, ',
    292. 

  292.      +                          nvars)
    293. 

  293.                 else if (ngatts .ne. ngatts0) then
    294. 

  294.                     call errori(
    295. 

  295.      +                  'nf_inq in define mode: ngatts wrong, ', ngatts)
    296. 

  296.                 else if (recdim .ne. recdim0) then
    297. 

  297.                     call errori('nf_inq in define mode: recdim wrong, ',
    298. 

  298.      +                          recdim)
    299. 

  299.                 end if
    300. 

  300. 

  301. C               /* Add dim, var, global att */
    302. 

  302.                 err = nf_def_dim(ncid2, 'inqd', 1, did)
    303. 

  303.                 if (err .ne. 0)
    304. 

  304.      +              call errore('nf_def_dim: ', err)
    305. 

  305.                 err = nf_def_var(ncid2, 'inqv', NF_FLOAT, 0, 0, vid)
    306. 

  306.                 if (err .ne. 0)
    307. 

  307.      +              call errore('nf_def_var: ', err)
    308. 

  308. 

  309.                 err = nf_put_att_text(ncid2, NF_GLOBAL, 'inqa', 
    310. 

  310.      +                                len('stuff'), 'stuff')
    311. 

  311.                 if (err .ne. 0)
    312. 

  312.      +              call errore('nf_put_att_text: ', err)
    313. 

  313. 

  314. C               /* Make sure nf_inq sees the additions while in define mode */
    315. 

  315.                 err = nf_inq(ncid2, ndims, nvars, ngatts, recdim)
    316. 

  316.                 if (err .ne. 0) then
    317. 

  317.                     call errore('nf_inq in define mode: ', err)
    318. 

  318.                 else if (ndims .ne. ndims0 + 1) then
    319. 

  319.                     call errori('nf_inq in define mode: ndims wrong, ',
    320. 

  320.      +                          ndims)
    321. 

  321.                 else if (nvars .ne. nvars0 + 1) then
    322. 

  322.                     call errori('nf_inq in define mode: nvars wrong, ',
    323. 

  323.      +                          nvars)
    324. 

  324.                 else if (ngatts .ne. ngatts0 + 1) then
    325. 

  325.                     call errori('nf_inq in define mode: ngatts wrong, ',
    326. 

  326.      +                          ngatts)
    327. 

  327.                 end if
    328. 

  328.                 err = nf_enddef(ncid2)
    329. 

  329.                 if (err .ne. 0)
    330. 

  330.      +              call errore('nf_enddef: ', err)
    331. 

  331. 

  332. C               /* Make sure nf_inq stills sees additions in data mode */
    333. 

  333.                 err = nf_inq(ncid2, ndims, nvars, ngatts, recdim)
    334. 

  334.                 if (err .ne. 0) then
    335. 

  335.                     call errore('nf_inq failed in data mode: ',err)
    336. 

  336.                 else if (ndims .ne. ndims0 + 1) then
    337. 

  337.                     call errori('nf_inq in define mode: ndims wrong, ',
    338. 

  338.      +                          ndims)
    339. 

  339.                 else if (nvars .ne. nvars0 + 1) then
    340. 

  340.                     call errori('nf_inq in define mode: nvars wrong, ',
    341. 

  341.      +                          nvars)
    342. 

  342.                 else if (ngatts .ne. ngatts0 + 1) then
    343. 

  343.                     call errori('nf_inq in define mode: ngatts wrong, ',
    344. 

  344.      +                          ngatts)
    345. 

  345.                 end if
    346. 

  346.                 err = nf_close(ncid2)
    347. 

  347.                 err = nf_delete(scratch)
    348. 

  348.                 if (err .ne. 0)
    349. 

  349.      +              call errorc('delete of scratch file failed:', 
    350. 

  350.      +                  scratch)
    351. 

  351.             end if
    352. 

  352.         end if
    353. 

  353. 

  354.         err = nf_close(ncid)
    355. 

  355.         if (err .ne. 0)
    356. 

  356.      +      call errore('nf_close: ', err)
    357. 

  357.         end
    358. 

  358. 

  359. 

  360.         subroutine test_nf_inq_natts()
    361. 

  361.         implicit        none
    362. 

  362. #include "tests.inc"
    363. 

  363.         integer ncid
    364. 

  364.         integer ngatts                  !/* number of global attributes */
    365. 

  365.         integer err
    366. 

  366. 

  367.         err = nf_inq_natts(BAD_ID, ngatts)
    368. 

  368.         if (err .ne. NF_EBADID)
    369. 

  369.      +      call errore('bad ncid: ', err)
    370. 

  370.         err = nf_open(testfile, NF_NOWRITE, ncid)
    371. 

  371.         if (err .ne. 0)
    372. 

  372.      +      call errore('nf_open: ', err)
    373. 

  373.         err = nf_inq_natts(ncid, ngatts)
    374. 

  374.         if (err .ne. 0) then
    375. 

  375.             call errore('nf_inq_natts: ', err)
    376. 

  376.         else if (ngatts .ne. NGATTS) then
    377. 

  377.             call errori(
    378. 

  378.      +          'nf_inq_natts: wrong number of global atts returned, ',
    379. 

  379.      +          ngatts)
    380. 

  380.         end if
    381. 

  381.         err = nf_close(ncid)
    382. 

  382.         if (err .ne. 0)
    383. 

  383.      +      call errore('nf_close: ', err)
    384. 

  384.         end
    385. 

  385. 

  386. 

  387.         subroutine test_nf_inq_ndims()
    388. 

  388.         implicit        none
    389. 

  389. #include "tests.inc"
    390. 

  390.         integer ncid
    391. 

  391.         integer ndims
    392. 

  392.         integer err
    393. 

  393. 

  394.         err = nf_inq_ndims(BAD_ID, ndims)
    395. 

  395.         if (err .ne. NF_EBADID)
    396. 

  396.      +      call errore('bad ncid: ', err)
    397. 

  397.         err = nf_open(testfile, NF_NOWRITE, ncid)
    398. 

  398.         if (err .ne. 0)
    399. 

  399.      +      call errore('nf_open: ', err)
    400. 

  400.         err = nf_inq_ndims(ncid, ndims)
    401. 

  401.         if (err .ne. 0) then
    402. 

  402.             call errore('nf_inq_ndims: ', err)
    403. 

  403.         else if (ndims .ne. NDIMS) then
    404. 

  404.             call errori('nf_inq_ndims: wrong number returned, ', ndims)
    405. 

  405.         end if
    406. 

  406.         err = nf_close(ncid)
    407. 

  407.         if (err .ne. 0)
    408. 

  408.      +      call errore('nf_close: ', err)
    409. 

  409.         end
    410. 

  410. 

  411. 

  412.         subroutine test_nf_inq_nvars()
    413. 

  413.         implicit        none
    414. 

  414. #include "tests.inc"
    415. 

  415.         integer ncid
    416. 

  416.         integer nvars
    417. 

  417.         integer err
    418. 

  418. 

  419.         err = nf_inq_nvars(BAD_ID, nvars)
    420. 

  420.         if (err .ne. NF_EBADID)
    421. 

  421.      +      call errore('bad ncid: ', err)
    422. 

  422.         err = nf_open(testfile, NF_NOWRITE, ncid)
    423. 

  423.         if (err .ne. 0)
    424. 

  424.      +      call errore('nf_open: ', err)
    425. 

  425.         err = nf_inq_nvars(ncid, nvars)
    426. 

  426.         if (err .ne. 0) then
    427. 

  427.             call errore('nf_inq_nvars: ', err)
    428. 

  428.         else if (nvars .ne. NVARS) then
    429. 

  429.             call errori('nf_inq_nvars: wrong number returned, ', nvars)
    430. 

  430.         end if
    431. 

  431.         err = nf_close(ncid)
    432. 

  432.         if (err .ne. 0)
    433. 

  433.      +      call errore('nf_close: ', err)
    434. 

  434.         end
    435. 

  435. 

  436. 

  437.         subroutine test_nf_inq_unlimdim()
    438. 

  438.         implicit        none
    439. 

  439. #include "tests.inc"
    440. 

  440.         integer ncid
    441. 

  441.         integer unlimdim
    442. 

  442.         integer err
    443. 

  443. 

  444.         err = nf_inq_unlimdim(BAD_ID, unlimdim)
    445. 

  445.         if (err .ne. NF_EBADID)
    446. 

  446.      +      call errore('bad ncid: ', err)
    447. 

  447.         err = nf_open(testfile, NF_NOWRITE, ncid)
    448. 

  448.         if (err .ne. 0)
    449. 

  449.      +      call errore('nf_open: ', err)
    450. 

  450.         err = nf_inq_unlimdim(ncid, unlimdim)
    451. 

  451.         if (err .ne. 0) then
    452. 

  452.             call errore('nf_inq_unlimdim: ', err)
    453. 

  453.         else if (unlimdim .ne. RECDIM) then
    454. 

  454.             call errori('nf_inq_unlimdim: wrong number returned, ', 
    455. 

  455.      +                  unlimdim)
    456. 

  456.         end if
    457. 

  457.         err = nf_close(ncid)
    458. 

  458.         if (err .ne. 0)
    459. 

  459.      +      call errore('nf_close: ', err)
    460. 

  460.         end
    461. 

  461. 

  462. 

  463.         subroutine test_nf_inq_format()
    464. 

  464.         implicit        none
    465. 

  465. #include "tests.inc"
    466. 

  466.         integer ncid
    467. 

  467.         integer nformat
    468. 

  468.         integer err
    469. 

  469. 

  470.         err = nf_inq_format(BAD_ID, nformat)
    471. 

  471.         if (err .ne. NF_EBADID)
    472. 

  472.      +      call errore('bad ncid: ', err)
    473. 

  473.         err = nf_open(testfile, NF_NOWRITE, ncid)
    474. 

  474.         if (err .ne. 0)
    475. 

  475.      +      call errore('nf_open: ', err)
    476. 

  476.         err = nf_inq_format(ncid, nformat)
    477. 

  477.         if (err .ne. 0) then
    478. 

  478.             call errore('nf_inq_format: ', err)
    479. 

  479.         else if (nformat .ne. nf_format_classic .and.
    480. 

  480.      +           nformat .ne. nf_format_64bit .and.
    481. 

  481.      +           nformat .ne. nf_format_netcdf4 .and.
    482. 

  482.      +           nformat .ne. nf_format_netcdf4_classic) then
    483. 

  483.             call errori('nf_inq_format: wrong format number returned, ', 
    484. 

  484.      +                  nformat)
    485. 

  485.         end if
    486. 

  486.         err = nf_close(ncid)
    487. 

  487.         if (err .ne. 0)
    488. 

  488.      +      call errore('nf_close: ', err)
    489. 

  489.         end
    490. 

  490. 

  491. 

  492.         subroutine test_nf_inq_dimid()
    493. 

  493.         implicit        none
    494. 

  494. #include "tests.inc"
    495. 

  495.         integer ncid
    496. 

  496.         integer dimid
    497. 

  497.         integer i
    498. 

  498.         integer err
    499. 

  499. 

  500.         err = nf_open(testfile, NF_NOWRITE, ncid)
    501. 

  501.         if (err .ne. 0)
    502. 

  502.      +      call errore('nf_open: ', err)
    503. 

  503.         err = nf_inq_dimid(ncid, 'noSuch', dimid)
    504. 

  504.         if (err .ne. NF_EBADDIM)
    505. 

  505.      +      call errore('bad dim name: ', err)
    506. 

  506.         do 1, i = 1, NDIMS
    507. 

  507.             err = nf_inq_dimid(BAD_ID, dim_name(i), dimid)
    508. 

  508.             if (err .ne. NF_EBADID)
    509. 

  509.      +          call errore('bad ncid: ', err)
    510. 

  510.             err = nf_inq_dimid(ncid, dim_name(i), dimid)
    511. 

  511.             if (err .ne. 0) then
    512. 

  512.                 call errore('nf_inq_dimid: ', err)
    513. 

  513.             else if (dimid .ne. i) then
    514. 

  514.                 call errori('expected ', i)
    515. 

  515.                 call errori('got ', dimid)
    516. 

  516.             end if
    517. 

  517. 1       continue
    518. 

  518.         err = nf_close(ncid)
    519. 

  519.         if (err .ne. 0)
    520. 

  520.      +      call errore('nf_close: ', err)
    521. 

  521.         end
    522. 

  522. 

  523. 

  524.         subroutine test_nf_inq_dim()
    525. 

  525.         implicit        none
    526. 

  526. #include "tests.inc"
    527. 

  527.         integer ncid
    528. 

  528.         integer i
    529. 

  529.         integer err
    530. 

  530.         character*(NF_MAX_NAME) name
    531. 

  531.         integer length
    532. 

  532. 

  533.         err = nf_open(testfile, NF_NOWRITE, ncid)
    534. 

  534.         if (err .ne. 0)
    535. 

  535.      +      call errore('nf_open: ', err)
    536. 

  536.         do 1, i = 1, NDIMS
    537. 

  537.             err = nf_inq_dim(BAD_ID, i, name, length)
    538. 

  538.             if (err .ne. NF_EBADID)
    539. 

  539.      +          call errore('bad ncid: ', err)
    540. 

  540.             err = nf_inq_dim(ncid, BAD_DIMID, name, length)
    541. 

  541.             if (err .ne. NF_EBADDIM)
    542. 

  542.      +          call errore('bad dimid: ', err)
    543. 

  543.             err = nf_inq_dim(ncid, i, name, length)
    544. 

  544.             if (err .ne. 0) then
    545. 

  545.                 call errore('nf_inq_dim: ', err)
    546. 

  546.             else if (dim_name(i) .ne. name)  then
    547. 

  547.                 call errorc('name unexpected: ', name)
    548. 

  548.             else if (dim_len(i) .ne. length) then
    549. 

  549.                 call errori('size unexpected: ', length)
    550. 

  550.             end if
    551. 

  551. 1       continue
    552. 

  552.         err = nf_close(ncid)
    553. 

  553.         if (err .ne. 0)
    554. 

  554.      +      call errore('nf_close: ', err)
    555. 

  555.         end
    556. 

  556. 

  557. 

  558.         subroutine test_nf_inq_dimlen()
    559. 

  559.         implicit        none
    560. 

  560. #include "tests.inc"
    561. 

  561.         integer ncid
    562. 

  562.         integer i
    563. 

  563.         integer err
    564. 

  564.         integer length
    565. 

  565. 

  566.         err = nf_open(testfile, NF_NOWRITE, ncid)
    567. 

  567.         if (err .ne. 0)
    568. 

  568.      +      call errore('nf_open: ', err)
    569. 

  569.         do 1, i = 1, NDIMS
    570. 

  570.             err = nf_inq_dimlen(BAD_ID, i, length)
    571. 

  571.             if (err .ne. NF_EBADID)
    572. 

  572.      +          call errore('bad ncid: ', err)
    573. 

  573.             err = nf_inq_dimlen(ncid, BAD_DIMID, length)
    574. 

  574.             if (err .ne. NF_EBADDIM)
    575. 

  575.      +          call errore('bad dimid: ', err)
    576. 

  576.             err = nf_inq_dimlen(ncid, i, length)
    577. 

  577.             if (err .ne. 0) then
    578. 

  578.                 call errore('nf_inq_dimlen: ', err)
    579. 

  579.             else if (dim_len(i) .ne. length) then
    580. 

  580.                 call errori('size unexpected: ', length)
    581. 

  581.             end if
    582. 

  582. 1       continue
    583. 

  583.         err = nf_close(ncid)
    584. 

  584.         if (err .ne. 0)
    585. 

  585.      +      call errore('nf_close: ', err)
    586. 

  586.         end
    587. 

  587. 

  588. 

  589.         subroutine test_nf_inq_dimname()
    590. 

  590.         implicit        none
    591. 

  591. #include "tests.inc"
    592. 

  592.         integer ncid
    593. 

  593.         integer i
    594. 

  594.         integer err
    595. 

  595.         character*(NF_MAX_NAME)  name
    596. 

  596. 

  597.         err = nf_open(testfile, NF_NOWRITE, ncid)
    598. 

  598.         if (err .ne. 0)
    599. 

  599.      +      call errore('nf_open: ', err)
    600. 

  600.         do 1, i = 1, NDIMS
    601. 

  601.             err = nf_inq_dimname(BAD_ID, i, name)
    602. 

  602.             if (err .ne. NF_EBADID)
    603. 

  603.      +          call errore('bad ncid: ', err)
    604. 

  604.             err = nf_inq_dimname(ncid, BAD_DIMID, name)
    605. 

  605.             if (err .ne. NF_EBADDIM)
    606. 

  606.      +          call errore('bad dimid: ', err)
    607. 

  607.             err = nf_inq_dimname(ncid, i, name)
    608. 

  608.             if (err .ne. 0) then
    609. 

  609.                 call errore('nf_inq_dimname: ', err)
    610. 

  610.             else if (dim_name(i) .ne. name)  then
    611. 

  611.                 call errorc('name unexpected: ', name)
    612. 

  612.             end if
    613. 

  613. 1       continue
    614. 

  614.         err = nf_close(ncid)
    615. 

  615.         if (err .ne. 0)
    616. 

  616.      +      call errore('nf_close: ', err)
    617. 

  617.         end
    618. 

  618. 

  619. 

  620.         subroutine test_nf_inq_varid()
    621. 

  621.         implicit        none
    622. 

  622. #include "tests.inc"
    623. 

  623.         integer ncid
    624. 

  624.         integer vid
    625. 

  625.         integer i
    626. 

  626.         integer err
    627. 

  627. 

  628.         err = nf_open(testfile, NF_NOWRITE, ncid)
    629. 

  629.         if (err .ne. 0)
    630. 

  630.      +      call errore('nf_open: ', err)
    631. 

  631. 

  632.         err = nf_inq_varid(ncid, 'noSuch', vid)
    633. 

  633.         if (err .ne. NF_ENOTVAR)
    634. 

  634.      +      call errore('bad ncid: ', err)
    635. 

  635. 

  636.         do 1, i = 1, NVARS
    637. 

  637.             err = nf_inq_varid(BAD_ID, var_name(i), vid)
    638. 

  638.             if (err .ne. NF_EBADID)
    639. 

  639.      +          call errore('bad ncid: ', err)
    640. 

  640.             err = nf_inq_varid(ncid, var_name(i), vid)
    641. 

  641.             if (err .ne. 0) then
    642. 

  642.                 call errore('nf_inq_varid: ', err)
    643. 

  643.             else if (vid .ne. i) then
    644. 

  644.                 call errori('varid unexpected: ', vid)
    645. 

  645.             endif
    646. 

  646. 1       continue
    647. 

  647. 

  648.         err = nf_close(ncid)
    649. 

  649.         if (err .ne. 0)
    650. 

  650.      +      call errore('nf_close: ', err)
    651. 

  651.         end
    652. 

  652. 

  653. 

  654.         subroutine test_nf_inq_var()
    655. 

  655.         implicit        none
    656. 

  656. #include "tests.inc"
    657. 

  657.         integer ncid
    658. 

  658.         integer i
    659. 

  659.         integer err
    660. 

  660.         character*(NF_MAX_NAME) name
    661. 

  661.         integer datatype
    662. 

  662.         integer ndims
    663. 

  663.         integer dimids(MAX_RANK)
    664. 

  664.         integer na
    665. 

  665. 

  666.         err = nf_open(testfile, NF_NOWRITE, ncid)
    667. 

  667.         if (err .ne. 0)
    668. 

  668.      +      call errore('nf_open: ', err)
    669. 

  669.         do 1, i = 1, NVARS
    670. 

  670.             err = nf_inq_var(BAD_ID, i, name, datatype, ndims, dimids, 
    671. 

  671.      +                       na)
    672. 

  672.             if (err .ne. NF_EBADID)
    673. 

  673.      +          call errore('bad ncid: ', err)
    674. 

  674.             err = nf_inq_var(ncid,BAD_VARID,name,datatype,ndims,dimids,
    675. 

  675.      +                       na)
    676. 

  676.             if (err .ne. NF_ENOTVAR)
    677. 

  677.      +          call errore('bad var id: ', err)
    678. 

  678.             err = nf_inq_var(ncid, i, name, datatype, ndims, dimids, 
    679. 

  679.      +                       na)
    680. 

  680.             if (err .ne. 0) then
    681. 

  681.                 call errore('nf_inq_var: ', err)
    682. 

  682.             else if (var_name(i) .ne. name)  then
    683. 

  683.                 call errorc('name unexpected: ', name)
    684. 

  684.             else if (var_type(i) .ne. datatype) then
    685. 

  685.                 call errori('type unexpected: ', datatype)
    686. 

  686.             else if (var_rank(i) .ne. ndims) then
    687. 

  687.                 call errori('ndims expected: ', ndims)
    688. 

  688.             else if (.not.int_vec_eq(var_dimid(1,i),dimids,ndims)) then
    689. 

  689.                 call error('unexpected dimid')
    690. 

  690.             else if (var_natts(i) .ne. na) then
    691. 

  691.                 call errori('natts unexpected: ', na)
    692. 

  692.             end if
    693. 

  693. 1       continue
    694. 

  694.         err = nf_close(ncid)
    695. 

  695.         if (err .ne. 0)
    696. 

  696.      +      call errore('nf_close: ', err)
    697. 

  697.         end
    698. 

  698. 

  699. 

  700.         subroutine test_nf_inq_vardimid()
    701. 

  701.         implicit        none
    702. 

  702. #include "tests.inc"
    703. 

  703.         integer ncid
    704. 

  704.         integer i
    705. 

  705.         integer err
    706. 

  706.         integer dimids(MAX_RANK)
    707. 

  707. 

  708.         err = nf_open(testfile, NF_NOWRITE, ncid)
    709. 

  709.         if (err .ne. 0)
    710. 

  710.      +      call errore('nf_open: ', err)
    711. 

  711.         do 1, i = 1, NVARS
    712. 

  712.             err = nf_inq_vardimid(BAD_ID, i, dimids)
    713. 

  713.             if (err .ne. NF_EBADID)
    714. 

  714.      +          call errore('bad ncid: ', err)
    715. 

  715.             err = nf_inq_vardimid(ncid, BAD_VARID, dimids)
    716. 

  716.             if (err .ne. NF_ENOTVAR)
    717. 

  717.      +          call errore('bad var id: ', err)
    718. 

  718.             err = nf_inq_vardimid(ncid, i, dimids)
    719. 

  719.             if (err .ne. 0) then
    720. 

  720.                 call errore('nf_inq_vardimid: ', err)
    721. 

  721.             else if (.not.int_vec_eq(var_dimid(1,i), dimids, 
    722. 

  722.      +               var_rank(i))) then
    723. 

  723.                 call error('unexpected dimid')
    724. 

  724.             end if
    725. 

  725. 1       continue
    726. 

  726.         err = nf_close(ncid)
    727. 

  727.         if (err .ne. 0)
    728. 

  728.      +      call errore('nf_close: ', err)
    729. 

  729.         end
    730. 

  730. 

  731. 

  732.         subroutine test_nf_inq_varname()
    733. 

  733.         implicit        none
    734. 

  734. #include "tests.inc"
    735. 

  735.         integer ncid
    736. 

  736.         integer i
    737. 

  737.         integer err
    738. 

  738.         character*(NF_MAX_NAME) name
    739. 

  739. 

  740.         err = nf_open(testfile, NF_NOWRITE, ncid)
    741. 

  741.         if (err .ne. 0)
    742. 

  742.      +      call errore('nf_open: ', err)
    743. 

  743.         do 1, i = 1, NVARS
    744. 

  744.             err = nf_inq_varname(BAD_ID, i, name)
    745. 

  745.             if (err .ne. NF_EBADID)
    746. 

  746.      +          call errore('bad ncid: ', err)
    747. 

  747.             err = nf_inq_varname(ncid, BAD_VARID, name)
    748. 

  748.             if (err .ne. NF_ENOTVAR)
    749. 

  749.      +          call errore('bad var id: ', err)
    750. 

  750.             err = nf_inq_varname(ncid, i, name)
    751. 

  751.             if (err .ne. 0) then
    752. 

  752.                 call errore('nf_inq_varname: ', err)
    753. 

  753.             else if (var_name(i) .ne. name)  then
    754. 

  754.                 call errorc('name unexpected: ', name)
    755. 

  755.             end if
    756. 

  756. 1       continue
    757. 

  757.         err = nf_close(ncid)
    758. 

  758.         if (err .ne. 0)
    759. 

  759.      +      call errore('nf_close: ', err)
    760. 

  760.         end
    761. 

  761. 

  762. 

  763.         subroutine test_nf_inq_varnatts()
    764. 

  764.         implicit        none
    765. 

  765. #include "tests.inc"
    766. 

  766.         integer ncid
    767. 

  767.         integer i
    768. 

  768.         integer err
    769. 

  769.         integer na
    770. 

  770. 

  771.         err = nf_open(testfile, NF_NOWRITE, ncid)
    772. 

  772.         if (err .ne. 0)
    773. 

  773.      +      call errore('nf_open: ', err)
    774. 

  774.         do 1, i = 0, NVARS	! start with global attributes
    775. 

  775.             err = nf_inq_varnatts(BAD_ID, i, na)
    776. 

  776.             if (err .ne. NF_EBADID)
    777. 

  777.      +          call errore('bad ncid: ', err)
    778. 

  778.             err = nf_inq_varnatts(ncid, BAD_VARID, na)
    779. 

  779.             if (err .ne. NF_ENOTVAR)
    780. 

  780.      +          call errore('bad var id: ', err)
    781. 

  781.             err = nf_inq_varnatts(ncid, VARID(i), na)
    782. 

  782.             if (err .ne. 0) then
    783. 

  783.                 call errore('nf_inq_varnatts: ', err)
    784. 

  784.             else if (NATTS(i) .ne. na) then	! works for global attributes
    785. 

  785.                 call errori('natts unexpected: ', na)
    786. 

  786.             end if
    787. 

  787. 1       continue
    788. 

  788.         err = nf_close(ncid)
    789. 

  789.         if (err .ne. 0)
    790. 

  790.      +      call errore('nf_close: ', err)
    791. 

  791.         end
    792. 

  792. 

  793. 

  794.         subroutine test_nf_inq_varndims()
    795. 

  795.         implicit        none
    796. 

  796. #include "tests.inc"
    797. 

  797.         integer ncid
    798. 

  798.         integer i
    799. 

  799.         integer err
    800. 

  800.         integer ndims
    801. 

  801. 

  802.         err = nf_open(testfile, NF_NOWRITE, ncid)
    803. 

  803.         if (err .ne. 0)
    804. 

  804.      +      call errore('nf_open: ', err)
    805. 

  805.         do 1, i = 1, NVARS
    806. 

  806.             err = nf_inq_varndims(BAD_ID, i, ndims)
    807. 

  807.             if (err .ne. NF_EBADID)
    808. 

  808.      +          call errore('bad ncid: ', err)
    809. 

  809.             err = nf_inq_varndims(ncid, BAD_VARID, ndims)
    810. 

  810.             if (err .ne. NF_ENOTVAR)
    811. 

  811.      +          call errore('bad var id: ', err)
    812. 

  812.             err = nf_inq_varndims(ncid, i, ndims)
    813. 

  813.             if (err .ne. 0) then
    814. 

  814.                 call errore('nf_inq_varndims: ', err)
    815. 

  815.             else if (var_rank(i) .ne. ndims) then
    816. 

  816.                 call errori('ndims unexpected: ', ndims)
    817. 

  817.             end if
    818. 

  818. 1       continue
    819. 

  819.         err = nf_close(ncid)
    820. 

  820.         if (err .ne. 0)
    821. 

  821.      +      call errore('nf_close: ', err)
    822. 

  822.         end
    823. 

  823. 

  824. 

  825.         subroutine test_nf_inq_vartype()
    826. 

  826.         implicit        none
    827. 

  827. #include "tests.inc"
    828. 

  828.         integer ncid
    829. 

  829.         integer i
    830. 

  830.         integer err
    831. 

  831.         integer datatype
    832. 

  832. 

  833.         err = nf_open(testfile, NF_NOWRITE, ncid)
    834. 

  834.         if (err .ne. 0)
    835. 

  835.      +      call errore('nf_open: ', err)
    836. 

  836.         do 1, i = 1, NVARS
    837. 

  837.             err = nf_inq_vartype(BAD_ID, i, datatype)
    838. 

  838.             if (err .ne. NF_EBADID)
    839. 

  839.      +          call errore('bad ncid: ', err)
    840. 

  840.             err = nf_inq_vartype(ncid, BAD_VARID, datatype)
    841. 

  841.             if (err .ne. NF_ENOTVAR)
    842. 

  842.      +          call errore('bad var id: ', err)
    843. 

  843.             err = nf_inq_vartype(ncid, i, datatype)
    844. 

  844.             if (err .ne. 0) then
    845. 

  845.                 call errore('nf_inq_vartype: ', err)
    846. 

  846.             else if (var_type(i) .ne. datatype) then
    847. 

  847.                 call errori('type unexpected: ', datatype)
    848. 

  848.             end if
    849. 

  849. 1       continue
    850. 

  850.         err = nf_close(ncid)
    851. 

  851.         if (err .ne. 0)
    852. 

  852.      +      call errore('nf_close: ', err)
    853. 

  853.         end
    854. 

  854. 

  855. 

  856.         subroutine test_nf_inq_att()
    857. 

  857.         implicit        none
    858. 

  858. #include "tests.inc"
    859. 

  859.         integer ncid
    860. 

  860.         integer i
    861. 

  861.         integer j
    862. 

  862.         integer err
    863. 

  863.         integer t
    864. 

  864.         integer n
    865. 

  865. 

  866.         err = nf_open(testfile, NF_NOWRITE, ncid)
    867. 

  867.         if (err .ne. 0) 
    868. 

  868.      +      call errore('nf_open: ', err)
    869. 

  869. 

  870.         do 1, i = 0, NVARS
    871. 

  871.             do 2, j = 1, NATTS(i)
    872. 

  872.                 err = nf_inq_att(BAD_ID, i, ATT_NAME(j,i), t, n)
    873. 

  873.                 if (err .ne. NF_EBADID) 
    874. 

  874.      +              call errore('bad ncid: ', err)
    875. 

  875.                 err = nf_inq_att(ncid, BAD_VARID, ATT_NAME(j,i), t, n)
    876. 

  876.                 if (err .ne. NF_ENOTVAR) 
    877. 

  877.      +              call errore('bad var id: ', err)
    878. 

  878.                 err = nf_inq_att(ncid, i, 'noSuch', t, n)
    879. 

  879.                 if (err .ne. NF_ENOTATT) 
    880. 

  880.      +              call errore('Bad attribute name: ', err)
    881. 

  881.                 err = nf_inq_att(ncid, i, ATT_NAME(j,i), t, n)
    882. 

  882.                 if (err .ne. 0) then
    883. 

  883.                     call error(nf_strerror(err))
    884. 

  884.                 else
    885. 

  885.                     if (t .ne. ATT_TYPE(j,i))
    886. 

  886.      +                  call error('type not that expected')
    887. 

  887.                     if (n .ne. ATT_LEN(j,i)) 
    888. 

  888.      +                  call error('length not that expected')
    889. 

  889.                 end if
    890. 

  890. 2           continue
    891. 

  891. 1       continue
    892. 

  892. 

  893.         err = nf_close(ncid)
    894. 

  894.         if (err .ne. 0)
    895. 

  895.      +      call errore('nf_close: ', err)
    896. 

  896.         end
    897. 

  897. 

  898. 

  899.         subroutine test_nf_inq_attlen()
    900. 

  900.         implicit        none
    901. 

  901. #include "tests.inc"
    902. 

  902.         integer ncid
    903. 

  903.         integer i
    904. 

  904.         integer j
    905. 

  905.         integer err
    906. 

  906.         integer len
    907. 

  907. 

  908.         err = nf_open(testfile, NF_NOWRITE, ncid)
    909. 

  909.         if (err .ne. 0)
    910. 

  910.      +      call errore('nf_open: ', err)
    911. 

  911. 

  912.         do 1, i = 0, NVARS
    913. 

  913.             err = nf_inq_attlen(ncid, i, 'noSuch', len)
    914. 

  914.             if (err .ne. NF_ENOTATT)
    915. 

  915.      +          call errore('Bad attribute name: ', err)
    916. 

  916.             do 2, j = 1, NATTS(i)
    917. 

  917.                 err = nf_inq_attlen(BAD_ID, i, ATT_NAME(j,i), len)
    918. 

  918.                 if (err .ne. NF_EBADID)
    919. 

  919.      +              call errore('bad ncid: ', err)
    920. 

  920.                 err = nf_inq_attlen(ncid, BAD_VARID, ATT_NAME(j,i), len)
    921. 

  921.                 if (err .ne. NF_ENOTVAR)
    922. 

  922.      +              call errore('bad varid: ', err)
    923. 

  923.                 err = nf_inq_attlen(ncid, i, ATT_NAME(j,i), len)
    924. 

  924.                 if (err .ne. 0) then
    925. 

  925.                     call error(nf_strerror(err))
    926. 

  926.                 else
    927. 

  927.                     if (len .ne. ATT_LEN(j,i))
    928. 

  928.      +                  call error('len not that expected')
    929. 

  929.                 end if
    930. 

  930. 2           continue
    931. 

  931. 1       continue
    932. 

  932. 

  933.         err = nf_close(ncid)
    934. 

  934.         if (err .ne. 0)
    935. 

  935.      +      call errore('nf_close: ', err)
    936. 

  936.         end
    937. 

  937. 

  938. 

  939.         subroutine test_nf_inq_atttype()
    940. 

  940.         implicit        none
    941. 

  941. #include "tests.inc"
    942. 

  942.         integer ncid
    943. 

  943.         integer i
    944. 

  944.         integer j
    945. 

  945.         integer err
    946. 

  946.         integer datatype
    947. 

  947. 

  948.         err = nf_open(testfile, NF_NOWRITE, ncid)
    949. 

  949.         if (err .ne. 0)
    950. 

  950.      +      call errore('nf_open: ', err)
    951. 

  951. 

  952.         do 1, i = 0, NVARS
    953. 

  953.             err = nf_inq_atttype(ncid, i, 'noSuch', datatype)
    954. 

  954.             if (err .ne. NF_ENOTATT)
    955. 

  955.      +          call errore('Bad attribute name: ', err)
    956. 

  956.             do 2, j = 1, NATTS(i)
    957. 

  957.                 err = nf_inq_atttype(BAD_ID, i, ATT_NAME(j,i), datatype)
    958. 

  958.                 if (err .ne. NF_EBADID)
    959. 

  959.      +              call errore('bad ncid: ', err)
    960. 

  960.                 err = nf_inq_atttype(ncid, BAD_VARID, ATT_NAME(j,i), 
    961. 

  961.      +                               datatype)
    962. 

  962.                 if (err .ne. NF_ENOTVAR)
    963. 

  963.      +              call errore('bad varid: ', err)
    964. 

  964.                 err = nf_inq_atttype(ncid, i, ATT_NAME(j,i), datatype)
    965. 

  965.                 if (err .ne. 0) then
    966. 

  966.                     call error(nf_strerror(err))
    967. 

  967.                 else
    968. 

  968.                     if (datatype .ne. ATT_TYPE(j,i))
    969. 

  969.      +                  call error('type not that expected')
    970. 

  970.                 end if
    971. 

  971. 2           continue
    972. 

  972. 1       continue
    973. 

  973. 

  974.         err = nf_close(ncid)
    975. 

  975.         if (err .ne. 0)
    976. 

  976.      +      call errore('nf_close: ', err)
    977. 

  977.         end
    978. 

  978. 

  979. 

  980.         subroutine test_nf_inq_attname()
    981. 

  981.         implicit        none
    982. 

  982. #include "tests.inc"
    983. 

  983.         integer ncid
    984. 

  984.         integer i
    985. 

  985.         integer j
    986. 

  986.         integer err
    987. 

  987.         character*(NF_MAX_NAME) name
    988. 

  988. 

  989.         err = nf_open(testfile, NF_NOWRITE, ncid)
    990. 

  990.         if (err .ne. 0)
    991. 

  991.      +      call errore('nf_open: ', err)
    992. 

  992. 

  993.         do 1, i = 0, NVARS
    994. 

  994.             err = nf_inq_attname(ncid, i, BAD_ATTNUM, name)
    995. 

  995.             if (err .ne. NF_ENOTATT)
    996. 

  996.      +          call errore('Bad attribute number: ', err)
    997. 

  997.             err = nf_inq_attname(ncid, i, NATTS(i)+1, name)
    998. 

  998.             if (err .ne. NF_ENOTATT)
    999. 

  999.      +          call errore('Bad attribute number: ', err)
    1000. 

  1000.             do 2, j = 1, NATTS(i)
    1001. 

  1001.                 err = nf_inq_attname(BAD_ID, i, j, name)
    1002. 

  1002.                 if (err .ne. NF_EBADID)
    1003. 

  1003.      +              call errore('bad ncid: ', err)
    1004. 

  1004.                 err = nf_inq_attname(ncid, BAD_VARID, j, name)
    1005. 

  1005.                 if (err .ne. NF_ENOTVAR)
    1006. 

  1006.      +              call errore('bad var id: ', err)
    1007. 

  1007.                 err = nf_inq_attname(ncid, i, j, name)
    1008. 

  1008.                 if (err .ne. 0) then
    1009. 

  1009.                     call error(nf_strerror(err))
    1010. 

  1010.                 else
    1011. 

  1011.                     if (ATT_NAME(j,i) .ne. name)
    1012. 

  1012.      +                  call error('name not that expected')
    1013. 

  1013.                 end if
    1014. 

  1014. 2           continue
    1015. 

  1015. 1       continue
    1016. 

  1016. 

  1017.         err = nf_close(ncid)
    1018. 

  1018.         if (err .ne. 0)
    1019. 

  1019.      +      call errore('nf_close: ', err)
    1020. 

  1020.         end
    1021. 

  1021. 

  1022. 

  1023.         subroutine test_nf_inq_attid()
    1024. 

  1024.         implicit        none
    1025. 

  1025. #include "tests.inc"
    1026. 

  1026.         integer ncid
    1027. 

  1027.         integer i
    1028. 

  1028.         integer j
    1029. 

  1029.         integer err
    1030. 

  1030.         integer attnum
    1031. 

  1031. 

  1032.         err = nf_open(testfile, NF_NOWRITE, ncid)
    1033. 

  1033.         if (err .ne. 0)
    1034. 

  1034.      +      call errore('nf_open: ', err)
    1035. 

  1035. 

  1036.         do 1, i = 0, NVARS
    1037. 

  1037.             err = nf_inq_attid(ncid, i, 'noSuch', attnum)
    1038. 

  1038.             if (err .ne. NF_ENOTATT)
    1039. 

  1039.      +          call errore('Bad attribute name: ', err)
    1040. 

  1040.             do 2, j = 1, NATTS(i)
    1041. 

  1041.                 err = nf_inq_attid(BAD_ID, i, ATT_NAME(j,i), attnum)
    1042. 

  1042.                 if (err .ne. NF_EBADID)
    1043. 

  1043.      +              call errore('bad ncid: ', err)
    1044. 

  1044.                 err = nf_inq_attid(ncid, BAD_VARID, ATT_NAME(j,i), 
    1045. 

  1045.      +                             attnum)
    1046. 

  1046.                 if (err .ne. NF_ENOTVAR)
    1047. 

  1047.      +              call errore('bad varid: ', err)
    1048. 

  1048.                 err = nf_inq_attid(ncid, i, ATT_NAME(j,i), attnum)
    1049. 

  1049.                 if (err .ne. 0) then
    1050. 

  1050.                     call error(nf_strerror(err))
    1051. 

  1051.                 else
    1052. 

  1052.                     if (attnum .ne. j)
    1053. 

  1053.      +                  call error('attnum not that expected')
    1054. 

  1054.                 end if
    1055. 

  1055. 2           continue
    1056. 

  1056. 1       continue
    1057. 

  1057. 

  1058.         err = nf_close(ncid)
    1059. 

  1059.         if (err .ne. 0)
    1060. 

  1060.      +      call errore('nf_close: ', err)
    1061. 

  1061.         end
    1062. 

  1062.