Changeset 14086 for NEMO/trunk/src/TOP/PISCES/P4Z/p4zint.F90
- Timestamp:
- 2020-12-04T12:37:14+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/TOP/PISCES/P4Z/p4zint.F90
r12377 r14086 19 19 REAL(wp) :: xksilim = 16.5e-6_wp ! Half-saturation constant for the Si half-saturation constant computation 20 20 21 # include "do_loop_substitute.h90" 21 22 !!---------------------------------------------------------------------- 22 23 !! NEMO/TOP 4.0 , NEMO Consortium (2018) … … 49 50 ! Computation of the silicon dependant half saturation constant for silica uptake 50 51 ! --------------------------------------------------- 51 DO ji = 1, jpi 52 DO jj = 1, jpj 53 zvar = tr(ji,jj,1,jpsil,Kbb) * tr(ji,jj,1,jpsil,Kbb) 54 xksimax(ji,jj) = MAX( xksimax(ji,jj), ( 1.+ 7.* zvar / ( xksilim * xksilim + zvar ) ) * 1e-6 ) 55 END DO 56 END DO 52 DO_2D( 1, 1, 1, 1 ) 53 zvar = tr(ji,jj,1,jpsil,Kbb) * tr(ji,jj,1,jpsil,Kbb) 54 xksimax(ji,jj) = MAX( xksimax(ji,jj), ( 1.+ 7.* zvar / ( xksilim * xksilim + zvar ) ) * 1e-6 ) 55 END_2D 57 56 ! 58 57 IF( nday_year == nyear_len(1) ) THEN
Note: See TracChangeset
for help on using the changeset viewer.