- Timestamp:
- 2012-02-25T16:50:01+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/TRP/trczdf.F90
r3294 r3318 4 4 !! Ocean Passive tracers : vertical diffusive trends 5 5 !!===================================================================== 6 !! History : 9.0 ! 2005-11(G. Madec) Original code7 !! NEMO3.0 ! 2008-01 (C. Ethe, G. Madec) merge TRC-TRA6 !! History : 1.0 ! 2005-11 (G. Madec) Original code 7 !! 3.0 ! 2008-01 (C. Ethe, G. Madec) merge TRC-TRA 8 8 !!---------------------------------------------------------------------- 9 9 #if defined key_top … … 14 14 !! ldf_ctl : initialization, namelist read, and parameters control 15 15 !!---------------------------------------------------------------------- 16 USE oce_trc 17 USE trc 18 USE trcnam_trp 19 USE trazdf_exp 20 USE trazdf_imp 21 USE trd mod_oce22 USE trdtra 23 USE prtctl_trc 16 USE oce_trc ! ocean dynamics and active tracers 17 USE trc ! ocean passive tracers variables 18 USE trcnam_trp ! passive tracers transport namelist variables 19 USE trazdf_exp ! vertical diffusion: explicit (tra_zdf_exp routine) 20 USE trazdf_imp ! vertical diffusion: implicit (tra_zdf_imp routine) 21 USE trd_oce ! trends: ocean variables 22 USE trdtra ! trends: tracer manager 23 USE prtctl_trc ! Print control 24 24 25 25 IMPLICIT NONE 26 26 PRIVATE 27 27 28 PUBLIC trc_zdf 29 PUBLIC trc_zdf_alloc 28 PUBLIC trc_zdf ! called by step.F90 29 PUBLIC trc_zdf_alloc ! called by nemogcm.F90 30 30 31 31 INTEGER :: nzdf = 0 ! type vertical diffusion algorithm used … … 50 50 !!---------------------------------------------------------------------- 51 51 ALLOCATE( r2dt(jpk) , STAT=trc_zdf_alloc ) 52 !53 52 IF( trc_zdf_alloc /= 0 ) CALL ctl_warn('trc_zdf_alloc : failed to allocate array.') 54 !55 53 END FUNCTION trc_zdf_alloc 56 54 … … 62 60 !! ** Purpose : compute the vertical ocean tracer physics. 63 61 !!--------------------------------------------------------------------- 64 INTEGER, INTENT( in ) :: kt! ocean time-step index62 INTEGER, INTENT( in ) :: kt ! ocean time-step index 65 63 ! 66 INTEGER 67 CHARACTER (len=22) 64 INTEGER :: jk, jn 65 CHARACTER (len=22) :: charout 68 66 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztrtrd ! 4D workspace 69 67 !!--------------------------------------------------------------------- … … 106 104 ztrtrd(:,:,jk,jn) = ( ( tra(:,:,jk,jn) - trb(:,:,jk,jn) ) / r2dt(jk) ) - ztrtrd(:,:,jk,jn) 107 105 END DO 108 CALL trd_tra( kt, 'TRC', jn, jptra_ trd_zdf, ztrtrd(:,:,:,jn) )106 CALL trd_tra( kt, 'TRC', jn, jptra_zdf, ztrtrd(:,:,:,jn) ) 109 107 END DO 110 108 CALL wrk_dealloc( jpi, jpj, jpk, jptra, ztrtrd )
Note: See TracChangeset
for help on using the changeset viewer.