Changeset 9292 for branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/plankton.F90
- Timestamp:
- 2018-01-30T19:41:58+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/plankton.F90
r8441 r9292 36 36 fdpn, fdpn2, fdzme, fdzme2, & 37 37 fdzmi, fdzmi2, fsdiss, fsin, & 38 # if defined key_foam_medusa 39 fdep1, fprn, fprd, & 40 fgmepd, fgmepn, fgmipn, & 41 # endif 38 42 zphd, zphn, zpds, zzme, zzmi 43 # if defined key_foam_medusa 44 USE dom_oce, ONLY: e3t_0, e3t_n, gdepw_0, gdepw_n, tmask 45 USE par_kind, ONLY: wp 46 # else 39 47 USE dom_oce, ONLY: tmask 48 # endif 40 49 USE par_oce, ONLY: jpim1, jpjm1 41 50 USE phytoplankton_mod, ONLY: phytoplankton 42 51 USE sms_medusa, ONLY: jmpd, jmpn, jmzme, jmzmi, & 52 # if defined key_foam_medusa 53 pgrow_avg, ploss_avg, phyt_avg, & 54 # endif 43 55 xkphd, xkphn, xkzme, xkzmi, & 44 56 xmetapd, xmetapn, xmetazme, xmetazmi, & 45 57 xmpd, xmpn, xmzme, xmzmi, xsdiss 58 # if defined key_foam_medusa 59 USE zdfmxl, ONLY: hmld 60 # endif 46 61 USE zooplankton_mod, ONLY: zooplankton 62 63 # if defined key_foam_medusa 64 !!* Substitution 65 # include "domzgr_substitute.h90" 66 # endif 47 67 48 68 !! Level … … 50 70 51 71 INTEGER :: ji, jj 72 73 # if defined key_foam_medusa 74 REAL(wp) :: fq0 75 # endif 52 76 53 77 !!------------------------------------------------------------------- … … 188 212 ENDDO 189 213 214 # if defined key_foam_medusa 215 !! Mixed layer averages for ocean colour assimilation 216 !! 217 DO jj = 2,jpjm1 218 DO ji = 2,jpim1 219 IF (tmask(ji,jj,jk) == 1) THEN 220 if (fdep1(ji,jj).le.hmld(ji,jj)) then 221 !! this level is entirely in the mixed layer 222 fq0 = 1.0 223 elseif (fsdepw(ji,jj,jk).ge.hmld(ji,jj)) then 224 !! this level is entirely below the mixed layer 225 fq0 = 0.0 226 else 227 !! this level straddles the mixed layer 228 fq0 = (hmld(ji,jj) - fsdepw(ji,jj,jk)) / fse3t(ji,jj,jk) 229 endif 230 !! 231 pgrow_avg(ji,jj) = pgrow_avg(ji,jj) + & 232 (((fprn(ji,jj) * zphn(ji,jj)) + & 233 (fprd(ji,jj) * zphd(ji,jj)) ) * & 234 fse3t(ji,jj,jk) * fq0) 235 ploss_avg(ji,jj) = ploss_avg(ji,jj) + & 236 ((fgmepd(ji,jj) + fdpd(ji,jj) + & 237 fdpd2(ji,jj) + & 238 fgmepn(ji,jj) + fdpn(ji,jj) + & 239 fdpn2(ji,jj) + fgmipn(ji,jj) ) * & 240 fse3t(ji,jj,jk) * fq0) 241 phyt_avg(ji,jj) = phyt_avg(ji,jj) + & 242 ((zphn(ji,jj) + zphd(ji,jj)) * & 243 fse3t(ji,jj,jk) * fq0) 244 ENDIF 245 ENDDO 246 ENDDO 247 # endif 248 190 249 END SUBROUTINE plankton 191 250
Note: See TracChangeset
for help on using the changeset viewer.