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 10302 for branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/TOP_SRC/MEDUSA/phytoplankton.F90 – NEMO

Ignore:
Timestamp:
2018-11-13T18:21:16+01:00 (5 years ago)
Author:
dford
Message:

Merge in revisions 8447:10159 of dev_r5518_GO6_package.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/TOP_SRC/MEDUSA/phytoplankton.F90

    r8441 r10302  
    66   !! History : 
    77   !!   -   ! 2017-04 (M. Stringer)        Code taken from trcbio_medusa.F90 
     8   !!   -   ! 2017-08 (A. Yool)            Mean mixed layer chlorophyll 
    89   !!---------------------------------------------------------------------- 
    910#if defined key_medusa 
     
    4243                                   zchd, zchn, zdet, zdin, zdtc,         & 
    4344                                   zfer, zpds, zphd, zphn, zsil,         & 
    44                                    zzme, zzmi 
    45       USE dom_oce,           ONLY: e3t_0, e3t_n, gdepw_0, gdepw_n, tmask 
     45                                   zzme, zzmi, fchl_ml 
     46      USE dom_oce,           ONLY: e3t_0, gdepw_0, tmask 
     47#if defined key_vvl 
     48      USE dom_oce,           ONLY: e3t_n, gdepw_n 
     49#endif 
    4650      USE in_out_manager,    ONLY: lwp, numout 
    4751      USE oce,               ONLY: tsn 
     
    5559                                   xvpd, xvpn, xxi 
    5660      USE zdfmxl,            ONLY: hmld 
     61      USE lbclnk,            ONLY: lbc_lnk 
    5762 
    5863   !!* Substitution 
     
    373378               fprd_ml(ji,jj) = fprd_ml(ji,jj) + (fprd(ji,jj) * zphd(ji,jj) * & 
    374379                                                  fse3t(ji,jj,jk) * fq0) 
    375             ENDIF 
    376          ENDDO 
    377       ENDDO 
     380          !! AXY (16/08/17) 
     381          fchl_ml(ji,jj) = fchl_ml(ji,jj) + ((zchn(ji,jj) + zchd(ji,jj)) * & 
     382                                             (fse3t(ji,jj,jk) * fq0) / hmld(ji,jj)) 
     383            ENDIF 
     384         ENDDO 
     385      ENDDO 
     386      CALL lbc_lnk(fchl_ml(:,:),'T',1. ) 
    378387 
    379388      DO jj = 2,jpjm1 
Note: See TracChangeset for help on using the changeset viewer.