Changeset 13669 for NEMO/branches/UKMO/NEMO_4.0.1_implicit_ice_drag/src/ICE/iceupdate.F90

Timestamp:

2020-10-23T19:23:00+02:00 (3 years ago)

Author:

dancopsey

Message:

Merge in implicit sea ice drag code from NEMO/branches/UKMO/NEMO_4.0.1_FKOSM_m11715

File:

• NEMO/branches/UKMO/NEMO_4.0.1_implicit_ice_drag/src/ICE/iceupdate.F90 (modified) (4 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_implicit_ice_drag/src/ICE/iceupdate.F90

@@ -26 +26 @@
    USE icectl         ! sea-ice: control prints
    USE bdy_oce , ONLY : ln_bdy
+   USE zdfdrg , ONLY: ln_drgice_imp
    !
    USE in_out_manager ! I/O manager
@@ -323 +324 @@
       REAL(wp) ::   zat_u, zutau_ice, zu_t, zmodt   ! local scalar
       REAL(wp) ::   zat_v, zvtau_ice, zv_t, zrhoco  !   -      -
+      REAL(wp) ::   zflagi                          !   -      -
       !!---------------------------------------------------------------------
       IF( ln_timing )   CALL timing_start('ice_update_tau')
@@ -355 +357 @@
       !
       !                                      !==  every ocean time-step  ==!
+      IF ( ln_drgice_imp ) THEN
+         ! Save drag with right sign to update top drag in the ocean implicit
+         ! friction 
+         rCdU_ice(:,:) = -r1_rau0 * tmod_io(:,:) * at_i(:,:) * tmask(:,:,1) 
+         zflagi = 0._wp
+      ELSE
+         zflagi = 1._wp
+      ENDIF
       !
       DO jj = 2, jpjm1                                !* update the stress WITHOUT an ice-ocean rotation angle
@@ -364 +374 @@
                &     / MAX( 1.0_wp , tmask(ji,jj,1) + tmask(ji  ,jj+1,1) )
             !                                                   ! linearized quadratic drag formulation
-            zutau_ice   = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji+1,jj) ) * ( u_ice(ji,jj) - pu_oce(ji,jj) )
-            zvtau_ice   = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji,jj+1) ) * ( v_ice(ji,jj) - pv_oce(ji,jj) )
+            zutau_ice   = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji+1,jj) ) * ( u_ice(ji,jj) - zflagi * pu_oce(ji,jj) )
+            zvtau_ice   = 0.5_wp * ( tmod_io(ji,jj) + tmod_io(ji,jj+1) ) * ( v_ice(ji,jj) - zflagi * pv_oce(ji,jj) )
             !                                                   ! stresses at the ocean surface
             utau(ji,jj) = ( 1._wp - zat_u ) * utau_oce(ji,jj) + zat_u * zutau_ice