Changeset 12377 for NEMO/trunk/src/OCE/TRD/trdmxl.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/OCE/TRD/trdmxl.F90
r11536 r12377 68 68 INTEGER :: ionce, icount 69 69 70 !! * Substitutions 71 # include "do_loop_substitute.h90" 70 72 !!---------------------------------------------------------------------- 71 73 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 86 88 87 89 88 SUBROUTINE trd_tra_mxl( ptrdx, ptrdy, ktrd, kt, p2dt, kmxln )90 SUBROUTINE trd_tra_mxl( ptrdx, ptrdy, ktrd, kt, p2dt, kmxln, Kmm ) 89 91 !!---------------------------------------------------------------------- 90 92 !! *** ROUTINE trd_tra_mng *** … … 98 100 INTEGER , INTENT(in ) :: ktrd ! tracer trend index 99 101 INTEGER , INTENT(in ) :: kt ! time step index 102 INTEGER , INTENT(in ) :: Kmm ! time level index 100 103 REAL(wp) , INTENT(in ) :: p2dt ! time step [s] 101 104 REAL(wp), DIMENSION(:,:) , INTENT(in ) :: kmxln ! number of t-box for the vertical average … … 116 119 ! 117 120 wkx(:,:,:) = 0._wp !== now ML weights for vertical averaging ==! 118 DO jk = 1, jpktrd ! initialize wkx with vertical scale factor in mixed-layer 119 DO jj = 1,jpj 120 DO ji = 1,jpi 121 IF( jk - kmxln(ji,jj) < 0 ) wkx(ji,jj,jk) = e3t_n(ji,jj,jk) * tmask(ji,jj,jk) 122 END DO 123 END DO 124 END DO 121 DO_3D_11_11( 1, jpktrd ) 122 IF( jk - kmxln(ji,jj) < 0 ) wkx(ji,jj,jk) = e3t(ji,jj,jk,Kmm) * tmask(ji,jj,jk) 123 END_3D 125 124 hmxl(:,:) = 0._wp ! NOW mixed-layer depth 126 125 DO jk = 1, jpktrd … … 136 135 tml(:,:) = 0._wp ; sml(:,:) = 0._wp 137 136 DO jk = 1, jpktrd 138 tml(:,:) = tml(:,:) + wkx(:,:,jk) * ts n(:,:,jk,jp_tem)139 sml(:,:) = sml(:,:) + wkx(:,:,jk) * ts n(:,:,jk,jp_sal)137 tml(:,:) = tml(:,:) + wkx(:,:,jk) * ts(:,:,jk,jp_tem,Kmm) 138 sml(:,:) = sml(:,:) + wkx(:,:,jk) * ts(:,:,jk,jp_sal,Kmm) 140 139 END DO 141 140 ! … … 371 370 hmxlbn(:,:) = hmxl(:,:) 372 371 373 IF( ln_ctl ) THEN372 IF( sn_cfctl%l_prtctl ) THEN 374 373 WRITE(numout,*) ' we reach kt == nit000 + 1 = ', nit000+1 375 374 CALL prt_ctl(tab2d_1=tmlbb , clinfo1=' tmlbb - : ', mask1=tmask) … … 380 379 END IF 381 380 382 IF( ( ln_rstart ) .AND. ( kt == nit000 ) .AND. ( ln_ctl )) THEN381 IF( ( ln_rstart ) .AND. ( kt == nit000 ) .AND. sn_cfctl%l_prtctl ) THEN 383 382 IF( ln_trdmxl_instant ) THEN 384 383 WRITE(numout,*) ' restart from kt == nit000 = ', nit000 … … 548 547 hmxlbn (:,:) = hmxl (:,:) 549 548 550 IF( ln_ctl ) THEN549 IF( sn_cfctl%l_prtctl ) THEN 551 550 IF( ln_trdmxl_instant ) THEN 552 551 CALL prt_ctl(tab2d_1=tmlbb , clinfo1=' tmlbb - : ', mask1=tmask) … … 732 731 !!---------------------------------------------------------------------- 733 732 ! 734 REWIND( numnam_ref ) ! Namelist namtrd_mxl in reference namelist : mixed layer trends diagnostic735 733 READ ( numnam_ref, namtrd_mxl, IOSTAT = ios, ERR = 901 ) 736 734 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrd_mxl in reference namelist' ) 737 735 738 REWIND( numnam_cfg ) ! Namelist namtrd_mxl in configuration namelist : mixed layer trends diagnostic739 736 READ ( numnam_cfg, namtrd_mxl, IOSTAT = ios, ERR = 902 ) 740 737 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namtrd_mxl in configuration namelist' )
Note: See TracChangeset
for help on using the changeset viewer.