Ignore:
Timestamp:
2020-06-11T19:32:37+02:00 (5 months ago)
Author:
dford
Message:

Allow nitrogen balancing for both 2D and 3D chlorophyll increments.

File:
1 moved

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_phytobal3d/NEMOGCM/NEMO/OPA_SRC/ASM/asmphytobal_hadocc.F90

    r13096 r13097  
    1 MODULE asmphyto2dbal_hadocc 
     1MODULE asmphytobal_hadocc 
    22   !!====================================================================== 
    3    !!                       ***  MODULE asmphyto2dbal_hadocc  *** 
    4    !! Calculate increments to HadOCC based on surface phyto2d increments 
     3   !!                       ***  MODULE asmphytobal_hadocc  *** 
     4   !! Calculate increments to HadOCC based on surface phyto increments 
    55   !! 
    66   !! IMPORTANT NOTE: This calls the bioanalysis routine of Hemmings et al. 
     
    1717   !! 'key_hadocc'          : HadOCC model 
    1818   !!---------------------------------------------------------------------- 
    19    !! asm_phyto2d_bal_hadocc : routine to calculate increments to HadOCC 
     19   !! asm_phyto_bal_hadocc : routine to calculate increments to HadOCC 
    2020   !!---------------------------------------------------------------------- 
    2121   USE par_kind,      ONLY: wp             ! kind parameters 
     
    3232   PRIVATE                    
    3333 
    34    PUBLIC asm_phyto2d_bal_hadocc 
     34   PUBLIC asm_phyto_bal_hadocc 
    3535 
    3636   ! Default values for biological assimilation parameters 
     
    6767CONTAINS 
    6868 
    69    SUBROUTINE asm_phyto2d_bal_hadocc( ld_chltot,                      & 
    70       &                               pinc_chltot,                    & 
    71       &                               ld_phytot,                      & 
    72       &                               pinc_phytot,                    & 
    73       &                               pincper,                        & 
    74       &                               p_maxchlinc, ld_phytobal, pmld, & 
    75       &                               pgrow_avg_bkg, ploss_avg_bkg,   & 
    76       &                               phyt_avg_bkg, mld_max_bkg,      & 
    77       &                               cchl_p_bkg,                     & 
    78       &                               tracer_bkg, phyto2d_balinc ) 
     69   SUBROUTINE asm_phyto_bal_hadocc( kdeps,                              & 
     70      &                             ld_chltot,                          & 
     71      &                             pinc_chltot_3d,                     & 
     72      &                             ld_phytot,                          & 
     73      &                             pinc_phytot_3d,                     & 
     74      &                             pincper,                            & 
     75      &                             p_maxchlinc, ld_phytobal, pmld,     & 
     76      &                             pgrow_avg_bkg_3d, ploss_avg_bkg_3d, & 
     77      &                             phyt_avg_bkg_3d, mld_max_bkg,       & 
     78      &                             cchl_p_bkg_3d,                      & 
     79      &                             tracer_bkg, phyto_balinc ) 
    7980      !!--------------------------------------------------------------------------- 
    80       !!                    ***  ROUTINE asm_phyto2d_bal_hadocc  *** 
     81      !!                    ***  ROUTINE asm_phyto_bal_hadocc  *** 
    8182      !! 
    8283      !! ** Purpose :   calculate increments to HadOCC from 2d phytoplankton increments 
     
    8485      !! ** Method  :   call nitrogen balancing scheme 
    8586      !! 
    86       !! ** Action  :   populate phyto2d_balinc 
     87      !! ** Action  :   populate phyto_balinc 
    8788      !! 
    8889      !! References :   Hemmings et al., 2008, J. Mar. Res. 
     
    9091      !!--------------------------------------------------------------------------- 
    9192      !! 
    92       LOGICAL,  INTENT(in   )                               :: ld_chltot      ! Assim chltot y/n 
    93       REAL(wp), INTENT(inout), DIMENSION(jpi,jpj)           :: pinc_chltot    ! chltot increments 
    94       LOGICAL,  INTENT(in   )                               :: ld_phytot      ! Assim phytot y/n 
    95       REAL(wp), INTENT(inout), DIMENSION(jpi,jpj)           :: pinc_phytot    ! phytot increments 
    96       REAL(wp), INTENT(in   )                               :: pincper        ! Assimilation period 
    97       REAL(wp), INTENT(in   )                               :: p_maxchlinc    ! Max chl increment 
    98       LOGICAL,  INTENT(in   )                               :: ld_phytobal    ! Balancing y/n 
    99       REAL(wp), INTENT(inout), DIMENSION(jpi,jpj)           :: pmld           ! Mixed layer depth 
    100       REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: pgrow_avg_bkg  ! Avg phyto growth 
    101       REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: ploss_avg_bkg  ! Avg phyto loss 
    102       REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: phyt_avg_bkg   ! Avg phyto 
    103       REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: mld_max_bkg    ! Max MLD 
    104       REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: cchl_p_bkg     ! Surface C:Chl 
    105       REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,jpk,jptra) :: tracer_bkg     ! State variables 
    106       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk,jptra) :: phyto2d_balinc ! Balancing increments 
    107       !! 
    108       INTEGER                                               :: ji, jj, jk, jn ! Loop counters 
    109       INTEGER                                               :: jkmax          ! Loop index 
    110       INTEGER,                 DIMENSION(6)                 :: i_tracer       ! Tracer indices 
    111       REAL(wp),                DIMENSION(16)                :: modparm        ! Model parameters 
    112       REAL(wp),                DIMENSION(20)                :: assimparm      ! Assimilation parameters 
    113       REAL(wp),                DIMENSION(jpi,jpj,jpk,6)     :: bstate         ! Background state 
    114       REAL(wp),                DIMENSION(jpi,jpj,jpk,6)     :: outincs        ! Balancing increments 
    115       REAL(wp),                DIMENSION(jpi,jpj,22)        :: diag           ! Depth-indep diagnostics 
    116       REAL(wp),                DIMENSION(jpi,jpj,jpk,22)    :: diag_fulldepth ! Full-depth diagnostics 
     93      INTEGER,  INTENT(in   )                               :: kdeps            ! No. inc deps 1 or jpk 
     94      LOGICAL,  INTENT(in   )                               :: ld_chltot        ! Assim chltot y/n 
     95      REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,kdeps)     :: pinc_chltot_3d   ! chltot increments 
     96      LOGICAL,  INTENT(in   )                               :: ld_phytot        ! Assim phytot y/n 
     97      REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,kdeps)     :: pinc_phytot_3d   ! phytot increments 
     98      REAL(wp), INTENT(in   )                               :: pincper          ! Assimilation period 
     99      REAL(wp), INTENT(in   )                               :: p_maxchlinc      ! Max chl increment 
     100      LOGICAL,  INTENT(in   )                               :: ld_phytobal      ! Balancing y/n 
     101      REAL(wp), INTENT(inout), DIMENSION(jpi,jpj)           :: pmld             ! Mixed layer depth 
     102      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,kdeps)     :: pgrow_avg_bkg_3d ! Avg phyto growth 
     103      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,kdeps)     :: ploss_avg_bkg_3d ! Avg phyto loss 
     104      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,kdeps)     :: phyt_avg_bkg_3d  ! Avg phyto 
     105      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: mld_max_bkg      ! Max MLD 
     106      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,kdeps)     :: cchl_p_bkg_3d    ! C:Chl 
     107      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,jpk,jptra) :: tracer_bkg       ! State variables 
     108      REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk,jptra) :: phyto_balinc     ! Balancing increments 
     109      !! 
     110      INTEGER                                               :: ji, jj, jk, jn    ! Loop counters 
     111      INTEGER                                               :: jkmax             ! Loop index 
     112      INTEGER,                 DIMENSION(6)                 :: i_tracer          ! Tracer indices 
     113      REAL(wp),                DIMENSION(16)                :: modparm           ! Model parameters 
     114      REAL(wp),                DIMENSION(20)                :: assimparm         ! Assimilation parameters 
     115      REAL(wp),                DIMENSION(jpi,jpj,1,6)       :: bstate_2d         ! Background state (2D) 
     116      REAL(wp),                DIMENSION(jpi,jpj,jpk,6)     :: bstate_3d         ! Background state (3D) 
     117      REAL(wp),                DIMENSION(jpi,jpj,1,6)       :: outincs_2d        ! Balancing increments (2D) 
     118      REAL(wp),                DIMENSION(jpi,jpj,jpk,6)     :: outincs_3d        ! Balancing increments (3D) 
     119      REAL(wp),                DIMENSION(jpi,jpj,22)        :: diag              ! Depth-indep diagnostics 
     120      REAL(wp),                DIMENSION(jpi,jpj,1,22)      :: diag_fulldepth_2d ! Full-depth diagnostics (2D) 
     121      REAL(wp),                DIMENSION(jpi,jpj,jpk,22)    :: diag_fulldepth_3d ! Full-depth diagnostics (3D) 
     122      REAL(wp),                DIMENSION(jpi,jpj)           :: cchl_p_bkg_2d     ! C:Chl for total phy (2D) 
     123      REAL(wp),                DIMENSION(jpi,jpj,1)         :: tmask_2d          ! Single-level tmask 
     124      REAL(wp),                DIMENSION(jpi,jpj)           :: pinc_chltot_2d    ! chltot increments (2D) 
     125      REAL(wp),                DIMENSION(jpi,jpj)           :: pinc_phytot_2d    ! phytot increments (2D) 
     126      REAL(wp),                DIMENSION(jpi,jpj)           :: pgrow_avg_bkg_2d  ! Avg phyto growth (2D) 
     127      REAL(wp),                DIMENSION(jpi,jpj)           :: ploss_avg_bkg_2d  ! Avg phyto loss (2D) 
     128      REAL(wp),                DIMENSION(jpi,jpj)           :: phyt_avg_bkg_2d   ! Avg phyto (2D) 
    117129      !!--------------------------------------------------------------------------- 
    118130 
    119131      IF ( ( .NOT. ld_chltot ) .AND. ( .NOT. ld_phytot ) ) THEN 
    120          CALL ctl_stop( ' Trying to do phyto2d balancing but nothing to assimilate' ) 
     132         CALL ctl_stop( ' Trying to do phyto balancing but nothing to assimilate' ) 
    121133      ENDIF 
    122134       
     
    124136      IF ( p_maxchlinc > 0.0 ) THEN 
    125137         IF ( ld_chltot ) THEN 
    126             DO jj = 1, jpj 
    127                DO ji = 1, jpi 
    128                   pinc_chltot(ji,jj) = MAX( -1.0 * p_maxchlinc, MIN( pinc_chltot(ji,jj), p_maxchlinc ) ) 
     138            DO jk = 1, kdeps 
     139               DO jj = 1, jpj 
     140                  DO ji = 1, jpi 
     141                     pinc_chltot_3d(ji,jj,jk) = MAX( -1.0 * p_maxchlinc, MIN( pinc_chltot_3d(ji,jj,jk), p_maxchlinc ) ) 
     142                  END DO 
    129143               END DO 
    130144            END DO 
     
    187201 
    188202         ! Set background state 
    189          bstate(:,:,:,i_tracer(1)) = tracer_bkg(:,:,:,jp_had_nut) 
    190          bstate(:,:,:,i_tracer(2)) = tracer_bkg(:,:,:,jp_had_phy) 
    191          bstate(:,:,:,i_tracer(3)) = tracer_bkg(:,:,:,jp_had_zoo) 
    192          bstate(:,:,:,i_tracer(4)) = tracer_bkg(:,:,:,jp_had_pdn) 
    193          bstate(:,:,:,i_tracer(5)) = tracer_bkg(:,:,:,jp_had_dic) 
    194          bstate(:,:,:,i_tracer(6)) = tracer_bkg(:,:,:,jp_had_alk) 
     203         bstate_3d(:,:,:,i_tracer(1)) = tracer_bkg(:,:,:,jp_had_nut) 
     204         bstate_3d(:,:,:,i_tracer(2)) = tracer_bkg(:,:,:,jp_had_phy) 
     205         bstate_3d(:,:,:,i_tracer(3)) = tracer_bkg(:,:,:,jp_had_zoo) 
     206         bstate_3d(:,:,:,i_tracer(4)) = tracer_bkg(:,:,:,jp_had_pdn) 
     207         bstate_3d(:,:,:,i_tracer(5)) = tracer_bkg(:,:,:,jp_had_dic) 
     208         bstate_3d(:,:,:,i_tracer(6)) = tracer_bkg(:,:,:,jp_had_alk) 
    195209 
    196210         ! Call nitrogen balancing routine 
    197          CALL bio_analysis( jpi, jpj, jpk, ZDZ(:,:,:), i_tracer, modparm,               & 
    198             &               n2be_p, n2be_z, n2be_d, assimparm,                          & 
    199             &               INT(pincper), 1, kmt(:,:), tmask(:,:,:),                    & 
    200             &               pmld(:,:), mld_max_bkg(:,:), pinc_chltot(:,:), cchl_p_bkg(:,:), & 
    201             &               nbal_active, phyt_avg_bkg(:,:),                             & 
    202             &               gl_active, pgrow_avg_bkg(:,:), ploss_avg_bkg(:,:),          & 
    203             &               subsurf_active, deepneg_active,                             & 
    204             &               deeppos_active, nutprof_active,                             & 
    205             &               bstate, outincs,                                            & 
    206             &               diag_active, diag,                                          & 
    207             &               diag_fulldepth_active, diag_fulldepth ) 
     211         IF (kdeps == 1) THEN 
     212            pinc_chltot_2d(:,:)   = pinc_chltot_3d(:,:,1) 
     213            cchl_p_bkg_2d(:,:)    = cchl_p_bkg_3d(:,:,1) 
     214            phyt_avg_bkg_2d(:,:)  = phyt_avg_bkg_3d(:,:,1) 
     215            pgrow_avg_bkg_2d(:,:) = pgrow_avg_bkg_3d(:,:,1) 
     216            ploss_avg_bkg_2d(:,:) = ploss_avg_bkg_3d(:,:,1) 
     217             
     218            CALL bio_analysis( jpi, jpj, jpk, ZDZ(:,:,:), i_tracer, modparm,               & 
     219               &               n2be_p, n2be_z, n2be_d, assimparm,                          & 
     220               &               INT(pincper), 1, kmt(:,:), tmask(:,:,:),                    & 
     221               &               pmld(:,:), mld_max_bkg(:,:), pinc_chltot_2d(:,:), cchl_p_bkg_2d(:,:), & 
     222               &               nbal_active, phyt_avg_bkg_2d(:,:),                          & 
     223               &               gl_active, pgrow_avg_bkg_2d(:,:), ploss_avg_bkg_2d(:,:),    & 
     224               &               subsurf_active, deepneg_active,                             & 
     225               &               deeppos_active, nutprof_active,                             & 
     226               &               bstate_3d, outincs_3d,                                      & 
     227               &               diag_active, diag,                                          & 
     228               &               diag_fulldepth_active, diag_fulldepth_3d ) 
     229         ELSE 
     230            pmld(:,:) = 0.5 
     231             
     232            DO jk = 1, kdeps 
     233               pinc_chltot_2d(:,:)   = pinc_chltot_3d(:,:,jk) 
     234               cchl_p_bkg_2d(:,:)    = cchl_p_bkg_3d(:,:,jk) 
     235               phyt_avg_bkg_2d(:,:)  = phyt_avg_bkg_3d(:,:,jk) 
     236               pgrow_avg_bkg_2d(:,:) = pgrow_avg_bkg_3d(:,:,jk) 
     237               ploss_avg_bkg_2d(:,:) = ploss_avg_bkg_3d(:,:,jk) 
     238               tmask_2d(:,:,1)       = tmask(:,:,jk) 
     239               bstate_2d(:,:,1,:)    = bstate_3d(:,:,jk,:) 
     240               outincs_2d(:,:,:,:)   = 0.0 
     241 
     242               CALL bio_analysis( jpi, jpj, 1, gdepw_n(:,:,2), i_tracer, modparm,            & 
     243                  &               n2be_p, n2be_z, n2be_d, assimparm,                         & 
     244                  &               INT(pincper), 1, INT(SUM(tmask_2d,3)), tmask_2d(:,:,:),    & 
     245                  &               pmld(:,:), pmld(:,:), pinc_chltot_2d(:,:), cchl_p_bkg_2d(:,:), & 
     246                  &               nbal_active, phyt_avg_bkg_2d(:,:),                         & 
     247                  &               gl_active, pgrow_avg_bkg_2d(:,:), ploss_avg_bkg_2d(:,:),   & 
     248                  &               subsurf_active, deepneg_active,                            & 
     249                  &               deeppos_active, nutprof_active,                            & 
     250                  &               bstate_2d, outincs_2d,                                     & 
     251                  &               diag_active, diag,                                         & 
     252                  &               diag_fulldepth_active, diag_fulldepth_2d ) 
     253 
     254               outincs_3d(:,:,jk,:) = outincs_2d(:,:,1,:) 
     255            END DO 
     256         ENDIF 
    208257 
    209258         ! Save balancing increments 
    210          phyto2d_balinc(:,:,:,jp_had_nut) = outincs(:,:,:,i_tracer(1)) 
    211          phyto2d_balinc(:,:,:,jp_had_phy) = outincs(:,:,:,i_tracer(2)) 
    212          phyto2d_balinc(:,:,:,jp_had_zoo) = outincs(:,:,:,i_tracer(3)) 
    213          phyto2d_balinc(:,:,:,jp_had_pdn) = outincs(:,:,:,i_tracer(4)) 
    214          phyto2d_balinc(:,:,:,jp_had_dic) = outincs(:,:,:,i_tracer(5)) 
    215          phyto2d_balinc(:,:,:,jp_had_alk) = outincs(:,:,:,i_tracer(6)) 
     259         phyto_balinc(:,:,:,jp_had_nut) = outincs_3d(:,:,:,i_tracer(1)) 
     260         phyto_balinc(:,:,:,jp_had_phy) = outincs_3d(:,:,:,i_tracer(2)) 
     261         phyto_balinc(:,:,:,jp_had_zoo) = outincs_3d(:,:,:,i_tracer(3)) 
     262         phyto_balinc(:,:,:,jp_had_pdn) = outincs_3d(:,:,:,i_tracer(4)) 
     263         phyto_balinc(:,:,:,jp_had_dic) = outincs_3d(:,:,:,i_tracer(5)) 
     264         phyto_balinc(:,:,:,jp_had_alk) = outincs_3d(:,:,:,i_tracer(6)) 
    216265       
    217266      ELSE   ! No nitrogen balancing 
    218267       
    219268         ! Initialise phytoplankton increment to zero 
    220          phyto2d_balinc(:,:,:,jp_had_phy) = 0.0 
     269         phyto_balinc(:,:,:,jp_had_phy) = 0.0 
    221270          
    222271         ! Convert surface chlorophyll increment to phytoplankton nitrogen 
    223          phyto2d_balinc(:,:,1,jp_had_phy) = ( cchl_p_bkg(:,:) / (mw_carbon * c2n_p) ) * pinc_chltot(:,:) 
     272         DO jk = 1, kdeps 
     273            phyto_balinc(:,:,jk,jp_had_phy) = ( cchl_p_bkg_3d(:,:,jk) / (mw_carbon * c2n_p) ) * pinc_chltot_3d(:,:,jk) 
     274         END DO 
    224275          
    225          ! Propagate through mixed layer 
    226          DO jj = 1, jpj 
    227             DO ji = 1, jpi 
    228                ! 
    229                jkmax = jpk-1 
    230                DO jk = jpk-1, 1, -1 
    231                   IF ( ( pmld(ji,jj) >  gdepw_n(ji,jj,jk)   ) .AND. & 
    232                      & ( pmld(ji,jj) <= gdepw_n(ji,jj,jk+1) ) ) THEN 
    233                      pmld(ji,jj) = gdepw_n(ji,jj,jk+1) 
    234                      jkmax = jk 
    235                   ENDIF 
     276         IF (kdeps == 1) THEN 
     277            ! Propagate through mixed layer 
     278            DO jj = 1, jpj 
     279               DO ji = 1, jpi 
     280                  ! 
     281                  jkmax = jpk-1 
     282                  DO jk = jpk-1, 1, -1 
     283                     IF ( ( pmld(ji,jj) >  gdepw_n(ji,jj,jk)   ) .AND. & 
     284                        & ( pmld(ji,jj) <= gdepw_n(ji,jj,jk+1) ) ) THEN 
     285                        pmld(ji,jj) = gdepw_n(ji,jj,jk+1) 
     286                        jkmax = jk 
     287                     ENDIF 
     288                  END DO 
     289                  ! 
     290                  DO jk = 2, jkmax 
     291                     phyto_balinc(ji,jj,jk,jp_had_phy) = phyto_balinc(ji,jj,1,jp_had_phy) 
     292                  END DO 
     293                  ! 
    236294               END DO 
    237                ! 
    238                DO jk = 2, jkmax 
    239                   phyto2d_balinc(ji,jj,jk,jp_had_phy) = phyto2d_balinc(ji,jj,1,jp_had_phy) 
    240                END DO 
    241                ! 
    242295            END DO 
    243          END DO 
     296         ENDIF 
    244297 
    245298         ! Set other balancing increments to zero 
    246          phyto2d_balinc(:,:,:,jp_had_nut) = 0.0 
    247          phyto2d_balinc(:,:,:,jp_had_zoo) = 0.0 
    248          phyto2d_balinc(:,:,:,jp_had_pdn) = 0.0 
    249          phyto2d_balinc(:,:,:,jp_had_dic) = 0.0 
    250          phyto2d_balinc(:,:,:,jp_had_alk) = 0.0 
     299         phyto_balinc(:,:,:,jp_had_nut) = 0.0 
     300         phyto_balinc(:,:,:,jp_had_zoo) = 0.0 
     301         phyto_balinc(:,:,:,jp_had_pdn) = 0.0 
     302         phyto_balinc(:,:,:,jp_had_dic) = 0.0 
     303         phyto_balinc(:,:,:,jp_had_alk) = 0.0 
    251304       
    252305      ENDIF 
    253306 
    254    END SUBROUTINE asm_phyto2d_bal_hadocc 
     307   END SUBROUTINE asm_phyto_bal_hadocc 
    255308 
    256309#else 
     
    259312   !!---------------------------------------------------------------------- 
    260313CONTAINS 
    261    SUBROUTINE asm_phyto2d_bal_hadocc( ld_chltot,                      & 
    262       &                               pinc_chltot,                    & 
    263       &                               ld_phytot,                      & 
    264       &                               pinc_phytot,                    & 
    265       &                               pincper,                        & 
    266       &                               p_maxchlinc, ld_phytobal, pmld, & 
    267       &                               pgrow_avg_bkg, ploss_avg_bkg,   & 
    268       &                               phyt_avg_bkg, mld_max_bkg,      & 
    269       &                               cchl_p_bkg,                     & 
    270       &                               tracer_bkg, phyto2d_balinc ) 
     314   SUBROUTINE asm_phyto_bal_hadocc( kdeps,                              & 
     315      &                             ld_chltot,                          & 
     316      &                             pinc_chltot_3d,                     & 
     317      &                             ld_phytot,                          & 
     318      &                             pinc_phytot,                        & 
     319      &                             pincper,                            & 
     320      &                             p_maxchlinc, ld_phytobal, pmld,     & 
     321      &                             pgrow_avg_bkg_3d, ploss_avg_bkg_3d, & 
     322      &                             phyt_avg_bkg_3d, mld_max_bkg,       & 
     323      &                             cchl_p_bkg_3d,                      & 
     324      &                             tracer_bkg, phyto_balinc ) 
     325      INTEGER :: kdeps 
    271326      LOGICAL :: ld_chltot 
    272       REAL    :: pinc_chltot(:,:) 
     327      REAL    :: pinc_chltot_3d(:,:,:) 
    273328      LOGICAL :: ld_phytot 
    274       REAL    :: pinc_phytot(:,:) 
     329      REAL    :: pinc_phytot_3d(:,:,:) 
    275330      REAL    :: pincper 
    276331      REAL    :: p_maxchlinc 
    277332      LOGICAL :: ld_phytobal 
    278333      REAL    :: pmld(:,:) 
    279       REAL    :: pgrow_avg_bkg(:,:) 
    280       REAL    :: ploss_avg_bkg(:,:) 
    281       REAL    :: phyt_avg_bkg(:,:) 
     334      REAL    :: pgrow_avg_bkg_3d(:,:,:) 
     335      REAL    :: ploss_avg_bkg_3d(:,:,:) 
     336      REAL    :: phyt_avg_bkg_3d(:,:,:) 
    282337      REAL    :: mld_max_bkg(:,:) 
    283       REAL    :: cchl_p_bkg(:,:) 
     338      REAL    :: cchl_p_bkg_3d(:,:,:) 
    284339      REAL    :: tracer_bkg(:,:,:,:) 
    285       REAL    :: phyto2d_balinc(:,:,:,:) 
    286       WRITE(*,*) 'asm_phyto2d_bal_hadocc: You should not have seen this print! error?' 
    287    END SUBROUTINE asm_phyto2d_bal_hadocc 
     340      REAL    :: phyto_balinc(:,:,:,:) 
     341      WRITE(*,*) 'asm_phyto_bal_hadocc: You should not have seen this print! error?' 
     342   END SUBROUTINE asm_phyto_bal_hadocc 
    288343#endif 
    289344 
    290345   !!====================================================================== 
    291 END MODULE asmphyto2dbal_hadocc 
     346END MODULE asmphytobal_hadocc 
Note: See TracChangeset for help on using the changeset viewer.