Changeset 2162

Timestamp:

2010-10-05T21:21:11+02:00 (14 years ago)

Author:

cetlod

Message:

bug correction on PISCES

File:

• branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zrem.F90 (modified) (1 diff)

branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zrem.F90

@@ -273 +273 @@
                !     Remineralization rate of BSi depedant on T and saturation
                !     ---------------------------------------------------------
-               zsatur   = ( sio3eq(ji,jj,jk) - trn(ji,jj,jk,jpsil) ) / ( sio3eq(ji,jj,jk) + rtrn )
-               zsatur   = MAX( rtrn, zsatur )
-               ztem1    = ( 1. + tsn(ji,jj,jk,jp_tem) / 15.)
-               ztem2    = ( 1. + tsn(ji,jj,jk,jp_tem) / 400.)
-               zsatur1  = zsatur * ztem1
-               zsatur2  = zsatur * ztem2 * ztem2 * ztem2 * ztem2
-               zsatur22 = zsatur2 * zsatur2
-               znusil   = 0.225  * zsatur1 + 0.775 * zsatur22 * zsatur22 * zsatur22 * zsatur22 * zsatur2
+               zsatur  = ( sio3eq(ji,jj,jk) - trn(ji,jj,jk,jpsil) ) / ( sio3eq(ji,jj,jk) + rtrn )
+               zsatur  = MAX( rtrn, zsatur )
+               zsatur2 = zsatur * ( 1. + tn(ji,jj,jk) / 400.)**4
+               znusil  = 0.225  * ( 1. + tn(ji,jj,jk) / 15.) * zsatur + 0.775 * zsatur2**9
                zsiremin = xsirem * zstep * znusil
-               zosil    = zsiremin * trn(ji,jj,jk,jpdsi)
+               zosil = zsiremin * trn(ji,jj,jk,jpdsi)
 
                tra(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) - zosil