Changeset 16 for trunk/NEMO/OPA_SRC/TRD
- Timestamp:
- 2004-02-17T09:06:15+01:00 (20 years ago)
- Location:
- trunk/NEMO/OPA_SRC/TRD
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRD/trddyn.F90
r3 r16 34 34 35 35 !! * Shared module vaiables 36 LOGICAL, PUBLIC, PARAMETER :: lk_trddyn = .TRUE. !momentum trend flag36 LOGICAL, PUBLIC, PARAMETER :: lk_trddyn = .TRUE. !: momentum trend flag 37 37 38 38 !! * Substitutions … … 231 231 DO ji = 1, jpi 232 232 zhke(10) = zhke(10) & 233 &+ ub(ji,jj,1) * tautrd(ji,jj,1) * e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,1) &234 &+ vb(ji,jj,1) * tautrd(ji,jj,2) * e1v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,1)233 & + ub(ji,jj,1) * tautrd(ji,jj,1) * e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,1) & 234 & + vb(ji,jj,1) * tautrd(ji,jj,2) * e1v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,1) 235 235 END DO 236 236 END DO … … 240 240 DO jj = 1, jpj 241 241 DO ji = 1, jpi 242 zpeke = zpeke + zkepe(ji,jj,jk) * g * fsdept(ji,jj,jk) &243 &* e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk)244 END DO 245 END DO 246 END DO 247 248 # if defined key_mpp 249 CALL mpp_sum( zpeke )250 CALL mpp_sum( zumo , 11 )251 CALL mpp_sum( zvmo , 11 )252 CALL mpp_sum( zhke , 10 )253 # endif 242 zpeke = zpeke + zkepe(ji,jj,jk) * grav * fsdept(ji,jj,jk) & 243 & * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) 244 END DO 245 END DO 246 END DO 247 248 IF( lk_mpp ) THEN 249 CALL mpp_sum( zpeke ) 250 CALL mpp_sum( zumo , 11 ) 251 CALL mpp_sum( zvmo , 11 ) 252 CALL mpp_sum( zhke , 10 ) 253 ENDIF 254 254 255 255 … … 423 423 END DO 424 424 END DO 425 # if defined key_mpp 426 CALL mpp_sum( tvols ) 427 CALL mpp_sum( tvolu ) 428 CALL mpp_sum( tvolv ) 429 # endif 425 IF( lk_mpp ) CALL mpp_sum( tvols ) ! sums over the global domain 426 IF( lk_mpp ) CALL mpp_sum( tvolu ) 427 IF( lk_mpp ) CALL mpp_sum( tvolv ) 430 428 431 429 IF(lwp) THEN … … 446 444 !! Default option : NO mementum trend diagnostics 447 445 !!---------------------------------------------------------------------- 448 LOGICAL, PUBLIC, PARAMETER :: lk_trddyn = .FALSE. ! momentum trend flag446 LOGICAL, PUBLIC, PARAMETER :: lk_trddyn = .FALSE. !: momentum trend flag 449 447 CONTAINS 450 448 SUBROUTINE trd_dyn( kt ) ! Empty routine 451 WRITE(*,*) kt449 WRITE(*,*) 'trd_dyn: You should not have seen this print! error?', kt 452 450 END SUBROUTINE trd_dyn 453 451 SUBROUTINE trd_dyn_init ! Empty routine -
trunk/NEMO/OPA_SRC/TRD/trdmld.F90
r3 r16 37 37 38 38 !! * Shared module variables 39 LOGICAL, PUBLIC, PARAMETER :: lk_trdmld = .TRUE. !momentum trend flag39 LOGICAL, PUBLIC, PARAMETER :: lk_trdmld = .TRUE. !: momentum trend flag 40 40 41 41 !! * Module variables … … 609 609 !! Default option : Empty module 610 610 !!---------------------------------------------------------------------- 611 LOGICAL, PUBLIC, PARAMETER :: lk_trdmld = .FALSE. ! momentum trend flag611 LOGICAL, PUBLIC, PARAMETER :: lk_trdmld = .FALSE. !: momentum trend flag 612 612 CONTAINS 613 613 SUBROUTINE trd_mld( kt ) ! Empty routine 614 WRITE(*,*) kt614 WRITE(*,*) 'trd_mld: You should not have seen this print! error?', kt 615 615 END SUBROUTINE trd_mld 616 616 #endif -
trunk/NEMO/OPA_SRC/TRD/trdtra.F90
r3 r16 33 33 34 34 !! * Shared module variables 35 LOGICAL, PUBLIC, PARAMETER :: lk_trdtra = .TRUE. !momentum trend flag35 LOGICAL, PUBLIC, PARAMETER :: lk_trdtra = .TRUE. !: momentum trend flag 36 36 37 37 !! * Substitutions … … 95 95 END DO 96 96 END DO 97 #if defined key_mpp 98 CALL mpp_sum( tvolt ) 99 #endif 97 IF( lk_mpp ) CALL mpp_sum( tvolt ) ! sum over the global domain 98 100 99 IF(lwp) THEN 101 100 WRITE(numout,*) … … 208 207 END DO 209 208 210 #if defined key_mpp 211 CALL mpp_sum( ztmo, 10 )212 CALL mpp_sum( zsmo, 10 )213 CALL mpp_sum( zt2 , 10 )214 CALL mpp_sum( zs2 , 10 )215 #endif 209 IF( lk_mpp ) THEN 210 CALL mpp_sum( ztmo, 10 ) ! sums over the global domain 211 CALL mpp_sum( zsmo, 10 ) 212 CALL mpp_sum( zt2 , 10 ) 213 CALL mpp_sum( zs2 , 10 ) 214 ENDIF 216 215 217 216 ! 4. Print … … 358 357 END DO 359 358 END DO 360 #if defined key_mpp 361 CALL mpp_sum( tvolt ) 362 #endif 363 IF(lwp) THEN 364 WRITE(numout,*) ' total ocean volume at T-point tvolt = ',tvolt 365 ENDIF 359 IF( lk_mpp ) CALL mpp_sum( tvolt ) ! sum over the global domain 360 361 IF(lwp) WRITE(numout,*) ' total ocean volume at T-point tvolt = ',tvolt 366 362 367 363 END SUBROUTINE trd_tra_init … … 371 367 !! Default case : Empty module 372 368 !!---------------------------------------------------------------------- 373 LOGICAL, PUBLIC, PARAMETER :: lk_trdtra = .FALSE. ! momentum trend flag369 LOGICAL, PUBLIC, PARAMETER :: lk_trdtra = .FALSE. !: momentum trend flag 374 370 CONTAINS 375 371 SUBROUTINE trd_tra( kt ) ! Empty routine 376 WRITE(*,*) kt372 WRITE(*,*) 'trd_tra: You should not have seen this print! error?', kt 377 373 END SUBROUTINE trd_tra 378 374 SUBROUTINE trd_tra_init ! Empty routine -
trunk/NEMO/OPA_SRC/TRD/trdtra_oce.F90
r3 r16 12 12 PUBLIC 13 13 14 INTEGER :: & !! !namdia : diagnostics on dynamics and/or tracer trends14 INTEGER :: & !!: namdia : diagnostics on dynamics and/or tracer trends 15 15 ntrd = 10 , & !: time step frequency dynamics and tracers trends 16 16 nctls = 0 !: control surface type for trends vertical integration … … 25 25 !! Trends diagnostics parameters 26 26 !!--------------------------------------------------------------------- 27 INTEGER, PARAMETER :: & 27 INTEGER, PARAMETER :: & !: 28 28 # if defined key_traldf_eiv 29 29 jptrdh = 4, & !: number of 3D horiz trends arrays … … 38 38 !! Trends diagnostics variables 39 39 !!--------------------------------------------------------------------- 40 REAL(wp) :: & 41 tvolt ! volume of the whole ocean computed at t-points42 REAL(wp), DIMENSION(jpi,jpj,jpk,7) :: & 40 REAL(wp) :: & !: 41 tvolt !: volume of the whole ocean computed at t-points 42 REAL(wp), DIMENSION(jpi,jpj,jpk,7) :: & !: 43 43 ttrd !: trends of the temperature tracer equations 44 44 ! ! ttrd(,,,1) : horizontal advection … … 49 49 ! ! ttrd(,,,6) : damping OR vertical EIV 50 50 ! ! ttrd(,,,7) : penetrative solar radiation (T only) 51 REAL(wp), DIMENSION(jpi,jpj,jpk,6) :: & 51 REAL(wp), DIMENSION(jpi,jpj,jpk,6) :: & !: 52 52 strd !: trends of the salinity tracer equations 53 53 ! ! same as ttrd() 54 REAL(wp), DIMENSION(jpi,jpj,jpk,jptrdh) :: & 54 REAL(wp), DIMENSION(jpi,jpj,jpk,jptrdh) :: & !: 55 55 ttrdh, strdh !: ttrdh(,,,1) : zonal advection 56 56 ! ! ttrdh(,,,2) : meridional advection 57 57 ! ! ttrdh(,,,3) : zonal EIV 58 58 ! ! ttrdh(,,,4) : meridional EIV 59 REAL(wp), DIMENSION(jpi,jpj,2) :: & 59 REAL(wp), DIMENSION(jpi,jpj,2) :: & !: 60 60 flxtrd, & !: tracer forcing trends 61 61 bbltrd !: tracer bottom boundary layer trends
Note: See TracChangeset
for help on using the changeset viewer.