Changeset 2512 for palm/trunk/SOURCE/netcdf_interface_mod.f90
- Timestamp:
- Oct 4, 2017 8:26:59 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/netcdf_interface_mod.f90
r2302 r2512 25 25 ! ----------------- 26 26 ! $Id$ 27 ! upper bounds of cross section and 3d output changed from nx+1,ny+1 to nx,ny 28 ! no output of ghost layer data any more 29 ! 30 ! 2302 2017-07-03 14:07:20Z suehring 27 31 ! Reading of 3D topography using NetCDF data type NC_BYTE 28 32 ! … … 1307 1311 ! 1308 1312 !-- Define x-axis (for scalar position) 1309 CALL netcdf_create_dim( id_set_3d(av), 'x', nx+ 2, id_dim_x_3d(av), &1313 CALL netcdf_create_dim( id_set_3d(av), 'x', nx+1, id_dim_x_3d(av), & 1310 1314 73 ) 1311 1315 CALL netcdf_create_var( id_set_3d(av), (/ id_dim_x_3d(av) /), 'x', & … … 1314 1318 ! 1315 1319 !-- Define x-axis (for u position) 1316 CALL netcdf_create_dim( id_set_3d(av), 'xu', nx+ 2, id_dim_xu_3d(av), &1320 CALL netcdf_create_dim( id_set_3d(av), 'xu', nx+1, id_dim_xu_3d(av), & 1317 1321 358 ) 1318 1322 CALL netcdf_create_var( id_set_3d(av), (/ id_dim_xu_3d(av) /), 'xu', & … … 1321 1325 ! 1322 1326 !-- Define y-axis (for scalar position) 1323 CALL netcdf_create_dim( id_set_3d(av), 'y', ny+ 2, id_dim_y_3d(av), &1327 CALL netcdf_create_dim( id_set_3d(av), 'y', ny+1, id_dim_y_3d(av), & 1324 1328 76 ) 1325 1329 CALL netcdf_create_var( id_set_3d(av), (/ id_dim_y_3d(av) /), 'y', & … … 1328 1332 ! 1329 1333 !-- Define y-axis (for v position) 1330 CALL netcdf_create_dim( id_set_3d(av), 'yv', ny+ 2, id_dim_yv_3d(av), &1334 CALL netcdf_create_dim( id_set_3d(av), 'yv', ny+1, id_dim_yv_3d(av), & 1331 1335 361 ) 1332 1336 CALL netcdf_create_var( id_set_3d(av), (/ id_dim_yv_3d(av) /), 'yv', & … … 1537 1541 ! 1538 1542 !-- Write data for x (shifted by +dx/2) and xu axis 1539 ALLOCATE( netcdf_data(0:nx +1) )1540 1541 DO i = 0, nx +11543 ALLOCATE( netcdf_data(0:nx) ) 1544 1545 DO i = 0, nx 1542 1546 netcdf_data(i) = ( i + 0.5 ) * dx 1543 1547 ENDDO … … 1545 1549 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_x_3d(av), & 1546 1550 netcdf_data, start = (/ 1 /), & 1547 count = (/ nx+ 2/) )1551 count = (/ nx+1 /) ) 1548 1552 CALL netcdf_handle_error( 'netcdf_define_header', 83 ) 1549 1553 1550 DO i = 0, nx +11554 DO i = 0, nx 1551 1555 netcdf_data(i) = i * dx 1552 1556 ENDDO … … 1554 1558 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_xu_3d(av), & 1555 1559 netcdf_data, start = (/ 1 /), & 1556 count = (/ nx+ 2/) )1560 count = (/ nx+1 /) ) 1557 1561 CALL netcdf_handle_error( 'netcdf_define_header', 385 ) 1558 1562 … … 1561 1565 ! 1562 1566 !-- Write data for y (shifted by +dy/2) and yv axis 1563 ALLOCATE( netcdf_data(0:ny +1) )1564 1565 DO i = 0, ny +11567 ALLOCATE( netcdf_data(0:ny) ) 1568 1569 DO i = 0, ny 1566 1570 netcdf_data(i) = ( i + 0.5_wp ) * dy 1567 1571 ENDDO … … 1569 1573 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_y_3d(av), & 1570 1574 netcdf_data, start = (/ 1 /), & 1571 count = (/ ny+ 2/) )1575 count = (/ ny+1 /) ) 1572 1576 CALL netcdf_handle_error( 'netcdf_define_header', 84 ) 1573 1577 1574 DO i = 0, ny +11578 DO i = 0, ny 1575 1579 netcdf_data(i) = i * dy 1576 1580 ENDDO … … 1578 1582 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_yv_3d(av), & 1579 1583 netcdf_data, start = (/ 1 /), & 1580 count = (/ ny+ 2/))1584 count = (/ ny+1 /)) 1581 1585 CALL netcdf_handle_error( 'netcdf_define_header', 387 ) 1582 1586 … … 1612 1616 netcdf_data_format > 4 ) THEN 1613 1617 1614 IF ( nxr == nx .AND. nyn /= ny ) THEN1615 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), &1616 zu_s_inner(nxl:nxr+1,nys:nyn), &1617 start = (/ nxl+1, nys+1 /), &1618 count = (/ nxr-nxl+2, nyn-nys+1 /) )1619 ELSEIF ( nxr /= nx .AND. nyn == ny ) THEN1620 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), &1621 zu_s_inner(nxl:nxr,nys:nyn+1), &1622 start = (/ nxl+1, nys+1 /), &1623 count = (/ nxr-nxl+1, nyn-nys+2 /) )1624 ELSEIF ( nxr == nx .AND. nyn == ny ) THEN1625 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), &1626 zu_s_inner(nxl:nxr+1,nys:nyn+1), &1627 start = (/ nxl+1, nys+1 /), &1628 count = (/ nxr-nxl+2, nyn-nys+2 /) )1629 ELSE1618 ! IF ( nxr == nx .AND. nyn /= ny ) THEN 1619 ! nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), & 1620 ! zu_s_inner(nxl:nxr+1,nys:nyn), & 1621 ! start = (/ nxl+1, nys+1 /), & 1622 ! count = (/ nxr-nxl+2, nyn-nys+1 /) ) 1623 ! ELSEIF ( nxr /= nx .AND. nyn == ny ) THEN 1624 ! nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), & 1625 ! zu_s_inner(nxl:nxr,nys:nyn+1), & 1626 ! start = (/ nxl+1, nys+1 /), & 1627 ! count = (/ nxr-nxl+1, nyn-nys+2 /) ) 1628 ! ELSEIF ( nxr == nx .AND. nyn == ny ) THEN 1629 ! nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), & 1630 ! zu_s_inner(nxl:nxr+1,nys:nyn+1), & 1631 ! start = (/ nxl+1, nys+1 /), & 1632 ! count = (/ nxr-nxl+2, nyn-nys+2 /) ) 1633 ! ELSE 1630 1634 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), & 1631 1635 zu_s_inner(nxl:nxr,nys:nyn), & 1632 1636 start = (/ nxl+1, nys+1 /), & 1633 1637 count = (/ nxr-nxl+1, nyn-nys+1 /) ) 1634 ENDIF1638 ! ENDIF 1635 1639 CALL netcdf_handle_error( 'netcdf_define_header', 419 ) 1636 1640 1637 IF ( nxr == nx .AND. nyn /= ny ) THEN1638 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), &1639 zw_w_inner(nxl:nxr+1,nys:nyn), &1640 start = (/ nxl+1, nys+1 /), &1641 count = (/ nxr-nxl+2, nyn-nys+1 /) )1642 ELSEIF ( nxr /= nx .AND. nyn == ny ) THEN1643 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), &1644 zw_w_inner(nxl:nxr,nys:nyn+1), &1645 start = (/ nxl+1, nys+1 /), &1646 count = (/ nxr-nxl+1, nyn-nys+2 /) )1647 ELSEIF ( nxr == nx .AND. nyn == ny ) THEN1648 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), &1649 zw_w_inner(nxl:nxr+1,nys:nyn+1), &1650 start = (/ nxl+1, nys+1 /), &1651 count = (/ nxr-nxl+2, nyn-nys+2 /) )1652 ELSE1641 ! IF ( nxr == nx .AND. nyn /= ny ) THEN 1642 ! nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), & 1643 ! zw_w_inner(nxl:nxr+1,nys:nyn), & 1644 ! start = (/ nxl+1, nys+1 /), & 1645 ! count = (/ nxr-nxl+2, nyn-nys+1 /) ) 1646 ! ELSEIF ( nxr /= nx .AND. nyn == ny ) THEN 1647 ! nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), & 1648 ! zw_w_inner(nxl:nxr,nys:nyn+1), & 1649 ! start = (/ nxl+1, nys+1 /), & 1650 ! count = (/ nxr-nxl+1, nyn-nys+2 /) ) 1651 ! ELSEIF ( nxr == nx .AND. nyn == ny ) THEN 1652 ! nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), & 1653 ! zw_w_inner(nxl:nxr+1,nys:nyn+1), & 1654 ! start = (/ nxl+1, nys+1 /), & 1655 ! count = (/ nxr-nxl+2, nyn-nys+2 /) ) 1656 ! ELSE 1653 1657 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), & 1654 1658 zw_w_inner(nxl:nxr,nys:nyn), & 1655 1659 start = (/ nxl+1, nys+1 /), & 1656 1660 count = (/ nxr-nxl+1, nyn-nys+1 /) ) 1657 ENDIF1661 ! ENDIF 1658 1662 CALL netcdf_handle_error( 'netcdf_define_header', 420 ) 1659 1663 … … 1961 1965 ! 1962 1966 !-- Define x-axis (for scalar position) 1963 CALL netcdf_create_dim( id_set_xy(av), 'x', nx+ 2, id_dim_x_xy(av), &1967 CALL netcdf_create_dim( id_set_xy(av), 'x', nx+1, id_dim_x_xy(av), & 1964 1968 113 ) 1965 1969 CALL netcdf_create_var( id_set_xy(av), (/ id_dim_x_xy(av) /), 'x', & … … 1968 1972 ! 1969 1973 !-- Define x-axis (for u position) 1970 CALL netcdf_create_dim( id_set_xy(av), 'xu', nx+ 2, &1974 CALL netcdf_create_dim( id_set_xy(av), 'xu', nx+1, & 1971 1975 id_dim_xu_xy(av), 388 ) 1972 1976 CALL netcdf_create_var( id_set_xy(av), (/ id_dim_xu_xy(av) /), 'xu', & … … 1975 1979 ! 1976 1980 !-- Define y-axis (for scalar position) 1977 CALL netcdf_create_dim( id_set_xy(av), 'y', ny+ 2, id_dim_y_xy(av), &1981 CALL netcdf_create_dim( id_set_xy(av), 'y', ny+1, id_dim_y_xy(av), & 1978 1982 116 ) 1979 1983 CALL netcdf_create_var( id_set_xy(av), (/ id_dim_y_xy(av) /), 'y', & … … 1982 1986 ! 1983 1987 !-- Define y-axis (for scalar position) 1984 CALL netcdf_create_dim( id_set_xy(av), 'yv', ny+ 2, &1988 CALL netcdf_create_dim( id_set_xy(av), 'yv', ny+1, & 1985 1989 id_dim_yv_xy(av), 364 ) 1986 1990 CALL netcdf_create_var( id_set_xy(av), (/ id_dim_yv_xy(av) /), 'yv', & … … 2255 2259 ! 2256 2260 !-- Write data for x (shifted by +dx/2) and xu axis 2257 ALLOCATE( netcdf_data(0:nx +1) )2258 2259 DO i = 0, nx +12261 ALLOCATE( netcdf_data(0:nx) ) 2262 2263 DO i = 0, nx 2260 2264 netcdf_data(i) = ( i + 0.5_wp ) * dx 2261 2265 ENDDO … … 2263 2267 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_x_xy(av), & 2264 2268 netcdf_data, start = (/ 1 /), & 2265 count = (/ nx+ 2/) )2269 count = (/ nx+1 /) ) 2266 2270 CALL netcdf_handle_error( 'netcdf_define_header', 127 ) 2267 2271 2268 DO i = 0, nx +12272 DO i = 0, nx 2269 2273 netcdf_data(i) = i * dx 2270 2274 ENDDO … … 2272 2276 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_xu_xy(av), & 2273 2277 netcdf_data, start = (/ 1 /), & 2274 count = (/ nx+ 2/) )2278 count = (/ nx+1 /) ) 2275 2279 CALL netcdf_handle_error( 'netcdf_define_header', 367 ) 2276 2280 … … 2281 2285 ALLOCATE( netcdf_data(0:ny+1) ) 2282 2286 2283 DO i = 0, ny +12287 DO i = 0, ny 2284 2288 netcdf_data(i) = ( i + 0.5_wp ) * dy 2285 2289 ENDDO … … 2287 2291 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_y_xy(av), & 2288 2292 netcdf_data, start = (/ 1 /), & 2289 count = (/ ny+ 2/))2293 count = (/ ny+1 /)) 2290 2294 CALL netcdf_handle_error( 'netcdf_define_header', 128 ) 2291 2295 2292 DO i = 0, ny +12296 DO i = 0, ny 2293 2297 netcdf_data(i) = i * dy 2294 2298 ENDDO … … 2296 2300 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_yv_xy(av), & 2297 2301 netcdf_data, start = (/ 1 /), & 2298 count = (/ ny+ 2/))2302 count = (/ ny+1 /)) 2299 2303 CALL netcdf_handle_error( 'netcdf_define_header', 368 ) 2300 2304 … … 2309 2313 netcdf_data_format > 4 ) THEN 2310 2314 2311 IF ( nxr == nx .AND. nyn /= ny ) THEN2312 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), &2313 zu_s_inner(nxl:nxr+1,nys:nyn), &2314 start = (/ nxl+1, nys+1 /), &2315 count = (/ nxr-nxl+2, nyn-nys+1 /) )2316 ELSEIF ( nxr /= nx .AND. nyn == ny ) THEN2317 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), &2318 zu_s_inner(nxl:nxr,nys:nyn+1), &2319 start = (/ nxl+1, nys+1 /), &2320 count = (/ nxr-nxl+1, nyn-nys+2 /) )2321 ELSEIF ( nxr == nx .AND. nyn == ny ) THEN2322 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), &2323 zu_s_inner(nxl:nxr+1,nys:nyn+1), &2324 start = (/ nxl+1, nys+1 /), &2325 count = (/ nxr-nxl+2, nyn-nys+2 /) )2326 ELSE2315 ! IF ( nxr == nx .AND. nyn /= ny ) THEN 2316 ! nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), & 2317 ! zu_s_inner(nxl:nxr+1,nys:nyn), & 2318 ! start = (/ nxl+1, nys+1 /), & 2319 ! count = (/ nxr-nxl+2, nyn-nys+1 /) ) 2320 ! ELSEIF ( nxr /= nx .AND. nyn == ny ) THEN 2321 ! nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), & 2322 ! zu_s_inner(nxl:nxr,nys:nyn+1), & 2323 ! start = (/ nxl+1, nys+1 /), & 2324 ! count = (/ nxr-nxl+1, nyn-nys+2 /) ) 2325 ! ELSEIF ( nxr == nx .AND. nyn == ny ) THEN 2326 ! nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), & 2327 ! zu_s_inner(nxl:nxr+1,nys:nyn+1), & 2328 ! start = (/ nxl+1, nys+1 /), & 2329 ! count = (/ nxr-nxl+2, nyn-nys+2 /) ) 2330 ! ELSE 2327 2331 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), & 2328 2332 zu_s_inner(nxl:nxr,nys:nyn), & 2329 2333 start = (/ nxl+1, nys+1 /), & 2330 2334 count = (/ nxr-nxl+1, nyn-nys+1 /) ) 2331 ENDIF2335 ! ENDIF 2332 2336 CALL netcdf_handle_error( 'netcdf_define_header', 427 ) 2333 2337 2334 IF ( nxr == nx .AND. nyn /= ny ) THEN2335 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), &2336 zw_w_inner(nxl:nxr+1,nys:nyn), &2337 start = (/ nxl+1, nys+1 /), &2338 count = (/ nxr-nxl+2, nyn-nys+1 /) )2339 ELSEIF ( nxr /= nx .AND. nyn == ny ) THEN2340 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), &2341 zw_w_inner(nxl:nxr,nys:nyn+1), &2342 start = (/ nxl+1, nys+1 /), &2343 count = (/ nxr-nxl+1, nyn-nys+2 /) )2344 ELSEIF ( nxr == nx .AND. nyn == ny ) THEN2345 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), &2346 zw_w_inner(nxl:nxr+1,nys:nyn+1), &2347 start = (/ nxl+1, nys+1 /), &2348 count = (/ nxr-nxl+2, nyn-nys+2 /) )2349 ELSE2338 ! IF ( nxr == nx .AND. nyn /= ny ) THEN 2339 ! nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), & 2340 ! zw_w_inner(nxl:nxr+1,nys:nyn), & 2341 ! start = (/ nxl+1, nys+1 /), & 2342 ! count = (/ nxr-nxl+2, nyn-nys+1 /) ) 2343 ! ELSEIF ( nxr /= nx .AND. nyn == ny ) THEN 2344 ! nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), & 2345 ! zw_w_inner(nxl:nxr,nys:nyn+1), & 2346 ! start = (/ nxl+1, nys+1 /), & 2347 ! count = (/ nxr-nxl+1, nyn-nys+2 /) ) 2348 ! ELSEIF ( nxr == nx .AND. nyn == ny ) THEN 2349 ! nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), & 2350 ! zw_w_inner(nxl:nxr+1,nys:nyn+1), & 2351 ! start = (/ nxl+1, nys+1 /), & 2352 ! count = (/ nxr-nxl+2, nyn-nys+2 /) ) 2353 ! ELSE 2350 2354 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), & 2351 2355 zw_w_inner(nxl:nxr,nys:nyn), & 2352 2356 start = (/ nxl+1, nys+1 /), & 2353 2357 count = (/ nxr-nxl+1, nyn-nys+1 /) ) 2354 ENDIF2358 ! ENDIF 2355 2359 CALL netcdf_handle_error( 'netcdf_define_header', 428 ) 2356 2360 … … 2684 2688 ! 2685 2689 !-- Define x-axis (for scalar position) 2686 CALL netcdf_create_dim( id_set_xz(av), 'x', nx+ 2, id_dim_x_xz(av), &2690 CALL netcdf_create_dim( id_set_xz(av), 'x', nx+1, id_dim_x_xz(av), & 2687 2691 150 ) 2688 2692 CALL netcdf_create_var( id_set_xz(av), (/ id_dim_x_xz(av) /), 'x', & … … 2691 2695 ! 2692 2696 !-- Define x-axis (for u position) 2693 CALL netcdf_create_dim( id_set_xz(av), 'xu', nx+ 2, id_dim_xu_xz(av), &2697 CALL netcdf_create_dim( id_set_xz(av), 'xu', nx+1, id_dim_xu_xz(av), & 2694 2698 372 ) 2695 2699 CALL netcdf_create_var( id_set_xz(av), (/ id_dim_xu_xz(av) /), 'xu', & … … 2940 2944 ! 2941 2945 !-- Write data for x (shifted by +dx/2) and xu axis 2942 ALLOCATE( netcdf_data(0:nx +1) )2943 2944 DO i = 0, nx +12946 ALLOCATE( netcdf_data(0:nx) ) 2947 2948 DO i = 0, nx 2945 2949 netcdf_data(i) = ( i + 0.5_wp ) * dx 2946 2950 ENDDO … … 2948 2952 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_x_xz(av), & 2949 2953 netcdf_data, start = (/ 1 /), & 2950 count = (/ nx+ 2/) )2954 count = (/ nx+1 /) ) 2951 2955 CALL netcdf_handle_error( 'netcdf_define_header', 165 ) 2952 2956 2953 DO i = 0, nx +12957 DO i = 0, nx 2954 2958 netcdf_data(i) = i * dx 2955 2959 ENDDO … … 2957 2961 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_xu_xz(av), & 2958 2962 netcdf_data, start = (/ 1 /), & 2959 count = (/ nx+ 2/) )2963 count = (/ nx+1 /) ) 2960 2964 CALL netcdf_handle_error( 'netcdf_define_header', 377 ) 2961 2965 … … 3337 3341 ! 3338 3342 !-- Define y-axis (for scalar position) 3339 CALL netcdf_create_dim( id_set_yz(av), 'y', ny+ 2, id_dim_y_yz(av), &3343 CALL netcdf_create_dim( id_set_yz(av), 'y', ny+1, id_dim_y_yz(av), & 3340 3344 189 ) 3341 3345 CALL netcdf_create_var( id_set_yz(av), (/ id_dim_y_yz(av) /), 'y', & … … 3344 3348 ! 3345 3349 !-- Define y-axis (for v position) 3346 CALL netcdf_create_dim( id_set_yz(av), 'yv', ny+ 2, id_dim_yv_yz(av), &3350 CALL netcdf_create_dim( id_set_yz(av), 'yv', ny+1, id_dim_yv_yz(av), & 3347 3351 380 ) 3348 3352 CALL netcdf_create_var( id_set_yz(av), (/ id_dim_yv_yz(av) /), 'yv', & … … 3592 3596 ! 3593 3597 !-- Write data for y (shifted by +dy/2) and yv axis 3594 ALLOCATE( netcdf_data(0:ny +1) )3595 3596 DO j = 0, ny +13598 ALLOCATE( netcdf_data(0:ny) ) 3599 3600 DO j = 0, ny 3597 3601 netcdf_data(j) = ( j + 0.5_wp ) * dy 3598 3602 ENDDO … … 3600 3604 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_y_yz(av), & 3601 3605 netcdf_data, start = (/ 1 /), & 3602 count = (/ ny+ 2/) )3606 count = (/ ny+1 /) ) 3603 3607 CALL netcdf_handle_error( 'netcdf_define_header', 204 ) 3604 3608 3605 DO j = 0, ny +13609 DO j = 0, ny 3606 3610 netcdf_data(j) = j * dy 3607 3611 ENDDO … … 3609 3613 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_yv_yz(av), & 3610 3614 netcdf_data, start = (/ 1 /), & 3611 count = (/ ny+ 2/) )3615 count = (/ ny+1 /) ) 3612 3616 CALL netcdf_handle_error( 'netcdf_define_header', 384 ) 3613 3617
Note: See TracChangeset
for help on using the changeset viewer.