Changeset 10667 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90
- Timestamp:
- 2019-02-13T11:46:09+01:00 (5 years ago)
- File:
-
- 1 edited
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
Note: See TracChangeset
for help on using the changeset viewer.