New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10141 for branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc_v3/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90 – NEMO

Ignore:
Timestamp:
2018-09-19T17:14:46+02:00 (6 years ago)
Author:
dford
Message:

Change cpp key key_foam_medusa to a logical switch ln_foam_medusa.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc_v3/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90

    r10077 r10141  
    1010   !! 'key_hadocc'          : HadOCC model 
    1111   !! 'key_medusa'          : MEDUSA model 
    12    !! 'key_foam_medusa'     : MEDUSA extras for FOAM OBS and ASM 
    1312   !! 'key_roam'            : MEDUSA extras for carbonate cycle 
    1413   !! 'key_karaml'          : Kara mixed layer depth 
     
    5655      & jptra 
    5756#endif 
    58 #if defined key_medusa && defined key_foam_medusa 
     57#if defined key_medusa 
    5958   USE asmphyto2dbal_medusa, ONLY: & ! phyto2d balancing for MEDUSA 
    6059      & asm_phyto2d_bal_medusa 
     
    161160   REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE :: pph_bkginc      ! pph inc 
    162161   REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE :: po2_bkginc      ! po2 inc 
    163 #if defined key_top 
    164162   REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: phyto2d_balinc  ! Balancing incs from ocean colour 
    165163   REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: phyto3d_balinc  ! Balancing incs from p(l)chltot 
    166164   REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: pco2_balinc     ! Balancing incs from spco2/sfco2 
    167165   REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: ph_balinc       ! Balancing incs from pph 
    168 #endif 
    169  
    170 #if defined key_hadocc || (defined key_medusa && defined key_foam_medusa) 
     166 
    171167   REAL(wp), DIMENSION(:,:),     ALLOCATABLE :: pgrow_avg_bkg   ! Background phyto growth 
    172168   REAL(wp), DIMENSION(:,:),     ALLOCATABLE :: ploss_avg_bkg   ! Background phyto loss 
     
    174170   REAL(wp), DIMENSION(:,:),     ALLOCATABLE :: mld_max_bkg     ! Background max MLD 
    175171   REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: tracer_bkg      ! Background tracer state 
    176 #endif 
    177 #if defined key_hadocc 
    178172   REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE :: chl_bkg         ! Background chl 
    179173   REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE :: cchl_p_bkg      ! Background c:chl 
    180 #endif 
    181174 
    182175#  include "domzgr_substitute.h90" 
     
    197190      !!------------------------------------------------------------------------ 
    198191 
    199 #if ! defined key_top || ( ! defined key_hadocc && ( ! defined key_medusa || ! defined key_foam_medusa ) ) 
     192#if ! defined key_top || ( ! defined key_hadocc && ! defined key_medusa ) 
    200193      CALL ctl_stop( ' Attempting to assimilate biogeochemical observations', & 
    201194         &           ' but no compatible biogeochemical model is available' ) 
     
    206199         & ln_slphynoninc .OR. ln_psi4inc .OR. ln_pphinc .OR. ln_po2inc ) THEN 
    207200         CALL ctl_stop( ' Cannot assimilate PFTs, Si4, pH or O2 into HadOCC' ) 
     201      ENDIF 
     202#endif 
     203 
     204#if defined key_medusa 
     205      IF ( .NOT. ln_foam_medusa ) THEN 
     206         CALL ctl_stop( ' MEDUSA data assimilation options not turned on: set ln_foam_medusa' ) 
    208207      ENDIF 
    209208#endif 
     
    517516      !!------------------------------------------------------------------------ 
    518517 
    519 #if defined key_top && ( defined key_hadocc || (defined key_medusa && defined key_foam_medusa) ) 
     518#if defined key_top && ( defined key_hadocc || defined key_medusa ) 
    520519      IF ( ln_slchltotinc .OR. ln_slchldiainc .OR. ln_slchlnoninc .OR. & 
    521520         & ln_schltotinc  .OR. ln_slphytotinc .OR. ln_slphydiainc .OR. & 
     
    552551         chl_bkg(:,:,:)    = chl_bkg(:,:,:)    * tmask(:,:,:) 
    553552         cchl_p_bkg(:,:,:) = cchl_p_bkg(:,:,:) * tmask(:,:,:) 
    554 #elif defined key_medusa && defined key_foam_medusa 
     553#elif defined key_medusa 
    555554         CALL iom_get( inum, jpdom_autoglo, 'medusa_chn', tracer_bkg(:,:,:,jpchn) ) 
    556555         CALL iom_get( inum, jpdom_autoglo, 'medusa_chd', tracer_bkg(:,:,:,jpchd) ) 
     
    578577            CALL iom_get( inum, jpdom_autoglo, 'hadocc_dic', tracer_bkg(:,:,:,jp_had_dic) ) 
    579578            CALL iom_get( inum, jpdom_autoglo, 'hadocc_alk', tracer_bkg(:,:,:,jp_had_alk) ) 
    580 #elif defined key_medusa && defined key_foam_medusa 
     579#elif defined key_medusa 
    581580            CALL iom_get( inum, jpdom_autoglo, 'medusa_zmi', tracer_bkg(:,:,:,jpzmi) ) 
    582581            CALL iom_get( inum, jpdom_autoglo, 'medusa_zme', tracer_bkg(:,:,:,jpzme) ) 
     
    594593            CALL iom_get( inum, jpdom_autoglo, 'hadocc_dic', tracer_bkg(:,:,:,jp_had_dic) ) 
    595594            CALL iom_get( inum, jpdom_autoglo, 'hadocc_alk', tracer_bkg(:,:,:,jp_had_alk) ) 
    596 #elif defined key_medusa && defined key_foam_medusa 
     595#elif defined key_medusa 
    597596            CALL iom_get( inum, jpdom_autoglo, 'medusa_dic', tracer_bkg(:,:,:,jpdic) ) 
    598597            CALL iom_get( inum, jpdom_autoglo, 'medusa_alk', tracer_bkg(:,:,:,jpalk) ) 
     
    620619         CALL iom_get( inum, jpdom_autoglo, 'hadocc_dic', tracer_bkg(:,:,:,jp_had_dic) ) 
    621620         CALL iom_get( inum, jpdom_autoglo, 'hadocc_alk', tracer_bkg(:,:,:,jp_had_alk) ) 
    622 #elif defined key_medusa && defined key_foam_medusa 
     621#elif defined key_medusa 
    623622         CALL iom_get( inum, jpdom_autoglo, 'medusa_dic', tracer_bkg(:,:,:,jpdic) ) 
    624623         CALL iom_get( inum, jpdom_autoglo, 'medusa_alk', tracer_bkg(:,:,:,jpalk) ) 
     
    690689            & ln_schltotinc  .OR. ln_slphytotinc .OR. ln_slphydiainc .OR. & 
    691690            & ln_slphynoninc ) THEN 
    692 #if defined key_medusa && defined key_foam_medusa 
     691#if defined key_medusa 
    693692            CALL iom_rstput( kt, kt, inum, 'phy2d_chn', phyto2d_balinc(:,:,:,jpchn) ) 
    694693            CALL iom_rstput( kt, kt, inum, 'phy2d_chd', phyto2d_balinc(:,:,:,jpchd) ) 
     
    721720 
    722721         IF ( ln_plchltotinc .OR. ln_pchltotinc ) THEN 
    723 #if defined key_medusa && defined key_foam_medusa 
     722#if defined key_medusa 
    724723            CALL iom_rstput( kt, kt, inum, 'phy3d_chn', phyto3d_balinc(:,:,:,jpchn) ) 
    725724            CALL iom_rstput( kt, kt, inum, 'phy3d_chd', phyto3d_balinc(:,:,:,jpchd) ) 
     
    733732 
    734733         IF ( ln_spco2inc ) THEN 
    735 #if defined key_medusa && defined key_foam_medusa 
     734#if defined key_medusa 
    736735            CALL iom_rstput( kt, kt, inum, 'pco2_dic', pco2_balinc(:,:,:,jpdic) ) 
    737736            CALL iom_rstput( kt, kt, inum, 'pco2_alk', pco2_balinc(:,:,:,jpalk) ) 
     
    741740#endif 
    742741         ELSE IF ( ln_sfco2inc ) THEN 
    743 #if defined key_medusa && defined key_foam_medusa 
     742#if defined key_medusa 
    744743            CALL iom_rstput( kt, kt, inum, 'fco2_dic', pco2_balinc(:,:,:,jpdic) ) 
    745744            CALL iom_rstput( kt, kt, inum, 'fco2_alk', pco2_balinc(:,:,:,jpalk) ) 
     
    751750 
    752751         IF ( ln_pphinc ) THEN 
    753 #if defined key_medusa && defined key_foam_medusa 
     752#if defined key_medusa 
    754753            CALL iom_rstput( kt, kt, inum, 'ph_dic', ph_balinc(:,:,:,jpdic) ) 
    755754            CALL iom_rstput( kt, kt, inum, 'ph_alk', ph_balinc(:,:,:,jpalk) ) 
     
    806805      CALL iom_rstput( kt, nitbkg_r, knum, 'hadocc_chl'  , HADOCC_CHL(:,:,:)     ) 
    807806      CALL iom_rstput( kt, nitbkg_r, knum, 'hadocc_cchl' , cchl_p(:,:,:)         ) 
    808 #elif defined key_medusa && defined key_foam_medusa 
     807#elif defined key_medusa 
    809808      CALL iom_rstput( kt, nitbkg_r, knum, 'pgrow_avg'   , pgrow_avg        ) 
    810809      CALL iom_rstput( kt, nitbkg_r, knum, 'ploss_avg'   , ploss_avg        ) 
     
    904903      REAL(wp), DIMENSION(jpi,jpj) :: zinc_phytot ! Local phytot incs 
    905904      REAL(wp), DIMENSION(jpi,jpj) :: zbkg_phytot ! Local phytot bkg 
    906 #if defined key_medusa && defined key_foam_medusa 
     905#if defined key_medusa 
    907906      REAL(wp), DIMENSION(jpi,jpj) :: zinc_chldia ! Local chldia incs 
    908907      REAL(wp), DIMENSION(jpi,jpj) :: zbkg_chldia ! Local chldia bkg 
     
    924923         ! Total chlorophyll 
    925924         IF ( ln_slchltotinc ) THEN 
    926 #if defined key_medusa && defined key_foam_medusa 
     925#if defined key_medusa 
    927926            zbkg_chltot(:,:) = tracer_bkg(:,:,1,jpchn) + tracer_bkg(:,:,1,jpchd) 
    928927#elif defined key_hadocc 
     
    936935         ENDIF 
    937936 
    938 #if defined key_medusa && defined key_foam_medusa 
     937#if defined key_medusa 
    939938         ! Diatom chlorophyll 
    940939         IF ( ln_slchldiainc ) THEN 
     
    946945#endif 
    947946 
    948 #if defined key_medusa && defined key_foam_medusa 
     947#if defined key_medusa 
    949948         ! Non-diatom chlorophyll 
    950949         IF ( ln_slchlnoninc ) THEN 
     
    958957         ! Total phytoplankton carbon 
    959958         IF ( ln_slphytotinc ) THEN 
    960 #if defined key_medusa && defined key_foam_medusa 
     959#if defined key_medusa 
    961960            zbkg_phytot(:,:) = (trn(:,:,1,jpphn) * xthetapn) + (trn(:,:,1,jpphd) * xthetapd) 
    962961#elif defined key_hadocc 
     
    968967         ENDIF 
    969968 
    970 #if defined key_medusa && defined key_foam_medusa 
     969#if defined key_medusa 
    971970         ! Diatom phytoplankton carbon 
    972971         IF ( ln_slphydiainc ) THEN 
     
    978977#endif 
    979978 
    980 #if defined key_medusa && defined key_foam_medusa 
     979#if defined key_medusa 
    981980         ! Non-diatom phytoplankton carbon 
    982981         IF ( ln_slphynoninc ) THEN 
     
    990989         ! Select mixed layer 
    991990         IF ( ll_asmdin ) THEN 
    992 #if defined key_top && ( defined key_hadocc || (defined key_medusa && defined key_foam_medusa) ) 
     991#if defined key_top && ( defined key_hadocc || defined key_medusa ) 
    993992            CALL ctl_warn( ' Doing direct initialisation with ocean colour assimilation', & 
    994993               &           ' Mixed layer depth taken to be background maximum mld_max_bkg' ) 
     
    10261025         zincper = (nitiaufin_r - nitiaustr_r + 1) * rdt 
    10271026 
    1028 #if defined key_medusa && defined key_foam_medusa 
     1027#if defined key_medusa 
    10291028         CALL asm_phyto2d_bal_medusa( (ln_slchltotinc .OR. ln_schltotinc), & 
    10301029            &                         zinc_chltot,                         & 
     
    10841083            ! Add directly to trn and trb, rather than to tra, because tra gets 
    10851084            ! reset to zero at the start of trc_stp, called after this routine 
    1086 #if defined key_medusa && defined key_foam_medusa 
     1085#if defined key_medusa 
    10871086            WHERE( phyto2d_balinc(:,:,:,jp_msa0:jp_msa1) > 0.0_wp .OR. & 
    10881087               &   trn(:,:,:,jp_msa0:jp_msa1) + phyto2d_balinc(:,:,:,jp_msa0:jp_msa1) * zincwgt > 0.0_wp ) 
     
    11201119            CALL ctl_warn( ' Doing direct initialisation with phyto2d assimilation', & 
    11211120               &           ' Background state is taken from model rather than background file' ) 
    1122 #if defined key_medusa && defined key_foam_medusa 
     1121#if defined key_medusa 
    11231122            WHERE( phyto2d_balinc(:,:,:,jp_msa0:jp_msa1) > 0.0_wp .OR. & 
    11241123               &   trn(:,:,:,jp_msa0:jp_msa1) + phyto2d_balinc(:,:,:,jp_msa0:jp_msa1) > 0.0_wp ) 
     
    11871186            ! 3) Subtract background from analysis to get chl incs 
    11881187            ! If rn_maxchlinc > 0 then cap total absolute chlorophyll increment at that value 
    1189 #if defined key_medusa && defined key_foam_medusa 
     1188#if defined key_medusa 
    11901189            bkg_chl(:,:,:) = tracer_bkg(:,:,:,jpchn) + tracer_bkg(:,:,:,jpchd) 
    11911190#elif defined key_hadocc 
     
    12201219         ENDIF 
    12211220 
    1222 #if defined key_medusa && defined key_foam_medusa 
     1221#if defined key_medusa 
    12231222         ! Loop over each grid point partioning the increments based on existing ratios 
    12241223         DO jk = 1, jpk 
     
    12711270            ! Add directly to trn and trb, rather than to tra, because tra gets 
    12721271            ! reset to zero at the start of trc_stp, called after this routine 
    1273 #if defined key_medusa && defined key_foam_medusa 
     1272#if defined key_medusa 
    12741273            WHERE( phyto3d_balinc(:,:,:,jp_msa0:jp_msa1) > 0.0_wp .OR. & 
    12751274               &   trn(:,:,:,jp_msa0:jp_msa1) + phyto3d_balinc(:,:,:,jp_msa0:jp_msa1) * zincwgt > 0.0_wp ) 
     
    13071306            CALL ctl_warn( ' Doing direct initialisation with phyto3d assimilation', & 
    13081307               &           ' Background state is taken from model rather than background file' ) 
    1309 #if defined key_medusa && defined key_foam_medusa 
     1308#if defined key_medusa 
    13101309            WHERE( phyto3d_balinc(:,:,:,jp_msa0:jp_msa1) > 0.0_wp .OR. & 
    13111310               &   trn(:,:,:,jp_msa0:jp_msa1) + phyto3d_balinc(:,:,:,jp_msa0:jp_msa1) > 0.0_wp ) 
     
    13881387 
    13891388         IF ( ln_sfco2inc ) THEN 
    1390 #if defined key_medusa && defined key_foam_medusa && defined key_roam 
     1389#if defined key_medusa && defined key_roam 
    13911390            ! If assimilating fCO2, then convert to pCO2 using MEDUSA MOCSY subroutine 
    13921391            patm(1) = 1.0 
     
    14211420         ENDIF 
    14221421 
    1423 #if defined key_medusa && defined key_foam_medusa 
     1422#if defined key_medusa 
    14241423         ! Account for phytoplankton balancing if required 
    14251424         IF ( ln_phytobal ) THEN 
     
    14561455         ! Select mixed layer 
    14571456         IF ( ll_asmdin ) THEN 
    1458 #if defined key_hadocc || (defined key_medusa && defined key_foam_medusa) 
     1457#if defined key_hadocc || defined key_medusa 
    14591458            CALL ctl_warn( ' Doing direct initialisation with pCO2 assimilation', & 
    14601459               &           ' Mixed layer depth taken to be background maximum mld_max_bkg' ) 
     
    15411540            ! Add directly to trn and trb, rather than to tra, because tra gets 
    15421541            ! reset to zero at the start of trc_stp, called after this routine 
    1543 #if defined key_medusa && defined key_foam_medusa 
     1542#if defined key_medusa 
    15441543            WHERE( pco2_balinc(:,:,:,jp_msa0:jp_msa1) > 0.0_wp .OR. & 
    15451544               &   trn(:,:,:,jp_msa0:jp_msa1) + pco2_balinc(:,:,:,jp_msa0:jp_msa1) * zincwgt > 0.0_wp ) 
     
    15781577            CALL ctl_warn( ' Doing direct initialisation with pCO2 assimilation', & 
    15791578               &           ' Background state is taken from model rather than background file' ) 
    1580 #if defined key_medusa && defined key_foam_medusa 
     1579#if defined key_medusa 
    15811580            WHERE( pco2_balinc(:,:,:,jp_msa0:jp_msa1) > 0.0_wp .OR. & 
    15821581               &   trn(:,:,:,jp_msa0:jp_msa1) + pco2_balinc(:,:,:,jp_msa0:jp_msa1) > 0.0_wp ) 
     
    16491648      !!------------------------------------------------------------------------ 
    16501649 
    1651 #if ! (defined key_medusa && defined key_foam_medusa) 
     1650#if ! defined key_medusa 
    16521651      CALL ctl_stop( ' pH balancing only implemented for MEDUSA' ) 
    16531652#else 
     
    18581857 
    18591858         IF ( ln_pno3inc ) THEN 
    1860 #if defined key_hadocc || ( defined key_medusa && defined key_foam_medusa ) 
     1859#if defined key_hadocc || defined key_medusa 
    18611860#if defined key_hadocc 
    18621861            it = jp_had_nut 
    1863 #elif defined key_medusa && defined key_foam_medusa 
     1862#elif defined key_medusa 
    18641863            it = jpdin 
    18651864#endif 
     
    18821881 
    18831882         IF ( ln_psi4inc ) THEN 
    1884 #if defined key_medusa && defined key_foam_medusa 
     1883#if defined key_medusa 
    18851884            it = jpsil 
    18861885            IF ( ln_phytobal ) THEN 
     
    19021901 
    19031902         IF ( ln_pdicinc ) THEN 
    1904 #if defined key_hadocc || ( defined key_medusa && defined key_foam_medusa ) 
     1903#if defined key_hadocc || defined key_medusa 
    19051904#if defined key_hadocc 
    19061905            it = jp_had_dic 
    1907 #elif defined key_medusa && defined key_foam_medusa 
     1906#elif defined key_medusa 
    19081907            it = jpdic 
    19091908#endif 
     
    19261925 
    19271926         IF ( ln_palkinc ) THEN 
    1928 #if defined key_hadocc || ( defined key_medusa && defined key_foam_medusa ) 
     1927#if defined key_hadocc || defined key_medusa 
    19291928#if defined key_hadocc 
    19301929            it = jp_had_alk 
    1931 #elif defined key_medusa && defined key_foam_medusa 
     1930#elif defined key_medusa 
    19321931            it = jpalk 
    19331932#endif 
     
    19501949 
    19511950         IF ( ln_po2inc ) THEN 
    1952 #if defined key_medusa && defined key_foam_medusa 
     1951#if defined key_medusa 
    19531952            it = jpoxy 
    19541953            IF ( ln_phytobal ) THEN 
     
    20022001                  trb(:,:,:,jp_had_nut) = trb(:,:,:,jp_had_nut) + pno3_bkginc(:,:,:) * zincwgt 
    20032002               END WHERE 
    2004 #elif defined key_medusa && defined key_foam_medusa 
     2003#elif defined key_medusa 
    20052004               WHERE( pno3_bkginc(:,:,:) > 0.0_wp .OR. & 
    20062005                  &   trn(:,:,:,jpdin) + pno3_bkginc(:,:,:) * zincwgt > 0.0_wp ) 
     
    20142013 
    20152014            IF ( ln_psi4inc ) THEN 
    2016 #if defined key_medusa && defined key_foam_medusa 
     2015#if defined key_medusa 
    20172016               WHERE( psi4_bkginc(:,:,:) > 0.0_wp .OR. & 
    20182017                  &   trn(:,:,:,jpsil) + psi4_bkginc(:,:,:) * zincwgt > 0.0_wp ) 
     
    20322031                  trb(:,:,:,jp_had_dic) = trb(:,:,:,jp_had_dic) + pdic_bkginc(:,:,:) * zincwgt 
    20332032               END WHERE 
    2034 #elif defined key_medusa && defined key_foam_medusa 
     2033#elif defined key_medusa 
    20352034               WHERE( pdic_bkginc(:,:,:) > 0.0_wp .OR. & 
    20362035                  &   trn(:,:,:,jpdic) + pdic_bkginc(:,:,:) * zincwgt > 0.0_wp ) 
     
    20502049                  trb(:,:,:,jp_had_alk) = trb(:,:,:,jp_had_alk) + palk_bkginc(:,:,:) * zincwgt 
    20512050               END WHERE 
    2052 #elif defined key_medusa && defined key_foam_medusa 
     2051#elif defined key_medusa 
    20532052               WHERE( palk_bkginc(:,:,:) > 0.0_wp .OR. & 
    20542053                  &   trn(:,:,:,jpalk) + palk_bkginc(:,:,:) * zincwgt > 0.0_wp ) 
     
    20622061 
    20632062            IF ( ln_po2inc ) THEN 
    2064 #if defined key_medusa && defined key_foam_medusa 
     2063#if defined key_medusa 
    20652064               WHERE( po2_bkginc(:,:,:) > 0.0_wp .OR. & 
    20662065                  &   trn(:,:,:,jpoxy) + po2_bkginc(:,:,:) * zincwgt > 0.0_wp ) 
     
    20982097            CALL ctl_warn( ' Doing direct initialisation of HadOCC with 3D BGC assimilation', & 
    20992098               &           ' Background state is taken from model rather than background file' ) 
    2100 #elif defined key_medusa && defined key_foam_medusa 
     2099#elif defined key_medusa 
    21012100            CALL ctl_warn( ' Doing direct initialisation of MEDUSA with 3D BGC assimilation', & 
    21022101               &           ' Background state is taken from model rather than background file' ) 
     
    21102109                  trb(:,:,:,jp_had_nut) = trn(:,:,:,jp_had_nut) 
    21112110               END WHERE 
    2112 #elif defined key_medusa && defined key_foam_medusa 
     2111#elif defined key_medusa 
    21132112               WHERE( pno3_bkginc(:,:,:) > 0.0_wp .OR. & 
    21142113                  &   trn(:,:,:,jpdin) + pno3_bkginc(:,:,:) > 0.0_wp ) 
     
    21222121 
    21232122            IF ( ln_psi4inc ) THEN 
    2124 #if defined key_medusa && defined key_foam_medusa 
     2123#if defined key_medusa 
    21252124               WHERE( psi4_bkginc(:,:,:) > 0.0_wp .OR. & 
    21262125                  &   trn(:,:,:,jpsil) + psi4_bkginc(:,:,:) > 0.0_wp ) 
     
    21402139                  trb(:,:,:,jp_had_dic) = trn(:,:,:,jp_had_dic) 
    21412140               END WHERE 
    2142 #elif defined key_medusa && defined key_foam_medusa 
     2141#elif defined key_medusa 
    21432142               WHERE( pdic_bkginc(:,:,:) > 0.0_wp .OR. & 
    21442143                  &   trn(:,:,:,jpdic) + pdic_bkginc(:,:,:) > 0.0_wp ) 
     
    21582157                  trb(:,:,:,jp_had_alk) = trn(:,:,:,jp_had_alk) 
    21592158               END WHERE 
    2160 #elif defined key_medusa && defined key_foam_medusa 
     2159#elif defined key_medusa 
    21612160               WHERE( palk_bkginc(:,:,:) > 0.0_wp .OR. & 
    21622161                  &   trn(:,:,:,jpalk) + palk_bkginc(:,:,:) > 0.0_wp ) 
     
    21702169 
    21712170            IF ( ln_po2inc ) THEN 
    2172 #if defined key_medusa && defined key_foam_medusa 
     2171#if defined key_medusa 
    21732172               WHERE( po2_bkginc(:,:,:) > 0.0_wp .OR. & 
    21742173                  &   trn(:,:,:,jpoxy) + po2_bkginc(:,:,:) > 0.0_wp ) 
Note: See TracChangeset for help on using the changeset viewer.