- Timestamp:
- 2019-08-22T11:49:08+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/dev_r10037_GPU/src/ICE/icedyn_adv_umx.F90
r10843 r11467 54 54 !! Reference : Leonard, B.P., 1991, Comput. Methods Appl. Mech. Eng., 88, 17-74. 55 55 !!---------------------------------------------------------------------- 56 USE scoce, ONLY : zudy => scr2D1, zvdx => scr2D2, zcu_box => scr2D3, zcv_box => scr2D4 56 57 INTEGER , INTENT(in ) :: k_order ! order of the scheme (1-5 or 20) 57 58 INTEGER , INTENT(in ) :: kt ! time step … … 72 73 INTEGER :: initad ! number of sub-timestep for the advection 73 74 REAL(wp) :: zcfl , zusnit, zdt ! - - 74 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zudy, zvdx, zcu_box, zcv_box75 75 !!---------------------------------------------------------------------- 76 76 ! 77 77 IF( kt == nit000 .AND. lwp ) WRITE(numout,*) '-- ice_dyn_adv_umx: Ultimate-Macho advection scheme' 78 78 ! 79 ALLOCATE( zudy(jpi,jpj) , zvdx(jpi,jpj) , zcu_box(jpi,jpj) , zcv_box(jpi,jpj) )80 79 ! 81 80 ! --- If ice drift field is too fast, use an appropriate time step for advection (CFL test for stability) --- ! … … 133 132 END DO 134 133 ! 135 DEALLOCATE( zudy, zvdx, zcu_box, zcv_box )136 !137 134 END SUBROUTINE ice_dyn_adv_umx 138 135 … … 151 148 !! ** Action : - pt the after advective tracer 152 149 !!---------------------------------------------------------------------- 150 USE scoce, ONLY : zfu_ups => scr2D5, zfu_ho => scr2D6, zt_u => scr2D7, zt_ups => scr2D8, & 151 zfv_ups => scr2D9, zfv_ho => scr2D10, zt_v => scr2D11, ztrd => scr2D12 153 152 INTEGER , INTENT(in ) :: k_order ! order of the ULTIMATE scheme 154 153 INTEGER , INTENT(in ) :: kt ! number of iteration … … 160 159 INTEGER :: ji, jj ! dummy loop indices 161 160 REAL(wp) :: ztra ! local scalar 162 REAL(wp), DIMENSION(jpi,jpj) :: zfu_ups, zfu_ho, zt_u, zt_ups163 REAL(wp), DIMENSION(jpi,jpj) :: zfv_ups, zfv_ho, zt_v, ztrd164 161 !!---------------------------------------------------------------------- 165 162 ! … … 245 242 !! Reference : Leonard, B.P., 1991, Comput. Methods Appl. Mech. Eng., 88, 17-74. 246 243 !!---------------------------------------------------------------------- 244 USE scoce, ONLY : zzt => scr2D13 247 245 INTEGER , INTENT(in ) :: k_order ! order of the ULTIMATE scheme 248 246 INTEGER , INTENT(in ) :: kt ! number of iteration … … 255 253 INTEGER :: ji, jj ! dummy loop indices 256 254 REAL(wp) :: zc_box ! - - 257 REAL(wp), DIMENSION(jpi,jpj) :: zzt258 255 !!---------------------------------------------------------------------- 259 256 ! … … 310 307 !! Reference : Leonard, B.P., 1991, Comput. Methods Appl. Mech. Eng., 88, 17-74. 311 308 !!---------------------------------------------------------------------- 309 USE scoce, ONLY : ztu1 => scr2D14, ztu2 => scr2D15, ztu3 => scr2D16, ztu4 => scr2D17 312 310 INTEGER , INTENT(in ) :: k_order ! ocean time-step index 313 311 REAL(wp) , INTENT(in ) :: pdt ! tracer time-step … … 318 316 INTEGER :: ji, jj ! dummy loop indices 319 317 REAL(wp) :: zcu, zdx2, zdx4 ! - - 320 REAL(wp), DIMENSION(jpi,jpj) :: ztu1, ztu2, ztu3, ztu4321 318 !!---------------------------------------------------------------------- 322 319 ! … … 428 425 !! Reference : Leonard, B.P., 1991, Comput. Methods Appl. Mech. Eng., 88, 17-74. 429 426 !!---------------------------------------------------------------------- 427 USE scoce, ONLY : ztv1 => scr2D14, ztv2 => scr2D15, ztv3 => scr2D16, ztv4 => scr2D17 430 428 INTEGER , INTENT(in ) :: k_order ! ocean time-step index 431 429 REAL(wp) , INTENT(in ) :: pdt ! tracer time-step … … 436 434 INTEGER :: ji, jj ! dummy loop indices 437 435 REAL(wp) :: zcv, zdy2, zdy4 ! - - 438 REAL(wp), DIMENSION(jpi,jpj) :: ztv1, ztv2, ztv3, ztv4439 436 !!---------------------------------------------------------------------- 440 437 ! … … 546 543 !! in-space based differencing for fluid 547 544 !!---------------------------------------------------------------------- 545 USE scoce, ONLY : zbetup => scr2D14, zbetdo => scr2D15, zbup => scr2D16, & 546 zbdo => scr2D17, zmsk => scr2D18, zdiv => scr2D19 548 547 REAL(wp) , INTENT(in ) :: pdt ! tracer time-step 549 548 REAL(wp), DIMENSION (jpi,jpj), INTENT(in ) :: pbef, paft ! before & after field … … 554 553 REAL(wp) :: zpos, zneg, zbt, za, zb, zc, zbig, zsml, z1_dt ! local scalars 555 554 REAL(wp) :: zau, zbu, zcu, zav, zbv, zcv, zup, zdo ! - - 556 REAL(wp), DIMENSION(jpi,jpj) :: zbetup, zbetdo, zbup, zbdo, zmsk, zdiv557 555 !!---------------------------------------------------------------------- 558 556 !
Note: See TracChangeset
for help on using the changeset viewer.