Changeset 9872
- Timestamp:
- 2018-07-04T10:48:01+02:00 (6 years ago)
- Location:
- NEMO/trunk/src
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/icestp.F90
r9866 r9872 137 137 ! !== AGRIF Parent to Child ==! 138 138 #if defined key_agrif 139 ! ! lim_nbstepranges from 1 to the nb of child ocean steps inside one parent ice step140 IF( .NOT. Agrif_Root() ) lim_nbstep = MOD( lim_nbstep, Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) + 1139 ! ! nbstep_ice ranges from 1 to the nb of child ocean steps inside one parent ice step 140 IF( .NOT. Agrif_Root() ) nbstep_ice = MOD( nbstep_ice, Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) + 1 141 141 ! ! these calls must remain here for restartability purposes 142 142 CALL agrif_interp_ice( 'T' ) -
NEMO/trunk/src/NST/agrif_ice.F90
r9656 r9872 16 16 17 17 INTEGER, PUBLIC :: u_ice_id, v_ice_id, tra_ice_id 18 INTEGER, PUBLIC :: lim_nbstep= 0 ! child time position in sea-ice model18 INTEGER, PUBLIC :: nbstep_ice = 0 ! child time position in sea-ice model 19 19 20 20 !!---------------------------------------------------------------------- -
NEMO/trunk/src/NST/agrif_ice_interp.F90
r9656 r9872 57 57 CASE('U','V') 58 58 IF( PRESENT( kiter ) ) THEN ! interpolation at the child ice sub-time step (only for ice rheology) 59 zbeta = ( REAL( lim_nbstep) - REAL(kitermax - kiter) / REAL(kitermax) ) / &59 zbeta = ( REAL(nbstep_ice) - REAL(kitermax - kiter) / REAL(kitermax) ) / & 60 60 & ( Agrif_Rhot() * REAL(Agrif_Parent(nn_fsbc)) / REAL(nn_fsbc) ) 61 61 ELSE ! interpolation at the child ice time step 62 zbeta = REAL( lim_nbstep) / ( Agrif_Rhot() * REAL(Agrif_Parent(nn_fsbc)) / REAL(nn_fsbc) )62 zbeta = REAL(nbstep_ice) / ( Agrif_Rhot() * REAL(Agrif_Parent(nn_fsbc)) / REAL(nn_fsbc) ) 63 63 ENDIF 64 64 CASE('T') 65 zbeta = REAL( lim_nbstep) / ( Agrif_Rhot() * REAL(Agrif_Parent(nn_fsbc)) / REAL(nn_fsbc) )65 zbeta = REAL(nbstep_ice) / ( Agrif_Rhot() * REAL(Agrif_Parent(nn_fsbc)) / REAL(nn_fsbc) ) 66 66 END SELECT 67 67 ! -
NEMO/trunk/src/NST/agrif_user.F90
r9788 r9872 483 483 CALL ctl_stop('rhot * nn_fsbc(parent) /= N * nn_fsbc(child), therefore nn_fsbc(child) should be set to 1 or nn_fsbc(parent)') 484 484 ENDIF 485 ! First Interpolations (using "after" ice subtime step => lim_nbstep=1)485 ! First Interpolations (using "after" ice subtime step => nbstep_ice=1) 486 486 !---------------------------------------------------------------------- 487 lim_nbstep= ( Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) ! clem: to have calledweight=1 in interp (otherwise the western border of the zoom is wrong)487 nbstep_ice = ( Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) ! clem: to have calledweight=1 in interp (otherwise the western border of the zoom is wrong) 488 488 CALL agrif_interp_ice('U') ! interpolation of ice velocities 489 489 CALL agrif_interp_ice('V') ! interpolation of ice velocities 490 490 CALL agrif_interp_ice('T') ! interpolation of ice tracers 491 lim_nbstep= 0491 nbstep_ice = 0 492 492 493 493 ! -
NEMO/trunk/src/OCE/SBC/sbccpl.F90
r9767 r9872 1352 1352 IF( srcv(jpr_ocx1)%laction ) THEN ! received by sas in case of opa <-> sas coupling 1353 1353 ssu_m(:,:) = frcv(jpr_ocx1)%z3(:,:,1) 1354 ub (:,:,1) = ssu_m(:,:) ! will be used in icestp in the call of lim_sbc_tau1354 ub (:,:,1) = ssu_m(:,:) ! will be used in icestp in the call of ice_forcing_tau 1355 1355 un (:,:,1) = ssu_m(:,:) ! will be used in sbc_cpl_snd if atmosphere coupling 1356 1356 CALL iom_put( 'ssu_m', ssu_m ) … … 1358 1358 IF( srcv(jpr_ocy1)%laction ) THEN 1359 1359 ssv_m(:,:) = frcv(jpr_ocy1)%z3(:,:,1) 1360 vb (:,:,1) = ssv_m(:,:) ! will be used in icestp in the call of lim_sbc_tau1360 vb (:,:,1) = ssv_m(:,:) ! will be used in icestp in the call of ice_forcing_tau 1361 1361 vn (:,:,1) = ssv_m(:,:) ! will be used in sbc_cpl_snd if atmosphere coupling 1362 1362 CALL iom_put( 'ssv_m', ssv_m )
Note: See TracChangeset
for help on using the changeset viewer.