Changeset 2789 for branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

Timestamp:

2011-06-27T13:18:25+02:00 (13 years ago)

Author:

cetlod

Message:

Implementation of the merge of TRA/TRP : first guess, see ticket #842

File:

• branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90 (modified) (3 diffs)

branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

@@ -131 +131 @@
       !!              coefficients using the GLS turbulent closure scheme.
       !!----------------------------------------------------------------------
-      USE oce,     z_elem_a  =>   ua   ! use ua as workspace
-      USE oce,     z_elem_b  =>   va   ! use va as workspace
-      USE oce,     z_elem_c  =>   ta   ! use ta as workspace
-      USE oce,     psi       =>   sa   ! use sa as workspace
+      USE oce     , ONLY z_elem_a  =>   ua   ! use ua as workspace
+      USE oce     , ONLY z_elem_b  =>   va   ! use va as workspace
+      USE oce     , ONLY tsa                 ! use tsa as workspace
       USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released
       USE wrk_nemo, ONLY: zdep  => wrk_2d_1
@@ -152 +151 @@
       REAL(wp) ::   prod, buoy, diss, zdiss, sm         !   -      -
       REAL(wp) ::   gh, gm, shr, dif, zsqen, zav        !   -      -
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: z_elem_c, psi
       !!--------------------------------------------------------------------
 
@@ -157 +157 @@
          CALL ctl_stop('zdf_gls: requested workspace arrays unavailable.')   ;   RETURN
       END IF
+      !
+      z_elem_c  => tsa(:,:,:,1)
+      psi       => tsa(:,:,:,2)
 
       ! Preliminary computing