Changeset 12377 for NEMO/trunk/src/TOP/PISCES/SED/trcdmp_sed.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/src/TOP/PISCES/SED/trcdmp_sed.F90
r10225 r12377 35 35 36 36 !! * Substitutions 37 # include " vectopt_loop_substitute.h90"37 # include "do_loop_substitute.h90" 38 38 !!---------------------------------------------------------------------- 39 39 !! NEMO/TOP 3.3 , NEMO Consortium (2010) … … 54 54 55 55 56 SUBROUTINE trc_dmp_sed( kt )56 SUBROUTINE trc_dmp_sed( kt, Kbb, Kmm, Krhs ) 57 57 !!---------------------------------------------------------------------- 58 58 !! *** ROUTINE trc_dmp_sed *** … … 64 64 !! ** Method : Newtonian damping towards trdta computed 65 65 !! and add to the general tracer trends: 66 !! tr n = tra + restotr * (trdta - trb)66 !! tr(Kmm) = tr(Krhs) + restotr * (trdta - tr(Kbb)) 67 67 !! The trend is computed either throughout the water column 68 68 !! (nlmdmptr=0) or in area of weak vertical mixing (nlmdmptr=1) or 69 69 !! below the well mixed layer (nlmdmptr=2) 70 70 !! 71 !! ** Action : - update the tracer trends tr awith the newtonian71 !! ** Action : - update the tracer trends tr(Krhs) with the newtonian 72 72 !! damping trends. 73 73 !! - save the trends ('key_trdmxl_trc') 74 74 !!---------------------------------------------------------------------- 75 INTEGER, INTENT(in) :: kt ! ocean time-step index 75 INTEGER, INTENT(in) :: kt ! ocean time-step index 76 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! time level index 76 77 ! 77 78 INTEGER :: ji, jj, jk, jn, jl, ikt ! dummy loop indices … … 90 91 ! 91 92 jl = n_trc_index(jn) 92 CALL trc_dta( kt, sf_trcdta(jl), rf_trfac(jl), ztrcdta ) ! read tracer data at nit00093 CALL trc_dta( kt, Kmm, sf_trcdta(jl), rf_trfac(jl), ztrcdta ) ! read tracer data at nit000 93 94 ! 94 DO jj = 1, jpj 95 DO ji = 1, jpi ! vector opt. 96 ikt = mbkt(ji,jj) 97 trb(ji,jj,ikt,jn) = ztrcdta(ji,jj,ikt) + ( trb(ji,jj,ikt,jn) - ztrcdta(ji,jj,ikt) ) & 98 & * exp( -restosed(ji,jj,ikt) * dtsed ) 99 END DO 100 END DO 95 DO_2D_11_11 96 ikt = mbkt(ji,jj) 97 tr(ji,jj,ikt,jn,Kbb) = ztrcdta(ji,jj,ikt) + ( tr(ji,jj,ikt,jn,Kbb) - ztrcdta(ji,jj,ikt) ) & 98 & * exp( -restosed(ji,jj,ikt) * dtsed ) 99 END_2D 101 100 ! 102 101 ENDIF … … 106 105 ! 107 106 ! ! print mean trends (used for debugging) 108 IF( ln_ctl) THEN107 IF( sn_cfctl%l_prttrc ) THEN 109 108 WRITE(charout, FMT="('dmp ')") 110 109 CALL prt_ctl_trc_info(charout) 111 CALL prt_ctl_trc( tab4d=tr a, mask=tmask, clinfo=ctrcnm, clinfo2='trd' )110 CALL prt_ctl_trc( tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 112 111 ENDIF 113 112 ! … … 148 147 !!---------------------------------------------------------------------- 149 148 CONTAINS 150 SUBROUTINE trc_dmp_sed( kt )! Empty routine149 SUBROUTINE trc_dmp_sed( kt, Kbb, Kmm, Krhs ) ! Empty routine 151 150 INTEGER, INTENT(in) :: kt 151 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs 152 152 WRITE(*,*) 'trc_dmp_sed: You should not have seen this print! error?', kt 153 153 END SUBROUTINE trc_dmp_sed
Note: See TracChangeset
for help on using the changeset viewer.