- Timestamp:
- 2015-07-10T13:28:53+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4765_CNRS_agrif/NEMOGCM/NEMO/TOP_SRC/TRP/trcdmp.F90
r4359 r5581 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 !! … … 125 126 DO jj = 2, jpjm1 126 127 DO ji = fs_2, fs_jpim1 ! vector opt. 127 IF( avt(ji,jj,jk) <= 5.e-4 ) THEN128 IF( avt(ji,jj,jk) <= 5.e-4_wp ) THEN 128 129 ztra = restotr(ji,jj,jk) * ( ztrcdta(ji,jj,jk) - trb(ji,jj,jk,jn) ) 129 130 tra(ji,jj,jk,jn) = tra(ji,jj,jk,jn) + ztra … … 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 ! ! =========== … … 184 185 INTEGER, INTENT( in ) :: kt ! ocean time-step index 185 186 ! 186 INTEGER :: ji, jj, jk, jn, jl, jc ! dummy loop indicesa 187 INTEGER :: ji , jj, jk, jn, jl, jc ! dummy loop indicesa 188 INTEGER :: isrow ! local index 187 189 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrcdta ! 3D workspace 188 190 … … 200 202 ! 201 203 SELECT CASE ( jp_cfg ) 204 ! ! ======================= 205 CASE ( 1 ) ! eORCA_R1 configuration 206 ! ! ======================= 207 isrow = 332 - jpjglo 208 ! 209 ! Caspian Sea 210 nctsi1(1) = 332 ; nctsj1(1) = 243 - isrow 211 nctsi2(1) = 344 ; nctsj2(1) = 275 - isrow 212 ! 202 213 ! ! ======================= 203 214 CASE ( 2 ) ! ORCA_R2 configuration … … 302 313 !!---------------------------------------------------------------------- 303 314 ! 315 INTEGER :: imask !local file handle 316 304 317 IF( nn_timing == 1 ) CALL timing_start('trc_dmp_init') 305 318 ! 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 319 314 320 IF( lzoom ) nn_zdmp_tr = 0 ! restoring to climatology at closed north or south boundaries … … 325 331 & CALL ctl_stop( 'passive trace damping need key_tradmp to compute damping coef.' ) 326 332 ! 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 ENDIF333 ! ! Read damping coefficients from file 334 !Read in mask from file 335 CALL iom_open ( cn_resto_tr, imask) 336 CALL iom_get ( imask, jpdom_autoglo, 'resto', restotr) 337 CALL iom_close( imask ) 332 338 ! 333 339 IF( nn_timing == 1 ) CALL timing_stop('trc_dmp_init')
Note: See TracChangeset
for help on using the changeset viewer.