Changeset ee13a7b in flexpart.git
- Timestamp:
- Apr 10, 2017, 1:55:02 AM (7 years ago)
- Branches:
- fp9.3.1-20161214-nc4
- Children:
- 030e3c1
- Parents:
- 8413bc9
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
flexpart_code/fpmetbinary_mod.F90
r8413bc9 ree13a7b 309 309 INTEGER :: nxmax_dimid, nymax_dimid, nzmax_dimid, nuvzmax_dimid, nwzmax_dimid, & 310 310 & maxspec_dimid, numclass_dimid, maxnests_dimid, nxmaxn_dimid, nymaxn_dimid, & 311 & zero_to_nzmax_dimid 311 & zero_to_nzmax_dimid, zero_to_maxnests_dimid, polemap_dimid 312 312 313 313 … … 326 326 & maxspec_dimname, numclass_dimname,& 327 327 & maxnests_dimname, nxmaxn_dimname, nymaxn_dimname, & 328 & zero_to_nzmax_dimname 328 & zero_to_nzmax_dimname, zero_to_maxnests_dimname, & 329 & polemap_dimname 329 330 330 331 ! These are temporary variables, used in the LOAD option, for … … 368 369 ncret = nf90_def_dim(ncid, 'zero_to_nzmax', nzmax+1, zero_to_nzmax_dimid) 369 370 call handle_nf90_err(ncret) 370 371 ncret = nf90_def_dim(ncid, 'zero_to_maxnests', maxnests+1, zero_to_maxnests_dimid) 372 call handle_nf90_err(ncret) 373 ncret = nf90_def_dim(ncid, 'polemap_dim', 9, polemap_dimid) 374 call handle_nf90_err(ncret) 371 375 372 376 ! Scalar values … … 1492 1496 WRITE(iounit) vdepn(:,:,:,cm_index,:) 1493 1497 1494 dim3dids = (/nxmax _dimid, nymax_dimid, maxnests_dimid/)1498 dim3dids = (/nxmaxn_dimid, nymaxn_dimid, maxnests_dimid/) 1495 1499 1496 1500 ncret = nf90_def_var(ncid, 'psn', NF90_FLOAT, & … … 1506 1510 call handle_nf90_err(ncret) 1507 1511 1508 1512 ncret = nf90_def_var(ncid, 'sdn', NF90_FLOAT, & 1513 & dim3dids, ncvarid) 1514 call handle_nf90_err(ncret) 1515 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1516 & shuffle=0, & 1517 & deflate=1, & 1518 & deflate_level=DEF_LEVEL) 1519 call handle_nf90_err(ncret) 1520 ncret = nf90_put_var(ncid, ncvarid, & 1521 & sdn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1522 call handle_nf90_err(ncret) 1523 1524 ncret = nf90_def_var(ncid, 'msln', NF90_FLOAT, & 1525 & dim3dids, ncvarid) 1526 call handle_nf90_err(ncret) 1527 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1528 & shuffle=0, & 1529 & deflate=1, & 1530 & deflate_level=DEF_LEVEL) 1531 call handle_nf90_err(ncret) 1532 ncret = nf90_put_var(ncid, ncvarid, & 1533 & msln(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1534 call handle_nf90_err(ncret) 1535 1536 ncret = nf90_def_var(ncid, 'tccn', NF90_FLOAT, & 1537 & dim3dids, ncvarid) 1538 call handle_nf90_err(ncret) 1539 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1540 & shuffle=0, & 1541 & deflate=1, & 1542 & deflate_level=DEF_LEVEL) 1543 call handle_nf90_err(ncret) 1544 ncret = nf90_put_var(ncid, ncvarid, & 1545 & tccn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1546 call handle_nf90_err(ncret) 1547 1548 ncret = nf90_def_var(ncid, 'u10n', NF90_FLOAT, & 1549 & dim3dids, ncvarid) 1550 call handle_nf90_err(ncret) 1551 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1552 & shuffle=0, & 1553 & deflate=1, & 1554 & deflate_level=DEF_LEVEL) 1555 call handle_nf90_err(ncret) 1556 ncret = nf90_put_var(ncid, ncvarid, & 1557 & u10n(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1558 call handle_nf90_err(ncret) 1559 1560 ncret = nf90_def_var(ncid, 'v10n', NF90_FLOAT, & 1561 & dim3dids, ncvarid) 1562 call handle_nf90_err(ncret) 1563 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1564 & shuffle=0, & 1565 & deflate=1, & 1566 & deflate_level=DEF_LEVEL) 1567 call handle_nf90_err(ncret) 1568 ncret = nf90_put_var(ncid, ncvarid, & 1569 & v10n(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1570 call handle_nf90_err(ncret) 1571 1572 ncret = nf90_def_var(ncid, 'tt2n', NF90_FLOAT, & 1573 & dim3dids, ncvarid) 1574 call handle_nf90_err(ncret) 1575 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1576 & shuffle=0, & 1577 & deflate=1, & 1578 & deflate_level=DEF_LEVEL) 1579 call handle_nf90_err(ncret) 1580 ncret = nf90_put_var(ncid, ncvarid, & 1581 & tt2n(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1582 call handle_nf90_err(ncret) 1583 1584 ncret = nf90_def_var(ncid, 'td2n', NF90_FLOAT, & 1585 & dim3dids, ncvarid) 1586 call handle_nf90_err(ncret) 1587 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1588 & shuffle=0, & 1589 & deflate=1, & 1590 & deflate_level=DEF_LEVEL) 1591 call handle_nf90_err(ncret) 1592 ncret = nf90_put_var(ncid, ncvarid, & 1593 & td2n(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1594 call handle_nf90_err(ncret) 1595 1596 ncret = nf90_def_var(ncid, 'lsprecn', NF90_FLOAT, & 1597 & dim3dids, ncvarid) 1598 call handle_nf90_err(ncret) 1599 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1600 & shuffle=0, & 1601 & deflate=1, & 1602 & deflate_level=DEF_LEVEL) 1603 call handle_nf90_err(ncret) 1604 ncret = nf90_put_var(ncid, ncvarid, & 1605 & lsprecn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1606 call handle_nf90_err(ncret) 1607 1608 ncret = nf90_def_var(ncid, 'convprecn', NF90_FLOAT, & 1609 & dim3dids, ncvarid) 1610 call handle_nf90_err(ncret) 1611 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1612 & shuffle=0, & 1613 & deflate=1, & 1614 & deflate_level=DEF_LEVEL) 1615 call handle_nf90_err(ncret) 1616 ncret = nf90_put_var(ncid, ncvarid, & 1617 & convprecn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1618 call handle_nf90_err(ncret) 1619 1620 ncret = nf90_def_var(ncid, 'sshfn', NF90_FLOAT, & 1621 & dim3dids, ncvarid) 1622 call handle_nf90_err(ncret) 1623 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1624 & shuffle=0, & 1625 & deflate=1, & 1626 & deflate_level=DEF_LEVEL) 1627 call handle_nf90_err(ncret) 1628 ncret = nf90_put_var(ncid, ncvarid, & 1629 & sshfn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1630 call handle_nf90_err(ncret) 1631 1632 ncret = nf90_def_var(ncid, 'ssrn', NF90_FLOAT, & 1633 & dim3dids, ncvarid) 1634 call handle_nf90_err(ncret) 1635 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1636 & shuffle=0, & 1637 & deflate=1, & 1638 & deflate_level=DEF_LEVEL) 1639 call handle_nf90_err(ncret) 1640 ncret = nf90_put_var(ncid, ncvarid, & 1641 & ssrn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1642 call handle_nf90_err(ncret) 1643 1644 ncret = nf90_def_var(ncid, 'surfstrn', NF90_FLOAT, & 1645 & dim3dids, ncvarid) 1646 call handle_nf90_err(ncret) 1647 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1648 & shuffle=0, & 1649 & deflate=1, & 1650 & deflate_level=DEF_LEVEL) 1651 call handle_nf90_err(ncret) 1652 ncret = nf90_put_var(ncid, ncvarid, & 1653 & surfstrn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1654 call handle_nf90_err(ncret) 1655 1656 ncret = nf90_def_var(ncid, 'ustarn', NF90_FLOAT, & 1657 & dim3dids, ncvarid) 1658 call handle_nf90_err(ncret) 1659 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1660 & shuffle=0, & 1661 & deflate=1, & 1662 & deflate_level=DEF_LEVEL) 1663 call handle_nf90_err(ncret) 1664 ncret = nf90_put_var(ncid, ncvarid, & 1665 & ustarn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1666 call handle_nf90_err(ncret) 1667 1668 ncret = nf90_def_var(ncid, 'wstarn', NF90_FLOAT, & 1669 & dim3dids, ncvarid) 1670 call handle_nf90_err(ncret) 1671 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1672 & shuffle=0, & 1673 & deflate=1, & 1674 & deflate_level=DEF_LEVEL) 1675 call handle_nf90_err(ncret) 1676 ncret = nf90_put_var(ncid, ncvarid, & 1677 & wstarn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1678 call handle_nf90_err(ncret) 1679 1680 ncret = nf90_def_var(ncid, 'hmixn', NF90_FLOAT, & 1681 & dim3dids, ncvarid) 1682 call handle_nf90_err(ncret) 1683 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1684 & shuffle=0, & 1685 & deflate=1, & 1686 & deflate_level=DEF_LEVEL) 1687 call handle_nf90_err(ncret) 1688 ncret = nf90_put_var(ncid, ncvarid, & 1689 & hmixn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1690 call handle_nf90_err(ncret) 1691 1692 ncret = nf90_def_var(ncid, 'tropopausen', NF90_FLOAT, & 1693 & dim3dids, ncvarid) 1694 call handle_nf90_err(ncret) 1695 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1696 & shuffle=0, & 1697 & deflate=1, & 1698 & deflate_level=DEF_LEVEL) 1699 call handle_nf90_err(ncret) 1700 ncret = nf90_put_var(ncid, ncvarid, & 1701 & tropopausen(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1702 call handle_nf90_err(ncret) 1703 1704 ncret = nf90_def_var(ncid, 'olin', NF90_FLOAT, & 1705 & dim3dids, ncvarid) 1706 call handle_nf90_err(ncret) 1707 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1708 & shuffle=0, & 1709 & deflate=1, & 1710 & deflate_level=DEF_LEVEL) 1711 call handle_nf90_err(ncret) 1712 ncret = nf90_put_var(ncid, ncvarid, & 1713 & olin(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1714 call handle_nf90_err(ncret) 1715 1716 ncret = nf90_def_var(ncid, 'diffkn', NF90_FLOAT, & 1717 & dim3dids, ncvarid) 1718 call handle_nf90_err(ncret) 1719 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1720 & shuffle=0, & 1721 & deflate=1, & 1722 & deflate_level=DEF_LEVEL) 1723 call handle_nf90_err(ncret) 1724 ncret = nf90_put_var(ncid, ncvarid, & 1725 & diffkn(0:nxmaxn-1, 0:nymaxn-1, 1, cm_index, 1:maxnests)) 1726 call handle_nf90_err(ncret) 1727 1728 dim4dids = (/nxmaxn_dimid, nymaxn_dimid, maxspec_dimid, maxnests_dimid/) 1729 1730 1731 1732 ncret = nf90_def_var(ncid, 'vdepn', NF90_FLOAT, & 1733 & dim4dids, ncvarid) 1734 call handle_nf90_err(ncret) 1735 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1736 & shuffle=0, & 1737 & deflate=1, & 1738 & deflate_level=DEF_LEVEL) 1739 call handle_nf90_err(ncret) 1740 ncret = nf90_put_var(ncid, ncvarid, & 1741 & vdepn(0:nxmaxn-1, 0:nymaxn-1, 1:maxspec, cm_index, 1:maxnests)) 1742 call handle_nf90_err(ncret) 1509 1743 1510 1744 … … 1512 1746 1513 1747 PRINT *, 'SUM(psn): ', SUM(psn(:,:,:,cm_index,:)) 1514 1515 1748 PRINT *, 'SUM(surfstrn): ', SUM(surfstrn(:,:,:,cm_index,:)) 1749 PRINT *, 'SUM(vdepn): ', SUM(vdepn(:,:,:,cm_index,:)) 1516 1750 1517 1751 … … 1527 1761 WRITE(iounit) yrn(:) 1528 1762 1763 dim1dids = (/zero_to_maxnests_dimid/) 1764 1765 ncret = nf90_def_var(ncid, 'xresoln', NF90_FLOAT, & 1766 & dim1dids, ncvarid) 1767 call handle_nf90_err(ncret) 1768 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1769 & shuffle=0, & 1770 & deflate=1, & 1771 & deflate_level=DEF_LEVEL) 1772 call handle_nf90_err(ncret) 1773 ncret = nf90_put_var(ncid, ncvarid, & 1774 & xresoln(0:maxnests)) 1775 call handle_nf90_err(ncret) 1776 1777 ncret = nf90_def_var(ncid, 'yresoln', NF90_FLOAT, & 1778 & dim1dids, ncvarid) 1779 call handle_nf90_err(ncret) 1780 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1781 & shuffle=0, & 1782 & deflate=1, & 1783 & deflate_level=DEF_LEVEL) 1784 call handle_nf90_err(ncret) 1785 ncret = nf90_put_var(ncid, ncvarid, & 1786 & yresoln(0:maxnests)) 1787 call handle_nf90_err(ncret) 1788 1789 dim1dids = (/maxnests_dimid/) 1790 1791 ncret = nf90_def_var(ncid, 'xln', NF90_FLOAT, & 1792 & dim1dids, ncvarid) 1793 call handle_nf90_err(ncret) 1794 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1795 & shuffle=0, & 1796 & deflate=1, & 1797 & deflate_level=DEF_LEVEL) 1798 call handle_nf90_err(ncret) 1799 ncret = nf90_put_var(ncid, ncvarid, & 1800 & xln(1:maxnests)) 1801 call handle_nf90_err(ncret) 1802 1803 ncret = nf90_def_var(ncid, 'yln', NF90_FLOAT, & 1804 & dim1dids, ncvarid) 1805 call handle_nf90_err(ncret) 1806 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1807 & shuffle=0, & 1808 & deflate=1, & 1809 & deflate_level=DEF_LEVEL) 1810 call handle_nf90_err(ncret) 1811 ncret = nf90_put_var(ncid, ncvarid, & 1812 & yln(1:maxnests)) 1813 call handle_nf90_err(ncret) 1814 1815 ncret = nf90_def_var(ncid, 'xrn', NF90_FLOAT, & 1816 & dim1dids, ncvarid) 1817 call handle_nf90_err(ncret) 1818 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1819 & shuffle=0, & 1820 & deflate=1, & 1821 & deflate_level=DEF_LEVEL) 1822 call handle_nf90_err(ncret) 1823 ncret = nf90_put_var(ncid, ncvarid, & 1824 & xrn(1:maxnests)) 1825 call handle_nf90_err(ncret) 1826 1827 ncret = nf90_def_var(ncid, 'yrn', NF90_FLOAT, & 1828 & dim1dids, ncvarid) 1829 call handle_nf90_err(ncret) 1830 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1831 & shuffle=0, & 1832 & deflate=1, & 1833 & deflate_level=DEF_LEVEL) 1834 call handle_nf90_err(ncret) 1835 ncret = nf90_put_var(ncid, ncvarid, & 1836 & yrn(1:maxnests)) 1837 call handle_nf90_err(ncret) 1838 1839 PRINT *, 'SUM(yresoln): ', SUM(yresoln) 1840 PRINT *, 'SUM(xrn): ', SUM(xrn) 1841 1842 1843 1529 1844 ! Variables for polar stereographic projection 1530 1845 WRITE(iounit) xglobal, sglobal, nglobal … … 1532 1847 WRITE(iounit) southpolemap(:) 1533 1848 WRITE(iounit) northpolemap(:) 1849 1850 1851 1852 dim1dids = (/polemap_dimid/) 1853 1854 ncret = nf90_def_var(ncid, 'southpolemap', NF90_FLOAT, & 1855 & dim1dids, ncvarid) 1856 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1857 & shuffle=0, & 1858 & deflate=1, & 1859 & deflate_level=DEF_LEVEL) 1860 ncret = nf90_put_var(ncid, ncvarid, & 1861 & southpolemap(:)) 1862 1863 ncret = nf90_def_var(ncid, 'northpolemap', NF90_FLOAT, & 1864 & dim1dids, ncvarid) 1865 ncret = nf90_def_var_deflate(ncid, ncvarid, & 1866 & shuffle=0, & 1867 & deflate=1, & 1868 & deflate_level=DEF_LEVEL) 1869 ncret = nf90_put_var(ncid, ncvarid, & 1870 & northpolemap(:)) 1871 1872 1873 1874 1875 ncret = nf90_def_var(ncid, 'xglobal', NF90_INT, ncvarid) 1876 call handle_nf90_err(ncret) 1877 ncret = nf90_put_var(ncid, ncvarid, xglobal) 1878 call handle_nf90_err(ncret) 1879 1880 ncret = nf90_def_var(ncid, 'sglobal', NF90_INT, ncvarid) 1881 call handle_nf90_err(ncret) 1882 ncret = nf90_put_var(ncid, ncvarid, sglobal) 1883 call handle_nf90_err(ncret) 1884 1885 ncret = nf90_def_var(ncid, 'nglobal', NF90_INT, ncvarid) 1886 call handle_nf90_err(ncret) 1887 ncret = nf90_put_var(ncid, ncvarid, nglobal) 1888 call handle_nf90_err(ncret) 1889 1890 ncret = nf90_def_var(ncid, 'switchnorthg', NF90_FLOAT, ncvarid) 1891 call handle_nf90_err(ncret) 1892 ncret = nf90_put_var(ncid, ncvarid, switchnorthg) 1893 call handle_nf90_err(ncret) 1894 1895 ncret = nf90_def_var(ncid, 'switchsouthg', NF90_FLOAT, ncvarid) 1896 call handle_nf90_err(ncret) 1897 ncret = nf90_put_var(ncid, ncvarid, switchsouthg) 1898 call handle_nf90_err(ncret) 1899 1900 1901 1902 PRINT *, 'SUM(northpolemap): ', SUM(northpolemap) 1903 PRINT *, 'xglobal: ', xglobal 1904 PRINT *, 'sglobal: ', sglobal 1905 PRINT *, 'nglobal: ', nglobal 1906 PRINT *, 'switchsouthg: ', switchsouthg 1534 1907 1535 1908 ! Variables declared in conv_mod (convection) … … 2305 2678 call handle_nf90_err(ncret) 2306 2679 2680 ncret = nf90_inq_varid(ncid, 'sdn', ncvarid) 2681 call handle_nf90_err(ncret) 2682 ncret = nf90_get_var(ncid, ncvarid, sdn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2683 call handle_nf90_err(ncret) 2684 2685 ncret = nf90_inq_varid(ncid, 'msln', ncvarid) 2686 call handle_nf90_err(ncret) 2687 ncret = nf90_get_var(ncid, ncvarid, msln(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2688 call handle_nf90_err(ncret) 2689 2690 ncret = nf90_inq_varid(ncid, 'tccn', ncvarid) 2691 call handle_nf90_err(ncret) 2692 ncret = nf90_get_var(ncid, ncvarid, tccn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2693 call handle_nf90_err(ncret) 2694 2695 ncret = nf90_inq_varid(ncid, 'u10n', ncvarid) 2696 call handle_nf90_err(ncret) 2697 ncret = nf90_get_var(ncid, ncvarid, u10n(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2698 call handle_nf90_err(ncret) 2699 2700 ncret = nf90_inq_varid(ncid, 'v10n', ncvarid) 2701 call handle_nf90_err(ncret) 2702 ncret = nf90_get_var(ncid, ncvarid, v10n(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2703 call handle_nf90_err(ncret) 2704 2705 ncret = nf90_inq_varid(ncid, 'tt2n', ncvarid) 2706 call handle_nf90_err(ncret) 2707 ncret = nf90_get_var(ncid, ncvarid, tt2n(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2708 call handle_nf90_err(ncret) 2709 2710 ncret = nf90_inq_varid(ncid, 'td2n', ncvarid) 2711 call handle_nf90_err(ncret) 2712 ncret = nf90_get_var(ncid, ncvarid, td2n(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2713 call handle_nf90_err(ncret) 2714 2715 ncret = nf90_inq_varid(ncid, 'lsprecn', ncvarid) 2716 call handle_nf90_err(ncret) 2717 ncret = nf90_get_var(ncid, ncvarid, lsprecn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2718 call handle_nf90_err(ncret) 2719 2720 ncret = nf90_inq_varid(ncid, 'convprecn', ncvarid) 2721 call handle_nf90_err(ncret) 2722 ncret = nf90_get_var(ncid, ncvarid, convprecn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2723 call handle_nf90_err(ncret) 2724 2725 ncret = nf90_inq_varid(ncid, 'sshfn', ncvarid) 2726 call handle_nf90_err(ncret) 2727 ncret = nf90_get_var(ncid, ncvarid, sshfn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2728 call handle_nf90_err(ncret) 2729 2730 ncret = nf90_inq_varid(ncid, 'ssrn', ncvarid) 2731 call handle_nf90_err(ncret) 2732 ncret = nf90_get_var(ncid, ncvarid, ssrn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2733 call handle_nf90_err(ncret) 2734 2735 ncret = nf90_inq_varid(ncid, 'surfstrn', ncvarid) 2736 call handle_nf90_err(ncret) 2737 ncret = nf90_get_var(ncid, ncvarid, surfstrn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2738 call handle_nf90_err(ncret) 2739 2740 ncret = nf90_inq_varid(ncid, 'ustarn', ncvarid) 2741 call handle_nf90_err(ncret) 2742 ncret = nf90_get_var(ncid, ncvarid, ustarn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2743 call handle_nf90_err(ncret) 2744 2745 ncret = nf90_inq_varid(ncid, 'wstarn', ncvarid) 2746 call handle_nf90_err(ncret) 2747 ncret = nf90_get_var(ncid, ncvarid, wstarn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2748 call handle_nf90_err(ncret) 2749 2750 ncret = nf90_inq_varid(ncid, 'hmixn', ncvarid) 2751 call handle_nf90_err(ncret) 2752 ncret = nf90_get_var(ncid, ncvarid, hmixn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2753 call handle_nf90_err(ncret) 2754 2755 ncret = nf90_inq_varid(ncid, 'tropopausen', ncvarid) 2756 call handle_nf90_err(ncret) 2757 ncret = nf90_get_var(ncid, ncvarid, tropopausen(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2758 call handle_nf90_err(ncret) 2759 2760 ncret = nf90_inq_varid(ncid, 'olin', ncvarid) 2761 call handle_nf90_err(ncret) 2762 ncret = nf90_get_var(ncid, ncvarid, olin(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2763 call handle_nf90_err(ncret) 2764 2765 ncret = nf90_inq_varid(ncid, 'diffkn', ncvarid) 2766 call handle_nf90_err(ncret) 2767 ncret = nf90_get_var(ncid, ncvarid, diffkn(0:nxmaxn-1,0:nymaxn-1,1,cm_index,1:maxnests)) 2768 call handle_nf90_err(ncret) 2769 2770 ncret = nf90_inq_varid(ncid, 'vdepn', ncvarid) 2771 call handle_nf90_err(ncret) 2772 ncret = nf90_get_var(ncid, ncvarid, vdepn(0:nxmaxn-1,0:nymaxn-1,1:maxspec,cm_index,1:maxnests)) 2773 call handle_nf90_err(ncret) 2774 2775 2307 2776 2308 2777 2309 2778 PRINT *, 'SUM(psn): ', SUM(psn(:,:,:,cm_index,:)) 2779 PRINT *, 'SUM(surfstrn): ', SUM(surfstrn(:,:,:,cm_index,:)) 2780 PRINT *, 'SUM(vdepn): ', SUM(vdepn(:,:,:,cm_index,:)) 2310 2781 2311 2782 … … 2319 2790 READ(iounit) yrn(:) 2320 2791 2792 ncret = nf90_inq_varid(ncid, 'xresoln', ncvarid) 2793 call handle_nf90_err(ncret) 2794 ncret = nf90_get_var(ncid, ncvarid, xresoln(0:maxnests)) 2795 call handle_nf90_err(ncret) 2796 2797 ncret = nf90_inq_varid(ncid, 'yresoln', ncvarid) 2798 call handle_nf90_err(ncret) 2799 ncret = nf90_get_var(ncid, ncvarid, yresoln(0:maxnests)) 2800 call handle_nf90_err(ncret) 2801 2802 ncret = nf90_inq_varid(ncid, 'xln', ncvarid) 2803 call handle_nf90_err(ncret) 2804 ncret = nf90_get_var(ncid, ncvarid, xln(1:maxnests)) 2805 call handle_nf90_err(ncret) 2806 2807 ncret = nf90_inq_varid(ncid, 'yln', ncvarid) 2808 call handle_nf90_err(ncret) 2809 ncret = nf90_get_var(ncid, ncvarid, yln(1:maxnests)) 2810 call handle_nf90_err(ncret) 2811 2812 ncret = nf90_inq_varid(ncid, 'xrn', ncvarid) 2813 call handle_nf90_err(ncret) 2814 ncret = nf90_get_var(ncid, ncvarid, xrn(1:maxnests)) 2815 call handle_nf90_err(ncret) 2816 2817 ncret = nf90_inq_varid(ncid, 'yrn', ncvarid) 2818 call handle_nf90_err(ncret) 2819 ncret = nf90_get_var(ncid, ncvarid, yrn(1:maxnests)) 2820 call handle_nf90_err(ncret) 2821 2822 2823 2824 PRINT *, 'SUM(yresoln): ', SUM(yresoln) 2825 PRINT *, 'SUM(xrn): ', SUM(xrn) 2826 2827 2828 2321 2829 ! Variables for polar stereographic projection 2322 2830 READ(iounit) xglobal, sglobal, nglobal … … 2324 2832 READ(iounit) southpolemap(:) 2325 2833 READ(iounit) northpolemap(:) 2834 2835 2836 ncret = nf90_inq_varid(ncid, 'southpolemap', ncvarid) 2837 call handle_nf90_err(ncret) 2838 ncret = nf90_get_var(ncid, ncvarid, southpolemap(:)) 2839 call handle_nf90_err(ncret) 2840 2841 ncret = nf90_inq_varid(ncid, 'northpolemap', ncvarid) 2842 call handle_nf90_err(ncret) 2843 ncret = nf90_get_var(ncid, ncvarid, northpolemap(:)) 2844 call handle_nf90_err(ncret) 2845 2846 ncret = nf90_inq_varid(ncid, 'xglobal', ncvarid) 2847 call handle_nf90_err(ncret) 2848 ncret = nf90_get_var(ncid, ncvarid, xglobal) 2849 call handle_nf90_err(ncret) 2850 2851 ncret = nf90_inq_varid(ncid, 'sglobal', ncvarid) 2852 call handle_nf90_err(ncret) 2853 ncret = nf90_get_var(ncid, ncvarid, sglobal) 2854 call handle_nf90_err(ncret) 2855 2856 ncret = nf90_inq_varid(ncid, 'nglobal', ncvarid) 2857 call handle_nf90_err(ncret) 2858 ncret = nf90_get_var(ncid, ncvarid, nglobal) 2859 call handle_nf90_err(ncret) 2860 2861 ncret = nf90_inq_varid(ncid, 'switchnorthg', ncvarid) 2862 call handle_nf90_err(ncret) 2863 ncret = nf90_get_var(ncid, ncvarid, switchnorthg) 2864 call handle_nf90_err(ncret) 2865 2866 ncret = nf90_inq_varid(ncid, 'switchsouthg', ncvarid) 2867 call handle_nf90_err(ncret) 2868 ncret = nf90_get_var(ncid, ncvarid, switchsouthg) 2869 call handle_nf90_err(ncret) 2870 2871 2872 PRINT *, 'SUM(northpolemap): ', SUM(northpolemap) 2873 PRINT *, 'xglobal: ', xglobal 2874 PRINT *, 'sglobal: ', sglobal 2875 PRINT *, 'nglobal: ', nglobal 2876 PRINT *, 'switchsouthg: ', switchsouthg 2877 2878 2879 2326 2880 2327 2881 ! Variables declared in conv_mod (convection)
Note: See TracChangeset
for help on using the changeset viewer.