Changeset 434 for trunk/NEMO/OPA_SRC/TRA

Timestamp:

2006-04-10T17:46:12+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_044 : CT : update the passive tracers TOP component and the standard GYRE configuration

File:

• trunk/NEMO/OPA_SRC/TRA/tradmp.F90 (modified) (7 diffs)

trunk/NEMO/OPA_SRC/TRA/tradmp.F90

@@ -24 +24 @@
    USE dtatem          ! temperature data
    USE dtasal          ! salinity data
+   USE dtasss          ! surface salinity data
    USE zdfmxl          ! mixed layer depth
    USE lib_mpp         ! distribued memory computing
@@ -47 +48 @@
    !! * Module variables
    INTEGER  ::             & !!! * newtonian damping namelist (mandmp) *
-      ndmp   =   -1 ,      &  ! = 0/-1/'latitude' for damping over T and S
+      ndmp   =   -1 ,      &  ! = 0/-1/-2/'latitude' for damping over T and S
       ndmpf  =    2 ,      &  ! = 1 create a damping.coeff NetCDF file 
-      nmldmp =    0           ! = 0/1/2 flag for damping in the mixed layer
+      nmldmp =    0           ! = 0/1/2/3 flag for damping in the mixed layer
    REAL(wp) ::             & !!!  * newtonian damping namelist *
       sdmp   =   50.,      &  ! surface time scale for internal damping (days)
@@ -177 +178 @@
          END DO
 
+      CASE( 3 )                ! newtonian damping of SSS only
+         DO jj = 2, jpjm1
+            DO ji = fs_2, fs_jpim1   ! vector opt.
+               zsa = resto(ji,jj,1) * ( sss(ji,jj) - sb(ji,jj,1) )
+               ! add the trends to the general tracer trends
+               sa(ji,jj,1) = sa(ji,jj,1) + zsa
+               ! save the salinity trend (used in flx to close the salt budget)
+               strdmp(ji,jj,1) = zsa
+            END DO
+         END DO
+
       END SELECT
 
@@ -237 +249 @@
       SELECT CASE ( ndmp )
 
+      CASE ( -2 )               ! GYRE: damping SSS only
+         IF(lwp) WRITE(numout,*) '          SSS damping in Gyre'
+
       CASE ( -1 )               ! ORCA: damping in Red & Med Seas only
          IF(lwp) WRITE(numout,*) '          tracer damping in the Med & Red seas only'
@@ -261 +276 @@
       CASE ( 2 )                ! no damping in the mixed layer 
          IF(lwp) WRITE(numout,*) '          no tracer damping in the mixed layer'
+
+      CASE ( 3 )                ! SSS damping only 
+         IF(lwp) WRITE(numout,*) '          SSS surface damping'
 
       CASE DEFAULT
@@ -275 +293 @@
          nstop = nstop + 1
       ENDIF
+
+      IF( lk_dtasss .AND. ( ( ndmp /= -2 ) .OR. (nmldmp /= 3) ) ) THEN
+         IF(lwp) WRITE(numout,cform_err)
+         IF(lwp) WRITE(numout,*) '          Check namelist for SSS restoring in Gyre '
+         IF(lwp) WRITE(numout,*) '          ndmp must be -2 and nmldmp must be 3'
+         nstop = nstop + 1
+      ENDIF
+
+
 
       strdmp(:,:,:) = 0.e0       ! internal damping salinity trend (used in ocesbc)
@@ -656 +683 @@
          resto(:,:,jpk) = 0.e0
 
+      ELSEIF( cp_cfg == "gyre" .AND. ( ndmp == -2 )) THEN
+
+         !                                         ! =========================
+         !                                         !  SSS damping
+         !                                         ! =========================
+         zsdmp = 1./(sdmp * rday)
+         IF(lwp)WRITE(numout,*)
+         IF(lwp)WRITE(numout,*) '              GYRE configuration: Damping coefficient'
+         resto(:,:, : ) = 0.
+         resto(:,:, 1 ) = zsdmp * tmask(:,:,1)
+
       ELSE
          !    ------------