Changeset 1308 for palm/trunk/SOURCE/netcdf.f90
- Timestamp:
- Mar 13, 2014 2:58:42 PM (11 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/netcdf.f90 (modified) (36 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/netcdf.f90
r1207 r1308 23 23 ! Current revisions: 24 24 ! ------------------ 25 ! 25 ! +ntime_count, oldmode 26 ! Adjust NF90_CREATE and NF90_OPEN statement for parallel output 27 ! To increase the performance for parallel output, the following is done: 28 ! - Limit time dimension 29 ! - Values of axis data are only written by PE0 30 ! - No fill is set for all variables 31 ! Check the number of output time levels for restart jobs 26 32 ! 27 33 ! Former revisions: … … 169 175 INTEGER :: av, cross_profiles_count, cross_profiles_maxi, delim, & 170 176 delim_old, file_id, i, id_last, id_x, id_y, id_z, j, & 171 k, kk, ns, ns_old, nz_old 177 k, kk, ns, ns_old, ntime_count, nz_old 178 179 INTEGER, SAVE :: oldmode 172 180 173 181 INTEGER, DIMENSION(1) :: id_dim_time_old, id_dim_x_yz_old, & … … 858 866 859 867 ! 860 !-- Define time coordinate for volume data (unlimited dimension) 861 nc_stat = NF90_DEF_DIM( id_set_3d(av), 'time', NF90_UNLIMITED, & 862 id_dim_time_3d(av) ) 863 CALL handle_netcdf_error( 'netcdf', 64 ) 868 !-- Define time coordinate for volume data. 869 !-- For parallel output the time dimensions has to be limited, otherwise 870 !-- the performance drops significantly. 871 IF ( netcdf_data_format < 5 ) THEN 872 nc_stat = NF90_DEF_DIM( id_set_3d(av), 'time', NF90_UNLIMITED, & 873 id_dim_time_3d(av) ) 874 CALL handle_netcdf_error( 'netcdf', 64 ) 875 ELSE 876 nc_stat = NF90_DEF_DIM( id_set_3d(av), 'time', ntdim_3d(av), & 877 id_dim_time_3d(av) ) 878 CALL handle_netcdf_error( 'netcdf', 523 ) 879 ENDIF 864 880 865 881 nc_stat = NF90_DEF_VAR( id_set_3d(av), 'time', NF90_DOUBLE, & … … 1088 1104 CALL handle_netcdf_error( 'netcdf', 357 ) 1089 1105 #if defined( __netcdf4_parallel ) 1090 !1091 !-- Set collective io operations for parallel io1092 1106 IF ( netcdf_data_format > 4 ) THEN 1107 ! 1108 !-- Set no fill for every variable to increase performance. 1109 nc_stat = NF90_DEF_VAR_FILL( id_set_3d(av), & 1110 id_var_do3d(av,i), & 1111 1, 0 ) 1112 CALL handle_netcdf_error( 'netcdf', 532 ) 1113 ! 1114 !-- Set collective io operations for parallel io 1093 1115 nc_stat = NF90_VAR_PAR_ACCESS( id_set_3d(av), & 1094 1116 id_var_do3d(av,i), & … … 1113 1135 1114 1136 ! 1137 !-- Set general no fill, otherwise the performance drops significantly for 1138 !-- parallel output. 1139 nc_stat = NF90_SET_FILL( id_set_3d(av), NF90_NOFILL, oldmode ) 1140 CALL handle_netcdf_error( 'netcdf', 528 ) 1141 1142 ! 1115 1143 !-- Leave netCDF define mode 1116 1144 nc_stat = NF90_ENDDEF( id_set_3d(av) ) … … 1118 1146 1119 1147 ! 1120 !-- Write data for x (shifted by +dx/2) and xu axis 1121 ALLOCATE( netcdf_data(0:nx+1) ) 1122 1123 DO i = 0, nx+1 1124 netcdf_data(i) = ( i + 0.5 ) * dx 1125 ENDDO 1126 1127 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_x_3d(av), netcdf_data, & 1128 start = (/ 1 /), count = (/ nx+2 /) ) 1129 CALL handle_netcdf_error( 'netcdf', 83 ) 1130 1131 DO i = 0, nx+1 1132 netcdf_data(i) = i * dx 1133 ENDDO 1134 1135 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_xu_3d(av), & 1136 netcdf_data, start = (/ 1 /), & 1137 count = (/ nx+2 /) ) 1138 CALL handle_netcdf_error( 'netcdf', 385 ) 1139 1140 DEALLOCATE( netcdf_data ) 1141 1142 ! 1143 !-- Write data for y (shifted by +dy/2) and yv axis 1144 ALLOCATE( netcdf_data(0:ny+1) ) 1145 1146 DO i = 0, ny+1 1147 netcdf_data(i) = ( i + 0.5 ) * dy 1148 ENDDO 1149 1150 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_y_3d(av), netcdf_data, & 1151 start = (/ 1 /), count = (/ ny+2 /)) 1152 CALL handle_netcdf_error( 'netcdf', 84 ) 1153 1154 DO i = 0, ny+1 1155 netcdf_data(i) = i * dy 1156 ENDDO 1157 1158 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_yv_3d(av), & 1159 netcdf_data, start = (/ 1 /), & 1160 count = (/ ny+2 /)) 1161 CALL handle_netcdf_error( 'netcdf', 387 ) 1162 1163 DEALLOCATE( netcdf_data ) 1164 1165 ! 1166 !-- Write zu and zw data (vertical axes) 1167 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zu_3d(av), & 1168 zu(nzb:nz_do3d), start = (/ 1 /), & 1169 count = (/ nz_do3d-nzb+1 /) ) 1170 CALL handle_netcdf_error( 'netcdf', 85 ) 1171 1172 1173 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zw_3d(av), & 1174 zw(nzb:nz_do3d), start = (/ 1 /), & 1175 count = (/ nz_do3d-nzb+1 /) ) 1176 CALL handle_netcdf_error( 'netcdf', 86 ) 1177 1178 1179 ! 1180 !-- In case of non-flat topography write height information 1181 IF ( TRIM( topography ) /= 'flat' ) THEN 1182 1183 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), & 1184 zu_s_inner(0:nx+1,0:ny+1), & 1185 start = (/ 1, 1 /), & 1186 count = (/ nx+2, ny+2 /) ) 1187 CALL handle_netcdf_error( 'netcdf', 419 ) 1188 1189 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), & 1190 zw_w_inner(0:nx+1,0:ny+1), & 1191 start = (/ 1, 1 /), & 1192 count = (/ nx+2, ny+2 /) ) 1193 CALL handle_netcdf_error( 'netcdf', 420 ) 1194 1148 !-- These data are only written by PE0 for parallel output to increase 1149 !-- the performance. 1150 IF ( myid == 0 .OR. netcdf_data_format < 5 ) THEN 1151 ! 1152 !-- Write data for x (shifted by +dx/2) and xu axis 1153 ALLOCATE( netcdf_data(0:nx+1) ) 1154 1155 DO i = 0, nx+1 1156 netcdf_data(i) = ( i + 0.5 ) * dx 1157 ENDDO 1158 1159 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_x_3d(av), & 1160 netcdf_data, start = (/ 1 /), & 1161 count = (/ nx+2 /) ) 1162 CALL handle_netcdf_error( 'netcdf', 83 ) 1163 1164 DO i = 0, nx+1 1165 netcdf_data(i) = i * dx 1166 ENDDO 1167 1168 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_xu_3d(av), & 1169 netcdf_data, start = (/ 1 /), & 1170 count = (/ nx+2 /) ) 1171 CALL handle_netcdf_error( 'netcdf', 385 ) 1172 1173 DEALLOCATE( netcdf_data ) 1174 1175 ! 1176 !-- Write data for y (shifted by +dy/2) and yv axis 1177 ALLOCATE( netcdf_data(0:ny+1) ) 1178 1179 DO i = 0, ny+1 1180 netcdf_data(i) = ( i + 0.5 ) * dy 1181 ENDDO 1182 1183 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_y_3d(av), & 1184 netcdf_data, start = (/ 1 /), & 1185 count = (/ ny+2 /) ) 1186 CALL handle_netcdf_error( 'netcdf', 84 ) 1187 1188 DO i = 0, ny+1 1189 netcdf_data(i) = i * dy 1190 ENDDO 1191 1192 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_yv_3d(av), & 1193 netcdf_data, start = (/ 1 /), & 1194 count = (/ ny+2 /)) 1195 CALL handle_netcdf_error( 'netcdf', 387 ) 1196 1197 DEALLOCATE( netcdf_data ) 1198 1199 ! 1200 !-- Write zu and zw data (vertical axes) 1201 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zu_3d(av), & 1202 zu(nzb:nz_do3d), start = (/ 1 /), & 1203 count = (/ nz_do3d-nzb+1 /) ) 1204 CALL handle_netcdf_error( 'netcdf', 85 ) 1205 1206 1207 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zw_3d(av), & 1208 zw(nzb:nz_do3d), start = (/ 1 /), & 1209 count = (/ nz_do3d-nzb+1 /) ) 1210 CALL handle_netcdf_error( 'netcdf', 86 ) 1211 1212 ! 1213 !-- In case of non-flat topography write height information 1214 IF ( TRIM( topography ) /= 'flat' ) THEN 1215 1216 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), & 1217 zu_s_inner(0:nx+1,0:ny+1), & 1218 start = (/ 1, 1 /), & 1219 count = (/ nx+2, ny+2 /) ) 1220 CALL handle_netcdf_error( 'netcdf', 419 ) 1221 1222 nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), & 1223 zw_w_inner(0:nx+1,0:ny+1), & 1224 start = (/ 1, 1 /), & 1225 count = (/ nx+2, ny+2 /) ) 1226 CALL handle_netcdf_error( 'netcdf', 420 ) 1227 1228 ENDIF 1195 1229 ENDIF 1196 1230 … … 1271 1305 1272 1306 nc_stat = NF90_INQUIRE_DIMENSION( id_set_3d(av), id_dim_time_3d(av), & 1273 len = do3d_time_count(av))1307 len = ntime_count ) 1274 1308 CALL handle_netcdf_error( 'netcdf', 93 ) 1309 1310 ! 1311 !-- For non-parallel output use the last output time level of the netcdf 1312 !-- file because the time dimension is unlimited. In case of parallel 1313 !-- output the variable ntime_count could get the value of 9*10E36 because 1314 !-- the time dimension is limited. 1315 IF ( netcdf_data_format < 5 ) do3d_time_count(av) = ntime_count 1275 1316 1276 1317 nc_stat = NF90_GET_VAR( id_set_3d(av), id_var_time_3d(av), & … … 1294 1335 ENDIF 1295 1336 1337 IF ( netcdf_data_format > 4 ) THEN 1338 IF ( ntdim_3d(av) > ntime_count ) THEN 1339 message_string = 'netCDF file for volume data ' // & 1340 TRIM( var ) // ' from previous run found,' // & 1341 '&but this file cannot be extended becaus' // & 1342 'e the number of output time levels&has b' // & 1343 'een increased compared to the previous s' // & 1344 'imulation.' // & 1345 '&New file is created instead.' 1346 CALL message( 'define_netcdf_header', 'PA0388', 0, 1, 0, 6, 0 ) 1347 do3d_time_count(av) = 0 1348 extend = .FALSE. 1349 RETURN 1350 ENDIF 1351 ENDIF 1352 1296 1353 ! 1297 1354 !-- Dataset seems to be extendable. … … 1367 1424 1368 1425 ! 1369 !-- Define time coordinate for xy sections (unlimited dimension) 1370 nc_stat = NF90_DEF_DIM( id_set_xy(av), 'time', NF90_UNLIMITED, & 1371 id_dim_time_xy(av) ) 1372 CALL handle_netcdf_error( 'netcdf', 99 ) 1426 !-- Define time coordinate for xy sections. 1427 !-- For parallel output the time dimensions has to be limited, otherwise 1428 !-- the performance drops significantly. 1429 IF ( netcdf_data_format < 5 ) THEN 1430 nc_stat = NF90_DEF_DIM( id_set_xy(av), 'time', NF90_UNLIMITED, & 1431 id_dim_time_xy(av) ) 1432 CALL handle_netcdf_error( 'netcdf', 99 ) 1433 ELSE 1434 nc_stat = NF90_DEF_DIM( id_set_xy(av), 'time', ntdim_2d_xy(av), & 1435 id_dim_time_xy(av) ) 1436 CALL handle_netcdf_error( 'netcdf', 524 ) 1437 ENDIF 1373 1438 1374 1439 nc_stat = NF90_DEF_VAR( id_set_xy(av), 'time', NF90_DOUBLE, & … … 1655 1720 CALL handle_netcdf_error( 'netcdf', 354 ) 1656 1721 #if defined( __netcdf4_parallel ) 1657 !1658 !-- Set collective io operations for parallel io1659 1722 IF ( netcdf_data_format > 4 ) THEN 1723 ! 1724 !-- Set no fill for every variable to increase performance. 1725 nc_stat = NF90_DEF_VAR_FILL( id_set_xy(av), & 1726 id_var_do2d(av,i), & 1727 1, 0 ) 1728 CALL handle_netcdf_error( 'netcdf', 533 ) 1729 ! 1730 !-- Set collective io operations for parallel io 1660 1731 nc_stat = NF90_VAR_PAR_ACCESS( id_set_xy(av), & 1661 1732 id_var_do2d(av,i), & … … 1682 1753 1683 1754 ! 1755 !-- Set general no fill, otherwise the performance drops significantly for 1756 !-- parallel output. 1757 nc_stat = NF90_SET_FILL( id_set_xy(av), NF90_NOFILL, oldmode ) 1758 CALL handle_netcdf_error( 'netcdf', 529 ) 1759 1760 ! 1684 1761 !-- Leave netCDF define mode 1685 1762 nc_stat = NF90_ENDDEF( id_set_xy(av) ) … … 1687 1764 1688 1765 ! 1689 !-- Write axis data: z_xy, x, y 1690 ALLOCATE( netcdf_data(1:ns) ) 1691 1692 ! 1693 !-- Write zu data 1694 DO i = 1, ns 1695 IF( section(i,1) == -1 ) THEN 1696 netcdf_data(i) = -1.0 ! section averaged along z 1697 ELSE 1698 netcdf_data(i) = zu( section(i,1) ) 1766 !-- These data are only written by PE0 for parallel output to increase 1767 !-- the performance. 1768 IF ( myid == 0 .OR. netcdf_data_format < 5 ) THEN 1769 1770 ! 1771 !-- Write axis data: z_xy, x, y 1772 ALLOCATE( netcdf_data(1:ns) ) 1773 1774 ! 1775 !-- Write zu data 1776 DO i = 1, ns 1777 IF( section(i,1) == -1 ) THEN 1778 netcdf_data(i) = -1.0 ! section averaged along z 1779 ELSE 1780 netcdf_data(i) = zu( section(i,1) ) 1781 ENDIF 1782 ENDDO 1783 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zu_xy(av), & 1784 netcdf_data, start = (/ 1 /), & 1785 count = (/ ns /) ) 1786 CALL handle_netcdf_error( 'netcdf', 123 ) 1787 1788 ! 1789 !-- Write zw data 1790 DO i = 1, ns 1791 IF( section(i,1) == -1 ) THEN 1792 netcdf_data(i) = -1.0 ! section averaged along z 1793 ELSE 1794 netcdf_data(i) = zw( section(i,1) ) 1795 ENDIF 1796 ENDDO 1797 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zw_xy(av), & 1798 netcdf_data, start = (/ 1 /), & 1799 count = (/ ns /) ) 1800 CALL handle_netcdf_error( 'netcdf', 124 ) 1801 1802 ! 1803 !-- Write gridpoint number data 1804 netcdf_data(1:ns) = section(1:ns,1) 1805 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_ind_z_xy(av), & 1806 netcdf_data, start = (/ 1 /), & 1807 count = (/ ns /) ) 1808 CALL handle_netcdf_error( 'netcdf', 125 ) 1809 1810 DEALLOCATE( netcdf_data ) 1811 1812 ! 1813 !-- Write the cross section height u*, t* 1814 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zu1_xy(av), & 1815 (/ zu(nzb+1) /), start = (/ 1 /), & 1816 count = (/ 1 /) ) 1817 CALL handle_netcdf_error( 'netcdf', 126 ) 1818 1819 ! 1820 !-- Write data for x (shifted by +dx/2) and xu axis 1821 ALLOCATE( netcdf_data(0:nx+1) ) 1822 1823 DO i = 0, nx+1 1824 netcdf_data(i) = ( i + 0.5 ) * dx 1825 ENDDO 1826 1827 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_x_xy(av), & 1828 netcdf_data, start = (/ 1 /), & 1829 count = (/ nx+2 /) ) 1830 CALL handle_netcdf_error( 'netcdf', 127 ) 1831 1832 DO i = 0, nx+1 1833 netcdf_data(i) = i * dx 1834 ENDDO 1835 1836 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_xu_xy(av), & 1837 netcdf_data, start = (/ 1 /), & 1838 count = (/ nx+2 /) ) 1839 CALL handle_netcdf_error( 'netcdf', 367 ) 1840 1841 DEALLOCATE( netcdf_data ) 1842 1843 ! 1844 !-- Write data for y (shifted by +dy/2) and yv axis 1845 ALLOCATE( netcdf_data(0:ny+1) ) 1846 1847 DO i = 0, ny+1 1848 netcdf_data(i) = ( i + 0.5 ) * dy 1849 ENDDO 1850 1851 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_y_xy(av), & 1852 netcdf_data, start = (/ 1 /), & 1853 count = (/ ny+2 /)) 1854 CALL handle_netcdf_error( 'netcdf', 128 ) 1855 1856 DO i = 0, ny+1 1857 netcdf_data(i) = i * dy 1858 ENDDO 1859 1860 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_yv_xy(av), & 1861 netcdf_data, start = (/ 1 /), & 1862 count = (/ ny+2 /)) 1863 CALL handle_netcdf_error( 'netcdf', 368 ) 1864 1865 DEALLOCATE( netcdf_data ) 1866 1867 ! 1868 !-- In case of non-flat topography write height information 1869 IF ( TRIM( topography ) /= 'flat' ) THEN 1870 1871 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), & 1872 zu_s_inner(0:nx+1,0:ny+1), & 1873 start = (/ 1, 1 /), & 1874 count = (/ nx+2, ny+2 /) ) 1875 CALL handle_netcdf_error( 'netcdf', 427 ) 1876 1877 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), & 1878 zw_w_inner(0:nx+1,0:ny+1), & 1879 start = (/ 1, 1 /), & 1880 count = (/ nx+2, ny+2 /) ) 1881 CALL handle_netcdf_error( 'netcdf', 428 ) 1882 1699 1883 ENDIF 1700 ENDDO 1701 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zu_xy(av), & 1702 netcdf_data, start = (/ 1 /), & 1703 count = (/ ns /) ) 1704 CALL handle_netcdf_error( 'netcdf', 123 ) 1705 1706 ! 1707 !-- Write zw data 1708 DO i = 1, ns 1709 IF( section(i,1) == -1 ) THEN 1710 netcdf_data(i) = -1.0 ! section averaged along z 1711 ELSE 1712 netcdf_data(i) = zw( section(i,1) ) 1713 ENDIF 1714 ENDDO 1715 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zw_xy(av), & 1716 netcdf_data, start = (/ 1 /), & 1717 count = (/ ns /) ) 1718 CALL handle_netcdf_error( 'netcdf', 124 ) 1719 1720 ! 1721 !-- Write gridpoint number data 1722 netcdf_data(1:ns) = section(1:ns,1) 1723 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_ind_z_xy(av), & 1724 netcdf_data, start = (/ 1 /), & 1725 count = (/ ns /) ) 1726 CALL handle_netcdf_error( 'netcdf', 125 ) 1727 1728 DEALLOCATE( netcdf_data ) 1729 1730 ! 1731 !-- Write the cross section height u*, t* 1732 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zu1_xy(av), & 1733 (/ zu(nzb+1) /), start = (/ 1 /), & 1734 count = (/ 1 /) ) 1735 CALL handle_netcdf_error( 'netcdf', 126 ) 1736 1737 ! 1738 !-- Write data for x (shifted by +dx/2) and xu axis 1739 ALLOCATE( netcdf_data(0:nx+1) ) 1740 1741 DO i = 0, nx+1 1742 netcdf_data(i) = ( i + 0.5 ) * dx 1743 ENDDO 1744 1745 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_x_xy(av), netcdf_data, & 1746 start = (/ 1 /), count = (/ nx+2 /) ) 1747 CALL handle_netcdf_error( 'netcdf', 127 ) 1748 1749 DO i = 0, nx+1 1750 netcdf_data(i) = i * dx 1751 ENDDO 1752 1753 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_xu_xy(av), & 1754 netcdf_data, start = (/ 1 /), & 1755 count = (/ nx+2 /) ) 1756 CALL handle_netcdf_error( 'netcdf', 367 ) 1757 1758 DEALLOCATE( netcdf_data ) 1759 1760 ! 1761 !-- Write data for y (shifted by +dy/2) and yv axis 1762 ALLOCATE( netcdf_data(0:ny+1) ) 1763 1764 DO i = 0, ny+1 1765 netcdf_data(i) = ( i + 0.5 ) * dy 1766 ENDDO 1767 1768 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_y_xy(av), netcdf_data, & 1769 start = (/ 1 /), count = (/ ny+2 /)) 1770 CALL handle_netcdf_error( 'netcdf', 128 ) 1771 1772 DO i = 0, ny+1 1773 netcdf_data(i) = i * dy 1774 ENDDO 1775 1776 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_yv_xy(av), & 1777 netcdf_data, start = (/ 1 /), & 1778 count = (/ ny+2 /)) 1779 CALL handle_netcdf_error( 'netcdf', 368 ) 1780 1781 DEALLOCATE( netcdf_data ) 1782 1783 ! 1784 !-- In case of non-flat topography write height information 1785 IF ( TRIM( topography ) /= 'flat' ) THEN 1786 1787 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), & 1788 zu_s_inner(0:nx+1,0:ny+1), & 1789 start = (/ 1, 1 /), & 1790 count = (/ nx+2, ny+2 /) ) 1791 CALL handle_netcdf_error( 'netcdf', 427 ) 1792 1793 nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), & 1794 zw_w_inner(0:nx+1,0:ny+1), & 1795 start = (/ 1, 1 /), & 1796 count = (/ nx+2, ny+2 /) ) 1797 CALL handle_netcdf_error( 'netcdf', 428 ) 1798 1799 ENDIF 1800 1884 ENDIF 1801 1885 1802 1886 CASE ( 'xy_ext' ) … … 1924 2008 1925 2009 nc_stat = NF90_INQUIRE_DIMENSION( id_set_xy(av), id_dim_time_xy(av), & 1926 len = do2d_xy_time_count(av))2010 len = ntime_count ) 1927 2011 CALL handle_netcdf_error( 'netcdf', 136 ) 2012 2013 ! 2014 !-- For non-parallel output use the last output time level of the netcdf 2015 !-- file because the time dimension is unlimited. In case of parallel 2016 !-- output the variable ntime_count could get the value of 9*10E36 because 2017 !-- the time dimension is limited. 2018 IF ( netcdf_data_format < 5 ) do2d_xy_time_count(av) = ntime_count 1928 2019 1929 2020 nc_stat = NF90_GET_VAR( id_set_xy(av), id_var_time_xy(av), & … … 1945 2036 extend = .FALSE. 1946 2037 RETURN 2038 ENDIF 2039 2040 IF ( netcdf_data_format > 4 ) THEN 2041 IF ( ntdim_2d_xy(av) > ntime_count ) THEN 2042 message_string = 'netCDF file for cross sections ' // & 2043 TRIM( var ) // ' from previous run found,' // & 2044 '&but this file cannot be extended becaus' // & 2045 'e the number of output time levels&has b' // & 2046 'een increased compared to the previous s' // & 2047 'imulation.' // & 2048 '&New file is created instead.' 2049 CALL message( 'define_netcdf_header', 'PA0389', 0, 1, 0, 6, 0 ) 2050 do2d_xy_time_count(av) = 0 2051 extend = .FALSE. 2052 RETURN 2053 ENDIF 1947 2054 ENDIF 1948 2055 … … 2023 2130 2024 2131 ! 2025 !-- Define time coordinate for xz sections (unlimited dimension) 2026 nc_stat = NF90_DEF_DIM( id_set_xz(av), 'time', NF90_UNLIMITED, & 2027 id_dim_time_xz(av) ) 2028 CALL handle_netcdf_error( 'netcdf', 142 ) 2132 !-- Define time coordinate for xz sections. 2133 !-- For parallel output the time dimensions has to be limited, otherwise 2134 !-- the performance drops significantly. 2135 IF ( netcdf_data_format < 5 ) THEN 2136 nc_stat = NF90_DEF_DIM( id_set_xz(av), 'time', NF90_UNLIMITED, & 2137 id_dim_time_xz(av) ) 2138 CALL handle_netcdf_error( 'netcdf', 142 ) 2139 ELSE 2140 nc_stat = NF90_DEF_DIM( id_set_xz(av), 'time', ntdim_2d_xz(av), & 2141 id_dim_time_xz(av) ) 2142 CALL handle_netcdf_error( 'netcdf', 525 ) 2143 ENDIF 2029 2144 2030 2145 nc_stat = NF90_DEF_VAR( id_set_xz(av), 'time', NF90_DOUBLE, & … … 2240 2355 CALL handle_netcdf_error( 'netcdf', 355 ) 2241 2356 #if defined( __netcdf4_parallel ) 2242 ! 2243 !-- Set independent io operations for parallel io. Collective io 2244 !-- is only allowed in case of a 1d-decomposition along x, because 2245 !-- otherwise, not all PEs have output data. 2357 2246 2358 IF ( netcdf_data_format > 4 ) THEN 2359 ! 2360 !-- Set no fill for every variable to increase performance. 2361 nc_stat = NF90_DEF_VAR_FILL( id_set_xz(av), & 2362 id_var_do2d(av,i), & 2363 1, 0 ) 2364 CALL handle_netcdf_error( 'netcdf', 534 ) 2365 ! 2366 !-- Set independent io operations for parallel io. Collective io 2367 !-- is only allowed in case of a 1d-decomposition along x, 2368 !-- because otherwise, not all PEs have output data. 2247 2369 IF ( npey == 1 ) THEN 2248 2370 nc_stat = NF90_VAR_PAR_ACCESS( id_set_xz(av), & … … 2251 2373 ELSE 2252 2374 ! 2253 !-- ATTENTION: Due to a probable bug in the netCDF42254 !-- installation, independet output does not work2255 !-- A workaround is used in data_output_2d on those2256 !-- PEs having no data2257 nc_stat = NF90_VAR_PAR_ACCESS( id_set_xz(av), &2258 2259 2375 !-- Test simulations showed that the output of cross sections 2376 !-- by all PEs in data_output_2d using NF90_COLLECTIVE is 2377 !-- faster than the output by the first row of PEs in 2378 !-- x-direction using NF90_INDEPENDENT. 2379 nc_stat = NF90_VAR_PAR_ACCESS( id_set_xz(av), & 2380 id_var_do2d(av,i), & 2381 NF90_COLLECTIVE ) 2260 2382 ! nc_stat = NF90_VAR_PAR_ACCESS( id_set_xz(av), & 2261 2383 ! id_var_do2d(av,i), & … … 2283 2405 2284 2406 ! 2407 !-- Set general no fill, otherwise the performance drops significantly for 2408 !-- parallel output. 2409 nc_stat = NF90_SET_FILL( id_set_xz(av), NF90_NOFILL, oldmode ) 2410 CALL handle_netcdf_error( 'netcdf', 530 ) 2411 2412 ! 2285 2413 !-- Leave netCDF define mode 2286 2414 nc_stat = NF90_ENDDEF( id_set_xz(av) ) … … 2288 2416 2289 2417 ! 2290 !-- Write axis data: y_xz, x, zu, zw 2291 ALLOCATE( netcdf_data(1:ns) ) 2292 2293 ! 2294 !-- Write y_xz data (shifted by +dy/2) 2295 DO i = 1, ns 2296 IF( section(i,2) == -1 ) THEN 2297 netcdf_data(i) = -1.0 ! section averaged along y 2298 ELSE 2299 netcdf_data(i) = ( section(i,2) + 0.5 ) * dy 2300 ENDIF 2301 ENDDO 2302 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_y_xz(av), netcdf_data, & 2303 start = (/ 1 /), count = (/ ns /) ) 2304 CALL handle_netcdf_error( 'netcdf', 163 ) 2305 2306 ! 2307 !-- Write yv_xz data 2308 DO i = 1, ns 2309 IF( section(i,2) == -1 ) THEN 2310 netcdf_data(i) = -1.0 ! section averaged along y 2311 ELSE 2312 netcdf_data(i) = section(i,2) * dy 2313 ENDIF 2314 ENDDO 2315 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_yv_xz(av), & 2316 netcdf_data, start = (/ 1 /), & 2317 count = (/ ns /) ) 2318 CALL handle_netcdf_error( 'netcdf', 375 ) 2319 2320 ! 2321 !-- Write gridpoint number data 2322 netcdf_data(1:ns) = section(1:ns,2) 2323 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_ind_y_xz(av), & 2324 netcdf_data, start = (/ 1 /), & 2325 count = (/ ns /) ) 2326 CALL handle_netcdf_error( 'netcdf', 164 ) 2327 2328 2329 DEALLOCATE( netcdf_data ) 2330 2331 ! 2332 !-- Write data for x (shifted by +dx/2) and xu axis 2333 ALLOCATE( netcdf_data(0:nx+1) ) 2334 2335 DO i = 0, nx+1 2336 netcdf_data(i) = ( i + 0.5 ) * dx 2337 ENDDO 2338 2339 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_x_xz(av), netcdf_data, & 2340 start = (/ 1 /), count = (/ nx+2 /) ) 2341 CALL handle_netcdf_error( 'netcdf', 165 ) 2342 2343 DO i = 0, nx+1 2344 netcdf_data(i) = i * dx 2345 ENDDO 2346 2347 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_xu_xz(av), & 2348 netcdf_data, start = (/ 1 /), & 2349 count = (/ nx+2 /) ) 2350 CALL handle_netcdf_error( 'netcdf', 377 ) 2351 2352 DEALLOCATE( netcdf_data ) 2353 2354 ! 2355 !-- Write zu and zw data (vertical axes) 2356 ALLOCATE( netcdf_data(0:nz+1) ) 2357 2358 netcdf_data(0:nz+1) = zu(nzb:nzt+1) 2359 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_zu_xz(av), & 2360 netcdf_data, start = (/ 1 /), & 2361 count = (/ nz+2 /) ) 2362 CALL handle_netcdf_error( 'netcdf', 166 ) 2363 2364 netcdf_data(0:nz+1) = zw(nzb:nzt+1) 2365 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_zw_xz(av), & 2366 netcdf_data, start = (/ 1 /), & 2367 count = (/ nz+2 /) ) 2368 CALL handle_netcdf_error( 'netcdf', 167 ) 2369 2370 DEALLOCATE( netcdf_data ) 2418 !-- These data are only written by PE0 for parallel output to increase 2419 !-- the performance. 2420 IF ( myid == 0 .OR. netcdf_data_format < 5 ) THEN 2421 2422 ! 2423 !-- Write axis data: y_xz, x, zu, zw 2424 ALLOCATE( netcdf_data(1:ns) ) 2425 2426 ! 2427 !-- Write y_xz data (shifted by +dy/2) 2428 DO i = 1, ns 2429 IF( section(i,2) == -1 ) THEN 2430 netcdf_data(i) = -1.0 ! section averaged along y 2431 ELSE 2432 netcdf_data(i) = ( section(i,2) + 0.5 ) * dy 2433 ENDIF 2434 ENDDO 2435 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_y_xz(av), & 2436 netcdf_data, start = (/ 1 /), & 2437 count = (/ ns /) ) 2438 CALL handle_netcdf_error( 'netcdf', 163 ) 2439 2440 ! 2441 !-- Write yv_xz data 2442 DO i = 1, ns 2443 IF( section(i,2) == -1 ) THEN 2444 netcdf_data(i) = -1.0 ! section averaged along y 2445 ELSE 2446 netcdf_data(i) = section(i,2) * dy 2447 ENDIF 2448 ENDDO 2449 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_yv_xz(av), & 2450 netcdf_data, start = (/ 1 /), & 2451 count = (/ ns /) ) 2452 CALL handle_netcdf_error( 'netcdf', 375 ) 2453 2454 ! 2455 !-- Write gridpoint number data 2456 netcdf_data(1:ns) = section(1:ns,2) 2457 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_ind_y_xz(av), & 2458 netcdf_data, start = (/ 1 /), & 2459 count = (/ ns /) ) 2460 CALL handle_netcdf_error( 'netcdf', 164 ) 2461 2462 2463 DEALLOCATE( netcdf_data ) 2464 2465 ! 2466 !-- Write data for x (shifted by +dx/2) and xu axis 2467 ALLOCATE( netcdf_data(0:nx+1) ) 2468 2469 DO i = 0, nx+1 2470 netcdf_data(i) = ( i + 0.5 ) * dx 2471 ENDDO 2472 2473 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_x_xz(av), & 2474 netcdf_data, start = (/ 1 /), & 2475 count = (/ nx+2 /) ) 2476 CALL handle_netcdf_error( 'netcdf', 165 ) 2477 2478 DO i = 0, nx+1 2479 netcdf_data(i) = i * dx 2480 ENDDO 2481 2482 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_xu_xz(av), & 2483 netcdf_data, start = (/ 1 /), & 2484 count = (/ nx+2 /) ) 2485 CALL handle_netcdf_error( 'netcdf', 377 ) 2486 2487 DEALLOCATE( netcdf_data ) 2488 2489 ! 2490 !-- Write zu and zw data (vertical axes) 2491 ALLOCATE( netcdf_data(0:nz+1) ) 2492 2493 netcdf_data(0:nz+1) = zu(nzb:nzt+1) 2494 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_zu_xz(av), & 2495 netcdf_data, start = (/ 1 /), & 2496 count = (/ nz+2 /) ) 2497 CALL handle_netcdf_error( 'netcdf', 166 ) 2498 2499 netcdf_data(0:nz+1) = zw(nzb:nzt+1) 2500 nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_zw_xz(av), & 2501 netcdf_data, start = (/ 1 /), & 2502 count = (/ nz+2 /) ) 2503 CALL handle_netcdf_error( 'netcdf', 167 ) 2504 2505 DEALLOCATE( netcdf_data ) 2506 2507 ENDIF 2371 2508 2372 2509 … … 2495 2632 2496 2633 nc_stat = NF90_INQUIRE_DIMENSION( id_set_xz(av), id_dim_time_xz(av), & 2497 len = do2d_xz_time_count(av))2634 len = ntime_count ) 2498 2635 CALL handle_netcdf_error( 'netcdf', 175 ) 2636 2637 ! 2638 !-- For non-parallel output use the last output time level of the netcdf 2639 !-- file because the time dimension is unlimited. In case of parallel 2640 !-- output the variable ntime_count could get the value of 9*10E36 because 2641 !-- the time dimension is limited. 2642 IF ( netcdf_data_format < 5 ) do2d_xz_time_count(av) = ntime_count 2499 2643 2500 2644 nc_stat = NF90_GET_VAR( id_set_xz(av), id_var_time_xz(av), & … … 2517 2661 RETURN 2518 2662 ENDIF 2663 2664 IF ( netcdf_data_format > 4 ) THEN 2665 IF ( ntdim_2d_xz(av) > ntime_count ) THEN 2666 message_string = 'netCDF file for cross sections ' // & 2667 TRIM( var ) // ' from previous run found,' // & 2668 '&but this file cannot be extended becaus' // & 2669 'e the number of output time levels&has b' // & 2670 'een increased compared to the previous s' // & 2671 'imulation.' // & 2672 '&New file is created instead.' 2673 CALL message( 'define_netcdf_header', 'PA0390', 0, 1, 0, 6, 0 ) 2674 do2d_xz_time_count(av) = 0 2675 extend = .FALSE. 2676 RETURN 2677 ENDIF 2678 ENDIF 2519 2679 2520 2680 ! … … 2539 2699 ELSE 2540 2700 ! 2541 !-- ATTENTION: Due to a probable bug in the netCDF42542 !-- installation, independet output does not work2543 !-- A workaround is used in data_output_2d on those2544 !-- PEs having no data2701 !-- Test simulations showed that the output of cross sections 2702 !-- by all PEs in data_output_2d using NF90_COLLECTIVE is 2703 !-- faster than the output by the first row of PEs in 2704 !-- x-direction using NF90_INDEPENDENT. 2545 2705 nc_stat = NF90_VAR_PAR_ACCESS( id_set_xz(av), & 2546 2706 id_var_do2d(av,i), & … … 2610 2770 2611 2771 ! 2612 !-- Define time coordinate for yz sections (unlimited dimension) 2613 nc_stat = NF90_DEF_DIM( id_set_yz(av), 'time', NF90_UNLIMITED, & 2614 id_dim_time_yz(av) ) 2615 CALL handle_netcdf_error( 'netcdf', 181 ) 2772 !-- Define time coordinate for yz sections. 2773 !-- For parallel output the time dimensions has to be limited, otherwise 2774 !-- the performance drops significantly. 2775 IF ( netcdf_data_format < 5 ) THEN 2776 nc_stat = NF90_DEF_DIM( id_set_yz(av), 'time', NF90_UNLIMITED, & 2777 id_dim_time_yz(av) ) 2778 CALL handle_netcdf_error( 'netcdf', 181 ) 2779 ELSE 2780 nc_stat = NF90_DEF_DIM( id_set_yz(av), 'time', ntdim_2d_yz(av), & 2781 id_dim_time_yz(av) ) 2782 CALL handle_netcdf_error( 'netcdf', 526 ) 2783 ENDIF 2616 2784 2617 2785 nc_stat = NF90_DEF_VAR( id_set_yz(av), 'time', NF90_DOUBLE, & … … 2827 2995 CALL handle_netcdf_error( 'netcdf', 356 ) 2828 2996 #if defined( __netcdf4_parallel ) 2829 !2830 !-- Set independent io operations for parallel io. Collective io2831 !-- is only allowed in case of a 1d-decomposition along y, because2832 !-- otherwise, not all PEs have output data.2833 2997 IF ( netcdf_data_format > 4 ) THEN 2998 ! 2999 !-- Set no fill for every variable to increase performance. 3000 nc_stat = NF90_DEF_VAR_FILL( id_set_yz(av), & 3001 id_var_do2d(av,i), & 3002 1, 0 ) 3003 CALL handle_netcdf_error( 'netcdf', 535 ) 3004 ! 3005 !-- Set independent io operations for parallel io. Collective io 3006 !-- is only allowed in case of a 1d-decomposition along y, 3007 !-- because otherwise, not all PEs have output data. 2834 3008 IF ( npex == 1 ) THEN 2835 3009 nc_stat = NF90_VAR_PAR_ACCESS( id_set_yz(av), & … … 2838 3012 ELSE 2839 3013 ! 2840 !-- ATTENTION: Due to a probable bug in the netCDF42841 !-- installation, independet output does not work2842 !-- A workaround is used in data_output_2d on those2843 !-- PEs having no data3014 !-- Test simulations showed that the output of cross sections 3015 !-- by all PEs in data_output_2d using NF90_COLLECTIVE is 3016 !-- faster than the output by the first row of PEs in 3017 !-- y-direction using NF90_INDEPENDENT. 2844 3018 nc_stat = NF90_VAR_PAR_ACCESS( id_set_yz(av), & 2845 3019 id_var_do2d(av,i), & … … 2870 3044 2871 3045 ! 3046 !-- Set general no fill, otherwise the performance drops significantly for 3047 !-- parallel output. 3048 nc_stat = NF90_SET_FILL( id_set_yz(av), NF90_NOFILL, oldmode ) 3049 CALL handle_netcdf_error( 'netcdf', 531 ) 3050 3051 ! 2872 3052 !-- Leave netCDF define mode 2873 3053 nc_stat = NF90_ENDDEF( id_set_yz(av) ) … … 2875 3055 2876 3056 ! 2877 !-- Write axis data: x_yz, y, zu, zw 2878 ALLOCATE( netcdf_data(1:ns) ) 2879 2880 ! 2881 !-- Write x_yz data (shifted by +dx/2) 2882 DO i = 1, ns 2883 IF( section(i,3) == -1 ) THEN 2884 netcdf_data(i) = -1.0 ! section averaged along x 2885 ELSE 2886 netcdf_data(i) = ( section(i,3) + 0.5 ) * dx 2887 ENDIF 2888 ENDDO 2889 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_x_yz(av), netcdf_data, & 2890 start = (/ 1 /), count = (/ ns /) ) 2891 CALL handle_netcdf_error( 'netcdf', 202 ) 2892 2893 ! 2894 !-- Write x_yz data (xu grid) 2895 DO i = 1, ns 2896 IF( section(i,3) == -1 ) THEN 2897 netcdf_data(i) = -1.0 ! section averaged along x 2898 ELSE 2899 netcdf_data(i) = section(i,3) * dx 2900 ENDIF 2901 ENDDO 2902 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_xu_yz(av), netcdf_data, & 2903 start = (/ 1 /), count = (/ ns /) ) 2904 CALL handle_netcdf_error( 'netcdf', 383 ) 2905 2906 ! 2907 !-- Write gridpoint number data 2908 netcdf_data(1:ns) = section(1:ns,3) 2909 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_ind_x_yz(av), & 2910 netcdf_data, start = (/ 1 /), & 2911 count = (/ ns /) ) 2912 CALL handle_netcdf_error( 'netcdf', 203 ) 2913 2914 DEALLOCATE( netcdf_data ) 2915 2916 ! 2917 !-- Write data for y (shifted by +dy/2) and yv axis 2918 ALLOCATE( netcdf_data(0:ny+1) ) 2919 2920 DO j = 0, ny+1 2921 netcdf_data(j) = ( j + 0.5 ) * dy 2922 ENDDO 2923 2924 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_y_yz(av), netcdf_data, & 2925 start = (/ 1 /), count = (/ ny+2 /) ) 2926 CALL handle_netcdf_error( 'netcdf', 204 ) 2927 2928 DO j = 0, ny+1 2929 netcdf_data(j) = j * dy 2930 ENDDO 2931 2932 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_yv_yz(av), & 2933 netcdf_data, start = (/ 1 /), & 2934 count = (/ ny+2 /) ) 2935 CALL handle_netcdf_error( 'netcdf', 384 ) 2936 2937 DEALLOCATE( netcdf_data ) 2938 2939 ! 2940 !-- Write zu and zw data (vertical axes) 2941 ALLOCATE( netcdf_data(0:nz+1) ) 2942 2943 netcdf_data(0:nz+1) = zu(nzb:nzt+1) 2944 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_zu_yz(av), & 2945 netcdf_data, start = (/ 1 /), & 2946 count = (/ nz+2 /) ) 2947 CALL handle_netcdf_error( 'netcdf', 205 ) 2948 2949 netcdf_data(0:nz+1) = zw(nzb:nzt+1) 2950 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_zw_yz(av), & 2951 netcdf_data, start = (/ 1 /), & 2952 count = (/ nz+2 /) ) 2953 CALL handle_netcdf_error( 'netcdf', 206 ) 2954 2955 DEALLOCATE( netcdf_data ) 3057 !-- These data are only written by PE0 for parallel output to increase 3058 !-- the performance. 3059 IF ( myid == 0 .OR. netcdf_data_format < 5 ) THEN 3060 3061 ! 3062 !-- Write axis data: x_yz, y, zu, zw 3063 ALLOCATE( netcdf_data(1:ns) ) 3064 3065 ! 3066 !-- Write x_yz data (shifted by +dx/2) 3067 DO i = 1, ns 3068 IF( section(i,3) == -1 ) THEN 3069 netcdf_data(i) = -1.0 ! section averaged along x 3070 ELSE 3071 netcdf_data(i) = ( section(i,3) + 0.5 ) * dx 3072 ENDIF 3073 ENDDO 3074 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_x_yz(av), & 3075 netcdf_data, start = (/ 1 /), & 3076 count = (/ ns /) ) 3077 CALL handle_netcdf_error( 'netcdf', 202 ) 3078 3079 ! 3080 !-- Write x_yz data (xu grid) 3081 DO i = 1, ns 3082 IF( section(i,3) == -1 ) THEN 3083 netcdf_data(i) = -1.0 ! section averaged along x 3084 ELSE 3085 netcdf_data(i) = section(i,3) * dx 3086 ENDIF 3087 ENDDO 3088 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_xu_yz(av), & 3089 netcdf_data, start = (/ 1 /), & 3090 count = (/ ns /) ) 3091 CALL handle_netcdf_error( 'netcdf', 383 ) 3092 3093 ! 3094 !-- Write gridpoint number data 3095 netcdf_data(1:ns) = section(1:ns,3) 3096 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_ind_x_yz(av), & 3097 netcdf_data, start = (/ 1 /), & 3098 count = (/ ns /) ) 3099 CALL handle_netcdf_error( 'netcdf', 203 ) 3100 3101 DEALLOCATE( netcdf_data ) 3102 3103 ! 3104 !-- Write data for y (shifted by +dy/2) and yv axis 3105 ALLOCATE( netcdf_data(0:ny+1) ) 3106 3107 DO j = 0, ny+1 3108 netcdf_data(j) = ( j + 0.5 ) * dy 3109 ENDDO 3110 3111 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_y_yz(av), & 3112 netcdf_data, start = (/ 1 /), & 3113 count = (/ ny+2 /) ) 3114 CALL handle_netcdf_error( 'netcdf', 204 ) 3115 3116 DO j = 0, ny+1 3117 netcdf_data(j) = j * dy 3118 ENDDO 3119 3120 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_yv_yz(av), & 3121 netcdf_data, start = (/ 1 /), & 3122 count = (/ ny+2 /) ) 3123 CALL handle_netcdf_error( 'netcdf', 384 ) 3124 3125 DEALLOCATE( netcdf_data ) 3126 3127 ! 3128 !-- Write zu and zw data (vertical axes) 3129 ALLOCATE( netcdf_data(0:nz+1) ) 3130 3131 netcdf_data(0:nz+1) = zu(nzb:nzt+1) 3132 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_zu_yz(av), & 3133 netcdf_data, start = (/ 1 /), & 3134 count = (/ nz+2 /) ) 3135 CALL handle_netcdf_error( 'netcdf', 205 ) 3136 3137 netcdf_data(0:nz+1) = zw(nzb:nzt+1) 3138 nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_zw_yz(av), & 3139 netcdf_data, start = (/ 1 /), & 3140 count = (/ nz+2 /) ) 3141 CALL handle_netcdf_error( 'netcdf', 206 ) 3142 3143 DEALLOCATE( netcdf_data ) 3144 3145 ENDIF 2956 3146 2957 3147 … … 3080 3270 3081 3271 nc_stat = NF90_INQUIRE_DIMENSION( id_set_yz(av), id_dim_time_yz(av), & 3082 len = do2d_yz_time_count(av))3272 len = ntime_count ) 3083 3273 CALL handle_netcdf_error( 'netcdf', 214 ) 3274 3275 ! 3276 !-- For non-parallel output use the last output time level of the netcdf 3277 !-- file because the time dimension is unlimited. In case of parallel 3278 !-- output the variable ntime_count could get the value of 9*10E36 because 3279 !-- the time dimension is limited. 3280 IF ( netcdf_data_format < 5 ) do2d_yz_time_count(av) = ntime_count 3084 3281 3085 3282 nc_stat = NF90_GET_VAR( id_set_yz(av), id_var_time_yz(av), & … … 3103 3300 ENDIF 3104 3301 3302 IF ( netcdf_data_format > 4 ) THEN 3303 IF ( ntdim_2d_yz(av) > ntime_count ) THEN 3304 message_string = 'netCDF file for cross sections ' // & 3305 TRIM( var ) // ' from previous run found,' // & 3306 '&but this file cannot be extended becaus' // & 3307 'e the number of output time levels&has b' // & 3308 'een increased compared to the previous s' // & 3309 'imulation.' // & 3310 '&New file is created instead.' 3311 CALL message( 'define_netcdf_header', 'PA0391', 0, 1, 0, 6, 0 ) 3312 do2d_yz_time_count(av) = 0 3313 extend = .FALSE. 3314 RETURN 3315 ENDIF 3316 ENDIF 3317 3105 3318 ! 3106 3319 !-- Dataset seems to be extendable. … … 3124 3337 ELSE 3125 3338 ! 3126 !-- ATTENTION: Due to a probable bug in the netCDF43127 !-- installation, independet output does not work3128 !-- A workaround is used in data_output_2d on those3129 !-- PEs having no data3339 !-- Test simulations showed that the output of cross sections 3340 !-- by all PEs in data_output_2d using NF90_COLLECTIVE is 3341 !-- faster than the output by the first row of PEs in 3342 !-- y-direction using NF90_INDEPENDENT. 3130 3343 nc_stat = NF90_VAR_PAR_ACCESS( id_set_yz(av), & 3131 3344 id_var_do2d(av,i), & … … 4651 4864 !-- 64bit-offset format 4652 4865 nc_stat = NF90_CREATE( filename, & 4653 OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ), id )4866 IOR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ), id ) 4654 4867 4655 4868 #if defined( __netcdf4 ) … … 4658 4871 ! 4659 4872 !-- netCDF4/HDF5 format 4660 nc_stat = NF90_CREATE( filename, OR( NF90_NOCLOBBER, NF90_NETCDF4 ), id )4873 nc_stat = NF90_CREATE( filename, IOR( NF90_NOCLOBBER, NF90_NETCDF4 ), id ) 4661 4874 4662 4875 ELSEIF ( netcdf_data_format == 4 .OR. & … … 4665 4878 !-- netCDF4/HDF5 format with classic model flag 4666 4879 nc_stat = NF90_CREATE( filename, & 4667 OR( NF90_NOCLOBBER, &4668 OR( NF90_CLASSIC_MODEL, NF90_HDF5 ) ), id )4880 IOR( NF90_NOCLOBBER, & 4881 IOR( NF90_CLASSIC_MODEL, NF90_HDF5 ) ), id ) 4669 4882 4670 4883 #if defined( __netcdf4_parallel ) … … 4672 4885 ! 4673 4886 !-- netCDF4/HDF5 format, parallel 4674 nc_stat = NF90_CREATE( filename, OR( NF90_NOCLOBBER, NF90_NETCDF4 ), & 4887 nc_stat = NF90_CREATE( filename, & 4888 IOR( NF90_NOCLOBBER, & 4889 IOR( NF90_NETCDF4, NF90_MPIIO ) ), & 4675 4890 id, COMM = comm2d, INFO = MPI_INFO_NULL ) 4676 4891 … … 4679 4894 !-- netCDF4/HDF5 format with classic model flag, parallel 4680 4895 nc_stat = NF90_CREATE( filename, & 4681 OR( NF90_NOCLOBBER, & 4682 OR( NF90_CLASSIC_MODEL, NF90_HDF5 ) ), & 4896 IOR( NF90_NOCLOBBER, & 4897 IOR( NF90_MPIIO, & 4898 IOR( NF90_CLASSIC_MODEL, NF90_HDF5 ) ) ), & 4683 4899 id, COMM = comm2d, INFO = MPI_INFO_NULL ) 4684 4900 … … 4720 4936 #if defined( __netcdf4_parallel ) 4721 4937 ELSEIF ( netcdf_data_format > 4 .AND. parallel ) THEN 4722 nc_stat = NF90_OPEN( filename, NF90_WRITE, id, COMM = comm2d,&4723 INFO = MPI_INFO_NULL )4938 nc_stat = NF90_OPEN( filename, IOR( NF90_WRITE, NF90_MPIIO ), id, & 4939 COMM = comm2d, INFO = MPI_INFO_NULL ) 4724 4940 #endif 4725 4941 #endif
Note: See TracChangeset
for help on using the changeset viewer.
