Changeset 334 for trunk/NEMO/TOP_SRC
- Timestamp:
- 2005-11-14T12:53:13+01:00 (18 years ago)
- Location:
- trunk/NEMO/TOP_SRC/TRP
- Files:
-
- 26 edited
-
trcadv_cen2.F90 (modified) (6 diffs)
-
trcadv_muscl.F90 (modified) (5 diffs)
-
trcadv_muscl2.F90 (modified) (5 diffs)
-
trcadv_smolar.F90 (modified) (4 diffs)
-
trcadv_tvd.F90 (modified) (6 diffs)
-
trcbbc.F90 (modified) (4 diffs)
-
trcbbl.F90 (modified) (10 diffs)
-
trcbbl_adv.h90 (modified) (3 diffs)
-
trcdmp.F90 (modified) (4 diffs)
-
trcldf_bilap.F90 (modified) (4 diffs)
-
trcldf_bilapg.F90 (modified) (4 diffs)
-
trcldf_iso.F90 (modified) (4 diffs)
-
trcldf_iso_zps.F90 (modified) (4 diffs)
-
trcldf_lap.F90 (modified) (4 diffs)
-
trcnxt.F90 (modified) (4 diffs)
-
trcrad.F90 (modified) (6 diffs)
-
trcsbc.F90 (modified) (4 diffs)
-
trcstp.F90 (modified) (3 diffs)
-
trctrp.F90 (modified) (3 diffs)
-
trctrp_ctl.F90 (modified) (2 diffs)
-
trctrp_lec.F90 (modified) (1 diff)
-
trczdf_exp.F90 (modified) (4 diffs)
-
trczdf_imp.F90 (modified) (5 diffs)
-
trczdf_iso.F90 (modified) (5 diffs)
-
trczdf_iso_vopt.F90 (modified) (8 diffs)
-
zpshde_trc.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/TRP/trcadv_cen2.F90
r274 r334 14 14 USE trc ! ocean passive tracers variables 15 15 USE trcbbl ! advective passive tracers in the BBL 16 USE prtctl_trc 16 17 17 18 IMPLICIT NONE … … 24 25 # include "passivetrc_substitute.h90" 25 26 !!---------------------------------------------------------------------- 26 !! TOP 1.0, LOCEAN-IPSL (2005) 27 !! $Header$ 28 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 27 !! OPA 9.0 , LODYC-IPSL (2003) 29 28 !!---------------------------------------------------------------------- 30 29 … … 133 132 zfui1, zfvj1 ! " " 134 133 #endif 134 CHARACTER (len=22) :: charout 135 135 !!---------------------------------------------------------------------- 136 136 … … 154 154 155 155 156 DO jn = 1, jptra157 158 156 ! Upstream / centered scheme indicator 159 157 ! ------------------------------------ 160 158 161 162 163 164 165 159 DO jk = 1, jpk 160 DO jj = 1, jpj 161 DO ji = 1, jpi 162 zind(ji,jj,jk) = MAX ( upsrnfh(ji,jj) * upsrnfz(jk), & ! changing advection scheme near runoff 163 & upsadv(ji,jj) & ! in the vicinity of some straits 166 164 #if defined key_ice_lim 167 & , tmask(ji,jj,jk) & ! half upstream tracer fluxes 168 & * MAX( 0., SIGN( 1., fzptn(ji,jj) & ! if tn < ("freezing"+0.1 ) 169 & +0.1-tn(ji,jj,jk) ) ) & 170 #endif 171 & ) 172 END DO 165 & , tmask(ji,jj,jk) & ! half upstream tracer fluxes 166 & * MAX( 0., SIGN( 1., fzptn(ji,jj) & ! if tn < ("freezing"+0.1 ) 167 & +0.1-tn(ji,jj,jk) ) ) & 168 #endif 169 & ) 173 170 END DO 174 171 END DO 175 176 172 END DO 173 174 175 DO jn = 1, jptra 177 176 ! I. Horizontal advective fluxes 178 177 ! ------------------------------ … … 256 255 END DO ! End of slab 257 256 ! ! =============== 258 259 IF(ln_ctl) THEN ! print mean trends (used for debugging) 260 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 261 WRITE(numout,*) ' trc/had - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 262 tra_ctl(jn) = ztra 263 ENDIF 264 265 ! II. Vertical advection 266 ! ---------------------- 257 ENDDO 258 259 IF(ln_ctl) THEN ! print mean trends (used for debugging) 260 WRITE(charout, FMT="('centered2 - had')") 261 CALL prt_ctl_trc_info(charout) 262 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 263 ENDIF 264 265 ! II. Vertical advection 266 ! ---------------------- 267 DO jn = 1, jptra 267 268 268 269 ! Bottom value : flux set to zero … … 320 321 END DO 321 322 322 IF(ln_ctl) THEN ! print mean trends (used for debugging)323 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )324 WRITE(numout,*) ' trc/zad - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' centered2'325 tra_ctl(jn) = ztra326 ENDIF327 328 323 END DO 329 324 325 IF(ln_ctl) THEN ! print mean trends (used for debugging) 326 WRITE(charout, FMT="('centered - zad')") 327 CALL prt_ctl_trc_info(charout) 328 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 329 ENDIF 330 330 331 331 END SUBROUTINE trc_adv_cen2 -
trunk/NEMO/TOP_SRC/TRP/trcadv_muscl.F90
r274 r334 15 15 USE trcbbl ! advective passive tracers in the BBL 16 16 USE lib_mpp 17 USE prtctl_trc ! Print control for debbuging 17 18 18 19 IMPLICIT NONE … … 25 26 # include "passivetrc_substitute.h90" 26 27 !!---------------------------------------------------------------------- 27 !! TOP 1.0, LOCEAN-IPSL (2005) 28 !! $Header$ 29 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 28 !! OPA 9.0 , LODYC-IPSL (2003) 30 29 !!---------------------------------------------------------------------- 31 30 … … 79 78 REAL(wp) :: zfui, zfvj 80 79 #endif 81 80 CHARACTER (len=22) :: charout 82 81 !!---------------------------------------------------------------------- 83 82 … … 231 230 END DO 232 231 END DO 233 234 IF(ln_ctl) THEN ! print mean trends (used for debugging) 235 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 236 WRITE(numout,*) ' trc/had - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' muscl' 237 tra_ctl(jn) = ztra 238 ENDIF 232 ENDDO 233 234 IF(ln_ctl) THEN ! print mean trends (used for debugging) 235 WRITE(charout, FMT="('muscl - had')") 236 CALL prt_ctl_trc_info(charout) 237 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 238 ENDIF 239 239 240 240 ! II. Vertical advective fluxes 241 241 ! ----------------------------- 242 242 243 DO jn = 1, jptra 243 244 ! First guess of the slope 244 245 ! interior values … … 317 318 END DO 318 319 319 IF(ln_ctl) THEN ! print mean trends (used for debugging)320 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )321 WRITE(numout,*) ' trc/zad - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' muscl'322 tra_ctl(jn) = ztra323 ENDIF324 325 320 END DO 321 322 IF(ln_ctl) THEN ! print mean trends (used for debugging) 323 WRITE(charout, FMT="('muscl - zad')") 324 CALL prt_ctl_trc_info(charout) 325 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 326 ENDIF 326 327 327 328 END SUBROUTINE trc_adv_muscl -
trunk/NEMO/TOP_SRC/TRP/trcadv_muscl2.F90
r274 r334 14 14 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 15 15 USE trcbbl ! advective passive tracers in the BBL 16 USE prtctl_trc ! Print control for debbuging 16 17 17 18 IMPLICIT NONE … … 28 29 # include "passivetrc_substitute.h90" 29 30 !!---------------------------------------------------------------------- 30 !! TOP 1.0, LOCEAN-IPSL (2005) 31 !! $Header$ 32 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 31 !! OPA 9.0 , LODYC-IPSL (2003) 33 32 !!---------------------------------------------------------------------- 34 33 … … 82 81 REAL(wp) :: zfui, zfvj 83 82 #endif 83 CHARACTER (len=22) :: charout 84 84 !!---------------------------------------------------------------------- 85 85 … … 279 279 END DO 280 280 END DO 281 282 IF(ln_ctl) THEN ! print mean trends (used for debugging) 283 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 284 WRITE(numout,*) ' trc/had - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' muscl' 285 tra_ctl(jn) = ztra 286 ENDIF 287 288 ! II. Vertical advective fluxes 289 ! ----------------------------- 281 ENDDO 282 283 IF(ln_ctl) THEN ! print mean trends (used for debugging) 284 WRITE(charout, FMT="('muscl2 - had')") 285 CALL prt_ctl_trc_info(charout) 286 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 287 ENDIF 288 289 ! II. Vertical advective fluxes 290 ! ----------------------------- 291 292 DO jn = 1, jptra 290 293 291 294 ! First guess of the slope … … 379 382 END DO 380 383 381 IF(ln_ctl) THEN ! print mean trends (used for debugging)382 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )383 WRITE(numout,*) ' trc/zad - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' muscl2'384 tra_ctl(jn) = ztra385 ENDIF386 387 384 END DO 385 386 IF(ln_ctl) THEN ! print mean trends (used for debugging) 387 WRITE(charout, FMT="('muscl2 - zad')") 388 CALL prt_ctl_trc_info(charout) 389 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 390 ENDIF 388 391 389 392 END SUBROUTINE trc_adv_muscl2 -
trunk/NEMO/TOP_SRC/TRP/trcadv_smolar.F90
r274 r334 15 15 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 16 16 USE trcbbl ! advective passive tracers in the BBL 17 USE prtctl_trc ! Print control for debbuging 17 18 18 19 IMPLICIT NONE … … 29 30 # include "passivetrc_substitute.h90" 30 31 !!---------------------------------------------------------------------- 31 !! TOP 1.0, LOCEAN-IPSL (2005) 32 !! $Header$ 33 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 32 !! OPA 9.0 , LODYC-IPSL (2003) 34 33 !!---------------------------------------------------------------------- 35 34 CONTAINS … … 96 95 zbuf 97 96 97 #if defined key_trc_diatrd 98 REAL(wp) :: zgm, zgz 99 #endif 100 98 101 REAL(wp) :: zbtr, ztra 99 102 REAL(wp) :: zfp_ui, zfp_vj, zfm_ui, zfm_vj, zfp_w, zfm_w 100 103 CHARACTER (len=22) :: charout 101 104 !!---------------------------------------------------------------------- 102 105 … … 496 499 ! Lateral boundary conditions on trtrd: 497 500 498 CALL lbc_lnk( trtrd(1,1,1,jn,1), 'T', 1. ) 499 CALL lbc_lnk( trtrd(1,1,1,jn,2), 'T', 1. ) 500 CALL lbc_lnk( trtrd(1,1,1,jn,3), 'T', 1. ) 501 #endif 502 503 IF(ln_ctl) THEN ! print mean trends (used for debugging) 504 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 505 WRITE(numout,*) ' trc/zad - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' smolar' 506 tra_ctl(jn) = ztra 507 ENDIF 501 CALL lbc_lnk( trtrd(:,:,:,jn,1), 'T', 1. ) 502 CALL lbc_lnk( trtrd(:,:,:,jn,2), 'T', 1. ) 503 CALL lbc_lnk( trtrd(:,:,:,jn,3), 'T', 1. ) 504 #endif 505 506 508 507 ! END of tracer loop 509 508 ! ================== 510 509 ENDDO 510 511 IF(ln_ctl) THEN ! print mean trends (used for debugging) 512 WRITE(charout, FMT="('smolar - adv')") 513 CALL prt_ctl_trc_info(charout) 514 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 515 ENDIF 511 516 512 517 END SUBROUTINE trc_adv_smolar -
trunk/NEMO/TOP_SRC/TRP/trcadv_tvd.F90
r274 r334 16 16 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 17 17 USE trcbbl ! advective passive tracers in the BBL 18 USE prtctl_trc ! Print control for debbuging 18 19 19 20 IMPLICIT NONE … … 26 27 # include "passivetrc_substitute.h90" 27 28 !!---------------------------------------------------------------------- 28 !! TOP 1.0, LOCEAN-IPSL (2005) 29 !! $Header$ 30 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 29 !! OPA 9.0 , LODYC-IPSL (2003) 31 30 !!---------------------------------------------------------------------- 32 31 … … 79 78 zfp_ui, zfp_vj, zfp_wk, & ! " " 80 79 zfm_ui, zfm_vj, zfm_wk ! " " 80 81 CHARACTER (len=22) :: charout 81 82 !!---------------------------------------------------------------------- 82 83 … … 234 235 END DO 235 236 236 IF(ln_ctl) THEN ! print mean trends (used for debugging)237 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 238 WRITE(numout,*) ' trc/zad - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn), ' tvd'239 tra_ctl(jn) = ztra240 ENDIF241 242 END DO237 END DO 238 239 IF(ln_ctl) THEN ! print mean trends (used for debugging) 240 WRITE(charout, FMT="('tvd - adv')") 241 CALL prt_ctl_trc_info(charout) 242 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 243 ENDIF 243 244 244 245 END SUBROUTINE trc_adv_tvd … … 289 290 ! -------------------- 290 291 ! large negative value (-zbig) inside land 291 WHERE( tmask(:,:,:) == 0. ) 292 pbef(:,:,:) = -zbig 293 paft(:,:,:) = -zbig 294 ENDWHERE 292 ! large negative value (-zbig) inside land 293 pbef(:,:,:) = pbef(:,:,:) * tmask(:,:,:) - zbig * ( 1.e0 - tmask(:,:,:) ) 294 paft(:,:,:) = paft(:,:,:) * tmask(:,:,:) - zbig * ( 1.e0 - tmask(:,:,:) ) 295 295 ! search maximum in neighbourhood 296 296 DO jk = 1, jpkm1 … … 309 309 END DO 310 310 ! large positive value (+zbig) inside land 311 WHERE( tmask(:,:,:) == 0. ) 312 pbef(:,:,:) = +zbig 313 paft(:,:,:) = +zbig 314 ENDWHERE 311 pbef(:,:,:) = pbef(:,:,:) * tmask(:,:,:) + zbig * ( 1.e0 - tmask(:,:,:) ) 312 paft(:,:,:) = paft(:,:,:) * tmask(:,:,:) + zbig * ( 1.e0 - tmask(:,:,:) ) 315 313 ! search minimum in neighbourhood 316 314 DO jk = 1, jpkm1 -
trunk/NEMO/TOP_SRC/TRP/trcbbc.F90
r274 r334 14 14 USE oce_trc ! ocean dynamics and active tracers variables 15 15 USE trc ! ocean passive tracers variables 16 USE prtctl_trc ! Print control for debbuging 16 17 17 18 IMPLICIT NONE … … 39 40 # include "passivetrc_substitute.h90" 40 41 !!---------------------------------------------------------------------- 41 !! TOP 1.0, LOCEAN-IPSL (2005) 42 !! $Header$ 43 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 42 !! OPA 9.0 , LODYC-IPSL (2003) 44 43 !!---------------------------------------------------------------------- 45 44 … … 83 82 #endif 84 83 REAL(wp) :: ztra ! temporary scalar 84 CHARACTER (len=22) :: charout 85 85 !!---------------------------------------------------------------------- 86 86 … … 106 106 END DO 107 107 #endif 108 IF(ln_ctl) THEN ! print mean trends (used for debugging)109 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 110 WRITE(numout,*) ' trc/bbc - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)111 tra_ctl(jn) = ztra112 ENDIF113 END DO114 108 END DO 109 110 IF(ln_ctl) THEN ! print mean trends (used for debugging) 111 WRITE(charout, FMT="('bbc')") 112 CALL prt_ctl_trc_info(charout) 113 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 114 ENDIF 115 115 END SELECT 116 116 -
trunk/NEMO/TOP_SRC/TRP/trcbbl.F90
r274 r334 18 18 USE oce_trc ! ocean dynamics and active tracers variables 19 19 USE trc ! ocean passive tracers variables 20 USE prtctl_trc ! Print control for debbuging 20 21 21 22 IMPLICIT NONE … … 48 49 # include "passivetrc_substitute.h90" 49 50 !!---------------------------------------------------------------------- 50 !! TOP 1.0, LOCEAN-IPSL (2005) 51 !! $Header$ 52 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 51 !! OPA 9.0 , LODYC-IPSL (2003) 53 52 !!---------------------------------------------------------------------- 54 53 … … 117 116 ztnb, zsnb, zdep, & 118 117 ztrb, zahu, zahv 118 CHARACTER (len=22) :: charout 119 119 REAL(wp) :: & 120 120 fsalbt, pft, pfs, pfh ! statement function … … 204 204 # endif 205 205 206 !! 207 !! OFFLINE VERSION OF DIFFUSIVE BBL 208 !! 209 #if defined key_off_tra 210 211 ! 2. Additional second order diffusive trends 212 ! ------------------------------------------- 213 214 DO jn = 1, jptra 215 ! first derivative (gradient) 216 217 # if defined key_vectopt_loop && ! defined key_autotasking 218 jj = 1 219 DO ji = 1, jpij ! vector opt. (forced unrolling) 220 # else 221 DO jj = 1, jpj 222 DO ji = 1, jpi 223 # endif 224 ik = mbkt(ji,jj) 225 ztrb(ji,jj) = trb(ji,jj,ik,jn) * tmask(ji,jj,1) 226 # if ! defined key_vectopt_loop || defined key_autotasking 227 END DO 228 # endif 229 END DO 230 231 # if defined key_vectopt_loop && ! defined key_autotasking 232 jj = 1 233 DO ji = 1, jpij-jpi ! vector opt. (forced unrolling) 234 # else 235 DO jj = 1, jpjm1 236 DO ji = 1, jpim1 237 # endif 238 zkx(ji,jj) = bblx(ji,jj) * zahu(ji,jj) * ( ztrb(ji+1,jj) - ztrb(ji,jj) ) 239 zky(ji,jj) = bbly(ji,jj) * zahv(ji,jj) * ( ztrb(ji,jj+1) - ztrb(ji,jj) ) 240 # if ! defined key_vectopt_loop || defined key_autotasking 241 END DO 242 # endif 243 END DO 244 !! 245 !! ONLINE VERSION OF DIFFUSIVE BBL 246 !! 247 #else 206 248 ! 1. Criteria of additional bottom diffusivity: grad(rho).grad(h)<0 207 249 ! -------------------------------------------- … … 257 299 END DO 258 300 259 260 301 ! 2. Additional second order diffusive trends 261 302 ! ------------------------------------------- … … 263 304 DO jn = 1, jptra 264 305 ! first derivative (gradient) 265 306 266 307 # if defined key_vectopt_loop && ! defined key_autotasking 267 308 jj = 1 … … 271 312 DO ji = 1, jpi 272 313 # endif 273 ik = mbkt(ji,jj) 314 ik = mbkt(ji,jj) 274 315 ztrb(ji,jj) = trb(ji,jj,ik,jn) * tmask(ji,jj,1) 275 316 # if ! defined key_vectopt_loop || defined key_autotasking … … 277 318 # endif 278 319 END DO 279 280 320 # if defined key_vectopt_loop && ! defined key_autotasking 281 321 jj = 1 … … 291 331 # endif 292 332 END DO 333 334 #endif 293 335 294 336 IF( cp_cfg == "orca" ) THEN … … 341 383 # endif 342 384 END DO 343 IF(ln_ctl) THEN ! print mean trends (used for debugging) 344 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 345 WRITE(numout,*) ' trc/bbl - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 346 tra_ctl(jn) = ztra 347 ENDIF 348 END DO 349 385 386 END DO 387 388 IF(ln_ctl) THEN ! print mean trends (used for debugging) 389 WRITE(charout, FMT="('bbl - dif')") 390 CALL prt_ctl_trc_info(charout) 391 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 392 ENDIF 350 393 351 394 END SUBROUTINE trc_bbl_dif -
trunk/NEMO/TOP_SRC/TRP/trcbbl_adv.h90
r274 r334 4 4 5 5 !!---------------------------------------------------------------------- 6 !! TOP 1.0, LOCEAN-IPSL (2005) 7 !! $Header$ 8 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 6 !! OPA 9.0 , LODYC-IPSL (2003) 9 7 !!---------------------------------------------------------------------- 10 8 … … 79 77 REAL(wp) :: & 80 78 fsalbt, pft, pfs, pfh ! statement function 79 CHARACTER (len=22) :: charout 81 80 !!---------------------------------------------------------------------- 82 81 ! ratio alpha/beta … … 384 383 #endif 385 384 END DO 386 IF(ln_ctl) THEN ! print mean trends (used for debugging) 387 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )388 WRITE(numout,*) ' trc/bbl - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 389 tra_ctl(jn) = ztra390 ENDIF391 392 END DO393 385 386 END DO 387 388 IF(ln_ctl) THEN ! print mean trends (used for debugging) 389 WRITE(charout, FMT="('bbl - adv')") 390 CALL prt_ctl_trc_info(charout) 391 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 392 ENDIF 394 393 ! 6. Vertical advection velocities 395 394 ! -------------------------------- -
trunk/NEMO/TOP_SRC/TRP/trcdmp.F90
r274 r334 19 19 USE trctrp_lec ! passive tracers transport 20 20 USE trcdta 21 USE prtctl_trc ! Print control for debbuging 21 22 22 23 IMPLICIT NONE … … 35 36 # include "passivetrc_substitute.h90" 36 37 !!---------------------------------------------------------------------- 37 !! TOP 1.0, LOCEAN-IPSL (2005) 38 !! $Header$ 39 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 38 !! OPA 9.0 , LODYC-IPSL (2003) 40 39 !!---------------------------------------------------------------------- 41 40 … … 74 73 INTEGER :: ji, jj, jk, jn ! dummy loop indices 75 74 REAL(wp) :: ztest, ztra, zdt ! temporary scalars 75 CHARACTER (len=22) :: charout 76 76 !!---------------------------------------------------------------------- 77 77 … … 157 157 ENDIF 158 158 159 IF(ln_ctl) THEN ! print mean trends (used for debugging)160 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )161 WRITE(numout,*) ' trc/dmp - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)162 tra_ctl(jn) = ztra163 ENDIF164 165 159 END DO 166 160 161 IF(ln_ctl) THEN ! print mean trends (used for debugging) 162 WRITE(charout, FMT="('dmp')") 163 CALL prt_ctl_trc_info(charout) 164 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 165 ENDIF 166 167 167 trb(:,:,:,:) = trn(:,:,:,:) 168 168 -
trunk/NEMO/TOP_SRC/TRP/trcldf_bilap.F90
r274 r334 13 13 USE trc ! ocean passive tracers variables 14 14 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 15 USE prtctl_trc ! Print control for debbuging 15 16 16 17 IMPLICIT NONE … … 23 24 # include "passivetrc_substitute.h90" 24 25 !!---------------------------------------------------------------------- 25 !! TOP 1.0, LOCEAN-IPSL (2005) 26 !! $Header$ 27 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 26 !! OPA 9.0 , LODYC-IPSL (2003) 28 27 !!---------------------------------------------------------------------- 29 28 … … 95 94 REAL(wp), DIMENSION(jpi,jpj,jpk) :: & 96 95 ztu, ztv ! workspace 96 CHARACTER (len=22) :: charout 97 97 !!---------------------------------------------------------------------- 98 98 … … 203 203 #if defined key_trc_diatrd 204 204 ! Lateral boundary conditions on the laplacian zlt (unchanged sgn) 205 CALL lbc_lnk( trtrd(1,1,1,jn,5), 'T', 1. ) 206 #endif 207 IF(ln_ctl) THEN ! print mean trends (used for debugging) 208 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 209 WRITE(numout,*) ' trc/ldf - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 210 tra_ctl(jn) = ztra 211 ENDIF 212 205 CALL lbc_lnk( trtrd(:,:,:,jn,5), 'T', 1. ) 206 #endif 213 207 END DO 208 209 IF(ln_ctl) THEN ! print mean trends (used for debugging) 210 WRITE(charout, FMT="('ldf - bilap')") 211 CALL prt_ctl_trc_info(charout) 212 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 213 ENDIF 214 214 215 215 END SUBROUTINE trc_ldf_bilap -
trunk/NEMO/TOP_SRC/TRP/trcldf_bilapg.F90
r274 r334 16 16 USE trc ! ocean passive tracers variables 17 17 USE lbclnk ! ocean lateral boundary condition (or mpp link) 18 USE prtctl_trc ! Print control for debbuging 18 19 19 20 IMPLICIT NONE … … 26 27 # include "passivetrc_substitute.h90" 27 28 !!---------------------------------------------------------------------- 28 !! TOP 1.0, LOCEAN-IPSL (2005) 29 !! $Header$ 30 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 29 !! OPA 9.0 , LODYC-IPSL (2003) 31 30 !!---------------------------------------------------------------------- 32 31 … … 71 70 wk1, wk2 ! work array used for rotated biharmonic 72 71 ! operator on tracers and/or momentum 72 CHARACTER (len=22) :: charout 73 73 !!---------------------------------------------------------------------- 74 74 … … 118 118 ! ! =============== 119 119 120 IF(ln_ctl) THEN ! print mean trends (used for debugging)121 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )122 WRITE(numout,*) ' trc/ldf - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)123 tra_ctl(jn) = ztra124 ENDIF125 126 120 END DO 121 122 IF(ln_ctl) THEN ! print mean trends (used for debugging) 123 WRITE(charout, FMT="('ldf - bilapg')") 124 CALL prt_ctl_trc_info(charout) 125 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 126 ENDIF 127 127 128 128 END SUBROUTINE trc_ldf_bilapg -
trunk/NEMO/TOP_SRC/TRP/trcldf_iso.F90
r274 r334 15 15 USE oce_trc ! ocean dynamics and tracers variables 16 16 USE trc ! ocean passive tracers variables 17 USE prtctl_trc ! Print control for debbuging 17 18 18 19 IMPLICIT NONE … … 24 25 !! * Substitutions 25 26 # include "passivetrc_substitute.h90" 26 !!----------------------------------------------------------------------27 !! TOP 1.0, LOCEAN-IPSL (2005)28 !! $Header$29 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt30 27 !!---------------------------------------------------------------------- 31 28 … … 98 95 REAL(wp), DIMENSION(jpi,jpj,jpk) :: & 99 96 zftv ! workspace 97 CHARACTER (len=22) :: charout 100 98 !!---------------------------------------------------------------------- 101 99 … … 203 201 END DO ! End of slab 204 202 ! ! =============== 205 IF(ln_ctl) THEN ! print mean trends (used for debugging)206 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )207 WRITE(numout,*) ' trc/ldf - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)208 tra_ctl(jn) = ztra209 ENDIF210 203 211 204 END DO 205 206 IF(ln_ctl) THEN ! print mean trends (used for debugging) 207 WRITE(charout, FMT="('ldf - iso')") 208 CALL prt_ctl_trc_info(charout) 209 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 210 ENDIF 212 211 213 212 END SUBROUTINE trc_ldf_iso -
trunk/NEMO/TOP_SRC/TRP/trcldf_iso_zps.F90
r274 r334 14 14 USE oce_trc ! ocean dynamics and active tracers variables 15 15 USE trc ! ocean passive tracers variables 16 USE prtctl_trc ! Print control for debbuging 16 17 17 18 IMPLICIT NONE … … 24 25 # include "passivetrc_substitute.h90" 25 26 !!---------------------------------------------------------------------- 26 !! TOP 1.0, LOCEAN-IPSL (2005) 27 !! $Header$ 28 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 27 !! OPA 9.0 , LODYC-IPSL (2003) 29 28 !!---------------------------------------------------------------------- 30 29 … … 102 101 zcg1,zcg2 103 102 #endif 104 103 CHARACTER (len=22) :: charout 105 104 !!---------------------------------------------------------------------- 106 105 … … 240 239 END DO ! End of slab 241 240 ! ! =============== 242 IF(ln_ctl) THEN ! print mean trends (used for debugging)243 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )244 WRITE(numout,*) ' trc/ldf - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)245 tra_ctl(jn) = ztra246 ENDIF247 248 241 END DO 242 243 IF(ln_ctl) THEN ! print mean trends (used for debugging) 244 WRITE(charout, FMT="('ldf - iso/zps')") 245 CALL prt_ctl_trc_info(charout) 246 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 247 ENDIF 249 248 250 249 END SUBROUTINE trc_ldf_iso_zps -
trunk/NEMO/TOP_SRC/TRP/trcldf_lap.F90
r274 r334 12 12 USE oce_trc ! ocean dynamics and active tracers variables 13 13 USE trc ! ocean passive tracers variables 14 USE prtctl_trc ! Print control for debbuging 14 15 15 16 IMPLICIT NONE … … 22 23 # include "passivetrc_substitute.h90" 23 24 !!---------------------------------------------------------------------- 24 !! TOP 1.0, LOCEAN-IPSL (2005) 25 !! $Header$ 26 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 25 !! OPA 9.0 , LODYC-IPSL (2003) 27 26 !!---------------------------------------------------------------------- 28 27 … … 79 78 REAL(wp) :: & 80 79 ztra, ztrax, ztray ! workspace 80 CHARACTER (len=22) :: charout 81 81 !!---------------------------------------------------------------------- 82 82 … … 139 139 ! ! ============= 140 140 141 IF(ln_ctl) THEN ! print mean trends (used for debugging) 142 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 143 WRITE(numout,*) ' trc/ldf - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 144 tra_ctl(jn) = ztra 145 ENDIF 141 END DO 146 142 147 END DO 143 IF(ln_ctl) THEN ! print mean trends (used for debugging) 144 WRITE(charout, FMT="('ldf - lap')") 145 CALL prt_ctl_trc_info(charout) 146 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 147 ENDIF 148 148 149 149 END SUBROUTINE trc_ldf_lap -
trunk/NEMO/TOP_SRC/TRP/trcnxt.F90
r274 r334 13 13 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 14 14 USE trctrp_lec ! pasive tracers transport 15 USE prtctl_trc ! Print control for debbuging 15 16 16 17 IMPLICIT NONE … … 20 21 PUBLIC trc_nxt ! routine called by step.F90 21 22 !!---------------------------------------------------------------------- 22 !! TOP 1.0, LOCEAN-IPSL (2005) 23 !! $Header$ 24 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 23 !! OPA 9.0 , LODYC-IPSL (2003) 25 24 !!---------------------------------------------------------------------- 26 25 … … 67 66 INTEGER :: ji, jj, jk,jn ! dummy loop indices 68 67 REAL(wp) :: zfact, ztra ! temporary scalar 68 CHARACTER (len=22) :: charout 69 69 !!---------------------------------------------------------------------- 70 70 … … 139 139 END DO ! End of slab 140 140 ! ! =============== 141 IF(ln_ctl) THEN ! print mean field (used for debugging) 142 ztra = SUM( trn(2:nictle,2:njctle,1:jpkm1,jn)*tmask(2:nictle,2:njctle,1:jpkm1) ) 143 WRITE(numout,*) ' trc/nxt - ',ctrcnm(jn),' : ', ztra 144 ENDIF 141 END DO 145 142 146 END DO 143 IF(ln_ctl) THEN ! print mean trends (used for debugging) 144 WRITE(charout, FMT="('nxt')") 145 CALL prt_ctl_trc_info(charout) 146 CALL prt_ctl_trc(tab4d=trn, mask=tmask, clinfo=ctrcnm) 147 ENDIF 148 147 149 148 150 END SUBROUTINE trc_nxt -
trunk/NEMO/TOP_SRC/TRP/trcrad.F90
r274 r334 12 12 USE trc ! ocean passive tracers variables 13 13 USE lib_mpp 14 USE prtctl_trc ! Print control for debbuging 14 15 15 16 IMPLICIT NONE … … 21 22 # include "passivetrc_substitute.h90" 22 23 !!---------------------------------------------------------------------- 23 !! TOP 1.0, LOCEAN-IPSL (2005) 24 !! $Header$ 25 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 24 !! OPA 9.0 , LODYC-IPSL (2003) 26 25 !!---------------------------------------------------------------------- 27 26 CONTAINS … … 47 46 !! * Local declarations 48 47 INTEGER :: ji, jj, jk, jn ! dummy loop indices 49 REAL(wp) :: ztra 50 #if defined key_trc_hamocc3 || defined key_trc_pisces 48 #if defined key_trc_pisces 51 49 REAL(wp) :: zvolk, trcorb, trmasb ,trcorn, trmasn 52 50 #endif 51 CHARACTER (len=22) :: charout 53 52 !!---------------------------------------------------------------------- 54 53 … … 59 58 ENDIF 60 59 61 #if defined key_trc_hamocc362 DO jn = 1, jptra63 trcorb = 0.64 trmasb = 0.65 trcorn = 0.66 trmasn = 0.67 DO jk = 1, jpkm168 DO jj = 2, jpjm169 DO ji = 2, jpim170 zvolk = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) &71 #if defined key_off_degrad72 & * facvol(ji,jj,jk) &73 #endif74 & * tmask(ji,jj,jk)75 60 76 trcorb = trcorb + MIN( 0., trb(ji,jj,jk,jn) ) * zvolk 77 trcorn = trcorn + MIN( 0., trn(ji,jj,jk,jn) ) * zvolk 78 79 trb(ji,jj,jk,jn) = MAX( 0. , trb(ji,jj,jk,jn) ) 80 trn(ji,jj,jk,jn) = MAX( 0. , trn(ji,jj,jk,jn) ) 81 82 trmasb = trmasb + trb(ji,jj,jk,jn) * zvolk 83 trmasn = trmasn + trn(ji,jj,jk,jn) * zvolk 84 END DO 85 END DO 86 END DO 87 88 IF( lk_mpp ) THEN 89 CALL mpp_sum( trcorb ) ! sum over the global domain 90 CALL mpp_sum( trcorn ) ! sum over the global domain 91 CALL mpp_sum( trmasb ) ! sum over the global domain 92 CALL mpp_sum( trmasn ) ! sum over the global domain 93 ENDIF 94 95 DO jk = 1, jpkm1 96 DO jj = 1, jpj 97 DO ji = 1, jpi 98 trb(ji,jj,jk,jn) = MAX( 0., trb(ji,jj,jk,jn) ) 99 trb(ji,jj,jk,jn) = trb(ji,jj,jk,jn) * ( 1. + trcorb/trmasb ) * tmask(ji,jj,jk) 100 trn(ji,jj,jk,jn) = MAX( 0., trn(ji,jj,jk,jn) ) 101 trn(ji,jj,jk,jn) = trn(ji,jj,jk,jn) * ( 1. + trcorb/trmasb ) * tmask(ji,jj,jk) 102 END DO 103 END DO 104 END DO 105 END DO 106 107 #elif defined key_trc_age || defined key_trc_lobster1 61 #if defined key_trc_lobster1 || defined key_cfc 108 62 DO jn = 1, jptra 109 63 DO jk = 1, jpkm1 … … 151 105 ENDIF 152 106 153 IF( trcorb /= 0 ) THEN107 IF( trcorb /= 0 ) THEN 154 108 DO jk = 1, jpkm1 155 109 DO jj = 1, jpj … … 176 130 177 131 #endif 178 DO jn = 1, jptra 179 IF(ln_ctl) THEN ! print mean field(used for debugging)180 ztra = SUM( trn(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )181 WRITE(numout,*) ' trc/rad - ',ctrcnm(jn),' : ', ztra182 ENDIF183 END DO132 133 IF(ln_ctl) THEN ! print mean trends (used for debugging) 134 WRITE(charout, FMT="('rad')") 135 CALL prt_ctl_trc_info(charout) 136 CALL prt_ctl_trc(tab4d=trn, mask=tmask, clinfo=ctrcnm) 137 ENDIF 184 138 185 139 -
trunk/NEMO/TOP_SRC/TRP/trcsbc.F90
r274 r334 11 11 USE oce_trc ! ocean dynamics and active tracers variables 12 12 USE trc ! ocean passive tracers variables 13 USE prtctl_trc ! Print control for debbuging 14 13 15 14 16 IMPLICIT NONE … … 21 23 # include "passivetrc_substitute.h90" 22 24 !!---------------------------------------------------------------------- 23 !! TOP 1.0, LOCEAN-IPSL (2005) 24 !! $Header$ 25 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 25 !! OPA 9.0 , LODYC-IPSL (2003) 26 26 !!---------------------------------------------------------------------- 27 27 … … 60 60 INTEGER :: ji, jj, jn ! dummy loop indices 61 61 REAL(wp) :: ztra, zsrau, zse3t ! temporary scalars 62 CHARACTER (len=22) :: charout 62 63 !!---------------------------------------------------------------------- 63 64 … … 89 90 END DO 90 91 91 IF(ln_ctl) THEN ! print mean trends (used for debugging)92 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )93 WRITE(numout,*) ' trc/sbc - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)94 tra_ctl(jn) = ztra95 ENDIF96 92 END DO 93 94 IF(ln_ctl) THEN ! print mean trends (used for debugging) 95 WRITE(charout, FMT="('sbc')") 96 CALL prt_ctl_trc_info(charout) 97 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 98 ENDIF 97 99 98 100 END SUBROUTINE trc_sbc -
trunk/NEMO/TOP_SRC/TRP/trcstp.F90
r276 r334 13 13 USE trctrp ! passive tracers transport 14 14 USE trcsms ! passive tracers sources and sinks 15 USE trcdia ! passive tracer diagnostics (trc_dia routine) 16 USE trcrst ! restart for passive tracers 15 USE prtctl_trc ! Print control for debbuging 16 USE trcdia 17 USE trcrst 17 18 18 19 IMPLICIT NONE … … 22 23 PUBLIC trc_stp ! called by step 23 24 !!---------------------------------------------------------------------- 24 !! TOP 1.0, LOCEAN-IPSL (2005) 25 !! $Header$ 26 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 25 !! OPA 9.0 , LODYC-IPSL (2003) 27 26 !!---------------------------------------------------------------------- 28 27 29 28 CONTAINS 30 29 31 SUBROUTINE trc_stp( kt )30 SUBROUTINE trc_stp( kt, kindic ) 32 31 !!------------------------------------------------------------------- 33 32 !! *** ROUTINE trc_stp *** … … 44 43 !! * Arguments 45 44 INTEGER, INTENT( in ) :: kt ! ocean time-step index 46 INTEGER :: jn47 REAL(wp) :: ztra45 INTEGER, INTENT( in ) :: kindic 46 CHARACTER (len=25) :: charout 48 47 49 !! * local declarations 50 INTEGER :: indic ! error indicator if < 0 51 !! --------------------------------------------------------------------- 48 IF( kt == nit000 ) CALL trc_dia( kt, kindic ) ! diagnostics init. 52 49 53 indic = 1 ! reset to no error condition 50 ! this ROUTINE is called only every ndttrc time step 51 IF( MOD( kt , ndttrc ) /= 0 ) RETURN 52 53 ! tracers: sink and source 54 IF(ln_ctl) THEN 55 WRITE(charout,FMT="('kt =', I4,' d/m/y =',I2,I2,I4)") kt, nday, nmonth, nyear 56 CALL prt_ctl_trc_info(charout) 57 ENDIF 54 58 55 59 60 CALL trc_sms( kt ) 56 61 57 ! diagnostics init 58 59 IF( kt == nit000 ) CALL trc_dia( kt, indic ) 60 61 ! these ROUTINES are called only every ndttrc time step 62 IF( MOD( kt , ndttrc ) == 0 ) THEN 63 64 ! tracers: sink and source 65 66 CALL trc_sms( kt ) 67 68 DO jn = 1, jptra 69 IF(ln_ctl) THEN ! print mean trends (used for debugging) 70 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 71 WRITE(numout,*) ' trc/sms - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 72 tra_ctl(jn) = ztra 73 ENDIF 74 END DO 75 62 IF(ln_ctl) THEN ! print mean trends (used for debugging) 63 WRITE(charout, FMT="('sms')") 64 CALL prt_ctl_trc_info(charout) 65 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm) 66 ENDIF 76 67 77 68 78 69 ! transport of passive tracers 79 70 80 71 CALL trc_trp( kt ) 81 72 82 DO jn = 1, jptra 83 IF(ln_ctl) THEN ! print mean trends (used for debugging) 84 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 85 WRITE(numout,*) ' trc/trp - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 86 tra_ctl(jn) = ztra 87 ENDIF 88 END DO 73 CALL trc_wri( kt ) ! outputs 89 74 90 ENDIF 91 92 93 94 ! outputs 95 96 CALL trc_wri( kt ) 97 98 99 100 ! diagnostics 101 102 CALL trc_dia( kt, indic ) 103 75 CALL trc_dia( kt, kindic ) ! diagnostics 104 76 105 77 -
trunk/NEMO/TOP_SRC/TRP/trctrp.F90
r274 r334 24 24 25 25 USE trcnxt ! time-stepping (trc_nxt routine) 26 !USE trcrad ! positivity (trc_rad routine)26 USE trcrad ! positivity (trc_rad routine) 27 27 28 28 USE trcadv_cen2 ! 2nd order centered advection (trc_adv_cen2 routine) … … 50 50 # include "domzgr_substitute.h90" 51 51 !!---------------------------------------------------------------------- 52 !! TOP 1.0, LOCEAN-IPSL (2005) 53 !! $Header$ 54 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 52 !! OPA 9.0 , LODYC-IPSL (2003) 55 53 !!---------------------------------------------------------------------- 56 54 … … 125 123 CALL trc_nxt( kt ) ! tracer fields at next time step 126 124 127 !!CALL trc_rad( kt ) ! Correct artificial negative concentrations125 CALL trc_rad( kt ) ! Correct artificial negative concentrations 128 126 ! ! for isopycnal scheme 129 127 130 !IF( lk_zps ) CALL zps_hde_trc( kt, trb, gtru, gtrv ) ! Partial steps: now horizontal gradient128 IF( lk_zps ) CALL zps_hde_trc( kt, trb, gtru, gtrv ) ! Partial steps: now horizontal gradient 131 129 ! ! of passive tracers at the bottom ocean level 132 130 -
trunk/NEMO/TOP_SRC/TRP/trctrp_ctl.F90
r281 r334 27 27 28 28 !!---------------------------------------------------------------------- 29 !! TOP 1.0, LOCEAN-IPSL (2005) 30 !! $Header$ 31 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 29 !! OPA 9.0 , LODYC-IPSL (2003) 32 30 !!---------------------------------------------------------------------- 33 31 … … 312 310 l_trczdf_imp = .TRUE. 313 311 #else 314 IF( l _trazdf_exp ) THEN312 IF( ln_zdfexp ) THEN 315 313 l_trczdf_exp = .TRUE. ! use explicit scheme 316 314 l_trczdf_imp = .FALSE. -
trunk/NEMO/TOP_SRC/TRP/trctrp_lec.F90
r274 r334 70 70 #endif 71 71 !!---------------------------------------------------------------------- 72 !! TOP 1.0, LOCEAN-IPSL (2005) 73 !! $Header$ 74 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 72 !! OPA 9.0 , LODYC-IPSL (2003) 75 73 !!---------------------------------------------------------------------- 76 74 -
trunk/NEMO/TOP_SRC/TRP/trczdf_exp.F90
r274 r334 14 14 USE trc ! ocean passive tracers variables 15 15 USE trctrp_lec ! passive tracers transport 16 USE prtctl_trc ! Print control for debbuging 16 17 17 18 IMPLICIT NONE … … 28 29 # include "passivetrc_substitute.h90" 29 30 !!---------------------------------------------------------------------- 30 !! TOP 1.0, LOCEAN-IPSL (2005) 31 !! $Header$ 32 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 31 !! OPA 9.0, LODYC-IPSL (2003) 33 32 !!---------------------------------------------------------------------- 34 33 … … 81 80 REAL(wp), DIMENSION(jpi,jpk) :: & 82 81 zwx, zwy 82 CHARACTER (len=22) :: charout 83 83 !!--------------------------------------------------------------------- 84 84 … … 161 161 END DO 162 162 163 IF(ln_ctl) THEN ! print mean trends (used for debugging) 164 WRITE(charout, FMT="('zdf - exp')") 165 CALL prt_ctl_trc_info(charout) 166 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 167 ENDIF 168 163 169 END SUBROUTINE trc_zdf_exp 164 170 -
trunk/NEMO/TOP_SRC/TRP/trczdf_imp.F90
r274 r334 13 13 USE trc ! ocean passive tracers variables 14 14 USE trctrp_lec ! passive tracers transport 15 USE prtctl_trc 15 16 16 17 IMPLICIT NONE … … 27 28 # include "passivetrc_substitute.h90" 28 29 !!---------------------------------------------------------------------- 29 !! TOP 1.0, LOCEAN-IPSL (2005) 30 !! $Header$ 31 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 30 !! OPA 9.0 , LODYC-IPSL (2003) 32 31 !!---------------------------------------------------------------------- 33 32 … … 81 80 REAL(wp), DIMENSION(jpi,jpj,jpk,jptra) :: & 82 81 ztrd 82 CHARACTER (len=22) :: charout 83 83 !!--------------------------------------------------------------------- 84 84 … … 183 183 END DO 184 184 185 IF(ln_ctl) THEN 185 IF(ln_ctl) THEN ! print mean trends (used for debugging) 186 ztrd(:,:,:,:) = 0. 186 187 DO jk = 1, jpkm1 187 188 DO jj = 2, jpjm1 … … 192 193 END DO 193 194 ENDIF 194 195 IF(ln_ctl) THEN ! print mean trends (used for debugging) 196 ztra = SUM( ztrd(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) ) 197 WRITE(numout,*) ' trc/zdf - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn) 198 tra_ctl(jn) = ztra 199 ENDIF 200 195 201 196 END DO 197 198 IF(ln_ctl) THEN ! print mean trends (used for debugging) 199 WRITE(charout, FMT="('zdf - imp')") 200 CALL prt_ctl_trc_info(charout) 201 CALL prt_ctl_trc(tab4d=ztrd, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 202 ENDIF 202 203 203 204 END SUBROUTINE trc_zdf_imp -
trunk/NEMO/TOP_SRC/TRP/trczdf_iso.F90
r274 r334 17 17 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 18 18 USE trctrp_lec ! passive tracers transport 19 USE prtctl_trc ! Print control for debbuging 19 20 20 21 IMPLICIT NONE … … 30 31 !! * Substitutions 31 32 # include "passivetrc_substitute.h90" 32 !!----------------------------------------------------------------------33 !! TOP 1.0, LOCEAN-IPSL (2005)34 !! $Header$35 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt36 33 !!---------------------------------------------------------------------- 37 34 … … 148 145 zuwki, zvwki 149 146 #endif 147 CHARACTER (len=22) :: charout 150 148 !!--------------------------------------------------------------------- 151 149 !! OPA 8.5, LODYC-IPSL (2002) … … 349 347 ! ! =============== 350 348 351 IF(ln_ctl) THEN ! print mean trends (used for debugging)352 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )353 WRITE(numout,*) ' trc/zdf 1 - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)354 tra_ctl(jn) = ztra355 ENDIF356 357 349 END DO 358 350 351 IF(ln_ctl) THEN ! print mean trends (used for debugging) 352 WRITE(charout, FMT="('zdf - 1')") 353 CALL prt_ctl_trc_info(charout) 354 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 355 ENDIF 359 356 360 357 DO jn = 1, jptra … … 439 436 ! ! =============== 440 437 441 442 IF(ln_ctl) THEN ! print mean trends (used for debugging)443 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )444 WRITE(numout,*) ' trc/zdf 2 - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)445 tra_ctl(jn) = ztra446 ENDIF447 448 438 END DO 439 440 IF(ln_ctl) THEN ! print mean trends (used for debugging) 441 WRITE(charout, FMT="('zdf - 2')") 442 CALL prt_ctl_trc_info(charout) 443 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 444 ENDIF 449 445 450 446 END SUBROUTINE trc_zdf_iso -
trunk/NEMO/TOP_SRC/TRP/trczdf_iso_vopt.F90
r274 r334 20 20 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 21 21 USE trctrp_lec ! passive tracers transport 22 USE prtctl_trc ! Print control for debbuging 22 23 23 24 IMPLICIT NONE … … 33 34 !! * Substitutions 34 35 # include "passivetrc_substitute.h90" 35 !!----------------------------------------------------------------------36 !! TOP 1.0, LOCEAN-IPSL (2005)37 !! $Header$38 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt39 36 !!---------------------------------------------------------------------- 40 37 … … 55 52 !! * Arguments 56 53 INTEGER, INTENT( in ) :: kt ! ocean time-step index 54 CHARACTER (len=22) :: charout 57 55 !!--------------------------------------------------------------------- 58 56 !! OPA 8.5, LODYC-IPSL (2002) … … 74 72 CALL trc_zdf_iso 75 73 74 IF(ln_ctl) THEN ! print mean trends (used for debugging) 75 WRITE(charout, FMT="('zdf - 1')") 76 CALL prt_ctl_trc_info(charout) 77 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 78 ENDIF 76 79 77 80 ! II. vertical diffusion (including the vertical diagonal part of the rotated tensor) … … 80 83 CALL trc_zdf_zdf( kt ) 81 84 85 IF(ln_ctl) THEN ! print mean trends (used for debugging) 86 WRITE(charout, FMT="('zdf - 2')") 87 CALL prt_ctl_trc_info(charout) 88 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 89 ENDIF 82 90 83 91 END SUBROUTINE trc_zdf_iso_vopt … … 318 326 # endif 319 327 #endif 320 IF(ln_ctl) THEN ! print mean trends (used for debugging)321 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )322 WRITE(numout,*) ' trc/zdf 1 - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)323 tra_ctl(jn) = ztra324 ENDIF325 328 326 329 END DO … … 405 408 REAL(wp), DIMENSION(jpi,jpj,jpk) :: & 406 409 zwz, zwt, ztfw ! temporary workspace arrays 407 REAL(wp) :: ztra !temporary scalars 410 408 411 !!--------------------------------------------------------------------- 409 412 !! OPA 8.5, LODYC-IPSL (2002) … … 552 555 END DO 553 556 554 IF(ln_ctl) THEN ! print mean trends (used for debugging)555 ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )556 WRITE(numout,*) ' trc/zdf 2 - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)557 tra_ctl(jn) = ztra558 ENDIF559 560 557 END DO 561 558 -
trunk/NEMO/TOP_SRC/TRP/zpshde_trc.F90
r274 r334 28 28 !! * Substitutions 29 29 # include "passivetrc_substitute.h90" 30 !!----------------------------------------------------------------------31 !! TOP 1.0, LOCEAN-IPSL (2005)32 !! $Header$33 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt34 30 !!---------------------------------------------------------------------- 35 31
Note: See TracChangeset
for help on using the changeset viewer.
