Changeset 5682 for branches/2015/dev_r5072_UKMO2_OBS_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90
- Timestamp:
- 2015-08-12T17:46:45+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5072_UKMO2_OBS_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90
r4990 r5682 27 27 USE dom_oce ! ocean space and time domain variables 28 28 USE sbc_oce ! surface boundary condition: ocean 29 USE sbcrnf ! river runoffs 30 USE sbcisf ! ice shelf melting/freezing 29 31 USE zdf_oce ! ocean vertical mixing 30 32 USE domvvl ! variable volume … … 45 47 USE timing ! Timing 46 48 #if defined key_agrif 47 USE agrif_opa_update48 49 USE agrif_opa_interp 49 50 #endif … … 109 110 ! Update after tracer on domain lateral boundaries 110 111 ! 112 #if defined key_agrif 113 CALL Agrif_tra ! AGRIF zoom boundaries 114 #endif 115 ! 111 116 CALL lbc_lnk( tsa(:,:,:,jp_tem), 'T', 1._wp ) ! local domain boundaries (T-point, unchanged sign) 112 117 CALL lbc_lnk( tsa(:,:,:,jp_sal), 'T', 1._wp ) … … 114 119 #if defined key_bdy 115 120 IF( lk_bdy ) CALL bdy_tra( kt ) ! BDY open boundaries 116 #endif117 #if defined key_agrif118 CALL Agrif_tra ! AGRIF zoom boundaries119 121 #endif 120 122 … … 143 145 ELSE ! Leap-Frog + Asselin filter time stepping 144 146 ! 145 IF( lk_vvl ) THEN ; CALL tra_nxt_vvl( kt, nit000, 'TRA', tsb, tsn, tsa, jpts ) ! variable volume level (vvl) 146 ELSE ; CALL tra_nxt_fix( kt, nit000, 'TRA', tsb, tsn, tsa, jpts ) ! fixed volume level 147 IF( lk_vvl ) THEN ; CALL tra_nxt_vvl( kt, nit000, rdttra, 'TRA', tsb, tsn, tsa, & 148 & sbc_tsc, sbc_tsc_b, jpts ) ! variable volume level (vvl) 149 ELSE ; CALL tra_nxt_fix( kt, nit000, 'TRA', tsb, tsn, tsa, jpts ) ! fixed volume level 147 150 ENDIF 148 ENDIF 149 ! 150 #if defined key_agrif 151 ! Update tracer at AGRIF zoom boundaries 152 IF( .NOT.Agrif_Root() ) CALL Agrif_Update_Tra( kt ) ! children only 153 #endif 151 ENDIF 154 152 ! 155 153 ! trends computation … … 241 239 242 240 243 SUBROUTINE tra_nxt_vvl( kt, kit000, cdtype, ptb, ptn, pta, kjpt )241 SUBROUTINE tra_nxt_vvl( kt, kit000, p2dt, cdtype, ptb, ptn, pta, psbc_tc, psbc_tc_b, kjpt ) 244 242 !!---------------------------------------------------------------------- 245 243 !! *** ROUTINE tra_nxt_vvl *** … … 265 263 !! - (ta,sa) time averaged (t,s) (ln_dynhpg_imp = T) 266 264 !!---------------------------------------------------------------------- 267 INTEGER , INTENT(in ) :: kt ! ocean time-step index 268 INTEGER , INTENT(in ) :: kit000 ! first time step index 269 CHARACTER(len=3), INTENT(in ) :: cdtype ! =TRA or TRC (tracer indicator) 270 INTEGER , INTENT(in ) :: kjpt ! number of tracers 271 REAL(wp) , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt) :: ptb ! before tracer fields 272 REAL(wp) , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt) :: ptn ! now tracer fields 273 REAL(wp) , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt) :: pta ! tracer trend 265 INTEGER , INTENT(in ) :: kt ! ocean time-step index 266 INTEGER , INTENT(in ) :: kit000 ! first time step index 267 REAL(wp) , INTENT(in ), DIMENSION(jpk) :: p2dt ! time-step 268 CHARACTER(len=3), INTENT(in ) :: cdtype ! =TRA or TRC (tracer indicator) 269 INTEGER , INTENT(in ) :: kjpt ! number of tracers 270 REAL(wp) , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt) :: ptb ! before tracer fields 271 REAL(wp) , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt) :: ptn ! now tracer fields 272 REAL(wp) , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt) :: pta ! tracer trend 273 REAL(wp) , INTENT(in ), DIMENSION(jpi,jpj,kjpt) :: psbc_tc ! surface tracer content 274 REAL(wp) , INTENT(in ), DIMENSION(jpi,jpj,kjpt) :: psbc_tc_b ! before surface tracer content 275 274 276 !! 275 LOGICAL :: ll_tra , ll_tra_hpg, ll_traqsr! local logical277 LOGICAL :: ll_tra_hpg, ll_traqsr, ll_rnf, ll_isf ! local logical 276 278 INTEGER :: ji, jj, jk, jn ! dummy loop indices 277 279 REAL(wp) :: zfact1, ztc_a , ztc_n , ztc_b , ztc_f , ztc_d ! local scalar … … 286 288 ! 287 289 IF( cdtype == 'TRA' ) THEN 288 ll_tra = .TRUE. ! active tracers case289 290 ll_tra_hpg = ln_dynhpg_imp ! active tracers case and semi-implicit hpg 290 291 ll_traqsr = ln_traqsr ! active tracers case and solar penetration 292 ll_rnf = ln_rnf ! active tracers case and river runoffs 293 IF (nn_isf .GE. 1) THEN 294 ll_isf = .TRUE. ! active tracers case and ice shelf melting/freezing 295 ELSE 296 ll_isf = .FALSE. 297 END IF 291 298 ELSE 292 ll_tra = .FALSE. ! passive tracers case293 299 ll_tra_hpg = .FALSE. ! passive tracers case or NO semi-implicit hpg 294 300 ll_traqsr = .FALSE. ! active tracers case and NO solar penetration 301 ll_rnf = .FALSE. ! passive tracers or NO river runoffs 302 ll_isf = .FALSE. ! passive tracers or NO ice shelf melting/freezing 295 303 ENDIF 296 304 ! 297 305 DO jn = 1, kjpt 298 306 DO jk = 1, jpkm1 299 zfact1 = atfp * rdttra(jk)307 zfact1 = atfp * p2dt(jk) 300 308 zfact2 = zfact1 / rau0 301 309 DO jj = 1, jpj … … 315 323 ztc_f = ztc_n + atfp * ztc_d 316 324 ! 317 IF( ll_tra .AND. jk == 1 ) THEN ! first level only for T & S 318 ze3t_f = ze3t_f - zfact2 * ( emp_b(ji,jj) - emp(ji,jj) ) 319 ztc_f = ztc_f - zfact1 * ( sbc_tsc(ji,jj,jn) - sbc_tsc_b(ji,jj,jn) ) 325 IF( jk == mikt(ji,jj) ) THEN ! first level 326 ze3t_f = ze3t_f - zfact2 * ( (emp_b(ji,jj) - emp(ji,jj) ) & 327 & - (rnf_b(ji,jj) - rnf(ji,jj) ) & 328 & + (fwfisf_b(ji,jj) - fwfisf(ji,jj)) ) 329 ztc_f = ztc_f - zfact1 * ( psbc_tc(ji,jj,jn) - psbc_tc_b(ji,jj,jn) ) 320 330 ENDIF 321 IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr ) & ! solar penetration (temperature only) 331 332 ! solar penetration (temperature only) 333 IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr ) & 322 334 & ztc_f = ztc_f - zfact1 * ( qsr_hc(ji,jj,jk) - qsr_hc_b(ji,jj,jk) ) 323 335 324 ze3t_f = 1.e0 / ze3t_f 325 ptb(ji,jj,jk,jn) = ztc_f * ze3t_f ! ptb <-- ptn filtered 326 ptn(ji,jj,jk,jn) = pta(ji,jj,jk,jn) ! ptn <-- pta 327 ! 328 IF( ll_tra_hpg ) THEN ! semi-implicit hpg (T & S only) 329 ze3t_d = 1.e0 / ( ze3t_n + rbcp * ze3t_d ) 330 pta(ji,jj,jk,jn) = ze3t_d * ( ztc_n + rbcp * ztc_d ) ! ta <-- Brown & Campana average 331 ENDIF 336 ! river runoff 337 IF( ll_rnf .AND. jk <= nk_rnf(ji,jj) ) & 338 & ztc_f = ztc_f - zfact1 * ( rnf_tsc(ji,jj,jn) - rnf_tsc_b(ji,jj,jn) ) & 339 & * fse3t_n(ji,jj,jk) / h_rnf(ji,jj) 340 341 ! ice shelf 342 IF( ll_isf ) THEN 343 ! level fully include in the Losch_2008 ice shelf boundary layer 344 IF ( jk >= misfkt(ji,jj) .AND. jk < misfkb(ji,jj) ) & 345 ztc_f = ztc_f - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) ) & 346 & * fse3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj) 347 ! level partially include in Losch_2008 ice shelf boundary layer 348 IF ( jk == misfkb(ji,jj) ) & 349 ztc_f = ztc_f - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) ) & 350 & * fse3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj) * ralpha(ji,jj) 351 END IF 352 353 ze3t_f = 1.e0 / ze3t_f 354 ptb(ji,jj,jk,jn) = ztc_f * ze3t_f ! ptb <-- ptn filtered 355 ptn(ji,jj,jk,jn) = pta(ji,jj,jk,jn) ! ptn <-- pta 356 ! 357 IF( ll_tra_hpg ) THEN ! semi-implicit hpg (T & S only) 358 ze3t_d = 1.e0 / ( ze3t_n + rbcp * ze3t_d ) 359 pta(ji,jj,jk,jn) = ze3t_d * ( ztc_n + rbcp * ztc_d ) ! ta <-- Brown & Campana average 360 ENDIF 332 361 END DO 333 362 END DO
Note: See TracChangeset
for help on using the changeset viewer.
