Ignore:
Timestamp:
2012-11-16T15:49:18+01:00 (8 years ago)
Author:
cetlod
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

    r3294 r3580  
    5050      !!            tra = tra + emp * trn / e3t   for k=1 
    5151      !!         where emp, the surface freshwater budget (evaporation minus 
    52       !!         precipitation minus runoff) given in kg/m2/s is divided 
     52      !!         precipitation ) given in kg/m2/s is divided 
    5353      !!         by 1035 kg/m3 (density of ocean water) to obtain m/s. 
    5454      !! 
     
    7979      ENDIF 
    8080 
     81      ! Coupling online : river runoff is added to the horizontal divergence (hdivn) in the subroutine sbc_rnf_div  
     82      ! one only consider the concentration/dilution effect due to evaporation minus precipitation + freezing/melting of sea-ice 
    8183 
    82       IF( lk_offline ) THEN          ! emps in dynamical files contains emps - rnf 
    83          zemps(:,:) = emps(:,:)   
    84       ELSE                           ! Concentration dilution effect on tracer due to evaporation, precipitation, and river runoff 
    85          IF( lk_vvl ) THEN                      ! volume variable 
    86             zemps(:,:) = emps(:,:) - emp(:,:)    
    87 !!ch         zemps(:,:) = 0. 
    88          ELSE                                   ! linear free surface 
    89             IF( ln_rnf ) THEN  ;  zemps(:,:) = emps(:,:) - rnf(:,:)   !  E-P-R 
    90             ELSE               ;  zemps(:,:) = emps(:,:) 
    91             ENDIF  
    92          ENDIF  
     84      ! Coupling in offline, hdivn is computed from ocean horizontal velocities only ; the runoff are not included. 
     85      ! emps in dynamical files contains (emps - rnf) 
     86      IF( .NOT. lk_offline .AND. lk_vvl ) THEN  ! online coupling + volume variable  
     87         zemps(:,:) = emps(:,:) - emp(:,:)    
     88      ELSE 
     89         zemps(:,:) = emps(:,:) 
    9390      ENDIF  
    9491 
Note: See TracChangeset for help on using the changeset viewer.