Changeset 2162 for branches/DEV_r2006_merge_TRA_TRC/NEMO
- Timestamp:
- 2010-10-05T21:21:11+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zrem.F90
r2124 r2162 273 273 ! Remineralization rate of BSi depedant on T and saturation 274 274 ! --------------------------------------------------------- 275 zsatur = ( sio3eq(ji,jj,jk) - trn(ji,jj,jk,jpsil) ) / ( sio3eq(ji,jj,jk) + rtrn ) 276 zsatur = MAX( rtrn, zsatur ) 277 ztem1 = ( 1. + tsn(ji,jj,jk,jp_tem) / 15.) 278 ztem2 = ( 1. + tsn(ji,jj,jk,jp_tem) / 400.) 279 zsatur1 = zsatur * ztem1 280 zsatur2 = zsatur * ztem2 * ztem2 * ztem2 * ztem2 281 zsatur22 = zsatur2 * zsatur2 282 znusil = 0.225 * zsatur1 + 0.775 * zsatur22 * zsatur22 * zsatur22 * zsatur22 * zsatur2 275 zsatur = ( sio3eq(ji,jj,jk) - trn(ji,jj,jk,jpsil) ) / ( sio3eq(ji,jj,jk) + rtrn ) 276 zsatur = MAX( rtrn, zsatur ) 277 zsatur2 = zsatur * ( 1. + tn(ji,jj,jk) / 400.)**4 278 znusil = 0.225 * ( 1. + tn(ji,jj,jk) / 15.) * zsatur + 0.775 * zsatur2**9 283 279 zsiremin = xsirem * zstep * znusil 284 zosil 280 zosil = zsiremin * trn(ji,jj,jk,jpdsi) 285 281 286 282 tra(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) - zosil
Note: See TracChangeset
for help on using the changeset viewer.