- Timestamp:
- 2018-10-29T15:55:40+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90
r10251 r10253 28 28 USE sbc_oce ! surface boundary condition: ocean 29 29 USE sbcrnf ! river runoffs 30 USE sbcisf ! ice shelf melting/freezing 30 31 USE zdf_oce ! ocean vertical mixing 31 32 USE domvvl ! variable volume … … 46 47 USE timing ! Timing 47 48 #if defined key_agrif 48 USE agrif_opa_update49 49 USE agrif_opa_interp 50 50 #endif … … 58 58 59 59 REAL(wp) :: rbcp ! Brown & Campana parameters for semi-implicit hpg 60 INTEGER :: warn_1, warn_2 ! indicators for warning statement 60 61 61 62 !! * Substitutions … … 93 94 INTEGER, INTENT(in) :: kt ! ocean time-step index 94 95 !! 95 INTEGER :: jk, jn ! dummy loop indices96 REAL(wp) :: zfact ! local scalars96 INTEGER :: jk, jn, ji, jj ! dummy loop indices 97 REAL(wp) :: zfact, zfreeze ! local scalars 97 98 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdt, ztrds 98 99 !!---------------------------------------------------------------------- … … 110 111 ! Update after tracer on domain lateral boundaries 111 112 ! 113 #if defined key_agrif 114 CALL Agrif_tra ! AGRIF zoom boundaries 115 #endif 116 ! 112 117 CALL lbc_lnk( tsa(:,:,:,jp_tem), 'T', 1._wp ) ! local domain boundaries (T-point, unchanged sign) 113 118 CALL lbc_lnk( tsa(:,:,:,jp_sal), 'T', 1._wp ) … … 115 120 #if defined key_bdy 116 121 IF( lk_bdy ) CALL bdy_tra( kt ) ! BDY open boundaries 117 #endif118 #if defined key_agrif119 CALL Agrif_tra ! AGRIF zoom boundaries120 122 #endif 121 123 … … 124 126 ELSEIF( kt <= nit000 + 1 ) THEN ; r2dtra(:) = 2._wp* rdttra(:) ! at nit000 or nit000+1 (Leapfrog) 125 127 ENDIF 128 129 #if ( ! defined key_lim3 && ! defined key_lim2 && ! key_cice ) 130 IF ( kt == nit000 ) warn_1=0 131 warn_2=0 132 DO jk = 1, jpkm1 133 DO jj = 1, jpj 134 DO ji = 1, jpi 135 IF ( tsa(ji,jj,jk,jp_tem) .lt. 0.0 ) THEN 136 ! calculate freezing point 137 zfreeze = ( -0.0575_wp + 1.710523E-3 * Sqrt(Abs(tsn(ji,jj,jk,jp_sal))) & 138 - 2.154996E-4 * tsn(ji,jj,jk,jp_sal) ) * tsn(ji,jj,jk,jp_sal) - 7.53E-4 * ( 10.0_wp + fsdept(ji,jj,jk) ) 139 IF ( tsa(ji,jj,jk,jp_tem) .lt. zfreeze ) THEN 140 tsa(ji,jj,jk,jp_tem)=zfreeze 141 warn_2=1 142 ENDIF 143 ENDIF 144 END DO 145 END DO 146 END DO 147 CALL mpp_max(warn_1) 148 CALL mpp_max(warn_2) 149 IF ( (warn_1 == 0) .and. (warn_2 /= 0) ) THEN 150 IF(lwp) THEN 151 CALL ctl_warn( ' Temperatures dropping below freezing point, ', & 152 & ' being forced to freezing point, no longer conservative' ) 153 ENDIF 154 warn_1=1 155 ENDIF 156 #endif 126 157 127 158 ! trends computation initialisation … … 148 179 ELSE ; CALL tra_nxt_fix( kt, nit000, 'TRA', tsb, tsn, tsa, jpts ) ! fixed volume level 149 180 ENDIF 150 ENDIF 151 ! 152 #if defined key_agrif 153 ! Update tracer at AGRIF zoom boundaries 154 IF( .NOT.Agrif_Root() ) CALL Agrif_Update_Tra( kt ) ! children only 155 #endif 156 ! 157 ! trends computation 181 ENDIF 182 ! 183 ! trends computation 158 184 IF( l_trdtra ) THEN ! trend of the Asselin filter (tb filtered - tb)/dt 159 185 DO jk = 1, jpkm1 … … 279 305 280 306 !! 281 LOGICAL :: ll_tra_hpg, ll_traqsr, ll_rnf ! local logical307 LOGICAL :: ll_tra_hpg, ll_traqsr, ll_rnf, ll_isf ! local logical 282 308 INTEGER :: ji, jj, jk, jn ! dummy loop indices 283 309 REAL(wp) :: zfact1, ztc_a , ztc_n , ztc_b , ztc_f , ztc_d ! local scalar … … 295 321 ll_traqsr = ln_traqsr ! active tracers case and solar penetration 296 322 ll_rnf = ln_rnf ! active tracers case and river runoffs 323 IF (nn_isf .GE. 1) THEN 324 ll_isf = .TRUE. ! active tracers case and ice shelf melting/freezing 325 ELSE 326 ll_isf = .FALSE. 327 END IF 297 328 ELSE 298 329 ll_tra_hpg = .FALSE. ! passive tracers case or NO semi-implicit hpg 299 330 ll_traqsr = .FALSE. ! active tracers case and NO solar penetration 300 331 ll_rnf = .FALSE. ! passive tracers or NO river runoffs 332 ll_isf = .FALSE. ! passive tracers or NO ice shelf melting/freezing 301 333 ENDIF 302 334 ! … … 321 353 ztc_f = ztc_n + atfp * ztc_d 322 354 ! 323 IF( jk == 1 ) THEN ! first level 324 ze3t_f = ze3t_f - zfact2 * ( emp_b(ji,jj) - emp(ji,jj) + rnf(ji,jj) - rnf_b(ji,jj) ) 355 IF( jk == mikt(ji,jj) ) THEN ! first level 356 ze3t_f = ze3t_f - zfact2 * ( (emp_b(ji,jj) - emp(ji,jj) ) & 357 & - (rnf_b(ji,jj) - rnf(ji,jj) ) & 358 & + (fwfisf_b(ji,jj) - fwfisf(ji,jj)) ) 325 359 ztc_f = ztc_f - zfact1 * ( psbc_tc(ji,jj,jn) - psbc_tc_b(ji,jj,jn) ) 326 360 ENDIF 327 361 328 IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr ) & ! solar penetration (temperature only) 362 ! solar penetration (temperature only) 363 IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr ) & 329 364 & ztc_f = ztc_f - zfact1 * ( qsr_hc(ji,jj,jk) - qsr_hc_b(ji,jj,jk) ) 330 365 331 IF( ll_rnf .AND. jk <= nk_rnf(ji,jj) ) & ! river runoffs 366 ! river runoff 367 IF( ll_rnf .AND. jk <= nk_rnf(ji,jj) ) & 332 368 & ztc_f = ztc_f - zfact1 * ( rnf_tsc(ji,jj,jn) - rnf_tsc_b(ji,jj,jn) ) & 333 369 & * fse3t_n(ji,jj,jk) / h_rnf(ji,jj) 370 371 ! ice shelf 372 IF( ll_isf ) THEN 373 ! level fully include in the Losch_2008 ice shelf boundary layer 374 IF ( jk >= misfkt(ji,jj) .AND. jk < misfkb(ji,jj) ) & 375 ztc_f = ztc_f - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) ) & 376 & * fse3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj) 377 ! level partially include in Losch_2008 ice shelf boundary layer 378 IF ( jk == misfkb(ji,jj) ) & 379 ztc_f = ztc_f - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) ) & 380 & * fse3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj) * ralpha(ji,jj) 381 END IF 334 382 335 383 ze3t_f = 1.e0 / ze3t_f
Note: See TracChangeset
for help on using the changeset viewer.