# Changeset 10267

Ignore:
Timestamp:
2018-10-31T19:09:43+01:00 (2 years ago)
Message:

Location:
Files:
2 edited

### Legend:

Unmodified
Removed

 r9604 INTEGER, INTENT(in) ::   kt     ! ice time step !! INTEGER ::   ji, jj, jl         ! dummy loop indices REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zhmax INTEGER  ::   ji, jj, jl        ! dummy loop indices REAL(wp) ::   zcoefu, zcoefv REAL(wp),              DIMENSION(jpi,jpj,jpl) ::   zhmax REAL(wp), ALLOCATABLE, DIMENSION(:,:)         ::   zdivu_i !!-------------------------------------------------------------------- ! CASE ( np_dynADV )           !==  pure advection ==!   (prescribed velocities) ALLOCATE( zdivu_i(jpi,jpj) ) u_ice(:,:) = rn_uice * umask(:,:,1) v_ice(:,:) = rn_vice * vmask(:,:,1) !!CALL RANDOM_NUMBER(u_ice(:,:)) !!CALL RANDOM_NUMBER(v_ice(:,:)) !CALL RANDOM_NUMBER(u_ice(:,:)) ; u_ice(:,:) = u_ice(:,:) * 0.1 + rn_uice * 0.9 * umask(:,:,1) !CALL RANDOM_NUMBER(v_ice(:,:)) ; v_ice(:,:) = v_ice(:,:) * 0.1 + rn_vice * 0.9 * vmask(:,:,1) ! --- monotonicity test from Lipscomb et al 2004 --- ! ! CFL = 0.5 at a distance from the bound of 1/6 of the basin length ! Then for dx = 2m and dt = 1s => rn_uice = u (1/6th) = 1m/s !DO jj = 1, jpj !   DO ji = 1, jpi !      zcoefu = ( REAL(jpiglo+1)*0.5 - REAL(ji+nimpp-1) ) / ( REAL(jpiglo+1)*0.5 - 1. ) !      zcoefv = ( REAL(jpjglo+1)*0.5 - REAL(jj+njmpp-1) ) / ( REAL(jpjglo+1)*0.5 - 1. ) !      u_ice(ji,jj) = rn_uice * 1.5 * SIGN( 1., zcoefu ) * ABS( zcoefu ) * umask(ji,jj,1) !      v_ice(ji,jj) = rn_vice * 1.5 * SIGN( 1., zcoefv ) * ABS( zcoefv ) * vmask(ji,jj,1) !   END DO !END DO ! --- CALL ice_dyn_adv   ( kt )                            ! -- advection of ice ! diagnostics: divergence at T points DO jj = 2, jpjm1 DO ji = 2, jpim1 zdivu_i(ji,jj) = ( e2u(ji,jj) * u_ice(ji,jj) - e2u(ji-1,jj) * u_ice(ji-1,jj)   & &             + e1v(ji,jj) * v_ice(ji,jj) - e1v(ji,jj-1) * v_ice(ji,jj-1) ) * r1_e1e2t(ji,jj) END DO END DO CALL lbc_lnk( zdivu_i, 'T', 1. ) IF( iom_use('icediv') )   CALL iom_put( "icediv" , zdivu_i(:,:) ) DEALLOCATE( zdivu_i ) END SELECT