- Timestamp:
- 2016-11-03T16:56:39+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zprod.F90
r7180 r7184 46 46 REAL(wp), PUBLIC :: chlcmin !: 47 47 REAL(wp), PUBLIC :: grosip !: 48 REAL(wp), PUBLIC :: zlimmxln !:49 REAL(wp), PUBLIC :: zlimmxlp !:50 REAL(wp), PUBLIC :: zlimmxld !:51 48 52 49 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: prmaxn !: optimal production = f(temperature) … … 161 158 IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN 162 159 zval = MAX( 1., zstrn(ji,jj) ) 163 zmxl_fac(ji,jj,jk) = zval 160 IF( fsdept(ji,jj,jk) <= hmld(ji,jj) ) THEN 161 zval = MIN(1., heup_01(ji,jj) / ( hmld(ji,jj) + rtrn )) 162 ENDIF 164 163 zmxl_chl(ji,jj,jk) = zval / 24. 165 IF( fsdepw(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 166 zval = MIN(1., heup_01(ji,jj) / ( hmld(ji,jj) + rtrn )) 167 zmxl_fac(ji,jj,jk) = zmxl_fac(ji,jj,jk) * zval 168 zmxl_chl(ji,jj,jk) = zmxl_chl(ji,jj,jk) * zval 169 ENDIF 170 zmxl_fac(ji,jj,jk) = ( 1. - exp( -0.2 * zmxl_fac(ji,jj,jk) ) ) 171 zmxl_chl(ji,jj,jk) = zmxl_chl(ji,jj,jk) 164 zmxl_fac(ji,jj,jk) = 1.5 * zval / ( 12. + zval ) 172 165 ENDIF 173 166 END DO … … 183 176 zdaylen(:,:) = MAX(1., zstrn(:,:)) / 24. 184 177 WHERE( zstrn(:,:) < 1.e0 ) zstrn(:,:) = 24. 185 zstrn(:,:) = 24. / zstrn(:,:)186 178 187 179 DO jk = 1, jpkm1 … … 202 194 & /( trb(ji,jj,jk,jpdia) * 12. + rtrn) 203 195 ! 204 zpislopen = zpislopeadn(ji,jj,jk) / ( prmaxn(ji,jj,jk) * rday * xlimphy(ji,jj,jk) + rtrn )205 zpislopep = zpislopeadp(ji,jj,jk) / ( prmaxp(ji,jj,jk) * rday * xlimpic(ji,jj,jk) + rtrn )206 zpisloped = zpislopeadd(ji,jj,jk) / ( prmaxd(ji,jj,jk) * rday * xlimdia(ji,jj,jk) + rtrn )196 zpislopen = zpislopeadn(ji,jj,jk) / ( zprbio(ji,jj,jk) * rday * xlimphy(ji,jj,jk) + rtrn ) 197 zpislopep = zpislopeadp(ji,jj,jk) / ( zprpic(ji,jj,jk) * rday * xlimpic(ji,jj,jk) + rtrn ) 198 zpisloped = zpislopeadd(ji,jj,jk) / ( zprdia(ji,jj,jk) * rday * xlimdia(ji,jj,jk) + rtrn ) 207 199 208 200 ! Computation of production function for Carbon … … 556 548 ! 557 549 NAMELIST/nampisprod/ pislopen, pislopep, pisloped, xadap, bresp, excretn, excretp, excretd, & 558 & thetannm, thetanpm, thetandm, chlcmin, grosip, zlimmxln, & 559 & zlimmxlp, zlimmxld 550 & thetannm, thetanpm, thetandm, chlcmin, grosip & 560 551 561 552 INTEGER :: ios ! Local integer output status for namelist read … … 588 579 WRITE(numout,*) ' Minimum Chl/N in picophytoplankton thetanpm =', thetanpm 589 580 WRITE(numout,*) ' Minimum Chl/N in diatoms thetandm =', thetandm 590 WRITE(numout,*) ' Critical time scale for mixing (nano) zlimmxln =', zlimmxln591 WRITE(numout,*) ' Critical time scale for mixing (pico) zlimmxlp =', zlimmxlp592 WRITE(numout,*) ' Critical time scale for mixing (diatoms) zlimmxld =', zlimmxld593 581 ENDIF 594 582 !
Note: See TracChangeset
for help on using the changeset viewer.