Changeset 1289
- Timestamp:
- 2009-02-03T15:18:40+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/PISCES/p4zbio.F90
r1264 r1289 59 59 INTEGER, INTENT(in) :: kt, jnt 60 60 INTEGER :: ji, jj, jk, jn 61 REAL(wp) :: ztra 61 62 #if defined key_kriest 62 63 REAL(wp) :: zcoef1, zcoef2 … … 69 70 ! OF PHYTOPLANKTON AND DETRITUS 70 71 71 xdiss(:,:,:) = 0.01 72 72 xdiss(:,:,:) = 1. 73 73 !!gm the use of nmld should be better here? 74 DO jk = 1, jpkm174 DO jk = 2, jpkm1 75 75 DO jj = 1, jpj 76 76 DO ji = 1, jpi 77 IF( fsdepw(ji,jj,jk+1) .le. hmld(ji,jj) ) xdiss(ji,jj,jk) = 1.e077 IF( fsdepw(ji,jj,jk+1) > hmld(ji,jj) ) xdiss(ji,jj,jk) = 0.01 78 78 END DO 79 79 END DO … … 92 92 CALL p4z_meso ( kt, jnt ) ! mesozooplankton 93 93 94 95 ! ! test if tracers concentrations fall below 0. 94 ! ! test if tracers concentrations fall below 0. 96 95 xnegtr(:,:,:) = 1.e0 97 96 DO jn = jp_pcs0, jp_pcs1 … … 99 98 DO jj = 1, jpj 100 99 DO ji = 1, jpi 101 IF( ( trn(ji,jj,jk,jn) + tra(ji,jj,jk,jn) ) < 0.e0 ) xnegtr(ji,jj,jk) = 0.e0 102 END DO 100 IF( ( trn(ji,jj,jk,jn) + tra(ji,jj,jk,jn) ) < 0.e0 ) THEN 101 ztra = ABS( ( trn(ji,jj,jk,jn) - rtrn ) & 102 / ( tra(ji,jj,jk,jn) + rtrn ) ) 103 xnegtr(ji,jj,jk) = MIN( xnegtr(ji,jj,jk), ztra ) 104 ENDIF 105 END DO 103 106 END DO 104 107 END DO 105 108 END DO 106 109 ! ! where at least 1 tracer concentration becomes negative 107 ! ! all tracer tendancy are set to zero (i.e. trn = trb)110 ! ! 108 111 DO jn = jp_pcs0, jp_pcs1 109 112 trn(:,:,:,jn) = trn(:,:,:,jn) + xnegtr(:,:,:) * tra(:,:,:,jn)
Note: See TracChangeset
for help on using the changeset viewer.