Changeset 5989 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OPA_SRC/CRS/crsfld.F90
- Timestamp:
- 2015-12-03T09:10:32+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OPA_SRC/CRS/crsfld.F90
r5260 r5989 11 11 !! other variables needed to be passed to TOP 12 12 !!---------------------------------------------------------------------- 13 USE oce ! ocean dynamics and tracers14 USE dom_oce ! ocean space and time domain15 USE ldftra_oce ! ocean active tracers: lateral physics16 USE sbc_oce ! Surface boundary condition: ocean fields17 USE zdf_oce ! vertical physics: ocean fields18 USE zdfddm ! vertical physics: double diffusion19 USE lbclnk ! ocean lateral boundary conditions (or mpp link)20 USE in_out_manager ! I/O manager21 USE timing ! preformance summary22 USE wrk_nemo ! working array23 13 USE crs 24 14 USE crsdom 25 15 USE crslbclnk 26 USE iom 16 USE oce ! ocean dynamics and tracers 17 USE dom_oce ! ocean space and time domain 18 USE sbc_oce ! Surface boundary condition: ocean fields 19 USE zdf_oce ! vertical physics: ocean fields 20 USE ldftra ! ocean active tracers: lateral diffusivity & EIV coefficients 21 USE zdfddm ! vertical physics: double diffusion 22 ! 23 USE in_out_manager ! I/O manager 24 USE iom ! 25 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 26 USE timing ! preformance summary 27 USE wrk_nemo ! working array 27 28 28 29 IMPLICIT NONE … … 30 31 31 32 PUBLIC crs_fld ! routines called by step.F90 32 33 33 34 34 !! * Substitutions … … 37 37 # include "vectopt_loop_substitute.h90" 38 38 !!---------------------------------------------------------------------- 39 !! NEMO/OPA 3. 3 , NEMO Consortium (2010)39 !! NEMO/OPA 3.7 , NEMO Consortium (2015) 40 40 !! $Id$ 41 41 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) … … 56 56 !! ** Method : 57 57 !!---------------------------------------------------------------------- 58 !! 59 60 INTEGER, INTENT( in ) :: kt ! ocean time-step index 61 !! 62 INTEGER :: ji, jj, jk ! dummy loop indices 63 !! 64 REAL(wp), POINTER, DIMENSION(:,:,:) :: zfse3t, zfse3u, zfse3v, zfse3w ! 3D workspace for e3 65 REAL(wp), POINTER, DIMENSION(:,:,:) :: zt, zs 66 REAL(wp), POINTER, DIMENSION(:,:,:) :: zt_crs, zs_crs ! 67 REAL(wp) :: z2dcrsu, z2dcrsv 68 !! 69 !!---------------------------------------------------------------------- 58 INTEGER, INTENT(in) :: kt ! ocean time-step index 59 ! 60 INTEGER :: ji, jj, jk ! dummy loop indices 61 REAL(wp) :: z2dcrsu, z2dcrsv ! local scalars 62 ! 63 REAL(wp), POINTER, DIMENSION(:,:,:) :: zfse3t, zfse3u, zfse3v, zfse3w ! 3D workspace for e3 64 REAL(wp), POINTER, DIMENSION(:,:,:) :: zt, zt_crs 65 REAL(wp), POINTER, DIMENSION(:,:,:) :: zs, zs_crs 66 !!---------------------------------------------------------------------- 70 67 ! 71 72 68 IF( nn_timing == 1 ) CALL timing_start('crs_fld') 73 69 74 70 ! Initialize arrays 75 CALL wrk_alloc( jpi, jpj, jpk,zfse3t, zfse3w )76 CALL wrk_alloc( jpi, jpj, jpk,zfse3u, zfse3v )77 CALL wrk_alloc( jpi, jpj, jpk, zt, zs)71 CALL wrk_alloc( jpi,jpj,jpk, zfse3t, zfse3w ) 72 CALL wrk_alloc( jpi,jpj,jpk, zfse3u, zfse3v ) 73 CALL wrk_alloc( jpi,jpj,jpk, zt , zs ) 78 74 ! 79 75 CALL wrk_alloc( jpi_crs, jpj_crs, jpk, zt_crs, zs_crs ) … … 169 165 CALL iom_put( "eken", rke_crs ) 170 166 171 ! Horizontal divergence ( following OPA_SRC/DYN/div cur.F90 )167 ! Horizontal divergence ( following OPA_SRC/DYN/divhor.F90 ) 172 168 DO jk = 1, jpkm1 173 169 DO ji = 2, jpi_crsm1
Note: See TracChangeset
for help on using the changeset viewer.