Changeset 7367 for branches/UKMO/CO5_package_branch/NEMOGCM/NEMO/OPA_SRC/LDF/ldfeiv.F90

Timestamp:

2016-11-29T11:52:31+01:00 (8 years ago)

Author:

deazer

Message:

Contains merged code for CO5 reference.

File:

• branches/UKMO/CO5_package_branch/NEMOGCM/NEMO/OPA_SRC/LDF/ldfeiv.F90 (modified) (3 diffs)

branches/UKMO/CO5_package_branch/NEMOGCM/NEMO/OPA_SRC/LDF/ldfeiv.F90

@@ -60 +60 @@
       REAL(wp) ::   zfw, ze3w, zn2, zf20, zaht, zaht_min      ! temporary scalars
       REAL(wp), DIMENSION(:,:), POINTER ::   zn, zah, zhw, zross   ! 2D workspace
+      REAL(wp), DIMENSION(:,:,:), POINTER ::   uslp, vslp, wslpi, wslpj ! 3D workspace
       !!----------------------------------------------------------------------
       !
@@ -65 +66 @@
       !
       CALL wrk_alloc( jpi,jpj, zn, zah, zhw, zross )
+      CALL wrk_alloc( jpi,jpj,jpk, uslp, vslp, wslpi, wslpj )
+ 
+      IF ( ln_traldf_iso ) THEN  
+         uslp  = uslp_iso  
+         vslp  = vslp_iso  
+         wslpi = wslpi_iso  
+         wslpj = wslpj_iso  
+      ELSEIF ( ln_traldf_hor ) THEN  
+         uslp  = uslp_hor  
+         vslp  = vslp_hor  
+         wslpi = wslpi_hor  
+         wslpj = wslpj_hor  
+      ENDIF  
 
       IF( kt == nit000 ) THEN
@@ -243 +257 @@
       !  
       CALL wrk_dealloc( jpi,jpj, zn, zah, zhw, zross )
+      CALL wrk_dealloc( jpi,jpj,jpk,  uslp, vslp, wslpi, wslpj )
       !
       IF( nn_timing == 1 )  CALL timing_stop('ldf_eiv')