Changeset 9811 for NEMO/releases/release-3.6/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

Timestamp:

2018-06-19T17:15:12+02:00 (6 years ago)

Author:

clem

Message:

debug bdy with sea ice if more than 1 boundary is used. See ticket #2063

File:

• NEMO/releases/release-3.6/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 (modified) (3 diffs)

NEMO/releases/release-3.6/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

@@ -58 +58 @@
    LOGICAL :: ll_bdylim3                  ! determine whether ice input is lim2 (F) or lim3 (T) type
    INTEGER :: jfld_hti, jfld_hts, jfld_ai ! indices of ice thickness, snow thickness and concentration in bf structure
+   INTEGER, DIMENSION(jp_bdy) :: jfld_htit, jfld_htst, jfld_ait ! indices of ice thickness, snow thickness and concentration in bf structure
 #endif
 
@@ -380 +381 @@
 #if defined key_lim3
                IF( .NOT. ll_bdylim3 .AND. cn_ice_lim(ib_bdy) /= 'none' .AND. nn_ice_lim_dta(ib_bdy) == 1 ) THEN ! bdy ice input (case input is lim2 type)
+                  jfld_hti = jfld_htit(ib_bdy)
+                  jfld_hts = jfld_htst(ib_bdy)
+                  jfld_ai  = jfld_ait(ib_bdy)
                 CALL lim_var_itd ( bf(jfld_hti)%fnow(:,1,1), bf(jfld_hts)%fnow(:,1,1), bf(jfld_ai)%fnow(:,1,1), &
                                   & dta_bdy(ib_bdy)%ht_i,     dta_bdy(ib_bdy)%ht_s,     dta_bdy(ib_bdy)%a_i     )
@@ -878 +882 @@
                   dta_bdy(ib_bdy)%ht_s => bf(jfld)%fnow(:,1,:)
                ELSE ! case input is lim2 type
-                  jfld_ai  = jfld + 1
-                  jfld_hti = jfld + 2
-                  jfld_hts = jfld + 3
+                  jfld_ait(ib_bdy)  = jfld + 1
+                  jfld_htit(ib_bdy) = jfld + 2
+                  jfld_htst(ib_bdy) = jfld + 3
                   jfld     = jfld + 3
                   ALLOCATE( dta_bdy(ib_bdy)%a_i (nblen(1),jpl) )