Changeset 13334 for NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_ice_interp.F90

Timestamp:

2020-07-22T16:20:32+02:00 (4 years ago)

Author:

jchanut

Message:

finish bypassing ocean/ice initialization with AGRIF, #2222, #2129

File:

• NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_ice_interp.F90 (modified) (3 diffs)

NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_ice_interp.F90

@@ -25 +25 @@
    USE agrif_oce
    USE phycst , ONLY: rt0
+   USE icevar
+   USE sbc_ice, ONLY : tn_ice
    
    IMPLICIT NONE
@@ -30 +32 @@
 
    PUBLIC   agrif_interp_ice   ! called by agrif_user.F90
-   PUBLIC   interp_tra_ice, interp_u_ice, interp_v_ice  ! called by iceistate.F90
+   PUBLIC   agrif_istate_ice   ! called by icerst.F90
 
    !!----------------------------------------------------------------------
@@ -39 +41 @@
 
 CONTAINS
+
+   SUBROUTINE agrif_istate_ice
+      !!-----------------------------------------------------------------------
+      !!                 *** ROUTINE agrif_istate_ice  ***
+      !!
+      !!  ** Method  : Set initial ice fields from parent grid
+      !!
+      !!-----------------------------------------------------------------------
+      IF(lwp) WRITE(numout,*) ' '
+      IF(lwp) WRITE(numout,*) 'Agrif_istate_ice : interp child ice initial state from parent'
+      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~~'
+      IF(lwp) WRITE(numout,*) ' '
+
+      ! Set a_i, v_i, v_s, sv_i, oa_i, a_ip, v_ip, t_su, e_s, e_i:
+      Agrif_SpecialValue    = -9999.
+      Agrif_UseSpecialValue = .TRUE.
+      CALL Agrif_init_variable(tra_iceini_id,procname=interp_tra_ice)
+      !
+      ! Set u_ice, v_ice:
+      use_sign_north = .TRUE.
+      sign_north = -1.
+      Agrif_UseSpecialValue = .TRUE.
+      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 ????
+      ! JC: do we really need the 3 lines below ?
+      CALL ice_var_glo2eqv
+      CALL ice_var_zapsmall
+      CALL ice_var_agg(2)
+
+      ! Melt ponds
+      WHERE( a_i > epsi10 )
+         a_ip_frac(:,:,:) = a_ip(:,:,:) / a_i(:,:,:)
+      ELSEWHERE
+         a_ip_frac(:,:,:) = 0._wp
+      END WHERE
+      WHERE( a_ip > 0._wp )       ! ???????    
+         h_ip(:,:,:) = v_ip(:,:,:) / a_ip(:,:,:)
+      ELSEWHERE
+         h_ip(:,:,:) = 0._wp
+      END WHERE   
+
+      tn_ice(:,:,:) = t_su(:,:,:)
+      t1_ice(:,:,:) = t_i (:,:,1,:) 
+
+   END SUBROUTINE agrif_istate_ice
 
    SUBROUTINE agrif_interp_ice( cd_type, kiter, kitermax )