Ignore:
Timestamp:
2017-08-14T15:22:09+02:00 (3 years ago)
Author:
dford
Message:

Implement initial version of surface chlorophyll assimilation for MEDUSA.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_hadocc.F90

    r8428 r8436  
    2525   USE iom                                 ! i/o 
    2626   USE trc,           ONLY: trn, trb,    & ! HadOCC variables 
    27       &                     HADOCC_CHL,  & 
    28       &                     pgrow_avg,   & 
    29       &                     ploss_avg,   & 
    30       &                     phyt_avg,    & 
    31       &                     mld_max 
     27      &                     HADOCC_CHL 
    3228   USE par_hadocc                          ! HadOCC parameters 
    3329   USE had_bgc_stnd,  ONLY: kmt            ! HadOCC parameters 
     
    7571 
    7672   SUBROUTINE asm_logchl_bal_hadocc( logchl_bkginc, aincper, mld_choice_bgc, & 
    77       &                              k_maxchlinc, ld_logchlbal, logchl_balinc ) 
     73      &                              k_maxchlinc, ld_logchlbal,              & 
     74      &                              pgrow_avg_bkg, ploss_avg_bkg,           & 
     75      &                              phyt_avg_bkg, mld_max_bkg,              & 
     76      &                              logchl_balinc ) 
    7877      !!--------------------------------------------------------------------------- 
    7978      !!                    ***  ROUTINE asm_logchl_bal_hadocc  *** 
     
    9594      REAL(wp), INTENT(in   )                               :: k_maxchlinc    ! Max chl increment 
    9695      LOGICAL,  INTENT(in   )                               :: ld_logchlbal   ! Balancing y/n 
     96      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: pgrow_avg_bkg  ! Avg phyto growth 
     97      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: ploss_avg_bkg  ! Avg phyto loss 
     98      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: phyt_avg_bkg   ! Avg phyto 
     99      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj)           :: mld_max_bkg    ! Max MLD 
    97100      REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk,jptra) :: logchl_balinc  ! Balancing increments 
    98101      !! 
    99       INTEGER                                               :: ji, jj, jk    ! Loop counters 
     102      INTEGER                                               :: ji, jj, jk, jn ! Loop counters 
    100103      INTEGER                                               :: jkmax          ! Loop index 
    101104      INTEGER,                 DIMENSION(6)                 :: i_tracer       ! Tracer indices 
     
    215218 
    216219         ! Call nitrogen balancing routine 
    217          CALL bio_analysis( jpi, jpj, jpk, ZDZ(:,:,:), i_tracer, modparm,          & 
    218             &               n2be_p, n2be_z, n2be_d, assimparm,                     & 
    219             &               INT(aincper), 1, kmt(:,:), tmask(:,:,:),               & 
    220             &               zmld(:,:), mld_max(:,:), chl_inc(:,:), cchl_p(:,:,1), & 
    221             &               nbal_active, phyt_avg(:,:),                            & 
    222             &               gl_active, pgrow_avg(:,:), ploss_avg(:,:),             & 
    223             &               subsurf_active, deepneg_active,                        & 
    224             &               deeppos_active, nutprof_active,                        & 
    225             &               bstate, outincs,                                       & 
    226             &               diag_active, diag,                                     & 
     220         CALL bio_analysis( jpi, jpj, jpk, ZDZ(:,:,:), i_tracer, modparm,             & 
     221            &               n2be_p, n2be_z, n2be_d, assimparm,                        & 
     222            &               INT(aincper), 1, kmt(:,:), tmask(:,:,:),                  & 
     223            &               zmld(:,:), mld_max_bkg(:,:), chl_inc(:,:), cchl_p(:,:,1), & 
     224            &               nbal_active, phyt_avg_bkg(:,:),                           & 
     225            &               gl_active, pgrow_avg_bkg(:,:), ploss_avg_bkg(:,:),        & 
     226            &               subsurf_active, deepneg_active,                           & 
     227            &               deeppos_active, nutprof_active,                           & 
     228            &               bstate, outincs,                                          & 
     229            &               diag_active, diag,                                        & 
    227230            &               diag_fulldepth_active, diag_fulldepth ) 
    228231 
Note: See TracChangeset for help on using the changeset viewer.