Changeset 10057
- Timestamp:
- 2018-08-20T17:56:23+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc_v3/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90
r10055 r10057 990 990 ! Select mixed layer 991 991 IF ( ll_asmdin ) THEN 992 #if defined key_top 992 993 CALL ctl_warn( ' Doing direct initialisation with ocean colour assimilation', & 993 994 & ' Mixed layer depth taken to be background maximum mld_max_bkg' ) 994 995 zmld(:,:) = mld_max_bkg(:,:) 996 #else 997 CALL ctl_stop( ' Should not be doing biogeochemical assimilation without key_top' ) 998 #endif 995 999 ELSE 996 1000 SELECT CASE( mld_choice_bgc ) … … 1854 1858 1855 1859 IF ( ln_pno3inc ) THEN 1860 #if defined key_hadocc || ( defined key_medusa && defined key_foam_medusa ) 1856 1861 #if defined key_hadocc 1857 1862 it = jp_had_nut 1858 #elif defined key_medusa && defined key_foam_medusa1863 #elif defined key_medusa 1859 1864 it = jpdin 1860 #else1861 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' )1862 1865 #endif 1863 1866 IF ( ln_phytobal ) THEN … … 1873 1876 pno3_bkginc(:,:,:) = pno3_bkginc(:,:,:) - ph_balinc(:,:,:,it) 1874 1877 ENDIF 1878 #else 1879 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 1880 #endif 1875 1881 ENDIF 1876 1882 … … 1878 1884 #if defined key_medusa && defined key_foam_medusa 1879 1885 it = jpsil 1880 #else1881 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' )1882 #endif1883 1886 IF ( ln_phytobal ) THEN 1884 1887 psi4_bkginc(:,:,:) = psi4_bkginc(:,:,:) - phyto2d_balinc(:,:,:,it) … … 1893 1896 psi4_bkginc(:,:,:) = psi4_bkginc(:,:,:) - ph_balinc(:,:,:,it) 1894 1897 ENDIF 1898 #else 1899 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 1900 #endif 1895 1901 ENDIF 1896 1902 1897 1903 IF ( ln_pdicinc ) THEN 1904 #if defined key_hadocc || ( defined key_medusa && defined key_foam_medusa ) 1898 1905 #if defined key_hadocc 1899 1906 it = jp_had_dic 1900 #elif defined key_medusa && defined key_foam_medusa1907 #elif defined key_medusa 1901 1908 it = jpdic 1902 #else1903 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' )1904 1909 #endif 1905 1910 IF ( ln_phytobal ) THEN … … 1915 1920 pdic_bkginc(:,:,:) = pdic_bkginc(:,:,:) - ph_balinc(:,:,:,it) 1916 1921 ENDIF 1922 #else 1923 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 1924 #endif 1917 1925 ENDIF 1918 1926 1919 1927 IF ( ln_palkinc ) THEN 1928 #if defined key_hadocc || ( defined key_medusa && defined key_foam_medusa ) 1920 1929 #if defined key_hadocc 1921 1930 it = jp_had_alk 1922 #elif defined key_medusa && defined key_foam_medusa1931 #elif defined key_medusa 1923 1932 it = jpalk 1924 #else1925 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' )1926 1933 #endif 1927 1934 IF ( ln_phytobal ) THEN … … 1937 1944 palk_bkginc(:,:,:) = palk_bkginc(:,:,:) - ph_balinc(:,:,:,it) 1938 1945 ENDIF 1946 #else 1947 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 1948 #endif 1939 1949 ENDIF 1940 1950 … … 1942 1952 #if defined key_medusa && defined key_foam_medusa 1943 1953 it = jpoxy 1944 #else1945 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' )1946 #endif1947 1954 IF ( ln_phytobal ) THEN 1948 1955 po2_bkginc(:,:,:) = po2_bkginc(:,:,:) - phyto2d_balinc(:,:,:,it) … … 1957 1964 po2_bkginc(:,:,:) = po2_bkginc(:,:,:) - ph_balinc(:,:,:,it) 1958 1965 ENDIF 1966 #else 1967 CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 1968 #endif 1959 1969 ENDIF 1960 1970 … … 2093 2103 #endif 2094 2104 2095 trn(:,:,:,jp_msa0:jp_msa1) = trn(:,:,:,jp_msa0:jp_msa1) + &2096 & pco2_balinc(:,:,:,jp_msa0:jp_msa1)2097 trb(:,:,:,jp_msa0:jp_msa1) = trn(:,:,:,jp_msa0:jp_msa1)2098 2099 2105 IF ( ln_pno3inc ) THEN 2100 2106 #if defined key_hadocc
Note: See TracChangeset
for help on using the changeset viewer.