Changeset 491 for trunk/NEMO/TOP_SRC

Timestamp:

2006-09-01T15:54:05+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_bugfix_054:CE+RB: correct numerical scheme in CFC forcing

File:

• trunk/NEMO/TOP_SRC/SMS/trcfreons.F90 (modified) (4 diffs)

trunk/NEMO/TOP_SRC/SMS/trcfreons.F90

@@ -162 +162 @@
       ENDDO
 
-
+  
       !-------------------------------
       ! Computation of speed transfert
@@ -183 +183 @@
       ENDDO
 
+
       !----------------------------------------------------------------
       ! Input function  : speed *( conc. at equil - concen at surface )
@@ -191 +192 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               qtr(ji,jj,jn) = -zak_cfc(ji,jj,jn) * ( trn(ji,jj,1,jn) - zca_cfc(ji,jj,jn) )   &
+               qtr(ji,jj,jn) = -zak_cfc(ji,jj,jn) * ( trb(ji,jj,1,jn) - zca_cfc(ji,jj,jn) )   &
+#if defined key_off_degrad
+                  &                               * facvol(ji,jj,1)                           &
+#endif
                   &                               * tmask(ji,jj,1) * ( 1. - freeze(ji,jj) )
             END DO
          END DO
       ENDDO
+
 
       ! ---------------------
@@ -278 +283 @@
       ENDDO
 
+      DO jn = 1, jptra
+         WRITE(numout,*) 'coefficient for solubility of tracer',ctrcnm(jn)
+         WRITE(numout,*) soa1(jn), soa2(jn),soa3(jn), soa4(jn), &
+            &           sob1(jn), sob2(jn),sob3(jn)
+         WRITE(numout,*) 
+         WRITE(numout,*) 'coefficient for schmidt of tracer',ctrcnm(jn)
+         WRITE(numout,*) sca1(jn), sca2(jn),sca3(jn), sca4(jn)
+      ENDDO
+
    END SUBROUTINE trc_freons_cst
 #else