Changeset 10267 for NEMO/branches/2018/dev_r9947_SI3_advection/src/ICE/icedyn.F90

Timestamp:

2018-10-31T19:09:43+01:00 (5 years ago)

Author:

clem

Message:

improve ice thickness after advection (advection of H using u*a from advection of concentration)

File:

• NEMO/branches/2018/dev_r9947_SI3_advection/src/ICE/icedyn.F90 (modified) (2 diffs)

NEMO/branches/2018/dev_r9947_SI3_advection/src/ICE/icedyn.F90

@@ -71 +71 @@
       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
       !!--------------------------------------------------------------------
       !
@@ -118 +120 @@
 
       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