Ignore:
Timestamp:
Oct 19, 2018, 6:16:41 PM (7 years ago)
Author:
mlefevre
Message:

Photochemistry architecture for Venus mesoscale and LES model

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/Venus/readmeteo.F90

    r1756 r2021  
    252252allocate(aps(altlen),bps(altlen),levels(altlen))
    253253#ifdef PHOTOCHEM
    254 nchemtrac = 14
     254nchemtrac = 34
    255255allocate(wtnom(nchemtrac))
    256 wtnom(1)  = "c_co2"
    257 wtnom(2)  = "c_co"
    258 wtnom(3)  = "c_o"
    259 wtnom(4)  = "c_o1d"
    260 wtnom(5)  = "c_o2"
    261 wtnom(6)  = "c_o3"
    262 wtnom(7)  = "c_h"
    263 wtnom(8)  = "c_h2"
    264 wtnom(9)  = "c_oh"
    265 wtnom(10) = "c_ho2"
    266 wtnom(11) = "c_h2o2"
    267 wtnom(12) = "c_ch4"
    268 wtnom(13) = "c_n2"
    269 wtnom(14) = "c_ar"
     256print*,'PHOTOCHEM2.1'
     257wtnom(1)  = "co2"
     258wtnom(2)  = "co"
     259wtnom(3)  = "h2"
     260wtnom(4) = "h2o"
     261wtnom(5)  = "o1d"
     262wtnom(6)  = "o"
     263wtnom(7)  = "o2"
     264wtnom(8)  = "o2dg"
     265wtnom(9)  = "o3"
     266wtnom(10)  = "h"
     267wtnom(11)  = "oh"
     268wtnom(12)  = "ho2"
     269wtnom(13) = "h2o2"
     270wtnom(14)  = "cl"
     271wtnom(15)  = "clo"
     272wtnom(16)  = "cl2"
     273wtnom(17)  = "hcl"
     274wtnom(18)  = "hocl"
     275wtnom(19)  = "clco"
     276wtnom(20)  = "clco3"
     277wtnom(21)  = "cocl2"
     278wtnom(22)  = "s"
     279wtnom(23)  = "so"
     280wtnom(24)  = "so2"
     281wtnom(25)  = "so3"
     282wtnom(26)  = "s2o2"
     283wtnom(27)  = "ocs"
     284wtnom(28)  = "hso3"
     285wtnom(29)  = "h2so4"
     286wtnom(30)  = "s2"
     287wtnom(31)  = "clso2"
     288wtnom(32)  = "oscl"
     289wtnom(33)  = "h2oliq"
     290wtnom(34)  = "h2so4liq"
    270291allocate(chemtrac(lonlen,latlen,altlen,timelen,nchemtrac))
    271292chemtrac(:,:,:,:,:)=0
    … …  
    772793  vfile(j,i,:,:)=vfile(lonlen+1-j,latlen+1-i,:,:)
    773794  tsoilfile(j,i,:,:)=tsoilfile(lonlen+1-j,latlen+1-i,:,:)
     795#ifdef PHOTOCHEM
     796 chemtrac(j,i,:,:,:)=chemtrac(lonlen+1-j,latlen+1-i,:,:,:)
     797#endif
    774798ENDDO
    775799ENDDO 
    … …  
    13431367!------------------------!
    13441368#ifdef PHOTOCHEM
    1345     DO i=1,nchemtrac
    1346        FIELD=wtnom(i)
    1347        UNITS='units'
    1348        DESC='desc'
    1349        XLVL=200100.
    1350        SLAB=chemtrac(:,:,1,time_out(l),i)
    1351        ! And now put everything in the destination file
    1352        ! ... Header
    1353        write(1) IFV
    1354        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
    1355        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
    1356        ! ... Data
    1357        write(1) SLAB
    1358     ENDDO
     1369        FIELD='CO2'
     1370        UNITS='kg/kg'
     1371        DESC='CO2 mixing ratio'
     1372        XLVL=levels(k)
     1373        SLAB=chemtrac(:,:,1,time_out(l),1)
     1374                ! And now put everything in the destination file
     1375                ! ... Header
     1376        write(1) IFV
     1377        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1378        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1379                ! ... Data
     1380        write(1) SLAB
     1381
     1382
     1383        FIELD='CO'
     1384        UNITS='kg/kg'
     1385        DESC='CO mixing ratio'
     1386        XLVL=levels(k)
     1387        SLAB=chemtrac(:,:,1,time_out(l),2)
     1388                ! And now put everything in the destination file
     1389                ! ... Header
     1390        write(1) IFV
     1391        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1392        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1393                ! ... Data
     1394        write(1) SLAB
     1395
     1396
     1397        FIELD='H2'
     1398        UNITS='kg/kg'
     1399        DESC='H2 mixing ratio'
     1400        XLVL=levels(k)
     1401        SLAB=chemtrac(:,:,1,time_out(l),3)
     1402                ! And now put everything in the destination file
     1403                ! ... Header
     1404        write(1) IFV
     1405        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1406        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1407                ! ... Data
     1408        write(1) SLAB
     1409
     1410
     1411        FIELD='H2O'
     1412        UNITS='kg/kg'
     1413        DESC='H2O mixing ratio'
     1414        XLVL=levels(k)
     1415        SLAB=chemtrac(:,:,1,time_out(l),4)
     1416                ! And now put everything in the destination file
     1417                ! ... Header
     1418        write(1) IFV
     1419        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1420        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1421                ! ... Data
     1422        write(1) SLAB
     1423
     1424
     1425        FIELD='O1D'
     1426        UNITS='kg/kg'
     1427        DESC='O1D mixing ratio'
     1428        XLVL=levels(k)
     1429        SLAB=chemtrac(:,:,1,time_out(l),5)
     1430                ! And now put everything in the destination file
     1431                ! ... Header
     1432        write(1) IFV
     1433        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1434        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1435                ! ... Data
     1436        write(1) SLAB
     1437
     1438
     1439        FIELD='O'
     1440        UNITS='kg/kg'
     1441        DESC='O mixing ratio'
     1442        XLVL=levels(k)
     1443        SLAB=chemtrac(:,:,1,time_out(l),6)
     1444                ! And now put everything in the destination file
     1445                ! ... Header
     1446        write(1) IFV
     1447        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1448        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1449                ! ... Data
     1450        write(1) SLAB
     1451
     1452
     1453        FIELD='O2'
     1454        UNITS='kg/kg'
     1455        DESC='O2 mixing ratio'
     1456        XLVL=levels(k)
     1457        SLAB=chemtrac(:,:,1,time_out(l),7)
     1458                ! And now put everything in the destination file
     1459                ! ... Header
     1460        write(1) IFV
     1461        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1462        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1463                ! ... Data
     1464        write(1) SLAB
     1465
     1466
     1467        FIELD='O2DG'
     1468        UNITS='kg/kg'
     1469        DESC='O2DG mixing ratio'
     1470        XLVL=levels(k)
     1471        SLAB=chemtrac(:,:,1,time_out(l),8)
     1472                ! And now put everything in the destination file
     1473                ! ... Header
     1474        write(1) IFV
     1475        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1476        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1477                ! ... Data
     1478        write(1) SLAB
     1479
     1480
     1481        FIELD='O3'
     1482        UNITS='kg/kg'
     1483        DESC='O3 mixing ratio'
     1484        XLVL=levels(k)
     1485        SLAB=chemtrac(:,:,1,time_out(l),9)
     1486                ! And now put everything in the destination file
     1487                ! ... Header
     1488        write(1) IFV
     1489        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1490        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1491                ! ... Data
     1492        write(1) SLAB
     1493
     1494
     1495        FIELD='H'
     1496        UNITS='kg/kg'
     1497        DESC='H mixing ratio'
     1498        XLVL=levels(k)
     1499        SLAB=chemtrac(:,:,1,time_out(l),10)
     1500                ! And now put everything in the destination file
     1501                ! ... Header
     1502        write(1) IFV
     1503        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1504        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1505                ! ... Data
     1506        write(1) SLAB
     1507
     1508
     1509        FIELD='OH'
     1510        UNITS='kg/kg'
     1511        DESC='OH mixing ratio'
     1512        XLVL=levels(k)
     1513        SLAB=chemtrac(:,:,1,time_out(l),11)
     1514                ! And now put everything in the destination file
     1515                ! ... Header
     1516        write(1) IFV
     1517        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1518        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1519                ! ... Data
     1520        write(1) SLAB
     1521
     1522
     1523        FIELD='HO2'
     1524        UNITS='kg/kg'
     1525        DESC='hO2 mixing ratio'
     1526        XLVL=levels(k)
     1527        SLAB=chemtrac(:,:,1,time_out(l),12)
     1528                ! And now put everything in the destination file
     1529                ! ... Header
     1530        write(1) IFV
     1531        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1532        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1533                ! ... Data
     1534        write(1) SLAB
     1535
     1536
     1537        FIELD='H2O2'
     1538        UNITS='kg/kg'
     1539        DESC='H2O2 mixing ratio'
     1540        XLVL=levels(k)
     1541        SLAB=chemtrac(:,:,1,time_out(l),13)
     1542                ! And now put everything in the destination file
     1543                ! ... Header
     1544        write(1) IFV
     1545        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1546        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1547                ! ... Data
     1548        write(1) SLAB
     1549
     1550
     1551        FIELD='Cl'
     1552        UNITS='kg/kg'
     1553        DESC='Cl mixing ratio'
     1554        XLVL=levels(k)
     1555        SLAB=chemtrac(:,:,1,time_out(l),14)
     1556                ! And now put everything in the destination file
     1557                ! ... Header
     1558        write(1) IFV
     1559        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1560        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1561                ! ... Data
     1562        write(1) SLAB
     1563
     1564
     1565        FIELD='ClO'
     1566        UNITS='kg/kg'
     1567        DESC='ClO mixing ratio'
     1568        XLVL=levels(k)
     1569        SLAB=chemtrac(:,:,1,time_out(l),15)
     1570                ! And now put everything in the destination file
     1571                ! ... Header
     1572        write(1) IFV
     1573        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1574        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1575                ! ... Data
     1576        write(1) SLAB
     1577
     1578
     1579        FIELD='Cl2'
     1580        UNITS='kg/kg'
     1581        DESC='Cl2 mixing ratio'
     1582        XLVL=levels(k)
     1583        SLAB=chemtrac(:,:,1,time_out(l),16)
     1584                ! And now put everything in the destination file
     1585                ! ... Header
     1586        write(1) IFV
     1587        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1588        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1589                ! ... Data
     1590        write(1) SLAB
     1591
     1592
     1593        FIELD='HCl'
     1594        UNITS='kg/kg'
     1595        DESC='HCl mixing ratio'
     1596        XLVL=levels(k)
     1597        SLAB=chemtrac(:,:,1,time_out(l),17)
     1598                ! And now put everything in the destination file
     1599                ! ... Header
     1600        write(1) IFV
     1601        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1602        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1603                ! ... Data
     1604        write(1) SLAB
     1605
     1606
     1607        FIELD='HOCl'
     1608        UNITS='kg/kg'
     1609        DESC='HOCl mixing ratio'
     1610        XLVL=levels(k)
     1611        SLAB=chemtrac(:,:,1,time_out(l),18)
     1612                ! And now put everything in the destination file
     1613                ! ... Header
     1614        write(1) IFV
     1615        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1616        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1617                ! ... Data
     1618        write(1) SLAB
     1619
     1620
     1621        FIELD='ClCO'
     1622        UNITS='kg/kg'
     1623        DESC='ClCO mixing ratio'
     1624        XLVL=levels(k)
     1625        SLAB=chemtrac(:,:,1,time_out(l),19)
     1626                ! And now put everything in the destination file
     1627                ! ... Header
     1628        write(1) IFV
     1629        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1630        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1631                ! ... Data
     1632        write(1) SLAB
     1633
     1634
     1635        FIELD='ClCO3'
     1636        UNITS='kg/kg'
     1637        DESC='ClCO3 mixing ratio'
     1638        XLVL=levels(k)
     1639        SLAB=chemtrac(:,:,1,time_out(l),20)
     1640                ! And now put everything in the destination file
     1641                ! ... Header
     1642        write(1) IFV
     1643        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1644        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1645                ! ... Data
     1646        write(1) SLAB
     1647
     1648
     1649        FIELD='COCl2'
     1650        UNITS='kg/kg'
     1651        DESC='COClC2 mixing ratio'
     1652        XLVL=levels(k)
     1653        SLAB=chemtrac(:,:,1,time_out(l),21)
     1654                ! And now put everything in the destination file
     1655                ! ... Header
     1656        write(1) IFV
     1657        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1658        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1659                ! ... Data
     1660        write(1) SLAB
     1661
     1662
     1663        FIELD='S'
     1664        UNITS='kg/kg'
     1665        DESC='S mixing ratio'
     1666        XLVL=levels(k)
     1667        SLAB=chemtrac(:,:,1,time_out(l),22)
     1668                ! And now put everything in the destination file
     1669                ! ... Header
     1670        write(1) IFV
     1671        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1672        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1673                ! ... Data
     1674        write(1) SLAB
     1675
     1676
     1677        FIELD='SO'
     1678        UNITS='kg/kg'
     1679        DESC='SO mixing ratio'
     1680        XLVL=levels(k)
     1681        SLAB=chemtrac(:,:,1,time_out(l),23)
     1682                ! And now put everything in the destination file
     1683                ! ... Header
     1684        write(1) IFV
     1685        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1686        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1687                ! ... Data
     1688        write(1) SLAB
     1689
     1690
     1691        FIELD='SO2'
     1692        UNITS='kg/kg'
     1693        DESC='SO2 mixing ratio'
     1694        XLVL=levels(k)
     1695        SLAB=chemtrac(:,:,1,time_out(l),24)
     1696                ! And now put everything in the destination file
     1697                ! ... Header
     1698        write(1) IFV
     1699        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1700        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1701                ! ... Data
     1702        write(1) SLAB
     1703
     1704
     1705        FIELD='SO3'
     1706        UNITS='kg/kg'
     1707        DESC='SO3 mixing ratio'
     1708        XLVL=levels(k)
     1709        SLAB=chemtrac(:,:,1,time_out(l),25)
     1710                ! And now put everything in the destination file
     1711                ! ... Header
     1712        write(1) IFV
     1713        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1714        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1715                ! ... Data
     1716        write(1) SLAB
     1717
     1718
     1719        FIELD='S2O2'
     1720        UNITS='kg/kg'
     1721        DESC='S2O2 mixing ratio'
     1722        XLVL=levels(k)
     1723        SLAB=chemtrac(:,:,1,time_out(l),26)
     1724                ! And now put everything in the destination file
     1725                ! ... Header
     1726        write(1) IFV
     1727        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1728        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1729                ! ... Data
     1730        write(1) SLAB
     1731
     1732
     1733        FIELD='OCS'
     1734        UNITS='kg/kg'
     1735        DESC='OCS mixing ratio'
     1736        XLVL=levels(k)
     1737        SLAB=chemtrac(:,:,1,time_out(l),27)
     1738                ! And now put everything in the destination file
     1739                ! ... Header
     1740        write(1) IFV
     1741        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1742        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1743                ! ... Data
     1744        write(1) SLAB
     1745
     1746
     1747        FIELD='HSO3'
     1748        UNITS='kg/kg'
     1749        DESC='HSO3 mixing ratio'
     1750        XLVL=levels(k)
     1751        SLAB=chemtrac(:,:,1,time_out(l),28)
     1752                ! And now put everything in the destination file
     1753                ! ... Header
     1754        write(1) IFV
     1755        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1756        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1757                ! ... Data
     1758        write(1) SLAB
     1759
     1760
     1761        FIELD='H2SO4'
     1762        UNITS='kg/kg'
     1763        DESC='H2SO4 mixing ratio'
     1764        XLVL=levels(k)
     1765        SLAB=chemtrac(:,:,1,time_out(l),29)
     1766                ! And now put everything in the destination file
     1767                ! ... Header
     1768        write(1) IFV
     1769        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1770        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1771                ! ... Data
     1772        write(1) SLAB
     1773
     1774
     1775        FIELD='S2'
     1776        UNITS='kg/kg'
     1777        DESC='S2 mixing ratio'
     1778        XLVL=levels(k)
     1779        SLAB=chemtrac(:,:,1,time_out(l),30)
     1780                ! And now put everything in the destination file
     1781                ! ... Header
     1782        write(1) IFV
     1783        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1784        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1785                ! ... Data
     1786        write(1) SLAB
     1787
     1788
     1789        FIELD='ClSO2'
     1790        UNITS='kg/kg'
     1791        DESC='ClSO2 mixing ratio'
     1792        XLVL=levels(k)
     1793        SLAB=chemtrac(:,:,1,time_out(l),31)
     1794                ! And now put everything in the destination file
     1795                ! ... Header
     1796        write(1) IFV
     1797        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1798        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1799                ! ... Data
     1800        write(1) SLAB
     1801
     1802
     1803        FIELD='OSCl'
     1804        UNITS='kg/kg'
     1805        DESC='OSCl mixing ratio'
     1806        XLVL=levels(k)
     1807        SLAB=chemtrac(:,:,1,time_out(l),32)
     1808                ! And now put everything in the destination file
     1809                ! ... Header
     1810        write(1) IFV
     1811        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1812        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1813                ! ... Data
     1814        write(1) SLAB
     1815
     1816
     1817        FIELD='H2Oliq'
     1818        UNITS='kg/kg'
     1819        DESC='H2Oliq mixing ratio'
     1820        XLVL=levels(k)
     1821        SLAB=chemtrac(:,:,1,time_out(l),33)
     1822                ! And now put everything in the destination file
     1823                ! ... Header
     1824        write(1) IFV
     1825        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1826        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1827                ! ... Data
     1828        write(1) SLAB
     1829
     1830
     1831        FIELD='H2SO4liq'
     1832        UNITS='kg/kg'
     1833        DESC='H2SO4liq mixing ratio'
     1834        XLVL=levels(k)
     1835        SLAB=chemtrac(:,:,1,time_out(l),34)
     1836                ! And now put everything in the destination file
     1837                ! ... Header
     1838        write(1) IFV
     1839        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1840        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1841                ! ... Data
     1842        write(1) SLAB
    13591843#endif
    13601844
    … …  
    16862170!------------------------!
    16872171#ifdef PHOTOCHEM
    1688     DO i=1,nchemtrac
    1689        FIELD=wtnom(i)
    1690        UNITS='units'
    1691        DESC='desc'
    1692        DO k = 1,altlen
    1693          XLVL=levels(k)
    1694          SLAB=chemtrac(:,:,k,time_out(l),i)
    1695          ! And now put everything in the destination file
    1696          ! ... Header
    1697          write(1) IFV
    1698          write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
    1699          write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
    1700          ! ... Data
    1701          write(1) SLAB
    1702        END DO
    1703     ENDDO
     2172FIELD='CO2'
     2173UNITS='kg/kg'
     2174DESC='CO2 mixing ratio'
     2175DO k = 1,altlen
     2176        XLVL=levels(k)
     2177        SLAB=chemtrac(:,:,k,time_out(l),1)
     2178                ! And now put everything in the destination file
     2179                ! ... Header
     2180        write(1) IFV
     2181        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2182        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2183                ! ... Data
     2184                write(1) SLAB
     2185END DO
     2186
     2187FIELD='CO'
     2188UNITS='kg/kg'
     2189DESC='CO mixing ratio'
     2190DO k = 1,altlen
     2191        XLVL=levels(k)
     2192        SLAB=chemtrac(:,:,k,time_out(l),2)
     2193                ! And now put everything in the destination file
     2194                ! ... Header
     2195        write(1) IFV
     2196        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2197        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2198                ! ... Data
     2199                write(1) SLAB
     2200END DO
     2201
     2202FIELD='H2'
     2203UNITS='kg/kg'
     2204DESC='H2 mixing ratio'
     2205DO k = 1,altlen
     2206        XLVL=levels(k)
     2207        SLAB=chemtrac(:,:,k,time_out(l),3)
     2208                ! And now put everything in the destination file
     2209                ! ... Header
     2210        write(1) IFV
     2211        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2212        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2213                ! ... Data
     2214                write(1) SLAB
     2215END DO
     2216
     2217FIELD='H2O'
     2218UNITS='kg/kg'
     2219DESC='H2O mixing ratio'
     2220DO k = 1,altlen
     2221        XLVL=levels(k)
     2222        SLAB=chemtrac(:,:,k,time_out(l),4)
     2223                ! And now put everything in the destination file
     2224                ! ... Header
     2225        write(1) IFV
     2226        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2227        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2228                ! ... Data
     2229                write(1) SLAB
     2230END DO
     2231
     2232FIELD='O1D'
     2233UNITS='kg/kg'
     2234DESC='O1D mixing ratio'
     2235DO k = 1,altlen
     2236        XLVL=levels(k)
     2237        SLAB=chemtrac(:,:,k,time_out(l),5)
     2238                ! And now put everything in the destination file
     2239                ! ... Header
     2240        write(1) IFV
     2241        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2242        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2243                ! ... Data
     2244                write(1) SLAB
     2245END DO
     2246
     2247FIELD='O'
     2248UNITS='kg/kg'
     2249DESC='O mixing ratio'
     2250DO k = 1,altlen
     2251        XLVL=levels(k)
     2252        SLAB=chemtrac(:,:,k,time_out(l),6)
     2253                ! And now put everything in the destination file
     2254                ! ... Header
     2255        write(1) IFV
     2256        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2257        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2258                ! ... Data
     2259                write(1) SLAB
     2260END DO
     2261
     2262FIELD='O2'
     2263UNITS='kg/kg'
     2264DESC='O2 mixing ratio'
     2265DO k = 1,altlen
     2266        XLVL=levels(k)
     2267        SLAB=chemtrac(:,:,k,time_out(l),7)
     2268                ! And now put everything in the destination file
     2269                ! ... Header
     2270        write(1) IFV
     2271        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2272        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2273                ! ... Data
     2274                write(1) SLAB
     2275END DO
     2276
     2277FIELD='O2DG'
     2278UNITS='kg/kg'
     2279DESC='O2DG mixing ratio'
     2280DO k = 1,altlen
     2281        XLVL=levels(k)
     2282        SLAB=chemtrac(:,:,k,time_out(l),8)
     2283                ! And now put everything in the destination file
     2284                ! ... Header
     2285        write(1) IFV
     2286        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2287        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2288                ! ... Data
     2289                write(1) SLAB
     2290END DO
     2291
     2292FIELD='O3'
     2293UNITS='kg/kg'
     2294DESC='O3 mixing ratio'
     2295DO k = 1,altlen
     2296        XLVL=levels(k)
     2297        SLAB=chemtrac(:,:,k,time_out(l),9)
     2298                ! And now put everything in the destination file
     2299                ! ... Header
     2300        write(1) IFV
     2301        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2302        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2303                ! ... Data
     2304                write(1) SLAB
     2305END DO
     2306
     2307FIELD='H'
     2308UNITS='kg/kg'
     2309DESC='H mixing ratio'
     2310DO k = 1,altlen
     2311        XLVL=levels(k)
     2312        SLAB=chemtrac(:,:,k,time_out(l),10)
     2313                ! And now put everything in the destination file
     2314                ! ... Header
     2315        write(1) IFV
     2316        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2317        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2318                ! ... Data
     2319                write(1) SLAB
     2320END DO
     2321
     2322FIELD='OH'
     2323UNITS='kg/kg'
     2324DESC='OH mixing ratio'
     2325DO k = 1,altlen
     2326        XLVL=levels(k)
     2327        SLAB=chemtrac(:,:,k,time_out(l),11)
     2328                ! And now put everything in the destination file
     2329                ! ... Header
     2330        write(1) IFV
     2331        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2332        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2333                ! ... Data
     2334                write(1) SLAB
     2335END DO
     2336
     2337FIELD='HO2'
     2338UNITS='kg/kg'
     2339DESC='hO2 mixing ratio'
     2340DO k = 1,altlen
     2341        XLVL=levels(k)
     2342        SLAB=chemtrac(:,:,k,time_out(l),12)
     2343                ! And now put everything in the destination file
     2344                ! ... Header
     2345        write(1) IFV
     2346        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2347        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2348                ! ... Data
     2349                write(1) SLAB
     2350END DO
     2351
     2352FIELD='H2O2'
     2353UNITS='kg/kg'
     2354DESC='H2O2 mixing ratio'
     2355DO k = 1,altlen
     2356        XLVL=levels(k)
     2357        SLAB=chemtrac(:,:,k,time_out(l),13)
     2358                ! And now put everything in the destination file
     2359                ! ... Header
     2360        write(1) IFV
     2361        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2362        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2363                ! ... Data
     2364                write(1) SLAB
     2365END DO
     2366
     2367FIELD='Cl'
     2368UNITS='kg/kg'
     2369DESC='Cl mixing ratio'
     2370DO k = 1,altlen
     2371        XLVL=levels(k)
     2372        SLAB=chemtrac(:,:,k,time_out(l),14)
     2373                ! And now put everything in the destination file
     2374                ! ... Header
     2375        write(1) IFV
     2376        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2377        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2378                ! ... Data
     2379                write(1) SLAB
     2380END DO
     2381
     2382FIELD='ClO'
     2383UNITS='kg/kg'
     2384DESC='ClO mixing ratio'
     2385DO k = 1,altlen
     2386        XLVL=levels(k)
     2387        SLAB=chemtrac(:,:,k,time_out(l),15)
     2388                ! And now put everything in the destination file
     2389                ! ... Header
     2390        write(1) IFV
     2391        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2392        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2393                ! ... Data
     2394                write(1) SLAB
     2395END DO
     2396
     2397FIELD='Cl2'
     2398UNITS='kg/kg'
     2399DESC='Cl2 mixing ratio'
     2400DO k = 1,altlen
     2401        XLVL=levels(k)
     2402        SLAB=chemtrac(:,:,k,time_out(l),16)
     2403                ! And now put everything in the destination file
     2404                ! ... Header
     2405        write(1) IFV
     2406        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2407        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2408                ! ... Data
     2409                write(1) SLAB
     2410END DO
     2411
     2412FIELD='HCl'
     2413UNITS='kg/kg'
     2414DESC='HCl mixing ratio'
     2415DO k = 1,altlen
     2416        XLVL=levels(k)
     2417        SLAB=chemtrac(:,:,k,time_out(l),17)
     2418                ! And now put everything in the destination file
     2419                ! ... Header
     2420        write(1) IFV
     2421        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2422        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2423                ! ... Data
     2424                write(1) SLAB
     2425END DO
     2426
     2427FIELD='HOCl'
     2428UNITS='kg/kg'
     2429DESC='HOCl mixing ratio'
     2430DO k = 1,altlen
     2431        XLVL=levels(k)
     2432        SLAB=chemtrac(:,:,k,time_out(l),18)
     2433                ! And now put everything in the destination file
     2434                ! ... Header
     2435        write(1) IFV
     2436        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2437        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2438                ! ... Data
     2439                write(1) SLAB
     2440END DO
     2441
     2442FIELD='ClCO'
     2443UNITS='kg/kg'
     2444DESC='ClCO mixing ratio'
     2445DO k = 1,altlen
     2446        XLVL=levels(k)
     2447        SLAB=chemtrac(:,:,k,time_out(l),19)
     2448                ! And now put everything in the destination file
     2449                ! ... Header
     2450        write(1) IFV
     2451        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2452        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2453                ! ... Data
     2454                write(1) SLAB
     2455END DO
     2456
     2457FIELD='ClCO3'
     2458UNITS='kg/kg'
     2459DESC='ClCO3 mixing ratio'
     2460DO k = 1,altlen
     2461        XLVL=levels(k)
     2462        SLAB=chemtrac(:,:,k,time_out(l),20)
     2463                ! And now put everything in the destination file
     2464                ! ... Header
     2465        write(1) IFV
     2466        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2467        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2468                ! ... Data
     2469                write(1) SLAB
     2470END DO
     2471
     2472FIELD='COCl2'
     2473UNITS='kg/kg'
     2474DESC='COClC2 mixing ratio'
     2475DO k = 1,altlen
     2476        XLVL=levels(k)
     2477        SLAB=chemtrac(:,:,k,time_out(l),21)
     2478                ! And now put everything in the destination file
     2479                ! ... Header
     2480        write(1) IFV
     2481        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2482        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2483                ! ... Data
     2484                write(1) SLAB
     2485END DO
     2486
     2487FIELD='S'
     2488UNITS='kg/kg'
     2489DESC='S mixing ratio'
     2490DO k = 1,altlen
     2491        XLVL=levels(k)
     2492        SLAB=chemtrac(:,:,k,time_out(l),22)
     2493                ! And now put everything in the destination file
     2494                ! ... Header
     2495        write(1) IFV
     2496        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2497        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2498                ! ... Data
     2499                write(1) SLAB
     2500END DO
     2501
     2502FIELD='SO'
     2503UNITS='kg/kg'
     2504DESC='SO mixing ratio'
     2505DO k = 1,altlen
     2506        XLVL=levels(k)
     2507        SLAB=chemtrac(:,:,k,time_out(l),23)
     2508                ! And now put everything in the destination file
     2509                ! ... Header
     2510        write(1) IFV
     2511        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2512        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2513                ! ... Data
     2514                write(1) SLAB
     2515END DO
     2516
     2517FIELD='SO2'
     2518UNITS='kg/kg'
     2519DESC='SO2 mixing ratio'
     2520DO k = 1,altlen
     2521        XLVL=levels(k)
     2522        SLAB=chemtrac(:,:,k,time_out(l),24)
     2523                ! And now put everything in the destination file
     2524                ! ... Header
     2525        write(1) IFV
     2526        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2527        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2528                ! ... Data
     2529                write(1) SLAB
     2530END DO
     2531
     2532FIELD='SO3'
     2533UNITS='kg/kg'
     2534DESC='SO3 mixing ratio'
     2535DO k = 1,altlen
     2536        XLVL=levels(k)
     2537        SLAB=chemtrac(:,:,k,time_out(l),25)
     2538                ! And now put everything in the destination file
     2539                ! ... Header
     2540        write(1) IFV
     2541        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2542        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2543                ! ... Data
     2544                write(1) SLAB
     2545END DO
     2546
     2547FIELD='S2O2'
     2548UNITS='kg/kg'
     2549DESC='S2O2 mixing ratio'
     2550DO k = 1,altlen
     2551        XLVL=levels(k)
     2552        SLAB=chemtrac(:,:,k,time_out(l),26)
     2553                ! And now put everything in the destination file
     2554                ! ... Header
     2555        write(1) IFV
     2556        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2557        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2558                ! ... Data
     2559                write(1) SLAB
     2560END DO
     2561
     2562FIELD='OCS'
     2563UNITS='kg/kg'
     2564DESC='OCS mixing ratio'
     2565DO k = 1,altlen
     2566        XLVL=levels(k)
     2567        SLAB=chemtrac(:,:,k,time_out(l),27)
     2568                ! And now put everything in the destination file
     2569                ! ... Header
     2570        write(1) IFV
     2571        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2572        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2573                ! ... Data
     2574                write(1) SLAB
     2575END DO
     2576
     2577FIELD='HSO3'
     2578UNITS='kg/kg'
     2579DESC='HSO3 mixing ratio'
     2580DO k = 1,altlen
     2581        XLVL=levels(k)
     2582        SLAB=chemtrac(:,:,k,time_out(l),28)
     2583                ! And now put everything in the destination file
     2584                ! ... Header
     2585        write(1) IFV
     2586        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2587        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2588                ! ... Data
     2589                write(1) SLAB
     2590END DO
     2591
     2592FIELD='H2SO4'
     2593UNITS='kg/kg'
     2594DESC='H2SO4 mixing ratio'
     2595DO k = 1,altlen
     2596        XLVL=levels(k)
     2597        SLAB=chemtrac(:,:,k,time_out(l),29)
     2598                ! And now put everything in the destination file
     2599                ! ... Header
     2600        write(1) IFV
     2601        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2602        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2603                ! ... Data
     2604                write(1) SLAB
     2605END DO
     2606
     2607FIELD='S2'
     2608UNITS='kg/kg'
     2609DESC='S2 mixing ratio'
     2610DO k = 1,altlen
     2611        XLVL=levels(k)
     2612        SLAB=chemtrac(:,:,k,time_out(l),30)
     2613                ! And now put everything in the destination file
     2614                ! ... Header
     2615        write(1) IFV
     2616        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2617        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2618                ! ... Data
     2619                write(1) SLAB
     2620END DO
     2621
     2622FIELD='ClSO2'
     2623UNITS='kg/kg'
     2624DESC='ClSO2 mixing ratio'
     2625DO k = 1,altlen
     2626        XLVL=levels(k)
     2627        SLAB=chemtrac(:,:,k,time_out(l),31)
     2628                ! And now put everything in the destination file
     2629                ! ... Header
     2630        write(1) IFV
     2631        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2632        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2633                ! ... Data
     2634                write(1) SLAB
     2635END DO
     2636
     2637FIELD='OSCl'
     2638UNITS='kg/kg'
     2639DESC='OSCl mixing ratio'
     2640DO k = 1,altlen
     2641        XLVL=levels(k)
     2642        SLAB=chemtrac(:,:,k,time_out(l),32)
     2643                ! And now put everything in the destination file
     2644                ! ... Header
     2645        write(1) IFV
     2646        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2647        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2648                ! ... Data
     2649                write(1) SLAB
     2650END DO
     2651
     2652FIELD='H2Oliq'
     2653UNITS='kg/kg'
     2654DESC='H2Oliq mixing ratio'
     2655DO k = 1,altlen
     2656        XLVL=levels(k)
     2657        SLAB=chemtrac(:,:,k,time_out(l),33)
     2658                ! And now put everything in the destination file
     2659                ! ... Header
     2660        write(1) IFV
     2661        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2662        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2663                ! ... Data
     2664                write(1) SLAB
     2665END DO
     2666
     2667FIELD='H2SO4liq'
     2668UNITS='kg/kg'
     2669DESC='H2SO4liq mixing ratio'
     2670DO k = 1,altlen
     2671        XLVL=levels(k)
     2672        SLAB=chemtrac(:,:,k,time_out(l),34)
     2673                ! And now put everything in the destination file
     2674                ! ... Header
     2675        write(1) IFV
     2676        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2677        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2678                ! ... Data
     2679                write(1) SLAB
     2680END DO
    17042681#endif
    17052682
Note: See TracChangeset for help on using the changeset viewer.