Changeset 10667
- Timestamp:
- 2019-02-13T11:46:09+01:00 (6 years ago)
- Location:
- branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO
- Files:
-
- 3 edited
-
OPA_SRC/ASM/asmbgc.F90 (modified) (31 diffs)
-
OPA_SRC/ASM/asmphyto2dbal_ersem.F90 (modified) (6 diffs)
-
TOP_SRC/FABM/trcsms_fabm.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90
r10661 r10667 641 641 CALL iom_get( inum, jpdom_autoglo, 'medusa_pds', tracer_bkg(:,:,:,jppds) ) 642 642 #elif defined key_fabm 643 CALL iom_get( inum, jpdom_autoglo, 'ersem_chl1', tracer_bkg(:,:,:,jp_fabm_ chl1) )644 CALL iom_get( inum, jpdom_autoglo, 'ersem_chl2', tracer_bkg(:,:,:,jp_fabm_ chl2) )645 CALL iom_get( inum, jpdom_autoglo, 'ersem_chl3', tracer_bkg(:,:,:,jp_fabm_ chl3) )646 CALL iom_get( inum, jpdom_autoglo, 'ersem_chl4', tracer_bkg(:,:,:,jp_fabm_ chl4) )647 CALL iom_get( inum, jpdom_autoglo, 'ersem_p1c', tracer_bkg(:,:,:,jp_fabm_ p1c) )648 CALL iom_get( inum, jpdom_autoglo, 'ersem_p1n', tracer_bkg(:,:,:,jp_fabm_ p1n) )649 CALL iom_get( inum, jpdom_autoglo, 'ersem_p1p', tracer_bkg(:,:,:,jp_fabm_ p1p) )650 CALL iom_get( inum, jpdom_autoglo, 'ersem_p1s', tracer_bkg(:,:,:,jp_fabm_ p1s) )651 CALL iom_get( inum, jpdom_autoglo, 'ersem_p2c', tracer_bkg(:,:,:,jp_fabm_ p2c) )652 CALL iom_get( inum, jpdom_autoglo, 'ersem_p2n', tracer_bkg(:,:,:,jp_fabm_ p2n) )653 CALL iom_get( inum, jpdom_autoglo, 'ersem_p2p', tracer_bkg(:,:,:,jp_fabm_ p2p) )654 CALL iom_get( inum, jpdom_autoglo, 'ersem_p3c', tracer_bkg(:,:,:,jp_fabm_ p3c) )655 CALL iom_get( inum, jpdom_autoglo, 'ersem_p3n', tracer_bkg(:,:,:,jp_fabm_ p3n) )656 CALL iom_get( inum, jpdom_autoglo, 'ersem_p3p', tracer_bkg(:,:,:,jp_fabm_ p3p) )657 CALL iom_get( inum, jpdom_autoglo, 'ersem_p4c', tracer_bkg(:,:,:,jp_fabm_ p4c) )658 CALL iom_get( inum, jpdom_autoglo, 'ersem_p4n', tracer_bkg(:,:,:,jp_fabm_ p4n) )659 CALL iom_get( inum, jpdom_autoglo, 'ersem_p4p', tracer_bkg(:,:,:,jp_fabm_ p4p) )643 CALL iom_get( inum, jpdom_autoglo, 'ersem_chl1', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_chl1) ) 644 CALL iom_get( inum, jpdom_autoglo, 'ersem_chl2', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_chl2) ) 645 CALL iom_get( inum, jpdom_autoglo, 'ersem_chl3', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_chl3) ) 646 CALL iom_get( inum, jpdom_autoglo, 'ersem_chl4', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_chl4) ) 647 CALL iom_get( inum, jpdom_autoglo, 'ersem_p1c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p1c) ) 648 CALL iom_get( inum, jpdom_autoglo, 'ersem_p1n', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p1n) ) 649 CALL iom_get( inum, jpdom_autoglo, 'ersem_p1p', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p1p) ) 650 CALL iom_get( inum, jpdom_autoglo, 'ersem_p1s', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p1s) ) 651 CALL iom_get( inum, jpdom_autoglo, 'ersem_p2c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p2c) ) 652 CALL iom_get( inum, jpdom_autoglo, 'ersem_p2n', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p2n) ) 653 CALL iom_get( inum, jpdom_autoglo, 'ersem_p2p', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p2p) ) 654 CALL iom_get( inum, jpdom_autoglo, 'ersem_p3c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p3c) ) 655 CALL iom_get( inum, jpdom_autoglo, 'ersem_p3n', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p3n) ) 656 CALL iom_get( inum, jpdom_autoglo, 'ersem_p3p', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p3p) ) 657 CALL iom_get( inum, jpdom_autoglo, 'ersem_p4c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p4c) ) 658 CALL iom_get( inum, jpdom_autoglo, 'ersem_p4n', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p4n) ) 659 CALL iom_get( inum, jpdom_autoglo, 'ersem_p4p', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p4p) ) 660 660 #endif 661 661 … … 690 690 CALL iom_get( inum, jpdom_autoglo, 'medusa_oxy', tracer_bkg(:,:,:,jpoxy) ) 691 691 #elif defined key_fabm 692 CALL iom_get( inum, jpdom_autoglo, 'ersem_z4c', tracer_bkg(:,:,:,jp_fabm_ z4c) )693 CALL iom_get( inum, jpdom_autoglo, 'ersem_z5c', tracer_bkg(:,:,:,jp_fabm_ z5c) )694 CALL iom_get( inum, jpdom_autoglo, 'ersem_z5n', tracer_bkg(:,:,:,jp_fabm_ z5n) )695 CALL iom_get( inum, jpdom_autoglo, 'ersem_z5p', tracer_bkg(:,:,:,jp_fabm_ z5p) )696 CALL iom_get( inum, jpdom_autoglo, 'ersem_z6c', tracer_bkg(:,:,:,jp_fabm_ z6c) )697 CALL iom_get( inum, jpdom_autoglo, 'ersem_z6n', tracer_bkg(:,:,:,jp_fabm_ z6n) )698 CALL iom_get( inum, jpdom_autoglo, 'ersem_z6p', tracer_bkg(:,:,:,jp_fabm_ z6p) )699 CALL iom_get( inum, jpdom_autoglo, 'ersem_n1p', tracer_bkg(:,:,:,jp_fabm_ n1p) )700 CALL iom_get( inum, jpdom_autoglo, 'ersem_n3n', tracer_bkg(:,:,:,jp_fabm_ n3n) )701 CALL iom_get( inum, jpdom_autoglo, 'ersem_n4n', tracer_bkg(:,:,:,jp_fabm_ n4n) )702 CALL iom_get( inum, jpdom_autoglo, 'ersem_n5s', tracer_bkg(:,:,:,jp_fabm_ n5s) )703 CALL iom_get( inum, jpdom_autoglo, 'ersem_o2o', tracer_bkg(:,:,:,jp_fabm_ o2o) )704 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3c', tracer_bkg(:,:,:,jp_fabm_ o3c) )705 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3ba', tracer_bkg(:,:,:,jp_fabm_ o3ba) )692 CALL iom_get( inum, jpdom_autoglo, 'ersem_z4c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z4c) ) 693 CALL iom_get( inum, jpdom_autoglo, 'ersem_z5c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z5c) ) 694 CALL iom_get( inum, jpdom_autoglo, 'ersem_z5n', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z5n) ) 695 CALL iom_get( inum, jpdom_autoglo, 'ersem_z5p', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z5p) ) 696 CALL iom_get( inum, jpdom_autoglo, 'ersem_z6c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z6c) ) 697 CALL iom_get( inum, jpdom_autoglo, 'ersem_z6n', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z6n) ) 698 CALL iom_get( inum, jpdom_autoglo, 'ersem_z6p', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z6p) ) 699 CALL iom_get( inum, jpdom_autoglo, 'ersem_n1p', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_n1p) ) 700 CALL iom_get( inum, jpdom_autoglo, 'ersem_n3n', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_n3n) ) 701 CALL iom_get( inum, jpdom_autoglo, 'ersem_n4n', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_n4n) ) 702 CALL iom_get( inum, jpdom_autoglo, 'ersem_n5s', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_n5s) ) 703 CALL iom_get( inum, jpdom_autoglo, 'ersem_o2o', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_o2o) ) 704 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_o3c) ) 705 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3ba', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_o3ba) ) 706 706 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3ta', totalk_bkg(:,:,:) ) 707 707 totalk_bkg(:,:,:) = totalk_bkg(:,:,:) * tmask(:,:,:) … … 715 715 CALL iom_get( inum, jpdom_autoglo, 'medusa_alk', tracer_bkg(:,:,:,jpalk) ) 716 716 #elif defined key_fabm 717 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3c', tracer_bkg(:,:,:,jp_fabm_ o3c) )718 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3ba', tracer_bkg(:,:,:,jp_fabm_ o3ba) )717 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_o3c) ) 718 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3ba', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_o3ba) ) 719 719 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3ta', totalk_bkg(:,:,:) ) 720 720 totalk_bkg(:,:,:) = totalk_bkg(:,:,:) * tmask(:,:,:) … … 746 746 CALL iom_get( inum, jpdom_autoglo, 'medusa_alk', tracer_bkg(:,:,:,jpalk) ) 747 747 #elif defined key_fabm 748 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3c', tracer_bkg(:,:,:,jp_fabm_ o3c) )749 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3ba', tracer_bkg(:,:,:,jp_fabm_ o3ba) )748 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3c', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_o3c) ) 749 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3ba', tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_o3ba) ) 750 750 CALL iom_get( inum, jpdom_autoglo, 'ersem_o3ta', totalk_bkg(:,:,:) ) 751 751 totalk_bkg(:,:,:) = totalk_bkg(:,:,:) * tmask(:,:,:) … … 846 846 ENDIF 847 847 #elif defined key_fabm 848 CALL iom_rstput( kt, kt, inum, 'phy2d_chl1', phyto2d_balinc(:,:,:,jp_fabm_ chl1) )849 CALL iom_rstput( kt, kt, inum, 'phy2d_chl2', phyto2d_balinc(:,:,:,jp_fabm_ chl2) )850 CALL iom_rstput( kt, kt, inum, 'phy2d_chl3', phyto2d_balinc(:,:,:,jp_fabm_ chl3) )851 CALL iom_rstput( kt, kt, inum, 'phy2d_chl4', phyto2d_balinc(:,:,:,jp_fabm_ chl4) )852 CALL iom_rstput( kt, kt, inum, 'phy2d_p1c', phyto2d_balinc(:,:,:,jp_fabm_ p1c) )853 CALL iom_rstput( kt, kt, inum, 'phy2d_p1n', phyto2d_balinc(:,:,:,jp_fabm_ p1n) )854 CALL iom_rstput( kt, kt, inum, 'phy2d_p1p', phyto2d_balinc(:,:,:,jp_fabm_ p1p) )855 CALL iom_rstput( kt, kt, inum, 'phy2d_p1s', phyto2d_balinc(:,:,:,jp_fabm_ p1s) )856 CALL iom_rstput( kt, kt, inum, 'phy2d_p2c', phyto2d_balinc(:,:,:,jp_fabm_ p2c) )857 CALL iom_rstput( kt, kt, inum, 'phy2d_p2n', phyto2d_balinc(:,:,:,jp_fabm_ p2n) )858 CALL iom_rstput( kt, kt, inum, 'phy2d_p2p', phyto2d_balinc(:,:,:,jp_fabm_ p2p) )859 CALL iom_rstput( kt, kt, inum, 'phy2d_p3c', phyto2d_balinc(:,:,:,jp_fabm_ p3c) )860 CALL iom_rstput( kt, kt, inum, 'phy2d_p3n', phyto2d_balinc(:,:,:,jp_fabm_ p3n) )861 CALL iom_rstput( kt, kt, inum, 'phy2d_p3p', phyto2d_balinc(:,:,:,jp_fabm_ p3p) )862 CALL iom_rstput( kt, kt, inum, 'phy2d_p4c', phyto2d_balinc(:,:,:,jp_fabm_ p4c) )863 CALL iom_rstput( kt, kt, inum, 'phy2d_p4n', phyto2d_balinc(:,:,:,jp_fabm_ p4n) )864 CALL iom_rstput( kt, kt, inum, 'phy2d_p4p', phyto2d_balinc(:,:,:,jp_fabm_ p4p) )848 CALL iom_rstput( kt, kt, inum, 'phy2d_chl1', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl1) ) 849 CALL iom_rstput( kt, kt, inum, 'phy2d_chl2', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl2) ) 850 CALL iom_rstput( kt, kt, inum, 'phy2d_chl3', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl3) ) 851 CALL iom_rstput( kt, kt, inum, 'phy2d_chl4', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl4) ) 852 CALL iom_rstput( kt, kt, inum, 'phy2d_p1c', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p1c) ) 853 CALL iom_rstput( kt, kt, inum, 'phy2d_p1n', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p1n) ) 854 CALL iom_rstput( kt, kt, inum, 'phy2d_p1p', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p1p) ) 855 CALL iom_rstput( kt, kt, inum, 'phy2d_p1s', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p1s) ) 856 CALL iom_rstput( kt, kt, inum, 'phy2d_p2c', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p2c) ) 857 CALL iom_rstput( kt, kt, inum, 'phy2d_p2n', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p2n) ) 858 CALL iom_rstput( kt, kt, inum, 'phy2d_p2p', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p2p) ) 859 CALL iom_rstput( kt, kt, inum, 'phy2d_p3c', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p3c) ) 860 CALL iom_rstput( kt, kt, inum, 'phy2d_p3n', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p3n) ) 861 CALL iom_rstput( kt, kt, inum, 'phy2d_p3p', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p3p) ) 862 CALL iom_rstput( kt, kt, inum, 'phy2d_p4c', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p4c) ) 863 CALL iom_rstput( kt, kt, inum, 'phy2d_p4n', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p4n) ) 864 CALL iom_rstput( kt, kt, inum, 'phy2d_p4p', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p4p) ) 865 865 IF ( ln_phytobal ) THEN 866 CALL iom_rstput( kt, kt, inum, 'phy2d_z4c', phyto2d_balinc(:,:,:,jp_fabm_ z4c) )867 CALL iom_rstput( kt, kt, inum, 'phy2d_z5c', phyto2d_balinc(:,:,:,jp_fabm_ z5c) )868 CALL iom_rstput( kt, kt, inum, 'phy2d_z5n', phyto2d_balinc(:,:,:,jp_fabm_ z5n) )869 CALL iom_rstput( kt, kt, inum, 'phy2d_z5p', phyto2d_balinc(:,:,:,jp_fabm_ z5p) )870 CALL iom_rstput( kt, kt, inum, 'phy2d_z6c', phyto2d_balinc(:,:,:,jp_fabm_ z6c) )871 CALL iom_rstput( kt, kt, inum, 'phy2d_z6n', phyto2d_balinc(:,:,:,jp_fabm_ z6n) )872 CALL iom_rstput( kt, kt, inum, 'phy2d_z6p', phyto2d_balinc(:,:,:,jp_fabm_ z6p) )873 CALL iom_rstput( kt, kt, inum, 'phy2d_n1p', phyto2d_balinc(:,:,:,jp_fabm_ n1p) )874 CALL iom_rstput( kt, kt, inum, 'phy2d_n3n', phyto2d_balinc(:,:,:,jp_fabm_ n3n) )875 CALL iom_rstput( kt, kt, inum, 'phy2d_n4n', phyto2d_balinc(:,:,:,jp_fabm_ n4n) )876 CALL iom_rstput( kt, kt, inum, 'phy2d_n5s', phyto2d_balinc(:,:,:,jp_fabm_ n5s) )877 CALL iom_rstput( kt, kt, inum, 'phy2d_o2o', phyto2d_balinc(:,:,:,jp_fabm_ o2o) )878 CALL iom_rstput( kt, kt, inum, 'phy2d_o3c', phyto2d_balinc(:,:,:,jp_fabm_ o3c) )879 CALL iom_rstput( kt, kt, inum, 'phy2d_o3ba', phyto2d_balinc(:,:,:,jp_fabm_ o3ba) )866 CALL iom_rstput( kt, kt, inum, 'phy2d_z4c', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_z4c) ) 867 CALL iom_rstput( kt, kt, inum, 'phy2d_z5c', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_z5c) ) 868 CALL iom_rstput( kt, kt, inum, 'phy2d_z5n', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_z5n) ) 869 CALL iom_rstput( kt, kt, inum, 'phy2d_z5p', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_z5p) ) 870 CALL iom_rstput( kt, kt, inum, 'phy2d_z6c', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_z6c) ) 871 CALL iom_rstput( kt, kt, inum, 'phy2d_z6n', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_z6n) ) 872 CALL iom_rstput( kt, kt, inum, 'phy2d_z6p', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_z6p) ) 873 CALL iom_rstput( kt, kt, inum, 'phy2d_n1p', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_n1p) ) 874 CALL iom_rstput( kt, kt, inum, 'phy2d_n3n', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_n3n) ) 875 CALL iom_rstput( kt, kt, inum, 'phy2d_n4n', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_n4n) ) 876 CALL iom_rstput( kt, kt, inum, 'phy2d_n5s', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_n5s) ) 877 CALL iom_rstput( kt, kt, inum, 'phy2d_o2o', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_o2o) ) 878 CALL iom_rstput( kt, kt, inum, 'phy2d_o3c', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_o3c) ) 879 CALL iom_rstput( kt, kt, inum, 'phy2d_o3ba', phyto2d_balinc(:,:,:,jp_fabm_m1+jp_fabm_o3ba) ) 880 880 ENDIF 881 881 #endif … … 892 892 CALL iom_rstput( kt, kt, inum, 'phy3d_phy', phyto3d_balinc(:,:,:,jp_had_phy) ) 893 893 #elif defined key_fabm 894 CALL iom_rstput( kt, kt, inum, 'phy3d_chl1', phyto3d_balinc(:,:,:,jp_fabm_ chl1) )895 CALL iom_rstput( kt, kt, inum, 'phy3d_chl2', phyto3d_balinc(:,:,:,jp_fabm_ chl2) )896 CALL iom_rstput( kt, kt, inum, 'phy3d_chl3', phyto3d_balinc(:,:,:,jp_fabm_ chl3) )897 CALL iom_rstput( kt, kt, inum, 'phy3d_chl4', phyto3d_balinc(:,:,:,jp_fabm_ chl4) )898 CALL iom_rstput( kt, kt, inum, 'phy3d_p1c', phyto3d_balinc(:,:,:,jp_fabm_ p1c) )899 CALL iom_rstput( kt, kt, inum, 'phy3d_p1n', phyto3d_balinc(:,:,:,jp_fabm_ p1n) )900 CALL iom_rstput( kt, kt, inum, 'phy3d_p1p', phyto3d_balinc(:,:,:,jp_fabm_ p1p) )901 CALL iom_rstput( kt, kt, inum, 'phy3d_p1s', phyto3d_balinc(:,:,:,jp_fabm_ p1s) )902 CALL iom_rstput( kt, kt, inum, 'phy3d_p2c', phyto3d_balinc(:,:,:,jp_fabm_ p2c) )903 CALL iom_rstput( kt, kt, inum, 'phy3d_p2n', phyto3d_balinc(:,:,:,jp_fabm_ p2n) )904 CALL iom_rstput( kt, kt, inum, 'phy3d_p2p', phyto3d_balinc(:,:,:,jp_fabm_ p2p) )905 CALL iom_rstput( kt, kt, inum, 'phy3d_p3c', phyto3d_balinc(:,:,:,jp_fabm_ p3c) )906 CALL iom_rstput( kt, kt, inum, 'phy3d_p3n', phyto3d_balinc(:,:,:,jp_fabm_ p3n) )907 CALL iom_rstput( kt, kt, inum, 'phy3d_p3p', phyto3d_balinc(:,:,:,jp_fabm_ p3p) )908 CALL iom_rstput( kt, kt, inum, 'phy3d_p4c', phyto3d_balinc(:,:,:,jp_fabm_ p4c) )909 CALL iom_rstput( kt, kt, inum, 'phy3d_p4n', phyto3d_balinc(:,:,:,jp_fabm_ p4n) )910 CALL iom_rstput( kt, kt, inum, 'phy3d_p4p', phyto3d_balinc(:,:,:,jp_fabm_ p4p) )894 CALL iom_rstput( kt, kt, inum, 'phy3d_chl1', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl1) ) 895 CALL iom_rstput( kt, kt, inum, 'phy3d_chl2', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl2) ) 896 CALL iom_rstput( kt, kt, inum, 'phy3d_chl3', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl3) ) 897 CALL iom_rstput( kt, kt, inum, 'phy3d_chl4', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl4) ) 898 CALL iom_rstput( kt, kt, inum, 'phy3d_p1c', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p1c) ) 899 CALL iom_rstput( kt, kt, inum, 'phy3d_p1n', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p1n) ) 900 CALL iom_rstput( kt, kt, inum, 'phy3d_p1p', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p1p) ) 901 CALL iom_rstput( kt, kt, inum, 'phy3d_p1s', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p1s) ) 902 CALL iom_rstput( kt, kt, inum, 'phy3d_p2c', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p2c) ) 903 CALL iom_rstput( kt, kt, inum, 'phy3d_p2n', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p2n) ) 904 CALL iom_rstput( kt, kt, inum, 'phy3d_p2p', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p2p) ) 905 CALL iom_rstput( kt, kt, inum, 'phy3d_p3c', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p3c) ) 906 CALL iom_rstput( kt, kt, inum, 'phy3d_p3n', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p3n) ) 907 CALL iom_rstput( kt, kt, inum, 'phy3d_p3p', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p3p) ) 908 CALL iom_rstput( kt, kt, inum, 'phy3d_p4c', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p4c) ) 909 CALL iom_rstput( kt, kt, inum, 'phy3d_p4n', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p4n) ) 910 CALL iom_rstput( kt, kt, inum, 'phy3d_p4p', phyto3d_balinc(:,:,:,jp_fabm_m1+jp_fabm_p4p) ) 911 911 #endif 912 912 ENDIF … … 920 920 CALL iom_rstput( kt, kt, inum, 'pco2_alk', pco2_balinc(:,:,:,jp_had_alk) ) 921 921 #elif defined key_fabm 922 CALL iom_rstput( kt, kt, inum, 'pco2_o3c', pco2_balinc(:,:,:,jp_fabm_ o3c) )923 CALL iom_rstput( kt, kt, inum, 'pco2_o3ba', pco2_balinc(:,:,:,jp_fabm_ o3ba) )922 CALL iom_rstput( kt, kt, inum, 'pco2_o3c', pco2_balinc(:,:,:,jp_fabm_m1+jp_fabm_o3c) ) 923 CALL iom_rstput( kt, kt, inum, 'pco2_o3ba', pco2_balinc(:,:,:,jp_fabm_m1+jp_fabm_o3ba) ) 924 924 #endif 925 925 ELSE IF ( ln_sfco2inc ) THEN … … 931 931 CALL iom_rstput( kt, kt, inum, 'fco2_alk', pco2_balinc(:,:,:,jp_had_alk) ) 932 932 #elif defined key_fabm 933 CALL iom_rstput( kt, kt, inum, 'fco2_o3c', pco2_balinc(:,:,:,jp_fabm_ o3c) )934 CALL iom_rstput( kt, kt, inum, 'fco2_o3ba', pco2_balinc(:,:,:,jp_fabm_ o3ba) )933 CALL iom_rstput( kt, kt, inum, 'fco2_o3c', pco2_balinc(:,:,:,jp_fabm_m1+jp_fabm_o3c) ) 934 CALL iom_rstput( kt, kt, inum, 'fco2_o3ba', pco2_balinc(:,:,:,jp_fabm_m1+jp_fabm_o3ba) ) 935 935 #endif 936 936 ENDIF … … 944 944 CALL iom_rstput( kt, kt, inum, 'ph_alk', ph_balinc(:,:,:,jp_had_alk) ) 945 945 #elif defined key_fabm 946 CALL iom_rstput( kt, kt, inum, 'ph_o3c', ph_balinc(:,:,:,jp_fabm_ o3c) )947 CALL iom_rstput( kt, kt, inum, 'ph_o3ba', ph_balinc(:,:,:,jp_fabm_ o3ba) )946 CALL iom_rstput( kt, kt, inum, 'ph_o3c', ph_balinc(:,:,:,jp_fabm_m1+jp_fabm_o3c) ) 947 CALL iom_rstput( kt, kt, inum, 'ph_o3ba', ph_balinc(:,:,:,jp_fabm_m1+jp_fabm_o3ba) ) 948 948 #endif 949 949 ENDIF … … 1159 1159 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'phyt_avg' , phyt_avg_tavg ) 1160 1160 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'mld_max' , mld_max ) 1161 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_chl1' , trn_tavg(:,:,:,jp_fabm_ chl1) )1162 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_chl2' , trn_tavg(:,:,:,jp_fabm_ chl2) )1163 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_chl3' , trn_tavg(:,:,:,jp_fabm_ chl3) )1164 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_chl4' , trn_tavg(:,:,:,jp_fabm_ chl4) )1165 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p1c' , trn_tavg(:,:,:,jp_fabm_ p1c) )1166 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p1n' , trn_tavg(:,:,:,jp_fabm_ p1n) )1167 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p1p' , trn_tavg(:,:,:,jp_fabm_ p1p) )1168 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p1s' , trn_tavg(:,:,:,jp_fabm_ p1s) )1169 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p2c' , trn_tavg(:,:,:,jp_fabm_ p2c) )1170 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p2n' , trn_tavg(:,:,:,jp_fabm_ p2n) )1171 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p2p' , trn_tavg(:,:,:,jp_fabm_ p2p) )1172 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p3c' , trn_tavg(:,:,:,jp_fabm_ p3c) )1173 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p3n' , trn_tavg(:,:,:,jp_fabm_ p3n) )1174 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p3p' , trn_tavg(:,:,:,jp_fabm_ p3p) )1175 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p4c' , trn_tavg(:,:,:,jp_fabm_ p4c) )1176 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p4n' , trn_tavg(:,:,:,jp_fabm_ p4n) )1177 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p4p' , trn_tavg(:,:,:,jp_fabm_ p4p) )1178 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z4c' , trn_tavg(:,:,:,jp_fabm_ z4c) )1179 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z5c' , trn_tavg(:,:,:,jp_fabm_ z5c) )1180 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z5n' , trn_tavg(:,:,:,jp_fabm_ z5n) )1181 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z5p' , trn_tavg(:,:,:,jp_fabm_ z5p) )1182 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z6c' , trn_tavg(:,:,:,jp_fabm_ z6c) )1183 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z6n' , trn_tavg(:,:,:,jp_fabm_ z6n) )1184 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z6p' , trn_tavg(:,:,:,jp_fabm_ z6p) )1185 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_n1p' , trn_tavg(:,:,:,jp_fabm_ n1p) )1186 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_n3n' , trn_tavg(:,:,:,jp_fabm_ n3n) )1187 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_n4n' , trn_tavg(:,:,:,jp_fabm_ n4n) )1188 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_n5s' , trn_tavg(:,:,:,jp_fabm_ n5s) )1189 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_o2o' , trn_tavg(:,:,:,jp_fabm_ o2o) )1190 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_o3c' , trn_tavg(:,:,:,jp_fabm_ o3c) )1191 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_o3ba' , trn_tavg(:,:,:,jp_fabm_ o3ba) )1161 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_chl1' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_chl1) ) 1162 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_chl2' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_chl2) ) 1163 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_chl3' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_chl3) ) 1164 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_chl4' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_chl4) ) 1165 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p1c' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p1c) ) 1166 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p1n' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p1n) ) 1167 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p1p' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p1p) ) 1168 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p1s' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p1s) ) 1169 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p2c' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p2c) ) 1170 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p2n' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p2n) ) 1171 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p2p' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p2p) ) 1172 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p3c' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p3c) ) 1173 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p3n' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p3n) ) 1174 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p3p' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p3p) ) 1175 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p4c' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p4c) ) 1176 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p4n' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p4n) ) 1177 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_p4p' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_p4p) ) 1178 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z4c' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_z4c) ) 1179 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z5c' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_z5c) ) 1180 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z5n' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_z5n) ) 1181 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z5p' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_z5p) ) 1182 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z6c' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_z6c) ) 1183 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z6n' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_z6n) ) 1184 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_z6p' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_z6p) ) 1185 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_n1p' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_n1p) ) 1186 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_n3n' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_n3n) ) 1187 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_n4n' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_n4n) ) 1188 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_n5s' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_n5s) ) 1189 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_o2o' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_o2o) ) 1190 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_o3c' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_o3c) ) 1191 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_o3ba' , trn_tavg(:,:,:,jp_fabm_m1+jp_fabm_o3ba) ) 1192 1192 CALL iom_rstput( kt, nitbgcbkg_r, knum, 'ersem_o3ta' , totalk_tavg ) 1193 1193 #endif … … 1302 1302 zbkg_chltot(:,:) = chl_bkg(:,:,1) 1303 1303 #elif defined key_fabm 1304 zbkg_chltot(:,:) = tracer_bkg(:,:,1,jp_fabm_ chl1) + &1305 & tracer_bkg(:,:,1,jp_fabm_ chl2) + &1306 & tracer_bkg(:,:,1,jp_fabm_ chl3) + &1307 & tracer_bkg(:,:,1,jp_fabm_ chl4)1304 zbkg_chltot(:,:) = tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_chl1) + & 1305 & tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_chl2) + & 1306 & tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_chl3) + & 1307 & tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_chl4) 1308 1308 #endif 1309 1309 CALL asm_bgc_unlog_2d( zbkg_chltot, slchltot_bkginc, zinc_chltot ) … … 1320 1320 zbkg_chldia(:,:) = tracer_bkg(:,:,1,jpchd) 1321 1321 #elif defined key_fabm 1322 zbkg_chldia(:,:) = tracer_bkg(:,:,1,jp_fabm_ chl1)1322 zbkg_chldia(:,:) = tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_chl1) 1323 1323 #endif 1324 1324 CALL asm_bgc_unlog_2d( zbkg_chldia, slchldia_bkginc, zinc_chldia ) … … 1341 1341 ! Nanophytoplankton chlorophyll 1342 1342 IF ( ln_slchlnaninc ) THEN 1343 zbkg_chlnan(:,:) = tracer_bkg(:,:,1,jp_fabm_ chl2)1343 zbkg_chlnan(:,:) = tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_chl2) 1344 1344 CALL asm_bgc_unlog_2d( zbkg_chlnan, slchlnan_bkginc, zinc_chlnan ) 1345 1345 ELSE … … 1349 1349 ! Picophytoplankton chlorophyll 1350 1350 IF ( ln_slchlpicinc ) THEN 1351 zbkg_chlpic(:,:) = tracer_bkg(:,:,1,jp_fabm_ chl3)1351 zbkg_chlpic(:,:) = tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_chl3) 1352 1352 CALL asm_bgc_unlog_2d( zbkg_chlpic, slchlpic_bkginc, zinc_chlpic ) 1353 1353 ELSE … … 1357 1357 ! Dinoflagellate chlorophyll 1358 1358 IF ( ln_slchldininc ) THEN 1359 zbkg_chldin(:,:) = tracer_bkg(:,:,1,jp_fabm_ chl4)1359 zbkg_chldin(:,:) = tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_chl4) 1360 1360 CALL asm_bgc_unlog_2d( zbkg_chldin, slchldin_bkginc, zinc_chldin ) 1361 1361 ELSE … … 1647 1647 bkg_chl(:,:,:) = chl_bkg(:,:,:) 1648 1648 #elif defined key_fabm 1649 bkg_chl(:,:,:) = tracer_bkg(:,:,:,jp_fabm_ chl1) + &1650 & tracer_bkg(:,:,:,jp_fabm_ chl2) + &1651 & tracer_bkg(:,:,:,jp_fabm_ chl3) + &1652 & tracer_bkg(:,:,:,jp_fabm_ chl4)1649 bkg_chl(:,:,:) = tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_chl1) + & 1650 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_chl2) + & 1651 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_chl3) + & 1652 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_chl4) 1653 1653 #endif 1654 1654 DO jk = 1, jpk … … 1707 1707 DO jj = 1, jpj 1708 1708 DO ji = 1, jpi 1709 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_ chl1) > 0.0 ) .AND. &1710 & ( tracer_bkg(ji,jj,jk,jp_fabm_ chl2) > 0.0 ) .AND. &1711 & ( tracer_bkg(ji,jj,jk,jp_fabm_ chl3) > 0.0 ) .AND. &1712 & ( tracer_bkg(ji,jj,jk,jp_fabm_ chl4) > 0.0 ) ) THEN1713 zfrac_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_ chl1) / bkg_chl(ji,jj,jk)1714 zfrac_chl2 = tracer_bkg(ji,jj,jk,jp_fabm_ chl2) / bkg_chl(ji,jj,jk)1715 zfrac_chl3 = tracer_bkg(ji,jj,jk,jp_fabm_ chl3) / bkg_chl(ji,jj,jk)1716 zfrac_chl4 = tracer_bkg(ji,jj,jk,jp_fabm_ chl4) / bkg_chl(ji,jj,jk)1717 phyto3d_balinc(ji,jj,jk,jp_fabm_ chl1) = chl_inc(ji,jj,jk) * zfrac_chl11718 phyto3d_balinc(ji,jj,jk,jp_fabm_ chl2) = chl_inc(ji,jj,jk) * zfrac_chl21719 phyto3d_balinc(ji,jj,jk,jp_fabm_ chl3) = chl_inc(ji,jj,jk) * zfrac_chl31720 phyto3d_balinc(ji,jj,jk,jp_fabm_ chl4) = chl_inc(ji,jj,jk) * zfrac_chl41721 zrat_p1c_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_ p1c) / tracer_bkg(ji,jj,jk,jp_fabm_chl1)1722 zrat_p1n_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_ p1n) / tracer_bkg(ji,jj,jk,jp_fabm_chl1)1723 zrat_p1p_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_ p1p) / tracer_bkg(ji,jj,jk,jp_fabm_chl1)1724 zrat_p1s_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_ p1s) / tracer_bkg(ji,jj,jk,jp_fabm_chl1)1725 zrat_p2c_chl2 = tracer_bkg(ji,jj,jk,jp_fabm_ p2c) / tracer_bkg(ji,jj,jk,jp_fabm_chl2)1726 zrat_p2n_chl2 = tracer_bkg(ji,jj,jk,jp_fabm_ p2n) / tracer_bkg(ji,jj,jk,jp_fabm_chl2)1727 zrat_p2p_chl2 = tracer_bkg(ji,jj,jk,jp_fabm_ p2p) / tracer_bkg(ji,jj,jk,jp_fabm_chl2)1728 zrat_p3c_chl3 = tracer_bkg(ji,jj,jk,jp_fabm_ p3c) / tracer_bkg(ji,jj,jk,jp_fabm_chl3)1729 zrat_p3n_chl3 = tracer_bkg(ji,jj,jk,jp_fabm_ p3n) / tracer_bkg(ji,jj,jk,jp_fabm_chl3)1730 zrat_p3p_chl3 = tracer_bkg(ji,jj,jk,jp_fabm_ p3p) / tracer_bkg(ji,jj,jk,jp_fabm_chl3)1731 zrat_p4c_chl4 = tracer_bkg(ji,jj,jk,jp_fabm_ p4c) / tracer_bkg(ji,jj,jk,jp_fabm_chl4)1732 zrat_p4n_chl4 = tracer_bkg(ji,jj,jk,jp_fabm_ p4n) / tracer_bkg(ji,jj,jk,jp_fabm_chl4)1733 zrat_p4p_chl4 = tracer_bkg(ji,jj,jk,jp_fabm_ p4p) / tracer_bkg(ji,jj,jk,jp_fabm_chl4)1734 phyto3d_balinc(ji,jj,jk,jp_fabm_ p1c) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl1) * zrat_p1c_chl11735 phyto3d_balinc(ji,jj,jk,jp_fabm_ p1n) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl1) * zrat_p1n_chl11736 phyto3d_balinc(ji,jj,jk,jp_fabm_ p1p) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl1) * zrat_p1p_chl11737 phyto3d_balinc(ji,jj,jk,jp_fabm_ p1s) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl1) * zrat_p1s_chl11738 phyto3d_balinc(ji,jj,jk,jp_fabm_ p2c) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl2) * zrat_p2c_chl21739 phyto3d_balinc(ji,jj,jk,jp_fabm_ p2n) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl2) * zrat_p2n_chl21740 phyto3d_balinc(ji,jj,jk,jp_fabm_ p2p) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl2) * zrat_p2p_chl21741 phyto3d_balinc(ji,jj,jk,jp_fabm_ p3c) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl3) * zrat_p3c_chl31742 phyto3d_balinc(ji,jj,jk,jp_fabm_ p3n) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl3) * zrat_p3n_chl31743 phyto3d_balinc(ji,jj,jk,jp_fabm_ p3p) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl3) * zrat_p3p_chl31744 phyto3d_balinc(ji,jj,jk,jp_fabm_ p4c) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl4) * zrat_p4c_chl41745 phyto3d_balinc(ji,jj,jk,jp_fabm_ p4n) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl4) * zrat_p4n_chl41746 phyto3d_balinc(ji,jj,jk,jp_fabm_ p4p) = phyto3d_balinc(ji,jj,jk,jp_fabm_chl4) * zrat_p4p_chl41709 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) > 0.0 ) .AND. & 1710 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) > 0.0 ) .AND. & 1711 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) > 0.0 ) .AND. & 1712 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) > 0.0 ) ) THEN 1713 zfrac_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) / bkg_chl(ji,jj,jk) 1714 zfrac_chl2 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) / bkg_chl(ji,jj,jk) 1715 zfrac_chl3 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) / bkg_chl(ji,jj,jk) 1716 zfrac_chl4 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) / bkg_chl(ji,jj,jk) 1717 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) = chl_inc(ji,jj,jk) * zfrac_chl1 1718 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) = chl_inc(ji,jj,jk) * zfrac_chl2 1719 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) = chl_inc(ji,jj,jk) * zfrac_chl3 1720 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) = chl_inc(ji,jj,jk) * zfrac_chl4 1721 zrat_p1c_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 1722 zrat_p1n_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 1723 zrat_p1p_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 1724 zrat_p1s_chl1 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 1725 zrat_p2c_chl2 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 1726 zrat_p2n_chl2 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 1727 zrat_p2p_chl2 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 1728 zrat_p3c_chl3 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 1729 zrat_p3n_chl3 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 1730 zrat_p3p_chl3 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 1731 zrat_p4c_chl4 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 1732 zrat_p4n_chl4 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 1733 zrat_p4p_chl4 = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 1734 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) * zrat_p1c_chl1 1735 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) * zrat_p1n_chl1 1736 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) * zrat_p1p_chl1 1737 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) * zrat_p1s_chl1 1738 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) * zrat_p2c_chl2 1739 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) * zrat_p2n_chl2 1740 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) * zrat_p2p_chl2 1741 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) * zrat_p3c_chl3 1742 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) * zrat_p3n_chl3 1743 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) * zrat_p3p_chl3 1744 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) * zrat_p4c_chl4 1745 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) * zrat_p4n_chl4 1746 phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) = phyto3d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) * zrat_p4p_chl4 1747 1747 ENDIF 1748 1748 END DO … … 1974 1974 ! Account for phytoplankton balancing if required 1975 1975 IF ( ln_phytobal ) THEN 1976 dic_bkg_temp(:,:) = tracer_bkg(:,:,1,jp_fabm_ o3c) + phyto2d_balinc(:,:,1,jp_fabm_o3c)1977 alk_bkg_temp(:,:) = totalk_bkg(:,:,1) + phyto2d_balinc(:,:,1,jp_fabm_ o3ba)1976 dic_bkg_temp(:,:) = tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_o3c) + phyto2d_balinc(:,:,1,jp_fabm_m1+jp_fabm_o3c) 1977 alk_bkg_temp(:,:) = totalk_bkg(:,:,1) + phyto2d_balinc(:,:,1,jp_fabm_m1+jp_fabm_o3ba) 1978 1978 ELSE 1979 dic_bkg_temp(:,:) = tracer_bkg(:,:,1,jp_fabm_ o3c)1979 dic_bkg_temp(:,:) = tracer_bkg(:,:,1,jp_fabm_m1+jp_fabm_o3c) 1980 1980 alk_bkg_temp(:,:) = totalk_bkg(:,:,1) 1981 1981 ENDIF … … 1983 1983 CALL asm_pco2_bal( pco2_bkginc_temp(:,:), dic_bkg_temp(:,:), alk_bkg_temp(:,:), & 1984 1984 & tem_bkg_temp(:,:), sal_bkg_temp(:,:), & 1985 & pco2_balinc(:,:,1,jp_fabm_ o3c), pco2_balinc(:,:,1,jp_fabm_o3ba) )1985 & pco2_balinc(:,:,1,jp_fabm_m1+jp_fabm_o3c), pco2_balinc(:,:,1,jp_fabm_m1+jp_fabm_o3ba) ) 1986 1986 1987 1987 #else … … 2592 2592 #elif defined key_fabm 2593 2593 WHERE( pno3_bkginc(:,:,:) > 0.0_wp .OR. & 2594 & trn(:,:,:,jp_fabm_ n3n) + pno3_bkginc(:,:,:) * zincwgt > 0.0_wp )2595 trn(:,:,:,jp_fabm_ n3n) = trn(:,:,:,jp_fabm_n3n) + pno3_bkginc(:,:,:) * zincwgt2596 trb(:,:,:,jp_fabm_ n3n) = trb(:,:,:,jp_fabm_n3n) + pno3_bkginc(:,:,:) * zincwgt2594 & trn(:,:,:,jp_fabm_m1+jp_fabm_n3n) + pno3_bkginc(:,:,:) * zincwgt > 0.0_wp ) 2595 trn(:,:,:,jp_fabm_m1+jp_fabm_n3n) = trn(:,:,:,jp_fabm_m1+jp_fabm_n3n) + pno3_bkginc(:,:,:) * zincwgt 2596 trb(:,:,:,jp_fabm_m1+jp_fabm_n3n) = trb(:,:,:,jp_fabm_m1+jp_fabm_n3n) + pno3_bkginc(:,:,:) * zincwgt 2597 2597 END WHERE 2598 2598 #else … … 2610 2610 #elif defined key_fabm 2611 2611 WHERE( psi4_bkginc(:,:,:) > 0.0_wp .OR. & 2612 & trn(:,:,:,jp_fabm_ n5s) + psi4_bkginc(:,:,:) * zincwgt > 0.0_wp )2613 trn(:,:,:,jp_fabm_ n5s) = trn(:,:,:,jp_fabm_n5s) + psi4_bkginc(:,:,:) * zincwgt2614 trb(:,:,:,jp_fabm_ n5s) = trb(:,:,:,jp_fabm_n5s) + psi4_bkginc(:,:,:) * zincwgt2612 & trn(:,:,:,jp_fabm_m1+jp_fabm_n5s) + psi4_bkginc(:,:,:) * zincwgt > 0.0_wp ) 2613 trn(:,:,:,jp_fabm_m1+jp_fabm_n5s) = trn(:,:,:,jp_fabm_m1+jp_fabm_n5s) + psi4_bkginc(:,:,:) * zincwgt 2614 trb(:,:,:,jp_fabm_m1+jp_fabm_n5s) = trb(:,:,:,jp_fabm_m1+jp_fabm_n5s) + psi4_bkginc(:,:,:) * zincwgt 2615 2615 END WHERE 2616 2616 #else … … 2622 2622 #if defined key_fabm 2623 2623 WHERE( ppo4_bkginc(:,:,:) > 0.0_wp .OR. & 2624 & trn(:,:,:,jp_fabm_ n1p) + ppo4_bkginc(:,:,:) * zincwgt > 0.0_wp )2625 trn(:,:,:,jp_fabm_ n1p) = trn(:,:,:,jp_fabm_n1p) + ppo4_bkginc(:,:,:) * zincwgt2626 trb(:,:,:,jp_fabm_ n1p) = trb(:,:,:,jp_fabm_n1p) + ppo4_bkginc(:,:,:) * zincwgt2624 & trn(:,:,:,jp_fabm_m1+jp_fabm_n1p) + ppo4_bkginc(:,:,:) * zincwgt > 0.0_wp ) 2625 trn(:,:,:,jp_fabm_m1+jp_fabm_n1p) = trn(:,:,:,jp_fabm_m1+jp_fabm_n1p) + ppo4_bkginc(:,:,:) * zincwgt 2626 trb(:,:,:,jp_fabm_m1+jp_fabm_n1p) = trb(:,:,:,jp_fabm_m1+jp_fabm_n1p) + ppo4_bkginc(:,:,:) * zincwgt 2627 2627 END WHERE 2628 2628 #else … … 2646 2646 #elif defined key_fabm 2647 2647 WHERE( pdic_bkginc(:,:,:) > 0.0_wp .OR. & 2648 & trn(:,:,:,jp_fabm_ o3c) + pdic_bkginc(:,:,:) * zincwgt > 0.0_wp )2649 trn(:,:,:,jp_fabm_ o3c) = trn(:,:,:,jp_fabm_o3c) + pdic_bkginc(:,:,:) * zincwgt2650 trb(:,:,:,jp_fabm_ o3c) = trb(:,:,:,jp_fabm_o3c) + pdic_bkginc(:,:,:) * zincwgt2648 & trn(:,:,:,jp_fabm_m1+jp_fabm_o3c) + pdic_bkginc(:,:,:) * zincwgt > 0.0_wp ) 2649 trn(:,:,:,jp_fabm_m1+jp_fabm_o3c) = trn(:,:,:,jp_fabm_m1+jp_fabm_o3c) + pdic_bkginc(:,:,:) * zincwgt 2650 trb(:,:,:,jp_fabm_m1+jp_fabm_o3c) = trb(:,:,:,jp_fabm_m1+jp_fabm_o3c) + pdic_bkginc(:,:,:) * zincwgt 2651 2651 END WHERE 2652 2652 #else … … 2670 2670 #elif defined key_fabm 2671 2671 WHERE( palk_bkginc(:,:,:) > 0.0_wp .OR. & 2672 & trn(:,:,:,jp_fabm_ o3ba) + palk_bkginc(:,:,:) * zincwgt > 0.0_wp )2673 trn(:,:,:,jp_fabm_ o3ba) = trn(:,:,:,jp_fabm_o3ba) + palk_bkginc(:,:,:) * zincwgt2674 trb(:,:,:,jp_fabm_ o3ba) = trb(:,:,:,jp_fabm_o3ba) + palk_bkginc(:,:,:) * zincwgt2672 & trn(:,:,:,jp_fabm_m1+jp_fabm_o3ba) + palk_bkginc(:,:,:) * zincwgt > 0.0_wp ) 2673 trn(:,:,:,jp_fabm_m1+jp_fabm_o3ba) = trn(:,:,:,jp_fabm_m1+jp_fabm_o3ba) + palk_bkginc(:,:,:) * zincwgt 2674 trb(:,:,:,jp_fabm_m1+jp_fabm_o3ba) = trb(:,:,:,jp_fabm_m1+jp_fabm_o3ba) + palk_bkginc(:,:,:) * zincwgt 2675 2675 END WHERE 2676 2676 #else … … 2688 2688 #elif defined key_fabm 2689 2689 WHERE( po2_bkginc(:,:,:) > 0.0_wp .OR. & 2690 & trn(:,:,:,jp_fabm_ o2o) + po2_bkginc(:,:,:) * zincwgt > 0.0_wp )2691 trn(:,:,:,jp_fabm_ o2o) = trn(:,:,:,jp_fabm_o2o) + po2_bkginc(:,:,:) * zincwgt2692 trb(:,:,:,jp_fabm_ o2o) = trb(:,:,:,jp_fabm_o2o) + po2_bkginc(:,:,:) * zincwgt2690 & trn(:,:,:,jp_fabm_m1+jp_fabm_o2o) + po2_bkginc(:,:,:) * zincwgt > 0.0_wp ) 2691 trn(:,:,:,jp_fabm_m1+jp_fabm_o2o) = trn(:,:,:,jp_fabm_m1+jp_fabm_o2o) + po2_bkginc(:,:,:) * zincwgt 2692 trb(:,:,:,jp_fabm_m1+jp_fabm_o2o) = trb(:,:,:,jp_fabm_m1+jp_fabm_o2o) + po2_bkginc(:,:,:) * zincwgt 2693 2693 END WHERE 2694 2694 #else … … 2738 2738 #elif defined key_fabm 2739 2739 WHERE( pno3_bkginc(:,:,:) > 0.0_wp .OR. & 2740 & trn(:,:,:,jp_fabm_ n3n) + pno3_bkginc(:,:,:) > 0.0_wp )2741 trn(:,:,:,jp_fabm_ n3n) = trn(:,:,:,jp_fabm_n3n) + pno3_bkginc(:,:,:)2742 trb(:,:,:,jp_fabm_ n3n) = trn(:,:,:,jp_fabm_n3n)2740 & trn(:,:,:,jp_fabm_m1+jp_fabm_n3n) + pno3_bkginc(:,:,:) > 0.0_wp ) 2741 trn(:,:,:,jp_fabm_m1+jp_fabm_n3n) = trn(:,:,:,jp_fabm_m1+jp_fabm_n3n) + pno3_bkginc(:,:,:) 2742 trb(:,:,:,jp_fabm_m1+jp_fabm_n3n) = trn(:,:,:,jp_fabm_m1+jp_fabm_n3n) 2743 2743 END WHERE 2744 2744 #else … … 2756 2756 #elif defined key_fabm 2757 2757 WHERE( psi4_bkginc(:,:,:) > 0.0_wp .OR. & 2758 & trn(:,:,:,jp_fabm_ n5s) + psi4_bkginc(:,:,:) > 0.0_wp )2759 trn(:,:,:,jp_fabm_ n5s) = trn(:,:,:,jp_fabm_n5s) + psi4_bkginc(:,:,:)2760 trb(:,:,:,jp_fabm_ n5s) = trn(:,:,:,jp_fabm_n5s)2758 & trn(:,:,:,jp_fabm_m1+jp_fabm_n5s) + psi4_bkginc(:,:,:) > 0.0_wp ) 2759 trn(:,:,:,jp_fabm_m1+jp_fabm_n5s) = trn(:,:,:,jp_fabm_m1+jp_fabm_n5s) + psi4_bkginc(:,:,:) 2760 trb(:,:,:,jp_fabm_m1+jp_fabm_n5s) = trn(:,:,:,jp_fabm_m1+jp_fabm_n5s) 2761 2761 END WHERE 2762 2762 #else … … 2768 2768 #if defined key_fabm 2769 2769 WHERE( ppo4_bkginc(:,:,:) > 0.0_wp .OR. & 2770 & trn(:,:,:,jp_fabm_ n1p) + ppo4_bkginc(:,:,:) > 0.0_wp )2771 trn(:,:,:,jp_fabm_ n1p) = trn(:,:,:,jp_fabm_n1p) + ppo4_bkginc(:,:,:)2772 trb(:,:,:,jp_fabm_ n1p) = trn(:,:,:,jp_fabm_n1p)2770 & trn(:,:,:,jp_fabm_m1+jp_fabm_n1p) + ppo4_bkginc(:,:,:) > 0.0_wp ) 2771 trn(:,:,:,jp_fabm_m1+jp_fabm_n1p) = trn(:,:,:,jp_fabm_m1+jp_fabm_n1p) + ppo4_bkginc(:,:,:) 2772 trb(:,:,:,jp_fabm_m1+jp_fabm_n1p) = trn(:,:,:,jp_fabm_m1+jp_fabm_n1p) 2773 2773 END WHERE 2774 2774 #else … … 2792 2792 #elif defined key_fabm 2793 2793 WHERE( pdic_bkginc(:,:,:) > 0.0_wp .OR. & 2794 & trn(:,:,:,jp_fabm_ o3c) + pdic_bkginc(:,:,:) > 0.0_wp )2795 trn(:,:,:,jp_fabm_ o3c) = trn(:,:,:,jp_fabm_o3c) + pdic_bkginc(:,:,:)2796 trb(:,:,:,jp_fabm_ o3c) = trn(:,:,:,jp_fabm_o3c)2794 & trn(:,:,:,jp_fabm_m1+jp_fabm_o3c) + pdic_bkginc(:,:,:) > 0.0_wp ) 2795 trn(:,:,:,jp_fabm_m1+jp_fabm_o3c) = trn(:,:,:,jp_fabm_m1+jp_fabm_o3c) + pdic_bkginc(:,:,:) 2796 trb(:,:,:,jp_fabm_m1+jp_fabm_o3c) = trn(:,:,:,jp_fabm_m1+jp_fabm_o3c) 2797 2797 END WHERE 2798 2798 #else … … 2816 2816 #elif defined key_fabm 2817 2817 WHERE( palk_bkginc(:,:,:) > 0.0_wp .OR. & 2818 & trn(:,:,:,jp_fabm_ o3ba) + palk_bkginc(:,:,:) > 0.0_wp )2819 trn(:,:,:,jp_fabm_ o3ba) = trn(:,:,:,jp_fabm_o3ba) + palk_bkginc(:,:,:)2820 trb(:,:,:,jp_fabm_ o3ba) = trn(:,:,:,jp_fabm_o3ba)2818 & trn(:,:,:,jp_fabm_m1+jp_fabm_o3ba) + palk_bkginc(:,:,:) > 0.0_wp ) 2819 trn(:,:,:,jp_fabm_m1+jp_fabm_o3ba) = trn(:,:,:,jp_fabm_m1+jp_fabm_o3ba) + palk_bkginc(:,:,:) 2820 trb(:,:,:,jp_fabm_m1+jp_fabm_o3ba) = trn(:,:,:,jp_fabm_m1+jp_fabm_o3ba) 2821 2821 END WHERE 2822 2822 #else … … 2834 2834 #elif defined key_fabm 2835 2835 WHERE( po2_bkginc(:,:,:) > 0.0_wp .OR. & 2836 & trn(:,:,:,jp_fabm_ o2o) + po2_bkginc(:,:,:) > 0.0_wp )2837 trn(:,:,:,jp_fabm_ o2o) = trn(:,:,:,jp_fabm_o2o) + po2_bkginc(:,:,:)2838 trb(:,:,:,jp_fabm_ o2o) = trn(:,:,:,jp_fabm_o2o)2836 & trn(:,:,:,jp_fabm_m1+jp_fabm_o2o) + po2_bkginc(:,:,:) > 0.0_wp ) 2837 trn(:,:,:,jp_fabm_m1+jp_fabm_o2o) = trn(:,:,:,jp_fabm_m1+jp_fabm_o2o) + po2_bkginc(:,:,:) 2838 trb(:,:,:,jp_fabm_m1+jp_fabm_o2o) = trn(:,:,:,jp_fabm_m1+jp_fabm_o2o) 2839 2839 END WHERE 2840 2840 #else -
branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/OPA_SRC/ASM/asmphyto2dbal_ersem.F90
r10664 r10667 339 339 340 340 ! Set background state 341 bstate(:,:,:,i_tracer(1)) = tracer_bkg(:,:,:,jp_fabm_ n3n) + &342 & tracer_bkg(:,:,:,jp_fabm_ n4n)343 bstate(:,:,:,i_tracer(2)) = tracer_bkg(:,:,:,jp_fabm_ p1n) + &344 & tracer_bkg(:,:,:,jp_fabm_ p2n) + &345 & tracer_bkg(:,:,:,jp_fabm_ p3n) + &346 & tracer_bkg(:,:,:,jp_fabm_ p4n)347 bstate(:,:,:,i_tracer(3)) = (tracer_bkg(:,:,:,jp_fabm_ z4c) * z4qnc) + &348 & tracer_bkg(:,:,:,jp_fabm_ z5n) + &349 & tracer_bkg(:,:,:,jp_fabm_ z6n)350 bstate(:,:,:,i_tracer(4)) = tracer_bkg(:,:,:,jp_fabm_ r4n) + &351 & tracer_bkg(:,:,:,jp_fabm_ r6n) + &352 & tracer_bkg(:,:,:,jp_fabm_ r8n)353 bstate(:,:,:,i_tracer(5)) = tracer_bkg(:,:,:,jp_fabm_ o3c)341 bstate(:,:,:,i_tracer(1)) = tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_n3n) + & 342 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_n4n) 343 bstate(:,:,:,i_tracer(2)) = tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p1n) + & 344 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p2n) + & 345 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p3n) + & 346 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_p4n) 347 bstate(:,:,:,i_tracer(3)) = (tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z4c) * z4qnc) + & 348 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z5n) + & 349 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_z6n) 350 bstate(:,:,:,i_tracer(4)) = tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_r4n) + & 351 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_r6n) + & 352 & tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_r8n) 353 bstate(:,:,:,i_tracer(5)) = tracer_bkg(:,:,:,jp_fabm_m1+jp_fabm_o3c) 354 354 bstate(:,:,:,i_tracer(6)) = totalk_bkg(:,:,:) 355 355 … … 359 359 DO jj = 1, jpj 360 360 DO ji = 1, jpi 361 IF ( ( tracer_bkg(ji,jj,1,jp_fabm_ chl1) + tracer_bkg(ji,jj,1,jp_fabm_chl2) + &362 & tracer_bkg(ji,jj,1,jp_fabm_ chl3) + tracer_bkg(ji,jj,1,jp_fabm_chl4) ) .GT. 0.0 ) THEN363 cchl_p(ji,jj) = zmassc * ( tracer_bkg(ji,jj,1,jp_fabm_ p1c) + &364 & tracer_bkg(ji,jj,1,jp_fabm_ p2c) + &365 & tracer_bkg(ji,jj,1,jp_fabm_ p3c) + &366 & tracer_bkg(ji,jj,1,jp_fabm_ p4c) ) / &367 & ( tracer_bkg(ji,jj,1,jp_fabm_ chl1) + &368 & tracer_bkg(ji,jj,1,jp_fabm_ chl2) + &369 & tracer_bkg(ji,jj,1,jp_fabm_ chl3) + &370 & tracer_bkg(ji,jj,1,jp_fabm_ chl4) )361 IF ( ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + & 362 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) ) .GT. 0.0 ) THEN 363 cchl_p(ji,jj) = zmassc * ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p1c) + & 364 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p2c) + & 365 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p3c) + & 366 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p4c) ) / & 367 & ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + & 368 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + & 369 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + & 370 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) ) 371 371 ENDIF 372 372 END DO … … 395 395 396 396 ! Phytoplankton 397 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_ p1n) > 0.0 ) .AND. &398 & ( tracer_bkg(ji,jj,jk,jp_fabm_ p2n) > 0.0 ) .AND. &399 & ( tracer_bkg(ji,jj,jk,jp_fabm_ p3n) > 0.0 ) .AND. &400 & ( tracer_bkg(ji,jj,jk,jp_fabm_ p4n) > 0.0 ) .AND. &397 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) > 0.0 ) .AND. & 398 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) > 0.0 ) .AND. & 399 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) > 0.0 ) .AND. & 400 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) > 0.0 ) .AND. & 401 401 & ( pinc_chltot(ji,jj) /= 0.0 ) ) THEN 402 402 IF ( ld_chltot ) THEN 403 403 ! Phytoplankton nitrogen split up based on existing ratios 404 zfrac_p1n = tracer_bkg(ji,jj,jk,jp_fabm_ p1n) / &405 & ( tracer_bkg(ji,jj,jk,jp_fabm_ p1n) + &406 & tracer_bkg(ji,jj,jk,jp_fabm_ p2n) + &407 & tracer_bkg(ji,jj,jk,jp_fabm_ p3n) + &408 & tracer_bkg(ji,jj,jk,jp_fabm_ p4n) )409 zfrac_p2n = tracer_bkg(ji,jj,jk,jp_fabm_ p2n) / &410 & ( tracer_bkg(ji,jj,jk,jp_fabm_ p1n) + &411 & tracer_bkg(ji,jj,jk,jp_fabm_ p2n) + &412 & tracer_bkg(ji,jj,jk,jp_fabm_ p3n) + &413 & tracer_bkg(ji,jj,jk,jp_fabm_ p4n) )414 zfrac_p3n = tracer_bkg(ji,jj,jk,jp_fabm_ p3n) / &415 & ( tracer_bkg(ji,jj,jk,jp_fabm_ p1n) + &416 & tracer_bkg(ji,jj,jk,jp_fabm_ p2n) + &417 & tracer_bkg(ji,jj,jk,jp_fabm_ p3n) + &418 & tracer_bkg(ji,jj,jk,jp_fabm_ p4n) )419 zfrac_p4n = tracer_bkg(ji,jj,jk,jp_fabm_ p4n) / &420 & ( tracer_bkg(ji,jj,jk,jp_fabm_ p1n) + &421 & tracer_bkg(ji,jj,jk,jp_fabm_ p2n) + &422 & tracer_bkg(ji,jj,jk,jp_fabm_ p3n) + &423 & tracer_bkg(ji,jj,jk,jp_fabm_ p4n) )404 zfrac_p1n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) / & 405 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) + & 406 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) + & 407 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) + & 408 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) ) 409 zfrac_p2n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) / & 410 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) + & 411 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) + & 412 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) + & 413 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) ) 414 zfrac_p3n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) / & 415 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) + & 416 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) + & 417 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) + & 418 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) ) 419 zfrac_p4n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) / & 420 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) + & 421 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) + & 422 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) + & 423 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) ) 424 424 ELSE 425 425 ! Phytoplankton nitrogen split up based on assimilation increments … … 431 431 432 432 ! Other phytoplankton variables split up based on existing ratios with nitrogen 433 zrat_chl1_p1n = tracer_bkg(ji,jj,jk,jp_fabm_ chl1) / tracer_bkg(ji,jj,jk,jp_fabm_p1n)434 zrat_p1c_p1n = tracer_bkg(ji,jj,jk,jp_fabm_ p1c) / tracer_bkg(ji,jj,jk,jp_fabm_p1n)435 zrat_p1p_p1n = tracer_bkg(ji,jj,jk,jp_fabm_ p1p) / tracer_bkg(ji,jj,jk,jp_fabm_p1n)436 zrat_p1s_p1n = tracer_bkg(ji,jj,jk,jp_fabm_ p1s) / tracer_bkg(ji,jj,jk,jp_fabm_p1n)437 zrat_chl2_p2n = tracer_bkg(ji,jj,jk,jp_fabm_ chl2) / tracer_bkg(ji,jj,jk,jp_fabm_p2n)438 zrat_p2c_p2n = tracer_bkg(ji,jj,jk,jp_fabm_ p2c) / tracer_bkg(ji,jj,jk,jp_fabm_p2n)439 zrat_p2p_p2n = tracer_bkg(ji,jj,jk,jp_fabm_ p2p) / tracer_bkg(ji,jj,jk,jp_fabm_p2n)440 zrat_chl3_p3n = tracer_bkg(ji,jj,jk,jp_fabm_ chl3) / tracer_bkg(ji,jj,jk,jp_fabm_p3n)441 zrat_p3c_p3n = tracer_bkg(ji,jj,jk,jp_fabm_ p3c) / tracer_bkg(ji,jj,jk,jp_fabm_p3n)442 zrat_p3p_p3n = tracer_bkg(ji,jj,jk,jp_fabm_ p3p) / tracer_bkg(ji,jj,jk,jp_fabm_p3n)443 zrat_chl4_p4n = tracer_bkg(ji,jj,jk,jp_fabm_ chl4) / tracer_bkg(ji,jj,jk,jp_fabm_p4n)444 zrat_p4c_p4n = tracer_bkg(ji,jj,jk,jp_fabm_ p4c) / tracer_bkg(ji,jj,jk,jp_fabm_p4n)445 zrat_p4p_p4n = tracer_bkg(ji,jj,jk,jp_fabm_ p4p) / tracer_bkg(ji,jj,jk,jp_fabm_p4n)433 zrat_chl1_p1n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) 434 zrat_p1c_p1n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) 435 zrat_p1p_p1n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) 436 zrat_p1s_p1n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) 437 zrat_chl2_p2n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) 438 zrat_p2c_p2n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) 439 zrat_p2p_p2n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) 440 zrat_chl3_p3n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) 441 zrat_p3c_p3n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) 442 zrat_p3p_p3n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) 443 zrat_chl4_p4n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) 444 zrat_p4c_p4n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) 445 zrat_p4p_p4n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) 446 446 447 phyto2d_balinc(ji,jj,jk,jp_fabm_ p1n) = outincs(ji,jj,jk,i_tracer(2)) * zfrac_p1n448 phyto2d_balinc(ji,jj,jk,jp_fabm_ p2n) = outincs(ji,jj,jk,i_tracer(2)) * zfrac_p2n449 phyto2d_balinc(ji,jj,jk,jp_fabm_ p3n) = outincs(ji,jj,jk,i_tracer(2)) * zfrac_p3n450 phyto2d_balinc(ji,jj,jk,jp_fabm_ p4n) = outincs(ji,jj,jk,i_tracer(2)) * zfrac_p4n451 phyto2d_balinc(ji,jj,jk,jp_fabm_ chl1) = phyto2d_balinc(ji,jj,jk,jp_fabm_p1n) * zrat_chl1_p1n452 phyto2d_balinc(ji,jj,jk,jp_fabm_ p1c) = phyto2d_balinc(ji,jj,jk,jp_fabm_p1n) * zrat_p1c_p1n453 phyto2d_balinc(ji,jj,jk,jp_fabm_ p1p) = phyto2d_balinc(ji,jj,jk,jp_fabm_p1n) * zrat_p1p_p1n454 phyto2d_balinc(ji,jj,jk,jp_fabm_ p1s) = phyto2d_balinc(ji,jj,jk,jp_fabm_p1n) * zrat_p1s_p1n455 phyto2d_balinc(ji,jj,jk,jp_fabm_ chl2) = phyto2d_balinc(ji,jj,jk,jp_fabm_p2n) * zrat_chl2_p2n456 phyto2d_balinc(ji,jj,jk,jp_fabm_ p2c) = phyto2d_balinc(ji,jj,jk,jp_fabm_p2n) * zrat_p2c_p2n457 phyto2d_balinc(ji,jj,jk,jp_fabm_ p2p) = phyto2d_balinc(ji,jj,jk,jp_fabm_p2n) * zrat_p2p_p2n458 phyto2d_balinc(ji,jj,jk,jp_fabm_ chl3) = phyto2d_balinc(ji,jj,jk,jp_fabm_p3n) * zrat_chl3_p3n459 phyto2d_balinc(ji,jj,jk,jp_fabm_ p3c) = phyto2d_balinc(ji,jj,jk,jp_fabm_p3n) * zrat_p3c_p3n460 phyto2d_balinc(ji,jj,jk,jp_fabm_ p3p) = phyto2d_balinc(ji,jj,jk,jp_fabm_p3n) * zrat_p3p_p3n461 phyto2d_balinc(ji,jj,jk,jp_fabm_ chl4) = phyto2d_balinc(ji,jj,jk,jp_fabm_p4n) * zrat_chl4_p4n462 phyto2d_balinc(ji,jj,jk,jp_fabm_ p4c) = phyto2d_balinc(ji,jj,jk,jp_fabm_p4n) * zrat_p4c_p4n463 phyto2d_balinc(ji,jj,jk,jp_fabm_ p4p) = phyto2d_balinc(ji,jj,jk,jp_fabm_p4n) * zrat_p4p_p4n447 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) = outincs(ji,jj,jk,i_tracer(2)) * zfrac_p1n 448 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) = outincs(ji,jj,jk,i_tracer(2)) * zfrac_p2n 449 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) = outincs(ji,jj,jk,i_tracer(2)) * zfrac_p3n 450 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) = outincs(ji,jj,jk,i_tracer(2)) * zfrac_p4n 451 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) * zrat_chl1_p1n 452 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) * zrat_p1c_p1n 453 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) * zrat_p1p_p1n 454 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) * zrat_p1s_p1n 455 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) * zrat_chl2_p2n 456 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) * zrat_p2c_p2n 457 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) * zrat_p2p_p2n 458 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) * zrat_chl3_p3n 459 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) * zrat_p3c_p3n 460 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) * zrat_p3p_p3n 461 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) * zrat_chl4_p4n 462 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) * zrat_p4c_p4n 463 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) * zrat_p4p_p4n 464 464 ENDIF 465 465 466 466 ! Zooplankton nitrogen split up based on existing ratios 467 467 ! Update carbon and phosphorus according to existing ratios 468 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_ z4c) > 0.0 ) .AND. &469 & ( tracer_bkg(ji,jj,jk,jp_fabm_ z5n) > 0.0 ) .AND. &470 & ( tracer_bkg(ji,jj,jk,jp_fabm_ z6n) > 0.0 ) ) THEN471 zfrac_z4n = ( tracer_bkg(ji,jj,jk,jp_fabm_ z4c) * z4qnc ) / &472 & ( ( tracer_bkg(ji,jj,jk,jp_fabm_ z4c) * z4qnc ) + &473 & tracer_bkg(ji,jj,jk,jp_fabm_ z5n) + &474 & tracer_bkg(ji,jj,jk,jp_fabm_ z6n) )475 zfrac_z5n = tracer_bkg(ji,jj,jk,jp_fabm_ z5n) / &476 & ( ( tracer_bkg(ji,jj,jk,jp_fabm_ z4c) * z4qnc ) + &477 & tracer_bkg(ji,jj,jk,jp_fabm_ z5n) + &478 & tracer_bkg(ji,jj,jk,jp_fabm_ z6n) )479 zfrac_z6n = tracer_bkg(ji,jj,jk,jp_fabm_ z6n) / &480 & ( ( tracer_bkg(ji,jj,jk,jp_fabm_ z4c) * z4qnc ) + &481 & tracer_bkg(ji,jj,jk,jp_fabm_ z5n) + &482 & tracer_bkg(ji,jj,jk,jp_fabm_ z6n) )468 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z4c) > 0.0 ) .AND. & 469 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) > 0.0 ) .AND. & 470 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) > 0.0 ) ) THEN 471 zfrac_z4n = ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z4c) * z4qnc ) / & 472 & ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z4c) * z4qnc ) + & 473 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) + & 474 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) ) 475 zfrac_z5n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) / & 476 & ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z4c) * z4qnc ) + & 477 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) + & 478 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) ) 479 zfrac_z6n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) / & 480 & ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z4c) * z4qnc ) + & 481 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) + & 482 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) ) 483 483 zrat_z4c_z4n = 1.0 / z4qnc 484 zrat_z5c_z5n = tracer_bkg(ji,jj,jk,jp_fabm_ z5c) / tracer_bkg(ji,jj,jk,jp_fabm_z5n)485 zrat_z5p_z5n = tracer_bkg(ji,jj,jk,jp_fabm_ z5p) / tracer_bkg(ji,jj,jk,jp_fabm_z5n)486 zrat_z6c_z6n = tracer_bkg(ji,jj,jk,jp_fabm_ z6c) / tracer_bkg(ji,jj,jk,jp_fabm_z6n)487 zrat_z6p_z6n = tracer_bkg(ji,jj,jk,jp_fabm_ z6p) / tracer_bkg(ji,jj,jk,jp_fabm_z6n)488 phyto2d_balinc(ji,jj,jk,jp_fabm_ z5n) = outincs(ji,jj,jk,i_tracer(3)) * zfrac_z5n489 phyto2d_balinc(ji,jj,jk,jp_fabm_ z6n) = outincs(ji,jj,jk,i_tracer(3)) * zfrac_z6n490 phyto2d_balinc(ji,jj,jk,jp_fabm_ z4c) = outincs(ji,jj,jk,i_tracer(3)) * zfrac_z4n * zrat_z4c_z4n491 phyto2d_balinc(ji,jj,jk,jp_fabm_ z5c) = phyto2d_balinc(ji,jj,jk,jp_fabm_z5n) * zrat_z5c_z5n492 phyto2d_balinc(ji,jj,jk,jp_fabm_ z6c) = phyto2d_balinc(ji,jj,jk,jp_fabm_z6n) * zrat_z6c_z6n493 phyto2d_balinc(ji,jj,jk,jp_fabm_ z5p) = phyto2d_balinc(ji,jj,jk,jp_fabm_z5n) * zrat_z5p_z5n494 phyto2d_balinc(ji,jj,jk,jp_fabm_ z6p) = phyto2d_balinc(ji,jj,jk,jp_fabm_z6n) * zrat_z6p_z6n484 zrat_z5c_z5n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z5c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) 485 zrat_z5p_z5n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z5p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) 486 zrat_z6c_z6n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z6c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) 487 zrat_z6p_z6n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z6p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) 488 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) = outincs(ji,jj,jk,i_tracer(3)) * zfrac_z5n 489 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) = outincs(ji,jj,jk,i_tracer(3)) * zfrac_z6n 490 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z4c) = outincs(ji,jj,jk,i_tracer(3)) * zfrac_z4n * zrat_z4c_z4n 491 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z5c) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) * zrat_z5c_z5n 492 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z6c) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) * zrat_z6c_z6n 493 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z5p) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z5n) * zrat_z5p_z5n 494 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z6p) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z6n) * zrat_z6p_z6n 495 495 ENDIF 496 496 497 497 ! Nitrogen nutrient split between nitrate and ammonium based on existing ratios 498 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_ n3n) > 0.0 ) .AND. &499 & ( tracer_bkg(ji,jj,jk,jp_fabm_ n4n) > 0.0 ) ) THEN500 zfrac_n3n = tracer_bkg(ji,jj,jk,jp_fabm_ n3n) / &501 & (tracer_bkg(ji,jj,jk,jp_fabm_ n3n) + tracer_bkg(ji,jj,jk,jp_fabm_n4n))502 zfrac_n4n = tracer_bkg(ji,jj,jk,jp_fabm_ n4n) / &503 & (tracer_bkg(ji,jj,jk,jp_fabm_ n3n) + tracer_bkg(ji,jj,jk,jp_fabm_n4n))504 phyto2d_balinc(ji,jj,jk,jp_fabm_ n3n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_n3n505 phyto2d_balinc(ji,jj,jk,jp_fabm_ n4n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_n4n498 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n3n) > 0.0 ) .AND. & 499 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n4n) > 0.0 ) ) THEN 500 zfrac_n3n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n3n) / & 501 & (tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n3n) + tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n4n)) 502 zfrac_n4n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n4n) / & 503 & (tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n3n) + tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n4n)) 504 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_n3n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_n3n 505 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_n4n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_n4n 506 506 ENDIF 507 507 508 508 ! Detritus nitrogen split up based on existing ratios 509 509 ! Update carbon and phosphorus according to existing ratios 510 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_ r4n) > 0.0 ) .AND. &511 & ( tracer_bkg(ji,jj,jk,jp_fabm_ r6n) > 0.0 ) .AND. &512 & ( tracer_bkg(ji,jj,jk,jp_fabm_ r8n) > 0.0 ) ) THEN513 zfrac_r4n = tracer_bkg(ji,jj,jk,jp_fabm_ r4n) / &514 & (tracer_bkg(ji,jj,jk,jp_fabm_ r4n) + &515 & tracer_bkg(ji,jj,jk,jp_fabm_ r6n) + &516 & tracer_bkg(ji,jj,jk,jp_fabm_ r8n))517 zfrac_r6n = tracer_bkg(ji,jj,jk,jp_fabm_ r6n) / &518 & (tracer_bkg(ji,jj,jk,jp_fabm_ r4n) + &519 & tracer_bkg(ji,jj,jk,jp_fabm_ r6n) + &520 & tracer_bkg(ji,jj,jk,jp_fabm_ r8n))521 zfrac_r8n = tracer_bkg(ji,jj,jk,jp_fabm_ r8n) / &522 & (tracer_bkg(ji,jj,jk,jp_fabm_ r4n) + &523 & tracer_bkg(ji,jj,jk,jp_fabm_ r6n) + &524 & tracer_bkg(ji,jj,jk,jp_fabm_ r8n))525 zrat_r4c_r4n = tracer_bkg(ji,jj,jk,jp_fabm_ r4c) / tracer_bkg(ji,jj,jk,jp_fabm_r4n)526 zrat_r4p_r4n = tracer_bkg(ji,jj,jk,jp_fabm_ r4p) / tracer_bkg(ji,jj,jk,jp_fabm_r4n)527 zrat_r6c_r6n = tracer_bkg(ji,jj,jk,jp_fabm_ r6c) / tracer_bkg(ji,jj,jk,jp_fabm_r6n)528 zrat_r6p_r6n = tracer_bkg(ji,jj,jk,jp_fabm_ r6p) / tracer_bkg(ji,jj,jk,jp_fabm_r6n)529 zrat_r6s_r6n = tracer_bkg(ji,jj,jk,jp_fabm_ r6s) / tracer_bkg(ji,jj,jk,jp_fabm_r6n)530 zrat_r8c_r8n = tracer_bkg(ji,jj,jk,jp_fabm_ r8c) / tracer_bkg(ji,jj,jk,jp_fabm_r8n)531 zrat_r8p_r8n = tracer_bkg(ji,jj,jk,jp_fabm_ r8p) / tracer_bkg(ji,jj,jk,jp_fabm_r8n)532 zrat_r8s_r8n = tracer_bkg(ji,jj,jk,jp_fabm_ r8s) / tracer_bkg(ji,jj,jk,jp_fabm_r8n)533 phyto2d_balinc(ji,jj,jk,jp_fabm_ r4n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_r4n534 phyto2d_balinc(ji,jj,jk,jp_fabm_ r6n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_r6n535 phyto2d_balinc(ji,jj,jk,jp_fabm_ r8n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_r8n536 phyto2d_balinc(ji,jj,jk,jp_fabm_ r4c) = phyto2d_balinc(ji,jj,jk,jp_fabm_r4n) * zrat_r4c_r4n537 phyto2d_balinc(ji,jj,jk,jp_fabm_ r4p) = phyto2d_balinc(ji,jj,jk,jp_fabm_r4n) * zrat_r4p_r4n538 phyto2d_balinc(ji,jj,jk,jp_fabm_ r6c) = phyto2d_balinc(ji,jj,jk,jp_fabm_r6n) * zrat_r6c_r6n539 phyto2d_balinc(ji,jj,jk,jp_fabm_ r6p) = phyto2d_balinc(ji,jj,jk,jp_fabm_r6n) * zrat_r6p_r6n540 phyto2d_balinc(ji,jj,jk,jp_fabm_ r6s) = phyto2d_balinc(ji,jj,jk,jp_fabm_r6n) * zrat_r6s_r6n541 phyto2d_balinc(ji,jj,jk,jp_fabm_ r8c) = phyto2d_balinc(ji,jj,jk,jp_fabm_r8n) * zrat_r8c_r8n542 phyto2d_balinc(ji,jj,jk,jp_fabm_ r8p) = phyto2d_balinc(ji,jj,jk,jp_fabm_r8n) * zrat_r8p_r8n543 phyto2d_balinc(ji,jj,jk,jp_fabm_ r8s) = phyto2d_balinc(ji,jj,jk,jp_fabm_r8n) * zrat_r8s_r8n510 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) > 0.0 ) .AND. & 511 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) > 0.0 ) .AND. & 512 & ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n) > 0.0 ) ) THEN 513 zfrac_r4n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) / & 514 & (tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) + & 515 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) + & 516 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n)) 517 zfrac_r6n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) / & 518 & (tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) + & 519 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) + & 520 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n)) 521 zfrac_r8n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n) / & 522 & (tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) + & 523 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) + & 524 & tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n)) 525 zrat_r4c_r4n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r4c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) 526 zrat_r4p_r4n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r4p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) 527 zrat_r6c_r6n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) 528 zrat_r6p_r6n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) 529 zrat_r6s_r6n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6s) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) 530 zrat_r8c_r8n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8c) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n) 531 zrat_r8p_r8n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8p) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n) 532 zrat_r8s_r8n = tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8s) / tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n) 533 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_r4n 534 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_r6n 535 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n) = outincs(ji,jj,jk,i_tracer(1)) * zfrac_r8n 536 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r4c) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) * zrat_r4c_r4n 537 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r4p) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r4n) * zrat_r4p_r4n 538 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r6c) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) * zrat_r6c_r6n 539 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r6p) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) * zrat_r6p_r6n 540 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r6s) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r6n) * zrat_r6s_r6n 541 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r8c) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n) * zrat_r8c_r8n 542 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r8p) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n) * zrat_r8p_r8n 543 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r8s) = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r8n) * zrat_r8s_r8n 544 544 ENDIF 545 545 546 546 ! DIC straight from balancing scheme 547 phyto2d_balinc(ji,jj,jk,jp_fabm_ o3c) = outincs(ji,jj,jk,i_tracer(5))547 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_o3c) = outincs(ji,jj,jk,i_tracer(5)) 548 548 549 549 ! Alkalinity straight from balancing scheme 550 phyto2d_balinc(ji,jj,jk,jp_fabm_ o3ba) = outincs(ji,jj,jk,i_tracer(6))550 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_o3ba) = outincs(ji,jj,jk,i_tracer(6)) 551 551 552 552 ! Remove P/R silicon increments from silicate to conserve mass 553 zfrac = phyto2d_balinc(ji,jj,jk,jp_fabm_ p1s) + &554 & phyto2d_balinc(ji,jj,jk,jp_fabm_ r6s) + &555 & phyto2d_balinc(ji,jj,jk,jp_fabm_ r8s)556 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_ n5s) - zfrac ) > 0.0 ) THEN557 phyto2d_balinc(ji,jj,jk,jp_fabm_ n5s) = zfrac * (-1.0)553 zfrac = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) + & 554 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r6s) + & 555 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r8s) 556 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n5s) - zfrac ) > 0.0 ) THEN 557 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_n5s) = zfrac * (-1.0) 558 558 ENDIF 559 559 560 560 ! Remove P/Z/R phosphorus increments from phosphate to conserve mass 561 zfrac = phyto2d_balinc(ji,jj,jk,jp_fabm_ p1p) + &562 & phyto2d_balinc(ji,jj,jk,jp_fabm_ p2p) + &563 & phyto2d_balinc(ji,jj,jk,jp_fabm_ p3p) + &564 & phyto2d_balinc(ji,jj,jk,jp_fabm_ p4p) + &565 & phyto2d_balinc(ji,jj,jk,jp_fabm_ z5p) + &566 & phyto2d_balinc(ji,jj,jk,jp_fabm_ z6p) + &567 & phyto2d_balinc(ji,jj,jk,jp_fabm_ r4p) + &568 & phyto2d_balinc(ji,jj,jk,jp_fabm_ r6p) + &569 & phyto2d_balinc(ji,jj,jk,jp_fabm_ r8p)570 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_ n1p) - zfrac ) > 0.0 ) THEN571 phyto2d_balinc(ji,jj,jk,jp_fabm_ n1p) = zfrac * (-1.0)561 zfrac = phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) + & 562 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) + & 563 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) + & 564 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) + & 565 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z5p) + & 566 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_z6p) + & 567 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r4p) + & 568 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r6p) + & 569 & phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_r8p) 570 IF ( ( tracer_bkg(ji,jj,jk,jp_fabm_m1+jp_fabm_n1p) - zfrac ) > 0.0 ) THEN 571 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_n1p) = zfrac * (-1.0) 572 572 ENDIF 573 573 … … 581 581 DO jj = 1, jpj 582 582 DO ji = 1, jpi 583 IF ( ( tracer_bkg(ji,jj,1,jp_fabm_ chl1) > 0.0 ) .AND. &584 & ( tracer_bkg(ji,jj,1,jp_fabm_ chl2) > 0.0 ) .AND. &585 & ( tracer_bkg(ji,jj,1,jp_fabm_ chl3) > 0.0 ) .AND. &586 & ( tracer_bkg(ji,jj,1,jp_fabm_ chl4) > 0.0 ) ) THEN583 IF ( ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) > 0.0 ) .AND. & 584 & ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) > 0.0 ) .AND. & 585 & ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) > 0.0 ) .AND. & 586 & ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) > 0.0 ) ) THEN 587 587 IF ( ld_chltot ) THEN 588 588 ! Chlorophyll split up based on existing ratios 589 zfrac_chl1 = tracer_bkg(ji,jj,1,jp_fabm_ chl1) / &590 & ( tracer_bkg(ji,jj,1,jp_fabm_ chl1) + &591 & tracer_bkg(ji,jj,1,jp_fabm_ chl2) + &592 & tracer_bkg(ji,jj,1,jp_fabm_ chl3) + &593 & tracer_bkg(ji,jj,1,jp_fabm_ chl4) )594 zfrac_chl2 = tracer_bkg(ji,jj,1,jp_fabm_ chl2) / &595 & ( tracer_bkg(ji,jj,1,jp_fabm_ chl1) + &596 & tracer_bkg(ji,jj,1,jp_fabm_ chl2) + &597 & tracer_bkg(ji,jj,1,jp_fabm_ chl3) + &598 & tracer_bkg(ji,jj,1,jp_fabm_ chl4) )599 zfrac_chl3 = tracer_bkg(ji,jj,1,jp_fabm_ chl3) / &600 & ( tracer_bkg(ji,jj,1,jp_fabm_ chl1) + &601 & tracer_bkg(ji,jj,1,jp_fabm_ chl2) + &602 & tracer_bkg(ji,jj,1,jp_fabm_ chl3) + &603 & tracer_bkg(ji,jj,1,jp_fabm_ chl4) )604 zfrac_chl4 = tracer_bkg(ji,jj,1,jp_fabm_ chl4) / &605 & ( tracer_bkg(ji,jj,1,jp_fabm_ chl1) + &606 & tracer_bkg(ji,jj,1,jp_fabm_ chl2) + &607 & tracer_bkg(ji,jj,1,jp_fabm_ chl3) + &608 & tracer_bkg(ji,jj,1,jp_fabm_ chl4) )609 phyto2d_balinc(ji,jj,1,jp_fabm_ chl1) = pinc_chltot(ji,jj) * zfrac_chl1610 phyto2d_balinc(ji,jj,1,jp_fabm_ chl2) = pinc_chltot(ji,jj) * zfrac_chl2611 phyto2d_balinc(ji,jj,1,jp_fabm_ chl3) = pinc_chltot(ji,jj) * zfrac_chl3612 phyto2d_balinc(ji,jj,1,jp_fabm_ chl4) = pinc_chltot(ji,jj) * zfrac_chl4589 zfrac_chl1 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) / & 590 & ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + & 591 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + & 592 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + & 593 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) ) 594 zfrac_chl2 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) / & 595 & ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + & 596 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + & 597 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + & 598 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) ) 599 zfrac_chl3 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) / & 600 & ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + & 601 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + & 602 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + & 603 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) ) 604 zfrac_chl4 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) / & 605 & ( tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + & 606 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + & 607 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + & 608 & tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) ) 609 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) = pinc_chltot(ji,jj) * zfrac_chl1 610 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) = pinc_chltot(ji,jj) * zfrac_chl2 611 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) = pinc_chltot(ji,jj) * zfrac_chl3 612 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) = pinc_chltot(ji,jj) * zfrac_chl4 613 613 ENDIF 614 614 IF( ld_chldia ) THEN 615 phyto2d_balinc(ji,jj,1,jp_fabm_ chl1) = pinc_chldia(ji,jj)615 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) = pinc_chldia(ji,jj) 616 616 ENDIF 617 617 IF( ld_chlnan ) THEN 618 phyto2d_balinc(ji,jj,1,jp_fabm_ chl2) = pinc_chlnan(ji,jj)618 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) = pinc_chlnan(ji,jj) 619 619 ENDIF 620 620 IF( ld_chlpic ) THEN 621 phyto2d_balinc(ji,jj,1,jp_fabm_ chl3) = pinc_chlpic(ji,jj)621 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) = pinc_chlpic(ji,jj) 622 622 ENDIF 623 623 IF( ld_chldin ) THEN 624 phyto2d_balinc(ji,jj,1,jp_fabm_ chl4) = pinc_chldin(ji,jj)624 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) = pinc_chldin(ji,jj) 625 625 ENDIF 626 626 627 627 ! Maintain stoichiometric ratios of carbon, nitrogen, phosphorus and silicon 628 628 IF ( ld_chltot .OR. ld_chldia ) THEN 629 zrat_p1c_chl1 = tracer_bkg(ji,jj,1,jp_fabm_ p1c) / tracer_bkg(ji,jj,1,jp_fabm_chl1)630 zrat_p1n_chl1 = tracer_bkg(ji,jj,1,jp_fabm_ p1n) / tracer_bkg(ji,jj,1,jp_fabm_chl1)631 zrat_p1p_chl1 = tracer_bkg(ji,jj,1,jp_fabm_ p1p) / tracer_bkg(ji,jj,1,jp_fabm_chl1)632 zrat_p1s_chl1 = tracer_bkg(ji,jj,1,jp_fabm_ p1s) / tracer_bkg(ji,jj,1,jp_fabm_chl1)633 phyto2d_balinc(ji,jj,1,jp_fabm_ p1c) = phyto2d_balinc(ji,jj,1,jp_fabm_chl1) * zrat_p1c_chl1634 phyto2d_balinc(ji,jj,1,jp_fabm_ p1n) = phyto2d_balinc(ji,jj,1,jp_fabm_chl1) * zrat_p1n_chl1635 phyto2d_balinc(ji,jj,1,jp_fabm_ p1p) = phyto2d_balinc(ji,jj,1,jp_fabm_chl1) * zrat_p1p_chl1636 phyto2d_balinc(ji,jj,1,jp_fabm_ p1s) = phyto2d_balinc(ji,jj,1,jp_fabm_chl1) * zrat_p1s_chl1629 zrat_p1c_chl1 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p1c) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) 630 zrat_p1n_chl1 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p1n) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) 631 zrat_p1p_chl1 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p1p) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) 632 zrat_p1s_chl1 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p1s) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) 633 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p1c) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) * zrat_p1c_chl1 634 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p1n) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) * zrat_p1n_chl1 635 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p1p) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) * zrat_p1p_chl1 636 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p1s) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) * zrat_p1s_chl1 637 637 ENDIF 638 638 IF ( ld_chltot .OR. ld_chlnan ) THEN 639 zrat_p2c_chl2 = tracer_bkg(ji,jj,1,jp_fabm_ p2c) / tracer_bkg(ji,jj,1,jp_fabm_chl2)640 zrat_p2n_chl2 = tracer_bkg(ji,jj,1,jp_fabm_ p2n) / tracer_bkg(ji,jj,1,jp_fabm_chl2)641 zrat_p2p_chl2 = tracer_bkg(ji,jj,1,jp_fabm_ p2p) / tracer_bkg(ji,jj,1,jp_fabm_chl2)642 phyto2d_balinc(ji,jj,1,jp_fabm_ p2c) = phyto2d_balinc(ji,jj,1,jp_fabm_chl2) * zrat_p2c_chl2643 phyto2d_balinc(ji,jj,1,jp_fabm_ p2n) = phyto2d_balinc(ji,jj,1,jp_fabm_chl2) * zrat_p2n_chl2644 phyto2d_balinc(ji,jj,1,jp_fabm_ p2p) = phyto2d_balinc(ji,jj,1,jp_fabm_chl2) * zrat_p2p_chl2639 zrat_p2c_chl2 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p2c) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) 640 zrat_p2n_chl2 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p2n) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) 641 zrat_p2p_chl2 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p2p) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) 642 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p2c) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) * zrat_p2c_chl2 643 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p2n) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) * zrat_p2n_chl2 644 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p2p) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) * zrat_p2p_chl2 645 645 ENDIF 646 646 IF ( ld_chltot .OR. ld_chlpic ) THEN 647 zrat_p3c_chl3 = tracer_bkg(ji,jj,1,jp_fabm_ p3c) / tracer_bkg(ji,jj,1,jp_fabm_chl3)648 zrat_p3n_chl3 = tracer_bkg(ji,jj,1,jp_fabm_ p3n) / tracer_bkg(ji,jj,1,jp_fabm_chl3)649 zrat_p3p_chl3 = tracer_bkg(ji,jj,1,jp_fabm_ p3p) / tracer_bkg(ji,jj,1,jp_fabm_chl3)650 phyto2d_balinc(ji,jj,1,jp_fabm_ p3c) = phyto2d_balinc(ji,jj,1,jp_fabm_chl3) * zrat_p3c_chl3651 phyto2d_balinc(ji,jj,1,jp_fabm_ p3n) = phyto2d_balinc(ji,jj,1,jp_fabm_chl3) * zrat_p3n_chl3652 phyto2d_balinc(ji,jj,1,jp_fabm_ p3p) = phyto2d_balinc(ji,jj,1,jp_fabm_chl3) * zrat_p3p_chl3647 zrat_p3c_chl3 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p3c) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) 648 zrat_p3n_chl3 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p3n) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) 649 zrat_p3p_chl3 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p3p) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) 650 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p3c) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) * zrat_p3c_chl3 651 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p3n) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) * zrat_p3n_chl3 652 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p3p) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) * zrat_p3p_chl3 653 653 ENDIF 654 654 IF ( ld_chltot .OR. ld_chldin ) THEN 655 zrat_p4c_chl4 = tracer_bkg(ji,jj,1,jp_fabm_ p4c) / tracer_bkg(ji,jj,1,jp_fabm_chl4)656 zrat_p4n_chl4 = tracer_bkg(ji,jj,1,jp_fabm_ p4n) / tracer_bkg(ji,jj,1,jp_fabm_chl4)657 zrat_p4p_chl4 = tracer_bkg(ji,jj,1,jp_fabm_ p4p) / tracer_bkg(ji,jj,1,jp_fabm_chl4)658 phyto2d_balinc(ji,jj,1,jp_fabm_ p4c) = phyto2d_balinc(ji,jj,1,jp_fabm_chl4) * zrat_p4c_chl4659 phyto2d_balinc(ji,jj,1,jp_fabm_ p4n) = phyto2d_balinc(ji,jj,1,jp_fabm_chl4) * zrat_p4n_chl4660 phyto2d_balinc(ji,jj,1,jp_fabm_ p4p) = phyto2d_balinc(ji,jj,1,jp_fabm_chl4) * zrat_p4p_chl4655 zrat_p4c_chl4 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p4c) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 656 zrat_p4n_chl4 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p4n) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 657 zrat_p4p_chl4 = tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_p4p) / tracer_bkg(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 658 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p4c) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) * zrat_p4c_chl4 659 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p4n) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) * zrat_p4n_chl4 660 phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p4p) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) * zrat_p4p_chl4 661 661 ENDIF 662 662 ENDIF … … 678 678 ! 679 679 DO jk = 2, jkmax 680 phyto2d_balinc(ji,jj,jk,jp_fabm_ chl1) = phyto2d_balinc(ji,jj,1,jp_fabm_chl1)681 phyto2d_balinc(ji,jj,jk,jp_fabm_ p1c) = phyto2d_balinc(ji,jj,1,jp_fabm_p1c)682 phyto2d_balinc(ji,jj,jk,jp_fabm_ p1n) = phyto2d_balinc(ji,jj,1,jp_fabm_p1n)683 phyto2d_balinc(ji,jj,jk,jp_fabm_ p1p) = phyto2d_balinc(ji,jj,1,jp_fabm_p1p)684 phyto2d_balinc(ji,jj,jk,jp_fabm_ p1s) = phyto2d_balinc(ji,jj,1,jp_fabm_p1s)685 phyto2d_balinc(ji,jj,jk,jp_fabm_ chl2) = phyto2d_balinc(ji,jj,1,jp_fabm_chl2)686 phyto2d_balinc(ji,jj,jk,jp_fabm_ p2c) = phyto2d_balinc(ji,jj,1,jp_fabm_p2c)687 phyto2d_balinc(ji,jj,jk,jp_fabm_ p2n) = phyto2d_balinc(ji,jj,1,jp_fabm_p2n)688 phyto2d_balinc(ji,jj,jk,jp_fabm_ p2p) = phyto2d_balinc(ji,jj,1,jp_fabm_p2p)689 phyto2d_balinc(ji,jj,jk,jp_fabm_ chl3) = phyto2d_balinc(ji,jj,1,jp_fabm_chl3)690 phyto2d_balinc(ji,jj,jk,jp_fabm_ p3c) = phyto2d_balinc(ji,jj,1,jp_fabm_p3c)691 phyto2d_balinc(ji,jj,jk,jp_fabm_ p3n) = phyto2d_balinc(ji,jj,1,jp_fabm_p3n)692 phyto2d_balinc(ji,jj,jk,jp_fabm_ p3p) = phyto2d_balinc(ji,jj,1,jp_fabm_p3p)693 phyto2d_balinc(ji,jj,jk,jp_fabm_ chl4) = phyto2d_balinc(ji,jj,1,jp_fabm_chl4)694 phyto2d_balinc(ji,jj,jk,jp_fabm_ p4c) = phyto2d_balinc(ji,jj,1,jp_fabm_p4c)695 phyto2d_balinc(ji,jj,jk,jp_fabm_ p4n) = phyto2d_balinc(ji,jj,1,jp_fabm_p4n)696 phyto2d_balinc(ji,jj,jk,jp_fabm_ p4p) = phyto2d_balinc(ji,jj,1,jp_fabm_p4p)680 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) 681 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p1c) 682 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p1n) 683 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p1p) 684 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p1s) 685 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) 686 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p2c) 687 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p2n) 688 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p2p) 689 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) 690 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p3c) 691 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p3n) 692 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p3p) 693 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 694 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p4c) 695 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p4n) 696 phyto2d_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) = phyto2d_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_p4p) 697 697 END DO 698 698 ! -
branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90
r10622 r10667 188 188 DO jk = 1, jkmax 189 189 PHYT_AVG(ji,jj) = PHYT_AVG(ji,jj) + & 190 & trn(ji,jj,jk,jp_fabm_ p1n) + &191 & trn(ji,jj,jk,jp_fabm_ p2n) + &192 & trn(ji,jj,jk,jp_fabm_ p3n) + &193 & trn(ji,jj,jk,jp_fabm_ p4n)190 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) + & 191 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) + & 192 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) + & 193 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) 194 194 IF ( pgrow_3d(ji,jj,jk) .GT. 0.0 ) THEN 195 195 PGROW_AVG(ji,jj) = PGROW_AVG(ji,jj) + &
Note: See TracChangeset
for help on using the changeset viewer.
