Changeset 463 for trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90

Timestamp:

2006-05-10T19:26:55+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_054:RB: take into account tracers and dynamics reorganization, change atsk to jki

File:

• trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90

@@ -16 +16 @@
    USE dom_oce         ! ocean space and time domain
    USE zdf_oce         ! ocean vertical physics
-   USE in_out_manager  ! I/O manager
-   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
    USE phycst          ! physical constants
    USE taumod          ! surface stress
@@ -23 +21 @@
    USE ocesbc          ! thermohaline fluxes
    USE zdfddm          ! double diffusion mixing
+   USE in_out_manager  ! I/O manager
+   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
    USE prtctl          ! Print control
 
@@ -30 +30 @@
    !! * Routine accessibility
    PUBLIC zdf_kpp   ! routine called by step.F90
+   PUBLIC tra_kpp   ! routine called by step.F90
 
    !! * Share Module variables
@@ -1262 +1263 @@
 
 
+   SUBROUTINE tra_kpp
+      !!----------------------------------------------------------------------
+      !!                  ***  ROUTINE tra_kpp  ***
+      !!
+      !! ** Purpose :   compute and add to the tracer trend the non-local
+      !!      tracer flux
+      !!
+      !! ** Method  :   ???
+      !!
+      !! history :
+      !!     9.0  ! 05-11 (G. Madec)  Original code
+      !!----------------------------------------------------------------------
+      !! * Modules used
+      USE oce, ONLY :    ztrdt => ua,       & ! use ua as 3D workspace
+                         ztrds => va          ! use va as 3D workspace
+      !!----------------------------------------------------------------------
+
+      IF( kt == nit000 ) THEN
+         IF(lwp) WRITE(numout,*)
+         IF(lwp) WRITE(numout,*) 'tra_kpp : KPP non-local tracer fluxes'
+         IF(lwp) WRITE(numout,*) '~~~~~~~   '
+      ENDIF
+
+
+      ! Save ta and sa trends
+      IF( l_trdtra )   THEN
+         ztrdt(:,:,:) = ta(:,:,:)
+         ztrds(:,:,:) = sa(:,:,:)
+      ENDIF
+
+      ! add non-local temperature and salinity flux ( in convective case only)
+      DO jk = 1, jpkm1
+         DO jj = 2, jpjm1
+            DO ji = fs_2, fs_jpim1
+               ta(ji,jj,jk) =  ta(ji,jj,jk)                           &
+                  &         - ( ghats(ji,jj,jk  ) * avt(ji,jj,jk  )   &
+                  &           - ghats(ji,jj,jk+1) * avt(ji,jj,jk+1) ) * wt0(ji,jj) / fse3t(ji,jj,jk)
+               sa(ji,jj,jk) = sa(ji,jj,jk)                              &
+                  &         - ( ghats(ji,jj,jk  ) * fsavs(ji,jj,jk  )   &
+                  &           - ghats(ji,jj,jk+1) * fsavs(ji,jj,jk+1) ) * ws0(ji,jj) / fse3t(ji,jj,jk)
+            END DO
+         END DO
+      END DO
+
+      ! save the non-local tracer flux trends for diagnostic
+      IF( l_trdtra )   THEN
+         ztrdt(:,:,:) = ta(:,:,:) - ztrdt(:,:,:)
+         ztrds(:,:,:) = sa(:,:,:) - ztrds(:,:,:)
+!!bug gm jpttdzdf ==> jpttkpp
+         CALL trd_mod(ztrdt, ztrds, jpttdzdf, 'TRA', kt)
+      ENDIF
+
+      IF(ln_ctl) THEN         ! print mean trends (used for debugging)
+         CALL prt_ctl(tab3d_1=ta, clinfo1=' kpp  - Ta: ', mask1=tmask, &
+            &         tab3d_2=sa, clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra')
+      ENDIF
+
+   END SUBROUTINE tra_kpp
+
 
    SUBROUTINE zdf_kpp_init
@@ -1547 +1607 @@
       WRITE(*,*) 'zdf_kpp: You should not have seen this print! error?', kt
    END SUBROUTINE zdf_kpp
+   SUBROUTINE tra_kpp( kt )          ! Empty routine
+      WRITE(*,*) 'tra_kpp: You should not have seen this print! error?', kt
+   END SUBROUTINE tra_kpp
 #endif