Changeset 5260 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/TOP_SRC/TRP/trcdmp.F90
- Timestamp:
- 2015-05-12T12:37:15+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/TOP_SRC/TRP/trcdmp.F90
r4359 r5260 23 23 USE prtctl_trc ! Print control for debbuging 24 24 USE trdtra 25 USE trdmod_oce 25 USE trd_oce 26 USE iom 26 27 27 28 IMPLICIT NONE … … 42 43 !!---------------------------------------------------------------------- 43 44 !! NEMO/TOP 3.3 , NEMO Consortium (2010) 44 !! $ Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trcdmp.F90,v 1.11 2006/09/01 14:03:49 opalod Exp$45 !! $Id$ 45 46 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 46 47 !!---------------------------------------------------------------------- … … 75 76 !! ** Action : - update the tracer trends tra with the newtonian 76 77 !! damping trends. 77 !! - save the trends ('key_trdm ld_trc')78 !! - save the trends ('key_trdmxl_trc') 78 79 !!---------------------------------------------------------------------- 79 80 !! … … 151 152 IF( l_trdtrc ) THEN 152 153 ztrtrd(:,:,:) = tra(:,:,:,jn) - ztrtrd(:,:,:) 153 CALL trd_tra( kt, 'TRC', jn, jptra_ trd_dmp, ztrtrd )154 CALL trd_tra( kt, 'TRC', jn, jptra_dmp, ztrtrd ) 154 155 END IF 155 156 ! ! =========== … … 302 303 !!---------------------------------------------------------------------- 303 304 ! 305 INTEGER :: imask !local file handle 306 304 307 IF( nn_timing == 1 ) CALL timing_start('trc_dmp_init') 305 308 ! 306 SELECT CASE ( nn_hdmp_tr )307 CASE ( -1 ) ; IF(lwp) WRITE(numout,*) ' tracer damping in the Med & Red seas only'308 CASE ( 1:90 ) ; IF(lwp) WRITE(numout,*) ' tracer damping poleward of', nn_hdmp_tr, ' degrees'309 CASE DEFAULT310 WRITE(ctmp1,*) ' bad flag value for nn_hdmp_tr = ', nn_hdmp_tr311 CALL ctl_stop(ctmp1)312 END SELECT313 309 314 310 IF( lzoom ) nn_zdmp_tr = 0 ! restoring to climatology at closed north or south boundaries … … 325 321 & CALL ctl_stop( 'passive trace damping need key_tradmp to compute damping coef.' ) 326 322 ! 327 ! ! Damping coefficients initialization328 IF( lzoom ) THEN ; CALL dtacof_zoom( restotr )329 ELSE ; CALL dtacof( nn_hdmp_tr, rn_surf_tr, rn_bot_tr, rn_dep_tr, &330 & nn_file_tr, 'TRC' , restotr)331 ENDIF323 ! ! Read damping coefficients from file 324 !Read in mask from file 325 CALL iom_open ( cn_resto_tr, imask) 326 CALL iom_get ( imask, jpdom_autoglo, 'resto', restotr) 327 CALL iom_close( imask ) 332 328 ! 333 329 IF( nn_timing == 1 ) CALL timing_stop('trc_dmp_init')
Note: See TracChangeset
for help on using the changeset viewer.