Changeset 9351
- Timestamp:
- 2018-02-23T14:18:47+01:00 (7 years ago)
- Location:
- branches/NERC/dev_r5518_GO6_rev9312_xCO2/NEMOGCM
- Files:
-
- 7 edited
- 1 copied
-
CONFIG/SHARED/field_def_bgc.xml (modified) (17 diffs)
-
CONFIG/SHARED/xco2.atm (copied) (copied from branches/NERC/dev_r5518_GO6_CO2_cmip/NEMOGCM/CONFIG/SHARED/xco2.atm)
-
NEMO/TOP_SRC/MEDUSA/air_sea.F90 (modified) (2 diffs)
-
NEMO/TOP_SRC/MEDUSA/carb_chem.F90 (modified) (1 diff)
-
NEMO/TOP_SRC/MEDUSA/sms_medusa.F90 (modified) (1 diff)
-
NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (15 diffs)
-
NEMO/TOP_SRC/MEDUSA/trcini_medusa.F90 (modified) (2 diffs)
-
NEMO/TOP_SRC/trcrst.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_GO6_rev9312_xCO2/NEMOGCM/CONFIG/SHARED/field_def_bgc.xml
r9257 r9351 86 86 <field id="PHN_E3T" long_name="PHN * E3T" unit="mmol-N/m2" > PHN * e3t </field > 87 87 <field id="PHN_INV" field_ref="PHN_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 88 <field id="SURF_PHN" field_ref="PHN" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 88 89 <field id="PHD" long_name="diatom phytoplankton" unit="mmol-N/m3" /> 89 90 <field id="PHD_E3T" long_name="PHD * E3T" unit="mmol-N/m2" > PHD * e3t </field > 90 91 <field id="PHD_INV" field_ref="PHD_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 92 <field id="SURF_PHD" field_ref="PHD" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 91 93 <field id="ZMI" long_name="micro zooplankton" unit="mmol-N/m3" /> 92 94 <field id="ZMI_E3T" long_name="ZMI * E3T" unit="mmol-N/m2" > ZMI * e3t </field > 93 95 <field id="ZMI_INV" field_ref="ZMI_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 96 <field id="SURF_ZMI" field_ref="ZMI" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 94 97 <field id="ZME" long_name="meso zooplankton" unit="mmol-N/m3" /> 95 98 <field id="ZME_E3T" long_name="ZME * E3T" unit="mmol-N/m2" > ZME * e3t </field > 96 99 <field id="ZME_INV" field_ref="ZME_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 100 <field id="SURF_ZME" field_ref="ZME" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 97 101 <field id="DIN" long_name="dissolved inorganic nitrogen" unit="mmol-N/m3" /> 98 102 <field id="DIN_E3T" long_name="DIN * E3T" unit="mmol-N/m2" > DIN * e3t </field > 99 103 <field id="DIN_INV" field_ref="DIN_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 104 <field id="SURF_DIN" field_ref="DIN" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 100 105 <field id="SIL" long_name="dissolved silicic acid" unit="mmolSi/m3" /> 101 106 <field id="SIL_E3T" long_name="SIL * E3T" unit="mmolSi/m2" > SIL * e3t </field > 102 107 <field id="SIL_INV" field_ref="SIL_E3T" unit="mmolSi/m2" grid_ref="vert_sum" /> 108 <field id="SURF_SIL" field_ref="SIL" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 103 109 <field id="FER" long_name="dissolved iron" unit="mmolFe/m3" /> 104 110 <field id="FER_E3T" long_name="FER * E3T" unit="mmolFe/m2" > FER * e3t </field > 105 111 <field id="FER_INV" field_ref="FER_E3T" unit="mmolFe/m2" grid_ref="vert_sum" /> 112 <field id="SURF_FER" field_ref="FER" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 106 113 <field id="DET" long_name="detrital nitrogen" unit="mmol-N/m3" /> 107 114 <field id="DET_E3T" long_name="DET * E3T" unit="mmol-N/m2" > DET * e3t </field > 108 115 <field id="DET_INV" field_ref="DET_E3T" unit="mmol-N/m2" grid_ref="vert_sum" /> 116 <field id="SURF_DET" field_ref="DET" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 109 117 <field id="PDS" long_name="biogenic silicon in diatom phytoplankton" unit="mmolSi/m3" /> 110 118 <field id="PDS_E3T" long_name="PDS * E3T" unit="mmolSi/m2" > PDS * e3t </field > 111 119 <field id="PDS_INV" field_ref="PDS_E3T" unit="mmolSi/m2" grid_ref="vert_sum" /> 120 <field id="SURF_PDS" field_ref="PDS" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 112 121 <field id="DTC" long_name="detrital carbon" unit="mmol-C/m3" /> 113 122 <field id="DTC_E3T" long_name="DTC * E3T" unit="mmol-C/m2" > DTC * e3t </field > 114 123 <field id="DTC_INV" field_ref="DTC_E3T" unit="mmol-C/m2" grid_ref="vert_sum" /> 124 <field id="SURF_DTC" field_ref="DTC" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 115 125 <field id="DiC" long_name="dissolved inorganic carbon" unit="mmol-C/m3" /> 116 126 <field id="DiC_E3T" long_name="DiC * E3T" unit="mmol-C/m2" > DiC * e3t </field > 117 127 <field id="DiC_INV" field_ref="DiC_E3T" unit="mmol-C/m2" grid_ref="vert_sum" /> 128 <field id="SURF_DiC" field_ref="DiC" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 118 129 <field id="ALK" long_name="alkalinity" unit="meq/m3" /> 119 130 <field id="ALK_E3T" long_name="ALK * E3T" unit="meq/m2" > ALK * e3t </field > 120 131 <field id="ALK_INV" field_ref="ALK_E3T" unit="meq/m2" grid_ref="vert_sum" /> 132 <field id="SURF_ALK" field_ref="ALK" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 121 133 <field id="OXY" long_name="dissolved oxygen" unit="mmolO2/m3" /> 122 134 <field id="OXY_E3T" long_name="OXY * E3T" unit="mmolO2/m2" > OXY * e3t </field > 123 135 <field id="OXY_INV" field_ref="OXY_E3T" unit="mmolO2/m2" grid_ref="vert_sum" /> 136 <field id="SURF_OXY" field_ref="OXY" unit="mg Chl/m3" grid_ref="grid_T_surface_extract" /> 124 137 125 138 <!-- Ideal traceurs : variables available with key_idtra --> … … 1646 1659 <field field_ref="CHN" name="CHN" /> 1647 1660 <field field_ref="CHN_E3T" name="CHN_E3T" /> 1648 <field field_ref="SURF_CHN" name="SURF_CHN" />1649 1661 <field field_ref="CHD" name="CHD" /> 1650 1662 <field field_ref="CHD_E3T" name="CHD_E3T" /> 1651 <field field_ref="SURF_CHD" name="SURF_CHD" />1652 1663 <field field_ref="PHN" name="PHN" /> 1653 1664 <field field_ref="PHN_E3T" name="PHN_E3T" /> … … 1676 1687 <field field_ref="OXY" name="OXY" /> 1677 1688 <field field_ref="OXY_E3T" name="OXY_E3T" /> 1678 </field_group> 1679 1689 </field_group> 1690 1691 <field_group id="groupMEDUSA_surf" > 1692 <field field_ref="SURF_CHN" name="SURF_CHN" /> 1693 <field field_ref="SURF_CHD" name="SURF_CHD" /> 1694 <field field_ref="SURF_PHN" name="SURF_PHN" /> 1695 <field field_ref="SURF_PHD" name="SURF_PHD" /> 1696 <field field_ref="SURF_ZMI" name="SURF_ZMI" /> 1697 <field field_ref="SURF_ZME" name="SURF_ZME" /> 1698 <field field_ref="SURF_DIN" name="SURF_DIN" /> 1699 <field field_ref="SURF_SIL" name="SURF_SIL" /> 1700 <field field_ref="SURF_FER" name="SURF_FER" /> 1701 <field field_ref="SURF_DET" name="SURF_DET" /> 1702 <field field_ref="SURF_PDS" name="SURF_PDS" /> 1703 <field field_ref="SURF_DTC" name="SURF_DTC" /> 1704 <field field_ref="SURF_DiC" name="SURF_DiC" /> 1705 <field field_ref="SURF_ALK" name="SURF_ALK" /> 1706 <field field_ref="SURF_OXY" name="SURF_OXY" /> 1707 </field_group> 1680 1708 1681 1709 <field_group id="groupMEDUSA_INV" > … … 1730 1758 1731 1759 1732 <field_group id="trendMEDUSA " >1760 <field_group id="trendMEDUSA_long" > 1733 1761 <!-- trd case: jptra_xad --> 1734 1762 <field field_ref="XAD_CHN" name="XAD_CHN" /> … … 1781 1809 <!-- trd case: jptra_totad --> 1782 1810 <field field_ref="TAD_CHN" name="TAD_CHN" /> 1783 <field field_ref="TAD_CHN_INV" name="TAD_CHN_INV" />1784 1811 <field field_ref="TAD_CHD" name="TAD_CHD" /> 1785 <field field_ref="TAD_CHD_INV" name="TAD_CHD_INV" />1786 1812 <field field_ref="TAD_PHN" name="TAD_PHN" /> 1787 <field field_ref="TAD_PHN_INV" name="TAD_PHN_INV" />1788 1813 <field field_ref="TAD_PHD" name="TAD_PHD" /> 1789 <field field_ref="TAD_PHD_INV" name="TAD_PHD_INV" />1790 1814 <field field_ref="TAD_ZMI" name="TAD_ZMI" /> 1791 <field field_ref="TAD_ZMI_INV" name="TAD_ZMI_INV" />1792 1815 <field field_ref="TAD_ZME" name="TAD_ZME" /> 1793 <field field_ref="TAD_ZME_INV" name="TAD_ZME_INV" />1794 1816 <field field_ref="TAD_DIN" name="TAD_DIN" /> 1795 <field field_ref="TAD_DIN_INV" name="TAD_DIN_INV" />1796 1817 <field field_ref="TAD_SIL" name="TAD_SIL" /> 1797 <field field_ref="TAD_SIL_INV" name="TAD_SIL_INV" />1798 1818 <field field_ref="TAD_FER" name="TAD_FER" /> 1799 <field field_ref="TAD_FER_INV" name="TAD_FER_INV" />1800 1819 <field field_ref="TAD_DET" name="TAD_DET" /> 1801 <field field_ref="TAD_DET_INV" name="TAD_DET_INV" />1802 1820 <field field_ref="TAD_PDS" name="TAD_PDS" /> 1803 <field field_ref="TAD_PDS_INV" name="TAD_PDS_INV" />1804 1821 <field field_ref="TAD_DTC" name="TAD_DTC" /> 1805 <field field_ref="TAD_DTC_INV" name="TAD_DTC_INV" />1806 1822 <field field_ref="TAD_DiC" name="TAD_DiC" /> 1807 <field field_ref="TAD_DiC_INV" name="TAD_DiC_INV" />1808 1823 <field field_ref="TAD_ALK" name="TAD_ALK" /> 1809 <field field_ref="TAD_ALK_INV" name="TAD_ALK_INV" />1810 1824 <field field_ref="TAD_OXY" name="TAD_OXY" /> 1811 <field field_ref="TAD_OXY_INV" name="TAD_OXY_INV" />1812 1825 <!-- trd case: jptra_ldf --> 1813 1826 <field field_ref="LDF_CHN" name="LDF_CHN" /> … … 1909 1922 <!-- trd case: jptra_sms --> 1910 1923 <field field_ref="SMS_CHN" name="SMS_CHN" /> 1911 <field field_ref="SMS_CHN_INV" name="SMS_CHN_INV" />1912 1924 <field field_ref="SMS_CHD" name="SMS_CHD" /> 1913 <field field_ref="SMS_CHD_INV" name="SMS_CHD_INV" />1914 1925 <field field_ref="SMS_PHN" name="SMS_PHN" /> 1915 <field field_ref="SMS_PHN_INV" name="SMS_PHN_INV" />1916 1926 <field field_ref="SMS_PHD" name="SMS_PHD" /> 1917 <field field_ref="SMS_PHD_INV" name="SMS_PHD_INV" />1918 1927 <field field_ref="SMS_ZMI" name="SMS_ZMI" /> 1919 <field field_ref="SMS_ZMI_INV" name="SMS_ZMI_INV" />1920 1928 <field field_ref="SMS_ZME" name="SMS_ZME" /> 1921 <field field_ref="SMS_ZME_INV" name="SMS_ZME_INV" />1922 1929 <field field_ref="SMS_DIN" name="SMS_DIN" /> 1923 <field field_ref="SMS_DIN_INV" name="SMS_DIN_INV" />1924 1930 <field field_ref="SMS_SIL" name="SMS_SIL" /> 1925 <field field_ref="SMS_SIL_INV" name="SMS_SIL_INV" />1926 1931 <field field_ref="SMS_FER" name="SMS_FER" /> 1927 <field field_ref="SMS_FER_INV" name="SMS_FER_INV" />1928 1932 <field field_ref="SMS_DET" name="SMS_DET" /> 1929 <field field_ref="SMS_DET_INV" name="SMS_DET_INV" />1930 1933 <field field_ref="SMS_PDS" name="SMS_PDS" /> 1931 <field field_ref="SMS_PDS_INV" name="SMS_PDS_INV" />1932 1934 <field field_ref="SMS_DTC" name="SMS_DTC" /> 1933 <field field_ref="SMS_DTC_INV" name="SMS_DTC_INV" />1934 1935 <field field_ref="SMS_DiC" name="SMS_DiC" /> 1935 <field field_ref="SMS_DiC_INV" name="SMS_DiC_INV" />1936 1936 <field field_ref="SMS_ALK" name="SMS_ALK" /> 1937 <field field_ref="SMS_ALK_INV" name="SMS_ALK_INV" />1938 1937 <field field_ref="SMS_OXY" name="SMS_OXY" /> 1939 <field field_ref="SMS_OXY_INV" name="SMS_OXY_INV" />1940 1938 <!-- trd case: jptra_atf --> 1941 1939 <field field_ref="ATF_CHN" name="ATF_CHN" /> … … 1956 1954 <!-- trd case: jptra_tot --> 1957 1955 <field field_ref="TOT_CHN" name="TOT_CHN" /> 1958 <field field_ref="TOT_CHN_INV" name="TOT_CHN_INV" />1959 1956 <field field_ref="TOT_CHD" name="TOT_CHD" /> 1960 <field field_ref="TOT_CHD_INV" name="TOT_CHD_INV" />1961 1957 <field field_ref="TOT_PHN" name="TOT_PHN" /> 1962 <field field_ref="TOT_PHN_INV" name="TOT_PHN_INV" />1963 1958 <field field_ref="TOT_PHD" name="TOT_PHD" /> 1964 <field field_ref="TOT_PHD_INV" name="TOT_PHD_INV" />1965 1959 <field field_ref="TOT_ZMI" name="TOT_ZMI" /> 1966 <field field_ref="TOT_ZMI_INV" name="TOT_ZMI_INV" />1967 1960 <field field_ref="TOT_ZME" name="TOT_ZME" /> 1968 <field field_ref="TOT_ZME_INV" name="TOT_ZME_INV" />1969 1961 <field field_ref="TOT_DIN" name="TOT_DIN" /> 1970 <field field_ref="TOT_DIN_INV" name="TOT_DIN_INV" />1971 1962 <field field_ref="TOT_SIL" name="TOT_SIL" /> 1972 <field field_ref="TOT_SIL_INV" name="TOT_SIL_INV" />1973 1963 <field field_ref="TOT_FER" name="TOT_FER" /> 1974 <field field_ref="TOT_FER_INV" name="TOT_FER_INV" />1975 1964 <field field_ref="TOT_DET" name="TOT_DET" /> 1976 <field field_ref="TOT_DET_INV" name="TOT_DET_INV" />1977 1965 <field field_ref="TOT_PDS" name="TOT_PDS" /> 1978 <field field_ref="TOT_PDS_INV" name="TOT_PDS_INV" />1979 1966 <field field_ref="TOT_DTC" name="TOT_DTC" /> 1980 <field field_ref="TOT_DTC_INV" name="TOT_DTC_INV" />1981 1967 <field field_ref="TOT_DiC" name="TOT_DiC" /> 1982 <field field_ref="TOT_DiC_INV" name="TOT_DiC_INV" />1983 1968 <field field_ref="TOT_ALK" name="TOT_ALK" /> 1984 <field field_ref="TOT_ALK_INV" name="TOT_ALK_INV" />1985 1969 <field field_ref="TOT_OXY" name="TOT_OXY" /> 1986 <field field_ref="TOT_OXY_INV" name="TOT_OXY_INV" />1987 1970 <!-- trd case: jptra_iso --> 1988 1971 <field field_ref="ISO_CHN" name="ISO_CHN" /> … … 2003 1986 </field_group> 2004 1987 2005 <field_group id="trendCFC11 " >1988 <field_group id="trendCFC11_long" > 2006 1989 <!-- cfc11 --> 2007 1990 <field field_ref="TOT_CFC11" name="TOT_CFC11" standard_name="CFC-11_Concentration" /> 2008 <field field_ref="TOT_CFC11_INV" name="TOT_CFC11_INV" standard_name="CFC-11_Concentration" />2009 1991 <field field_ref="XAD_CFC11" name="XAD_CFC11" standard_name="CFC-11_Concentration" /> 2010 1992 <field field_ref="YAD_CFC11" name="YAD_CFC11" standard_name="CFC-11_Concentration" /> 2011 1993 <field field_ref="ZAD_CFC11" name="ZAD_CFC11" standard_name="CFC-11_Concentration" /> 2012 1994 <field field_ref="TAD_CFC11" name="TAD_CFC11" standard_name="CFC-11_Concentration" /> 2013 <field field_ref="TAD_CFC11_INV" name="TAD_CFC11_INV" standard_name="CFC-11_Concentration" />2014 1995 <field field_ref="LDF_CFC11" name="LDF_CFC11" standard_name="CFC-11_Concentration" /> 2015 1996 <field field_ref="BBL_CFC11" name="BBL_CFC11" standard_name="CFC-11_Concentration" /> … … 2019 2000 <!-- <field field_ref="DMP_CFC11" name="DMP_CFC11" standard_name="CFC-11_Concentration" /> --> 2020 2001 <field field_ref="SMS_CFC11" name="SMS_CFC11" standard_name="CFC-11_Concentration" /> 2021 <field field_ref="SMS_CFC11_INV" name="SMS_CFC11_INV" standard_name="CFC-11_Concentration" />2022 2002 <field field_ref="ATF_CFC11" name="ATF_CFC11" standard_name="CFC-11_Concentration" /> 2023 2003 <field field_ref="ISO_CFC11" name="ISO_CFC11" standard_name="CFC-11_Concentration" /> 2024 2004 </field_group> 2025 2005 2026 <field_group id="trendCFC12 " >2006 <field_group id="trendCFC12_long" > 2027 2007 <!-- cfc12 --> 2028 2008 <field field_ref="TOT_CFC12" name="TOT_CFC12" standard_name="CFC-12_Concentration" /> 2029 <field field_ref="TOT_CFC12_INV" name="TOT_CFC12_INV" standard_name="CFC-12_Concentration" />2030 2009 <field field_ref="XAD_CFC12" name="XAD_CFC12" standard_name="CFC-12_Concentration" /> 2031 2010 <field field_ref="YAD_CFC12" name="YAD_CFC12" standard_name="CFC-12_Concentration" /> 2032 2011 <field field_ref="ZAD_CFC12" name="ZAD_CFC12" standard_name="CFC-12_Concentration" /> 2033 2012 <field field_ref="TAD_CFC12" name="TAD_CFC12" standard_name="CFC-12_Concentration" /> 2034 <field field_ref="TAD_CFC12_INV" name="TAD_CFC12_INV" standard_name="CFC-12_Concentration" />2035 2013 <field field_ref="LDF_CFC12" name="LDF_CFC12" standard_name="CFC-12_Concentration" /> 2036 2014 <field field_ref="BBL_CFC12" name="BBL_CFC12" standard_name="CFC-12_Concentration" /> … … 2040 2018 <!-- <field field_ref="DMP_CFC12" name="DMP_CFC12" standard_name="CFC-12_Concentration" /> --> 2041 2019 <field field_ref="SMS_CFC12" name="SMS_CFC12" standard_name="CFC-12_Concentration" /> 2042 <field field_ref="SMS_CFC12_INV" name="SMS_CFC12_INV" standard_name="CFC-12_Concentration" />2043 2020 <field field_ref="ATF_CFC12" name="ATF_CFC12" standard_name="CFC-12_Concentration" /> 2044 2021 <field field_ref="ISO_CFC12" name="ISO_CFC12" standard_name="CFC-12_Concentration" /> 2045 2022 </field_group> 2046 2023 2047 <field_group id="trendSF6 " >2024 <field_group id="trendSF6_long" > 2048 2025 <!-- sf6 --> 2049 2026 <field field_ref="TOT_SF6" name="TOT_SF6" standard_name="SF6_Concentration" /> 2050 <field field_ref="TOT_SF6_INV" name="TOT_SF6_INV" standard_name="SF6_Concentration" />2051 2027 <field field_ref="XAD_SF6" name="XAD_SF6" standard_name="SF6_Concentration" /> 2052 2028 <field field_ref="YAD_SF6" name="YAD_SF6" standard_name="SF6_Concentration" /> 2053 2029 <field field_ref="ZAD_SF6" name="ZAD_SF6" standard_name="SF6_Concentration" /> 2054 2030 <field field_ref="TAD_SF6" name="TAD_SF6" standard_name="SF6_Concentration" /> 2055 <field field_ref="TAD_SF6_INV" name="TAD_SF6_INV" standard_name="SF6_Concentration" />2056 2031 <field field_ref="LDF_SF6" name="LDF_SF6" standard_name="SF6_Concentration" /> 2057 2032 <field field_ref="BBL_SF6" name="BBL_SF6" standard_name="SF6_Concentration" /> … … 2061 2036 <!-- <field field_ref="DMP_SF6" name="DMP_SF6" standard_name="SF6_Concentration" /> --> 2062 2037 <field field_ref="SMS_SF6" name="SMS_SF6" standard_name="SF6_Concentration" /> 2063 <field field_ref="SMS_SF6_INV" name="SMS_SF6_INV" standard_name="SF6_Concentration" />2064 2038 <field field_ref="ATF_SF6" name="ATF_SF6" standard_name="SF6_Concentration" /> 2065 2039 <field field_ref="ISO_SF6" name="ISO_SF6" standard_name="SF6_Concentration" /> 2066 2040 </field_group> 2067 2041 2068 <field_group id="trendAGE " >2042 <field_group id="trendAGE_long" > 2069 2043 <!-- age --> 2070 2044 <field field_ref="TOT_Age" name="TOT_Age" standard_name="sea_water_age_since_surface_contact" /> 2071 <field field_ref="TOT_Age_INV" name="TOT_Age_INV" standard_name="sea_water_age_since_surface_contact" />2072 2045 <field field_ref="XAD_Age" name="XAD_Age" standard_name="sea_water_age_since_surface_contact" /> 2073 2046 <field field_ref="YAD_Age" name="YAD_Age" standard_name="sea_water_age_since_surface_contact" /> 2074 2047 <field field_ref="ZAD_Age" name="ZAD_Age" standard_name="sea_water_age_since_surface_contact" /> 2075 2048 <field field_ref="TAD_Age" name="TAD_Age" standard_name="sea_water_age_since_surface_contact" /> 2076 <field field_ref="TAD_Age_INV" name="TAD_Age_INV" standard_name="sea_water_age_since_surface_contact" />2077 2049 <field field_ref="LDF_Age" name="LDF_Age" standard_name="sea_water_age_since_surface_contact" /> 2078 2050 <field field_ref="BBL_Age" name="BBL_Age" standard_name="sea_water_age_since_surface_contact" /> … … 2082 2054 <!-- <field field_ref="DMP_Age" name="DMP_Age" standard_name="sea_water_age_since_surface_contact" /> --> 2083 2055 <field field_ref="SMS_Age" name="SMS_Age" standard_name="sea_water_age_since_surface_contact" /> 2084 <field field_ref="SMS_Age_INV" name="SMS_Age_INV" standard_name="sea_water_age_since_surface_contact" />2085 2056 <field field_ref="ATF_Age" name="ATF_Age" standard_name="sea_water_age_since_surface_contact" /> 2086 2057 <field field_ref="ISO_Age" name="ISO_Age" standard_name="sea_water_age_since_surface_contact" /> 2087 2058 </field_group> 2088 2059 2089 <field_group id="trendIDTRA " >2060 <field_group id="trendIDTRA_long" > 2090 2061 <!-- Ideal tracer --> 2091 2062 <field field_ref="TOT_IDTRA" name="TOT_IDTRA" standard_name="IDTRA_Concentration" /> 2092 <field field_ref="TOT_IDTRA_INV" name="TOT_IDTRA_INV" standard_name="IDTRA_Concentration" />2093 2063 <field field_ref="XAD_IDTRA" name="XAD_IDTRA" standard_name="IDTRA_Concentration" /> 2094 2064 <field field_ref="YAD_IDTRA" name="YAD_IDTRA" standard_name="IDTRA_Concentration" /> 2095 2065 <field field_ref="ZAD_IDTRA" name="ZAD_IDTRA" standard_name="IDTRA_Concentration" /> 2096 2066 <field field_ref="TAD_IDTRA" name="TAD_IDTRA" standard_name="IDTRA_Concentration" /> 2097 <field field_ref="TAD_IDTRA_INV" name="TAD_IDTRA_INV" standard_name="IDTRA_Concentration" />2098 2067 <field field_ref="LDF_IDTRA" name="LDF_IDTRA" standard_name="IDTRA_Concentration" /> 2099 2068 <field field_ref="BBL_IDTRA" name="BBL_IDTRA" standard_name="IDTRA_Concentration" /> … … 2103 2072 <!-- <field field_ref="DMP_IDTRA" name="DMP_IDTRA" standard_name="IDTRA_Concentration" /> --> 2104 2073 <field field_ref="SMS_IDTRA" name="SMS_IDTRA" standard_name="IDTRA_Concentration" /> 2105 <field field_ref="SMS_IDTRA_INV" name="SMS_IDTRA_INV" standard_name="IDTRA_Concentration" />2106 2074 <field field_ref="ATF_IDTRA" name="ATF_IDTRA" standard_name="IDTRA_Concentration" /> 2107 2075 <field field_ref="ISO_IDTRA" name="ISO_IDTRA" standard_name="IDTRA_Concentration" /> 2108 2076 </field_group> 2109 2077 2110 <field_group id="trendC14B " >2078 <field_group id="trendC14B_long" > 2111 2079 <!-- c14 bomb --> 2112 2080 <field field_ref="TOT_C14B" name="TOT_C14B" standard_name="C14b_Concentration" /> 2113 <field field_ref="TOT_C14B_INV" name="TOT_C14B_INV" standard_name="C14b_Concentration" />2114 2081 <field field_ref="XAD_C14B" name="XAD_C14B" standard_name="C14b_Concentration" /> 2115 2082 <field field_ref="YAD_C14B" name="YAD_C14B" standard_name="C14b_Concentration" /> 2116 2083 <field field_ref="ZAD_C14B" name="ZAD_C14B" standard_name="C14b_Concentration" /> 2117 2084 <field field_ref="TAD_C14B" name="TAD_C14B" standard_name="C14b_Concentration" /> 2118 <field field_ref="TAD_C14B_INV" name="TAD_C14B_INV" standard_name="C14b_Concentration" />2119 2085 <field field_ref="LDF_C14B" name="LDF_C14B" standard_name="C14b_Concentration" /> 2120 2086 <field field_ref="BBL_C14B" name="BBL_C14B" standard_name="C14b_Concentration" /> … … 2124 2090 <!-- <field field_ref="DMP_C14B" name="DMP_C14B" standard_name="C14b_Concentration" /> --> 2125 2091 <field field_ref="SMS_C14B" name="SMS_C14B" standard_name="C14b_Concentration" /> 2126 <field field_ref="SMS_C14B_INV" name="SMS_C14B_INV" standard_name="C14b_Concentration" />2127 2092 <field field_ref="ATF_C14B" name="ATF_C14B" standard_name="C14b_Concentration" /> 2128 2093 <field field_ref="ISO_C14B" name="ISO_C14B" standard_name="C14b_Concentration" /> … … 2161 2126 2162 2127 2163 <field_group id="trendCFC11_ short" >2128 <field_group id="trendCFC11_cmip6" > 2164 2129 <!-- cfc11 --> 2165 2130 <field field_ref="TOT_CFC11_E3T" name="TOT_CFC11_E3T" standard_name="CFC-11_Concentration" /> 2131 <field field_ref="SMS_CFC11_E3T" name="SMS_CFC11_E3T" standard_name="CFC-11_Concentration" /> 2132 </field_group> 2133 2134 <field_group id="trendCFC12_cmip6" > 2135 <!-- cfc12 --> 2136 <field field_ref="TOT_CFC12_E3T" name="TOT_CFC12_E3T" standard_name="CFC-12_Concentration" /> 2137 <field field_ref="SMS_CFC12_E3T" name="SMS_CFC12_E3T" standard_name="CFC-12_Concentration" /> 2138 </field_group> 2139 2140 <field_group id="trendSF6_cmip6" > 2141 <!-- sf6 --> 2142 <field field_ref="TOT_SF6_E3T" name="TOT_SF6_E3T" standard_name="SF6_Concentration" /> 2143 <field field_ref="SMS_SF6_E3T" name="SMS_SF6_E3T" standard_name="SF6_Concentration" /> 2144 </field_group> 2145 2146 <field_group id="trendAGE_cmip6" > 2147 <!-- age --> 2148 <field field_ref="TOT_Age_E3T" name="TOT_Age_E3T" standard_name="sea_water_age_since_surface_contact" /> 2149 <field field_ref="SMS_Age_E3T" name="SMS_Age_E3T" standard_name="sea_water_age_since_surface_contact" /> 2150 </field_group> 2151 2152 <field_group id="trendIDTRA_cmip6" > 2153 <!-- Ideal tracer --> 2154 <field field_ref="TOT_IDTRA_E3T" name="TOT_IDTRA_E3T" standard_name="IDTRA_Concentration" /> 2155 <field field_ref="SMS_IDTRA_E3T" name="SMS_IDTRA_E3T" standard_name="IDTRA_Concentration" /> 2156 </field_group> 2157 2158 <field_group id="trendC14B_cmip6" > 2159 <!-- c14 bomb --> 2160 <field field_ref="TOT_C14B_E3T" name="TOT_C14B_E3T" standard_name="C14b_Concentration" /> 2161 <field field_ref="SMS_C14B_E3T" name="SMS_C14B_E3T" standard_name="C14b_Concentration" /> 2162 </field_group> 2163 2164 <field_group id="trendCFC11_cmip6_inv" > 2165 <!-- cfc11 --> 2166 2166 <field field_ref="TOT_CFC11_INV" name="TOT_CFC11_INV" standard_name="CFC-11_Concentration" /> 2167 <field field_ref="SMS_CFC11_E3T" name="SMS_CFC11_E3T" standard_name="CFC-11_Concentration" />2168 2167 <field field_ref="SMS_CFC11_INV" name="SMS_CFC11_INV" standard_name="CFC-11_Concentration" /> 2169 2168 </field_group> 2170 2169 2171 <field_group id="trendCFC12_ short" >2170 <field_group id="trendCFC12_cmip6_inv" > 2172 2171 <!-- cfc12 --> 2173 <field field_ref="TOT_CFC12_E3T" name="TOT_CFC12_E3T" standard_name="CFC-12_Concentration" />2174 2172 <field field_ref="TOT_CFC12_INV" name="TOT_CFC12_INV" standard_name="CFC-12_Concentration" /> 2175 <field field_ref="SMS_CFC12_E3T" name="SMS_CFC12_E3T" standard_name="CFC-12_Concentration" />2176 2173 <field field_ref="SMS_CFC12_INV" name="SMS_CFC12_INV" standard_name="CFC-12_Concentration" /> 2177 2174 </field_group> 2178 2175 2179 <field_group id="trendSF6_ short" >2176 <field_group id="trendSF6_cmip6_inv" > 2180 2177 <!-- sf6 --> 2181 <field field_ref="TOT_SF6_E3T" name="TOT_SF6_E3T" standard_name="SF6_Concentration" />2182 2178 <field field_ref="TOT_SF6_INV" name="TOT_SF6_INV" standard_name="SF6_Concentration" /> 2183 <field field_ref="SMS_SF6_E3T" name="SMS_SF6_E3T" standard_name="SF6_Concentration" />2184 2179 <field field_ref="SMS_SF6_INV" name="SMS_SF6_INV" standard_name="SF6_Concentration" /> 2185 2180 </field_group> 2186 2181 2187 <field_group id="trendAGE_ short" >2182 <field_group id="trendAGE_cmip6_inv" > 2188 2183 <!-- age --> 2189 <field field_ref="TOT_Age_E3T" name="TOT_Age_E3T" standard_name="sea_water_age_since_surface_contact" />2190 2184 <field field_ref="TOT_Age_INV" name="TOT_Age_INV" standard_name="sea_water_age_since_surface_contact" /> 2191 <field field_ref="SMS_Age_E3T" name="SMS_Age_E3T" standard_name="sea_water_age_since_surface_contact" />2192 2185 <field field_ref="SMS_Age_INV" name="SMS_Age_INV" standard_name="sea_water_age_since_surface_contact" /> 2193 2186 </field_group> 2194 2187 2195 <field_group id="trendIDTRA_ short" >2188 <field_group id="trendIDTRA_cmip6_inv" > 2196 2189 <!-- Ideal tracer --> 2197 <field field_ref="TOT_IDTRA_E3T" name="TOT_IDTRA_E3T" standard_name="IDTRA_Concentration" />2198 2190 <field field_ref="TOT_IDTRA_INV" name="TOT_IDTRA_INV" standard_name="IDTRA_Concentration" /> 2199 <field field_ref="SMS_IDTRA_E3T" name="SMS_IDTRA_E3T" standard_name="IDTRA_Concentration" />2200 2191 <field field_ref="SMS_IDTRA_INV" name="SMS_IDTRA_INV" standard_name="IDTRA_Concentration" /> 2201 2192 </field_group> 2202 2193 2203 <field_group id="trendC14B_ short" >2194 <field_group id="trendC14B_cmip6_inv" > 2204 2195 <!-- c14 bomb --> 2205 <field field_ref="TOT_C14B_E3T" name="TOT_C14B_E3T" standard_name="C14b_Concentration" />2206 2196 <field field_ref="TOT_C14B_INV" name="TOT_C14B_INV" standard_name="C14b_Concentration" /> 2207 <field field_ref="SMS_C14B_E3T" name="SMS_C14B_E3T" standard_name="C14b_Concentration" />2208 2197 <field field_ref="SMS_C14B_INV" name="SMS_C14B_INV" standard_name="C14b_Concentration" /> 2209 2198 </field_group> … … 2421 2410 </field_group> 2422 2411 2412 <field_group id="groupMEDUSA_cmip6_duplicates" > 2413 <field field_ref= "FBDDTALK" name="FBDDTALK" /> 2414 <field field_ref= "FBDDTDIC" name="FBDDTDIC" /> 2415 <field field_ref= "FBDDTDIFE" name="FBDDTDIFE" /> 2416 <field field_ref= "FBDDTDIN" name="FBDDTDIN" /> 2417 <field field_ref= "FBDDTDISI" name="FBDDTDISI" /> 2418 <field field_ref= "BDDTALK3" name="BDDTALK3" /> 2419 <field field_ref= "BDDTDIC3" name="BDDTDIC3" /> 2420 <field field_ref= "BDDTDIFE3" name="BDDTDIFE3" /> 2421 <field field_ref= "BDDTDIN3" name="BDDTDIN3" /> 2422 <field field_ref= "BDDTDISI3" name="BDDTDISI3" /> 2423 </field_group> 2424 2423 2425 <field_group id="groupMEDUSA_cmip6" > 2424 2426 <field field_ref= "epC100" name="epC100" /> … … 2432 2434 <field field_ref= "O2min" name="O2min" /> 2433 2435 <field field_ref= "ZO2min" name="ZO2min" /> 2434 <field field_ref= "FBDDTALK" name="FBDDTALK" />2435 <field field_ref= "FBDDTDIC" name="FBDDTDIC" />2436 <field field_ref= "FBDDTDIFE" name="FBDDTDIFE" />2437 <field field_ref= "FBDDTDIN" name="FBDDTDIN" />2438 <field field_ref= "FBDDTDISI" name="FBDDTDISI" />2439 2436 <field field_ref= "TPPD3" name="TPPD3" /> 2440 <field field_ref= "BDDTALK3" name="BDDTALK3" />2441 <field field_ref= "BDDTDIC3" name="BDDTDIC3" />2442 <field field_ref= "BDDTDIFE3" name="BDDTDIFE3" />2443 <field field_ref= "BDDTDIN3" name="BDDTDIN3" />2444 <field field_ref= "BDDTDISI3" name="BDDTDISI3" />2445 2437 <field field_ref= "FD_NIT3" name="FD_NIT3" /> 2446 2438 <field field_ref= "FD_SIL3" name="FD_SIL3" /> -
branches/NERC/dev_r5518_GO6_rev9312_xCO2/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90
r9257 r9351 153 153 IF (tmask(ji,jj,1) == 1) then 154 154 IF (lk_oasis) THEN 155 !! use 2D atm xCO2 from atm coupling156 f_xco2a(ji,jj) = PCO2a_in_cpl(ji,jj)157 !!!158 155 !!! Jpalm test on atm xCO2 159 IF ( (f_xco2a(ji,jj) > 1500.0 ).OR.(f_xco2a(ji,jj) < 100.0 ) ) THEN156 IF ( (f_xco2a(ji,jj) .GT. 10000.0 ).OR.(f_xco2a(ji,jj) .LE. 0.0 ) ) THEN 160 157 IF(lwp) THEN 161 158 WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj), & 162 159 ' -- ji =', mig(ji),' jj = ', mjg(jj) 163 160 ENDIF 164 CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic atm xCO2 ' )161 CALL ctl_stop( 'MEDUSA - Air-Sea :', 'unrealistic coupled atm xCO2 ' ) 165 162 ENDIF 166 163 ENDIF … … 261 258 !! must be associated to air-sea flux or air xCO2... 262 259 !! Check MOCSY outputs 263 IF ( (f_pco2w(ji,jj) > 1.E4 ).OR.(f_pco2w(ji,jj) < 0.0 ) .OR. & 264 (f_fco2w(ji,jj) > 1.E4 ).OR.(f_fco2w(ji,jj) < 0.0 ) .OR. & 265 (f_fco2atm(ji,jj) > 1.E4 ).OR.(f_fco2atm(ji,jj) < 0.0 ) .OR. & 266 (f_co2flux(ji,jj) > 1.E-2 ).OR.(f_co2flux(ji,jj) < -1.E-2 ) .OR. & 267 (f_dpco2(ji,jj) > 1.E4 ).OR.(f_dpco2(ji,jj) < -1.E4 ) ) THEN 268 IF(lwp) THEN 269 WRITE(numout,*) ' surface T = ',ztmp(ji,jj) 270 WRITE(numout,*) ' surface S = ',zsal(ji,jj) 271 WRITE(numout,*) ' surface ALK = ',zalk(ji,jj) 272 WRITE(numout,*) ' surface DIC = ',zdic(ji,jj) 273 WRITE(numout,*) ' KW660 = ',f_kw660(ji,jj) 274 WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj) 275 WRITE(numout,*) ' surface pco2w = ',f_pco2w(ji,jj) 276 WRITE(numout,*) ' surface fco2w = ',f_fco2w(ji,jj) 277 WRITE(numout,*) ' surface fco2a = ',f_fco2atm(ji,jj) 278 WRITE(numout,*) ' surface co2flx = ',f_co2flux(ji,jj) 279 WRITE(numout,*) ' surface dpco2 = ',f_dpco2(ji,jj) 280 WRITE(numout,*) ' MOCSY output: ji =', mig(ji),' jj = ', mjg(jj), & 281 ' kt = ', kt 282 WRITE(numout,*) 'MEDUSA - Air-Sea OUTPUT: unrealistic surface Carb. Chemistry' 283 ENDIF 284 CALL ctl_stop( 'MEDUSA - Air-Sea OUTPUT: ', & 285 'unrealistic surface Carb. Chemistry -- OUTPUTS' ) 286 ENDIF 260 !!=================== 261 !! Jpalm -- 19-02-2018 -- remove the cap - only check MOCSY inputs 262 !! because of specific area in arabic sea where strangely 263 !! with core 2 forcing, ALK is lower than DIC and result in 264 !! Enormous dpco2 - even if all carb chem caract are OK. 265 !! and this check stops the model. 266 !! ==> let's run the model without it and see how it goes. 267 !!=================== 268 !!IF ( (f_pco2w(ji,jj) > 1.E4 ).OR.(f_pco2w(ji,jj) < 0.0 ) .OR. & 269 !! (f_fco2w(ji,jj) > 1.E4 ).OR.(f_fco2w(ji,jj) < 0.0 ) .OR. & 270 !! (f_fco2atm(ji,jj) > 1.E4 ).OR.(f_fco2atm(ji,jj) < 0.0 ) .OR. & 271 !! (f_co2flux(ji,jj) > 1.E-1 ).OR.(f_co2flux(ji,jj) < -1.E-1 ) .OR. & 272 !! (f_dpco2(ji,jj) > 1.E4 ).OR.(f_dpco2(ji,jj) < -1.E4 ) ) THEN 273 !! IF(lwp) THEN 274 !! WRITE(numout,*) ' surface T = ',ztmp(ji,jj) 275 !! WRITE(numout,*) ' surface S = ',zsal(ji,jj) 276 !! WRITE(numout,*) ' surface ALK = ',zalk(ji,jj) 277 !! WRITE(numout,*) ' surface DIC = ',zdic(ji,jj) 278 !! WRITE(numout,*) ' KW660 = ',f_kw660(ji,jj) 279 !! WRITE(numout,*) ' atm xCO2 = ',f_xco2a(ji,jj) 280 !! WRITE(numout,*) ' surface pco2w = ',f_pco2w(ji,jj) 281 !! WRITE(numout,*) ' surface fco2w = ',f_fco2w(ji,jj) 282 !! WRITE(numout,*) ' surface fco2a = ',f_fco2atm(ji,jj) 283 !! WRITE(numout,*) ' surface co2flx = ',f_co2flux(ji,jj) 284 !! WRITE(numout,*) ' surface dpco2 = ',f_dpco2(ji,jj) 285 !! WRITE(numout,*) ' MOCSY output: ji =', mig(ji),' jj = ', mjg(jj), & 286 !! ' kt = ', kt 287 !! WRITE(numout,*) 'MEDUSA - Air-Sea OUTPUT: unrealistic surface Carb. Chemistry' 288 !! ENDIF 289 !! CALL ctl_stop( 'MEDUSA - Air-Sea OUTPUT: ', & 290 !! 'unrealistic surface Carb. Chemistry -- OUTPUTS' ) 291 !!ENDIF 287 292 ENDIF 288 293 ENDDO -
branches/NERC/dev_r5518_GO6_rev9312_xCO2/NEMOGCM/NEMO/TOP_SRC/MEDUSA/carb_chem.F90
r8441 r9351 103 103 !! OPEN wet point IF..THEN loop 104 104 IF (tmask(ji,jj,jk).eq.1) THEN 105 IF (lk_oasis) THEN106 !! use 2D atm xCO2 from atm coupling107 f_xco2a(ji,jj) = PCO2a_in_cpl(ji,jj)108 ENDIF109 105 !! Do carbonate chemistry 110 106 !! -
branches/NERC/dev_r5518_GO6_rev9312_xCO2/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90
r9258 r9351 283 283 #if defined key_roam 284 284 !!---------------------------------------------------------------------- 285 !! Atmospheric pCO2 data (1859 to 2100 inclusive) 286 !!---------------------------------------------------------------------- 287 !! 288 REAL(wp), DIMENSION(242) :: hist_pco2 !: pCO2 289 290 # if defined key_rcp26 291 !! UKMO, run AJKKH + KAAEC, RCP 2.6, pCO2 time evolution 292 DATA hist_pco2 / 286.0230, 286.1730, 286.3230, 286.4480, 286.5730, & 293 & 286.7230, 286.8480, 286.9480, 287.0480, 287.1730, & 294 & 287.3230, 287.4730, 287.6480, 287.8480, 288.0730, & 295 & 288.3480, 288.6480, 288.9730, 289.3470, 289.7470, & 296 & 290.1730, 290.6470, 291.1470, 291.6220, 292.0720, & 297 & 292.5220, 292.9220, 293.2470, 293.5220, 293.7470, & 298 & 293.9470, 294.1220, 294.2720, 294.4220, 294.5470, & 299 & 294.6470, 294.7470, 294.8470, 294.9710, 295.1710, & 300 & 295.4460, 295.7470, 296.0720, 296.4210, 296.7710, & 301 & 297.1460, 297.5710, 298.0210, 298.4460, 298.8460, & 302 & 299.2460, 299.6450, 300.0210, 300.3710, 300.7200, & 303 & 301.0450, 301.3460, 301.6710, 302.0200, 302.3450, & 304 & 302.6450, 302.9700, 303.3450, 303.7200, 304.0700, & 305 & 304.4700, 304.9200, 305.3440, 305.7700, 306.2450, & 306 & 306.7190, 307.1700, 307.6440, 308.1190, 308.5440, & 307 & 308.9440, 309.3440, 309.6940, 309.9440, 310.1190, & 308 & 310.2440, 310.3190, 310.3190, 310.2440, 310.1440, & 309 & 310.0690, 310.0440, 310.0690, 310.1440, 310.2690, & 310 & 310.4440, 310.6940, 311.0430, 311.4440, 311.8690, & 311 & 312.3680, 312.9430, 313.5430, 314.1680, 314.7900, & 312 & 315.4430, 316.2150, 317.0170, 317.7370, 318.3400, & 313 & 318.8680, 319.5900, 320.5890, 321.5470, 322.5770, & 314 & 323.8440, 324.9260, 325.7960, 327.0810, 328.6180, & 315 & 329.6830, 330.5250, 331.6880, 333.2120, 334.7870, & 316 & 336.4640, 338.2990, 339.6660, 340.7310, 342.1360, & 317 & 343.7200, 345.2200, 346.7350, 348.5820, 350.6740, & 318 & 352.4230, 353.7910, 354.9530, 355.8210, 356.7130, & 319 & 358.0630, 359.7720, 361.3970, 363.0900, 365.2560, & 320 & 367.2810, 368.7980, 370.4000, 372.4550, 374.6920, & 321 & 376.7440, 378.7440, 380.7580, 382.7080, 384.7300, & 322 & 386.9310, 389.2150, 391.4910, 393.7710, 396.0460, & 323 & 398.3240, 400.6080, 402.8950, 405.1780, 407.4550, & 324 & 409.7260, 411.9930, 414.2500, 416.4410, 418.5280, & 325 & 420.5250, 422.4390, 424.2720, 426.0200, 427.6750, & 326 & 429.2360, 430.7050, 432.0850, 433.3580, 434.5140, & 327 & 435.5740, 436.5490, 437.4420, 438.2550, 438.9810, & 328 & 439.6110, 440.1430, 440.5770, 440.9450, 441.2660, & 329 & 441.5410, 441.7840, 442.0050, 442.2040, 442.3780, & 330 & 442.5210, 442.6200, 442.6720, 442.6810, 442.6540, & 331 & 442.5830, 442.4670, 442.3270, 442.1680, 441.9960, & 332 & 441.8060, 441.5930, 441.3440, 441.0540, 440.7230, & 333 & 440.3510, 439.9300, 439.4650, 438.9730, 438.4630, & 334 & 437.9400, 437.4020, 436.8400, 436.2640, 435.6850, & 335 & 435.1030, 434.5160, 433.9170, 433.3060, 432.7010, & 336 & 432.1110, 431.5380, 430.9810, 430.4320, 429.8860, & 337 & 429.3370, 428.7810, 428.2220, 427.6490, 427.0660, & 338 & 426.4890, 425.9270, 425.3840, 424.8610, 424.3540, & 339 & 423.8540, 423.3540, 422.8530, 422.3510, 421.8410, & 340 & 421.3250, 420.8190 / 341 # else 342 !! UKMO, run AJKKH + KAAEF, RCP 8.5, pCO2 time evolution 343 DATA hist_pco2 / 286.0230, 286.1730, 286.3230, 286.4480, 286.5730, & 344 & 286.7230, 286.8480, 286.9480, 287.0480, 287.1730, & 345 & 287.3230, 287.4730, 287.6480, 287.8480, 288.0730, & 346 & 288.3480, 288.6480, 288.9730, 289.3470, 289.7470, & 347 & 290.1730, 290.6470, 291.1470, 291.6220, 292.0720, & 348 & 292.5220, 292.9220, 293.2470, 293.5220, 293.7470, & 349 & 293.9470, 294.1220, 294.2720, 294.4220, 294.5470, & 350 & 294.6470, 294.7470, 294.8470, 294.9710, 295.1710, & 351 & 295.4460, 295.7470, 296.0720, 296.4210, 296.7710, & 352 & 297.1460, 297.5710, 298.0210, 298.4460, 298.8460, & 353 & 299.2460, 299.6450, 300.0210, 300.3710, 300.7200, & 354 & 301.0450, 301.3460, 301.6710, 302.0200, 302.3450, & 355 & 302.6450, 302.9700, 303.3450, 303.7200, 304.0700, & 356 & 304.4700, 304.9200, 305.3440, 305.7700, 306.2450, & 357 & 306.7190, 307.1700, 307.6440, 308.1190, 308.5440, & 358 & 308.9440, 309.3440, 309.6940, 309.9440, 310.1190, & 359 & 310.2440, 310.3190, 310.3190, 310.2440, 310.1440, & 360 & 310.0690, 310.0440, 310.0690, 310.1440, 310.2690, & 361 & 310.4440, 310.6940, 311.0430, 311.4440, 311.8690, & 362 & 312.3680, 312.9430, 313.5430, 314.1680, 314.7900, & 363 & 315.4430, 316.2150, 317.0170, 317.7370, 318.3400, & 364 & 318.8680, 319.5900, 320.5890, 321.5470, 322.5770, & 365 & 323.8440, 324.9260, 325.7960, 327.0810, 328.6180, & 366 & 329.6830, 330.5250, 331.6880, 333.2120, 334.7870, & 367 & 336.4640, 338.2990, 339.6660, 340.7310, 342.1360, & 368 & 343.7200, 345.2200, 346.7350, 348.5820, 350.6740, & 369 & 352.4230, 353.7910, 354.9530, 355.8210, 356.7130, & 370 & 358.0630, 359.7720, 361.3970, 363.0900, 365.2560, & 371 & 367.2810, 368.7980, 370.4000, 372.4550, 374.6920, & 372 & 376.7440, 378.7440, 380.7580, 382.7080, 384.7300, & 373 & 386.9420, 389.2540, 391.5670, 393.9370, 396.3920, & 374 & 398.9320, 401.5550, 404.2550, 407.0220, 409.8530, & 375 & 412.7470, 415.7050, 418.7210, 421.7880, 424.9180, & 376 & 428.1200, 431.3970, 434.7470, 438.1650, 441.6410, & 377 & 445.1700, 448.7530, 452.3920, 456.0950, 459.8810, & 378 & 463.7680, 467.7660, 471.8750, 476.0960, 480.4210, & 379 & 484.8390, 489.3470, 493.9430, 498.6400, 503.4380, & 380 & 508.3410, 513.3630, 518.5160, 523.8050, 529.2290, & 381 & 534.7780, 540.4450, 546.2230, 552.1120, 558.1110, & 382 & 564.2110, 570.4130, 576.7390, 583.1990, 589.7980, & 383 & 596.5390, 603.4110, 610.4060, 617.4940, 624.6500, & 384 & 631.8800, 639.1750, 646.5360, 653.9800, 661.5230, & 385 & 669.1840, 676.9570, 684.8290, 692.7790, 700.7690, & 386 & 708.8050, 716.8870, 725.0020, 733.1770, 741.3900, & 387 & 749.6700, 758.0480, 766.5050, 775.0350, 783.6110, & 388 & 792.2200, 800.8740, 809.5680, 818.2760, 827.0090, & 389 & 835.8020, 844.6550, 853.5730, 862.5690, 871.6190, & 390 & 880.7020, 889.8240, 898.9590, 908.1270, 917.3080, & 391 & 926.4960, 935.7040 / 392 # endif 285 !! JPALM -- change hist_pco2 init 286 !!---------------------------------------------------------------------- 287 !! 288 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) :: hist_pco2 !: pCO2 289 INTEGER :: co2_rec 290 REAL(wp) :: co2_yinit, co2_yend 291 #endif 292 293 !!---------------------------------------------------------------------- 294 !! JPALM -- PI CO2 key 295 !!---------------------------------------------------------------------- 296 !! 297 #if defined key_axy_pi_co2 298 LOGICAL , PUBLIC :: lk_pi_co2 = .TRUE. !: PI xCO2 used 299 #else 300 LOGICAL , PUBLIC :: lk_pi_co2 = .FALSE. !: PI xCO2 unused 393 301 #endif 394 302 -
branches/NERC/dev_r5518_GO6_rev9312_xCO2/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r9257 r9351 81 81 gdept_0, gdept_n, & 82 82 gdepw_0, gdepw_n, & 83 nday_year, nsec_day, nyear, & 83 nday_year, nsec_day, & 84 nyear, nyear_len, ndastp, & 85 nsec_month, & 84 86 rdt, tmask, mig, mjg, nimpp, & 85 87 njmpp … … 92 94 mpp_min, mpp_minloc, & 93 95 ctl_stop, ctl_warn, lk_mpp 94 USE oce, ONLY: tsb, tsn 96 USE oce, ONLY: tsb, tsn, PCO2a_in_cpl 95 97 USE par_kind, ONLY: wp 96 98 USE par_medusa, ONLY: jpalk, jpchd, jpchn, jpdet, & … … 102 104 !! JPALM (27-06-2016): add lk_oasis for CO2 and DMS coupling with atm 103 105 USE sbc_oce, ONLY: lk_oasis 104 USE sms_medusa, ONLY: hist_pco2 106 USE sms_medusa, ONLY: hist_pco2, co2_yinit, co2_yend, & 107 lk_pi_co2 105 108 USE trc, ONLY: ln_rsttr, nittrc000, trn 106 109 USE bio_medusa_init_mod, ONLY: bio_medusa_init … … 114 117 USE bio_medusa_diag_slice_mod, ONLY: bio_medusa_diag_slice 115 118 USE bio_medusa_fin_mod, ONLY: bio_medusa_fin 119 USE trcstat, ONLY: trc_rst_dia_stat 116 120 117 121 IMPLICIT NONE … … 181 185 !! 182 186 !! temporary variables 183 REAL(wp) :: fq0,fq1,fq2,fq3,fq4 187 REAL(wp) :: fq3,fq4 188 REAL(wp) :: this_y, this_d, this_s, fyear 184 189 !! 185 190 !! T and S check temporary variable … … 293 298 !!------------------------------------------------------------------ 294 299 !! 295 !! what's atmospheric pCO2 doing? (data start in 1859) 296 iyr1 = nyear - 1859 + 1 297 iyr2 = iyr1 + 1 298 if (iyr1 .le. 1) then 299 !! before 1860 300 f_xco2a(:,:) = hist_pco2(1) 301 elseif (iyr2 .ge. 242) then 302 !! after 2099 303 f_xco2a(:,:) = hist_pco2(242) 304 else 305 !! just right 306 fq0 = hist_pco2(iyr1) 307 fq1 = hist_pco2(iyr2) 308 fq2 = real(nsec_day) / (60.0 * 60.0 * 24.0) 309 !! AXY (14/06/12): tweaked to make more sense (and be correct) 310 # if defined key_bs_axy_yrlen 311 !! bugfix: for 360d year with HadGEM2-ES forcing 312 fq3 = (real(nday_year) - 1.0 + fq2) / 360.0 313 # else 314 !! original use of 365 days (not accounting for leap year or 315 !! 360d year) 316 fq3 = (real(nday_year) - 1.0 + fq2) / 365.0 317 # endif 318 fq4 = (fq0 * (1.0 - fq3)) + (fq1 * fq3) 319 f_xco2a(:,:) = fq4 320 endif 321 # if defined key_axy_pi_co2 322 !! OCMIP pre-industrial pCO2 323 !! f_xco2a(:,:) = 284.725 !! CMIP5 pre-industrial pCO2 324 f_xco2a = 284.317 !! CMIP6 pre-industrial pCO2 325 # endif 326 !! IF(lwp) WRITE(numout,*) ' MEDUSA nyear =', nyear 327 !! IF(lwp) WRITE(numout,*) ' MEDUSA nsec_day =', real(nsec_day) 328 !! IF(lwp) WRITE(numout,*) ' MEDUSA nday_year =', real(nday_year) 329 !! AXY (29/01/14): remove surplus diagnostics 330 !! IF(lwp) WRITE(numout,*) ' MEDUSA fq0 =', fq0 331 !! IF(lwp) WRITE(numout,*) ' MEDUSA fq1 =', fq1 332 !! IF(lwp) WRITE(numout,*) ' MEDUSA fq2 =', fq2 333 !! IF(lwp) WRITE(numout,*) ' MEDUSA fq3 =', fq3 334 IF(lwp) WRITE(numout,*) ' MEDUSA atm pCO2 =', f_xco2a(1,1) 300 IF (lk_oasis) THEN 301 !! xCO2 from coupled 302 IF ( ( kt == nittrc000 ) .AND. lwp ) & 303 WRITE(numout,*) '** MEDUSA Atm xCO2 given by the UM **' 304 f_xco2a(:,:) = PCO2a_in_cpl(:,:) 305 ELSEIF (lk_pi_co2) THEN 306 !! OCMIP pre-industrial xCO2 307 IF ( ( kt == nittrc000 ) .AND. lwp ) & 308 WRITE(numout,*) '** MEDUSA Atm xCO2 fixed to pre-industrial value **' 309 !! f_xco2a(:,:) = 284.725 !! CMIP5 pre-industrial pCO2 310 f_xco2a(:,:) = 284.317 !! CMIP6 pre-industrial pCO2 311 ELSE 312 !! xCO2 from file 313 !! AXY - JPALM new interpolation scheme usinf nyear_len 314 this_y = real(nyear) 315 this_d = real(nday_year) 316 this_s = real(nsec_day) 317 !! 318 fyear = this_y + ((this_d - 1) + (this_s / (60. * 60. * 24.))) / real(nyear_len(1)) 319 !! 320 IF ( ( kt == nittrc000 ) .AND. lwp ) THEN 321 WRITE(numout,*) '** MEDUSA Atm xCO2 from file **' 322 WRITE(numout,*) ' MEDUSA year =', this_y 323 WRITE(numout,*) ' Year length =', real(nyear_len(1)) 324 WRITE(numout,*) ' MEDUSA nday_year =', this_d 325 WRITE(numout,*) ' MEDUSA nsec_day =', this_s 326 ENDIF 327 !! 328 !! different case test 329 IF (fyear .LE. co2_yinit) THEN 330 !! before first record -- pre-industrial value 331 f_xco2a(:,:) = hist_pco2(1) 332 ELSEIF (fyear .GE. co2_yend) THEN 333 !! after last record - continue to use the last value 334 f_xco2a(:,:) = hist_pco2(int(co2_yend - co2_yinit + 1.) ) 335 ELSE 336 !! just right 337 iyr1 = int(fyear - co2_yinit) + 1 338 iyr2 = iyr1 + 1 339 fq3 = fyear - real(iyr1) - co2_yinit + 1. 340 fq4 = ((1 - fq3) * hist_pco2(iyr1)) + (fq3 * hist_pco2(iyr2)) 341 f_xco2a(:,:) = fq4 342 !! 343 IF ( ( kt == nittrc000 ) .AND. lwp ) THEN 344 WRITE(numout,*) ' MEDUSA year1 =', iyr1 345 WRITE(numout,*) ' MEDUSA year2 =', iyr2 346 WRITE(numout,*) ' xCO2 year1 =', hist_pco2(iyr1) 347 WRITE(numout,*) ' xCO2 year2 =', hist_pco2(iyr2) 348 WRITE(numout,*) ' Year2 weight =', fq3 349 ENDIF 350 ENDIF 351 ENDIF 352 353 !! Writing xCO2 in output on start and on the 1st tsp of each month 354 IF ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR. & 355 ( nsec_month .LE. INT(rdt) ) ) THEN 356 IF ( lwp ) WRITE(numout,*) ' *** Atm xCO2 *** -- kt:', kt, & 357 '; current date:', ndastp 358 call trc_rst_dia_stat(f_xco2a(:,:), 'atm xCO2') 359 ENDIF 335 360 # endif 336 361 … … 358 383 !! x * 30d + 1*rdt(i.e: mod = rdt) 359 384 !! ++ need to pass carb-chem output var through restarts 360 If ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR. & 361 ((86400*mod(nn_date0,100) + mod(kt*rdt,2592000.)) == rdt) ) THEN 385 !!If ( (kt == nitt8rc000 .AND. .NOT.ln_rsttr) .OR. & 386 !! ( (mod(kt*rdt,2592000.)) == rdt) THEN 387 !!============================= 388 !! (Jpalm -- updated for restartability issues) 389 !! We want this to be start of month or if starting afresh from 390 !! climatology - marc 20/6/17 391 !!If ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR. & 392 !! ((86400*mod(nn_date0,100) + mod(kt*rdt,2592000.)) == rdt) ) THEN 393 !!============================= 394 !! Jpalm -- 15-02-2018 -- need to change 3D carb-chem call freq again. 395 !! previous call did not work, probably the (86400*mod(nn_date0,100) part 396 !! should not be in... 397 !! now use the NEMO calendar tool : nsec_month to be sure to call 398 !! at the beginning of a new month . 399 IF ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR. & 400 ( nsec_month .LE. INT(rdt) ) ) THEN 401 IF ( lwp ) WRITE(numout,*) & 402 ' *** 3D carb chem call *** -- kt:', kt, & 403 '; current date:', ndastp 362 404 !!--------------------------------------------------------------- 363 405 !! Calculate the carbonate chemistry for the whole ocean on the first … … 502 544 !! Exceptionnal value did exist 503 545 !! 504 Call trc_bio_check(kt )546 Call trc_bio_check(kt, jk) 505 547 506 548 !!================================================================ … … 810 852 END SUBROUTINE trc_bio_exceptionnal_fix 811 853 812 SUBROUTINE trc_bio_check(kt )854 SUBROUTINE trc_bio_check(kt, jk) 813 855 !!----------------------------------- 814 856 !! JPALM -- 14-12-2017 -- Still dealing with this micro-boil/carb failure … … 824 866 INTEGER :: ii,ij ! temporary scalars 825 867 INTEGER, DIMENSION(2) :: ilocs ! 826 INTEGER, INTENT( in ) :: kt 868 INTEGER, INTENT( in ) :: kt, jk 827 869 !! 828 870 !!========================== … … 852 894 IF(lwp) THEN 853 895 WRITE(numout,*) 'trc_bio:tracer anomaly: ***** WARNING *****' 854 WRITE(numout,*) 'trc_bio:tracer anomaly: sea surface DIC> 4000 '855 WRITE(numout,9600) kt, zmax, ii, ij 896 WRITE(numout,*) 'trc_bio:tracer anomaly: DIC concentration > 4000 ' 897 WRITE(numout,9600) kt, zmax, ii, ij, jk 856 898 WRITE(numout,*) 'trc_bio:tracer anomaly: ***** END OF WARNING *****' 857 899 ENDIF … … 869 911 IF(lwp) THEN 870 912 WRITE(numout,*) 'trc_bio:tracer anomaly: ***** WARNING *****' 871 WRITE(numout,*) 'trc_bio:tracer anomaly: sea surface DIC<= 0 '872 WRITE(numout,9700) kt, zmin, ii, ij 913 WRITE(numout,*) 'trc_bio:tracer anomaly: DIC concentration <= 0 ' 914 WRITE(numout,9700) kt, zmin, ii, ij, jk 873 915 WRITE(numout,*) 'trc_bio:tracer anomaly: ***** END OF WARNING *****' 874 916 ENDIF … … 901 943 IF(lwp) THEN 902 944 WRITE(numout,*) 'trc_bio:tracer anomaly: ***** WARNING *****' 903 WRITE(numout,*) 'trc_bio:tracer anomaly: sea surface Alkalinity> 4000 '904 WRITE(numout,9800) kt, zmax, ii, ij 945 WRITE(numout,*) 'trc_bio:tracer anomaly: ALK concentration > 4000 ' 946 WRITE(numout,9800) kt, zmax, ii, ij, jk 905 947 WRITE(numout,*) 'trc_bio:tracer anomaly: ***** END OF WARNING *****' 906 948 ENDIF … … 918 960 IF(lwp) THEN 919 961 WRITE(numout,*) 'trc_bio:tracer anomaly: ***** WARNING *****' 920 WRITE(numout,*) 'trc_bio:tracer anomaly: sea surface Alkalinity<= 0 '921 WRITE(numout,9900) kt, zmin, ii, ij 962 WRITE(numout,*) 'trc_bio:tracer anomaly: ALK concentration <= 0 ' 963 WRITE(numout,9900) kt, zmin, ii, ij, jk 922 964 WRITE(numout,*) 'trc_bio:tracer anomaly: ***** END OF WARNING *****' 923 965 ENDIF … … 925 967 926 968 927 9600 FORMAT ('trc_bio:tracer anomaly: kt=',i6,' max DIC: ',f16.10,', i j : ',2i5)928 9700 FORMAT ('trc_bio:tracer anomaly: kt=',i6,' min DIC: ',f16.10,', i j : ',2i5)929 9800 FORMAT ('trc_bio:tracer anomaly: kt=',i6,' max ALK: ',f16.10,', i j : ',2i5)930 9900 FORMAT ('trc_bio:tracer anomaly: kt=',i6,' min ALK: ',f16.10,', i j : ',2i5)969 9600 FORMAT ('trc_bio:tracer anomaly: kt=',i6,' max DIC: ',f16.10,', i j k: ',3i5) 970 9700 FORMAT ('trc_bio:tracer anomaly: kt=',i6,' min DIC: ',f16.10,', i j k: ',3i5) 971 9800 FORMAT ('trc_bio:tracer anomaly: kt=',i6,' max ALK: ',f16.10,', i j k: ',3i5) 972 9900 FORMAT ('trc_bio:tracer anomaly: kt=',i6,' min ALK: ',f16.10,', i j k: ',3i5) 931 973 932 974 END SUBROUTINE trc_bio_check -
branches/NERC/dev_r5518_GO6_rev9312_xCO2/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcini_medusa.F90
r9114 r9351 324 324 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~' 325 325 IF(lwp) CALL flush(numout) 326 !! 327 !!---------------------------------------------------------------------- 328 !! JPALM (23-01-2017): new way to initialize CO2-atm for cmip6 329 !! initially done in trcsms_medusa 330 !!---------------------------------------------------------------------- 331 !! 332 IF( ( .NOT.lk_oasis ) .AND. ( .NOT.lk_pi_co2 ) ) THEN 333 IF(lwp) WRITE(numout,*) ' trc_ini_medusa: initialisating atm CO2 record' 334 CALL trc_ini_medusa_co2atm 335 ENDIF 326 336 327 337 END SUBROUTINE trc_ini_medusa … … 480 490 END SUBROUTINE trc_ini_medusa_river 481 491 492 SUBROUTINE trc_ini_medusa_co2atm 493 !!---------------------------------------------------------------------- 494 !! *** trc_ini_medusa_co2atm *** 495 !! 496 !! ** Purpose : initialization atmospheric co2 record 497 !! 498 !! ** Method : - Read the xco2 file 499 !!---------------------------------------------------------------------- 500 INTEGER :: jn, jm, io, ierr, inum 501 INTEGER :: iskip = 4 ! number of 1st descriptor lines 502 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) :: zyy !: xCO2 record years 503 CHARACTER (len=10) :: clname = 'xco2.atm' !! atm CO2 record file 504 !!---------------------------------------------------------------------- 505 506 IF(lwp) WRITE(numout,*) 507 IF(lwp) WRITE(numout,*) ' trc_ini_medusa_co2atm: initialisation of atm CO2 historical record' 508 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~' 509 510 511 IF(lwp) WRITE(numout,*) 'read of formatted file xco2.atm' 512 513 CALL ctl_opn( inum, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 514 REWIND(inum) 515 516 ! compute the number of year in the file 517 ! file starts in 1849 do jn represent the record number in the file. 518 ! the year is jn + yinit - 1 519 jn = 1 520 DO 521 READ(inum,'(1x)',END=100) 522 jn = jn + 1 523 END DO 524 100 co2_rec = jn - 1 - iskip 525 IF ( lwp) WRITE(numout,*) ' ', co2_rec ,' years read in the file' 526 ! ! Allocate CO2 hist arrays 527 ierr = 0 528 ALLOCATE( hist_pco2(co2_rec),zyy(co2_rec), STAT=ierr ) 529 IF( ierr > 0 ) THEN 530 CALL ctl_stop( 'trc_ini_medusa_co2atm: unable to allocate array' ) ; RETURN 531 ENDIF 532 533 REWIND(inum) 534 535 DO jm = 1, iskip ! Skip over 1st six descriptor lines 536 READ(inum,'(1x)') 537 END DO 538 ! file starts in 1931 do jn represent the year in the century.jhh 539 ! Read file till the end 540 ! allocate start and end year of the file 541 jn = 1 542 DO 543 READ(inum,*, IOSTAT=io) zyy(jn), hist_pco2(jn) 544 IF( io < 0 ) exit 545 IF(jn==1) co2_yinit = zyy(jn) 546 jn = jn + 1 547 END DO 548 co2_yend = co2_yinit + real(co2_rec) - 1. 549 550 IF(lwp) THEN ! Control print 551 WRITE(numout,*) 552 WRITE(numout,*) 'CO2 hist start year: ', co2_yinit 553 WRITE(numout,*) 'CO2 hist end year: ', co2_yend 554 WRITE(numout,*) ' Year xCO2 atm ' 555 DO jn = 1, co2_rec 556 WRITE(numout, '( 5F7.1, 6F9.2)') zyy(jn), hist_pco2(jn) 557 END DO 558 ENDIF 559 560 END SUBROUTINE trc_ini_medusa_co2atm 561 562 482 563 #else 483 564 !!---------------------------------------------------------------------- -
branches/NERC/dev_r5518_GO6_rev9312_xCO2/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r9262 r9351 378 378 !! 379 379 IF(lwp) WRITE(numout,*) ' CFC averaged properties absent - setting to zero ...' 380 qint_cfc(:,:,j n) = 0.0 !! CHN380 qint_cfc(:,:,jl) = 0.0 !! CHN 381 381 ENDIF 382 382 !!
Note: See TracChangeset
for help on using the changeset viewer.
