Changeset 3580 for trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

Timestamp:

2012-11-16T15:49:18+01:00 (10 years ago)

Author:

cetlod

Message:

Correction of surface boundary condition on passive tracers : see ticket #999

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90 (modified) (2 diffs)

trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

@@ -50 +50 @@
       !!            tra = tra + emp * trn / e3t   for k=1
       !!         where emp, the surface freshwater budget (evaporation minus
-      !!         precipitation minus runoff) given in kg/m2/s is divided
+      !!         precipitation ) given in kg/m2/s is divided
       !!         by 1035 kg/m3 (density of ocean water) to obtain m/s.
       !!
@@ -79 +79 @@
       ENDIF
 
+      ! Coupling online : river runoff is added to the horizontal divergence (hdivn) in the subroutine sbc_rnf_div 
+      ! one only consider the concentration/dilution effect due to evaporation minus precipitation + freezing/melting of sea-ice
 
-      IF( lk_offline ) THEN          ! emps in dynamical files contains emps - rnf
-         zemps(:,:) = emps(:,:)  
-      ELSE                           ! Concentration dilution effect on tracer due to evaporation, precipitation, and river runoff
-         IF( lk_vvl ) THEN                      ! volume variable
-            zemps(:,:) = emps(:,:) - emp(:,:)   
-!!ch         zemps(:,:) = 0.
-         ELSE                                   ! linear free surface
-            IF( ln_rnf ) THEN  ;  zemps(:,:) = emps(:,:) - rnf(:,:)   !  E-P-R
-            ELSE               ;  zemps(:,:) = emps(:,:)
-            ENDIF 
-         ENDIF 
+      ! Coupling in offline, hdivn is computed from ocean horizontal velocities only ; the runoff are not included.
+      ! emps in dynamical files contains (emps - rnf)
+      IF( .NOT. lk_offline .AND. lk_vvl ) THEN  ! online coupling + volume variable 
+         zemps(:,:) = emps(:,:) - emp(:,:)   
+      ELSE
+         zemps(:,:) = emps(:,:)
       ENDIF