- Timestamp:
- 2013-03-06T12:00:26+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zmeso.F90
r3295 r3830 74 74 REAL(wp) :: zfact , zstep, zfood, zfoodlim 75 75 REAL(wp) :: zepshert, zepsherv, zgrarsig, zgraztot, zgraztotf 76 REAL(wp) :: zgrarem2, zgrafer2, zgrapoc2, zprcaca, zmortz2, z grasrat76 REAL(wp) :: zgrarem2, zgrafer2, zgrapoc2, zprcaca, zmortz2, zmortzgoc, zgrasrat 77 77 #if defined key_kriest 78 78 REAL znumpoc … … 97 97 zstep = xstep 98 98 # endif 99 zfact = zstep * tgfunc (ji,jj,jk) * zcompam99 zfact = zstep * tgfunc2(ji,jj,jk) * zcompam 100 100 101 101 ! Respiration rates of both zooplankton … … 155 155 zepshert = epsher2 * MIN( 1., zncratio ) 156 156 zepsherv = zepshert * MIN( 1., zgrasrat / ferat3 ) 157 zgrarem2 = zgraztot * ( 1. - zepsherv - unass2 ) 158 zgrafer2 = zgraztot * MAX( 0. , ( 1. - unass2 ) * zgrasrat - ferat3 * zepshert ) 157 zgrarem2 = zgraztot * ( 1. - zepsherv - unass2 ) + zrespz2 & 158 & + ( 1. - zepsherv - unass2 ) /( 1. - zepsherv + rtrn) * ztortz2 159 zgrafer2 = zgraztot * MAX( 0. , ( 1. - unass2 ) * zgrasrat - ferat3 * zepsherv ) & 160 & + ferat3 * ( zrespz2 + ( 1. - zepsherv - unass2 ) /( 1. - zepsherv + rtrn) * ztortz2 ) 159 161 zgrapoc2 = zgraztot * unass2 160 162 … … 177 179 #endif 178 180 zmortz2 = ztortz2 + zrespz2 181 zmortzgoc = unass2 / ( 1. - zepsherv + rtrn ) * ztortz2 179 182 tra(ji,jj,jk,jpmes) = tra(ji,jj,jk,jpmes) - zmortz2 + zepsherv * zgraztot 180 183 tra(ji,jj,jk,jpdia) = tra(ji,jj,jk,jpdia) - zgrazd … … 198 201 #if defined key_kriest 199 202 znumpoc = trn(ji,jj,jk,jpnum) / ( trn(ji,jj,jk,jppoc) + rtrn ) 200 tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zmortz 2- zgrazpoc - zgrazffe203 tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zmortzgoc - zgrazpoc - zgrazffe 201 204 tra(ji,jj,jk,jpnum) = tra(ji,jj,jk,jpnum) - zgrazpoc * znumpoc & 202 & + zmortz 2* xkr_dmeso - zgrazffe * znumpoc * wsbio4(ji,jj,jk) / ( wsbio3(ji,jj,jk) + rtrn )205 & + zmortzgoc * xkr_dmeso - zgrazffe * znumpoc * wsbio4(ji,jj,jk) / ( wsbio3(ji,jj,jk) + rtrn ) 203 206 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + ferat3 * zmortz2 - zgrazfff - zgrazpof 204 207 #else 205 208 tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) - zgrazpoc 206 tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) + zmortz 2- zgrazffe209 tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) + zmortzgoc - zgrazffe 207 210 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) - zgrazpof 208 tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + ferat3 * zmortz 2- zgrazfff211 tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + ferat3 * zmortzgoc - zgrazfff 209 212 #endif 210 213
Note: See TracChangeset
for help on using the changeset viewer.