Changeset 15031 for NEMO/trunk/src/NST/agrif_ice_interp.F90

Timestamp:

2021-06-21T12:05:41+02:00 (3 years ago)

Author:

jchanut

Message:

#2638, add lbclnks for ice initial state if interpolated from parent

File:

• NEMO/trunk/src/NST/agrif_ice_interp.F90 (modified) (4 diffs)

NEMO/trunk/src/NST/agrif_ice_interp.F90

@@ -27 +27 @@
    USE icevar
    USE sbc_ice, ONLY : tn_ice
-   
+   USE lbclnk  
+ 
    IMPLICIT NONE
    PRIVATE
@@ -59 +60 @@
       CALL Agrif_init_variable(tra_iceini_id,procname=interp_tra_ice)
       !
+      CALL lbc_lnk( 'agrif_istate_ice',  a_i,'T',1._wp,  v_i,'T',1._wp, &
+               &         v_s,'T',1._wp, sv_i,'T',1._wp, oa_i,'T',1._wp, &
+               &        a_ip,'T',1._wp, v_ip,'T',1._wp, v_il,'T',1._wp )
+      CALL lbc_lnk( 'agrif_istate_ice', t_su,'T',1._wp )
+      CALL lbc_lnk( 'agrif_istate_ice',  e_s,'T',1._wp )
+      CALL lbc_lnk( 'agrif_istate_ice',  e_i,'T',1._wp )
+      !
       ! Set u_ice, v_ice:
       use_sign_north = .TRUE.
       sign_north = -1.
-      Agrif_UseSpecialValue = .TRUE.
+      ! JC: setting special value to -9999. with north Fold crossing
+      !     does not work probably because of the sign change.
+      !     it's likely that the same issue could occur at boundaries
+      !     but leave it as is for the time being
+      Agrif_SpecialValue = 0._wp
       CALL Agrif_init_variable(u_iceini_id  ,procname=interp_u_ice)
       CALL Agrif_init_variable(v_iceini_id  ,procname=interp_v_ice)
-      Agrif_SpecialValue = 0._wp
       use_sign_north = .FALSE.
       Agrif_UseSpecialValue = .FALSE.
-      ! lbc ????
+      ! 
+      CALL lbc_lnk( 'agrif_istate_ice', u_ice, 'U', -1._wp, v_ice, 'V', -1._wp )
       !
       CALL ice_var_glo2eqv
@@ -139 +151 @@
       !
       IF( before ) THEN  ! parent grid
-         ptab(:,:) = e2u(i1:i2,j1:j2) * u_ice(i1:i2,j1:j2)
+         ptab(i1:i2,j1:j2) = e2u(i1:i2,j1:j2) * u_ice(i1:i2,j1:j2)
          WHERE( umask(i1:i2,j1:j2,1) == 0. )   ptab(i1:i2,j1:j2) = Agrif_SpecialValue
       ELSE               ! child grid
@@ -166 +178 @@
       !
       IF( before ) THEN  ! parent grid
-         ptab(:,:) = e1v(i1:i2,j1:j2) * v_ice(i1:i2,j1:j2)
+         ptab(i1:i2,j1:j2) = e1v(i1:i2,j1:j2) * v_ice(i1:i2,j1:j2)
          WHERE( vmask(i1:i2,j1:j2,1) == 0. )   ptab(i1:i2,j1:j2) = Agrif_SpecialValue
       ELSE               ! child grid