Changeset 9810

Timestamp:

2018-06-19T17:09:31+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/trunk/src/OCE/BDY/bdydta.F90 (modified) (3 diffs)

NEMO/trunk/src/OCE/BDY/bdydta.F90

@@ -54 +54 @@
    INTEGER ::   nice_cat                      ! number of categories in the input file
    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
 #endif
 
@@ -351 +352 @@
 #if defined key_si3
                IF( cn_ice(jbdy) /= 'none' .AND. nn_ice_dta(jbdy) == 1 ) THEN
+                  jfld_hti = jfld_htit(jbdy)
+                  jfld_hts = jfld_htst(jbdy)
+                  jfld_ai  = jfld_ait(jbdy)
                   IF( nice_cat == 1 ) THEN ! case input cat = 1
                      CALL ice_var_itd ( bf(jfld_hti)%fnow(:,1,1), bf(jfld_hts)%fnow(:,1,1), bf(jfld_ai)%fnow(:,1,1), &
@@ -812 +816 @@
                   dta_bdy(jbdy)%h_s => bf(jfld)%fnow(:,1,:)
                ELSE                        ! case input cat = 1 OR (/=1 and /=jpl)
-                  jfld_ai  = jfld + 1
-                  jfld_hti = jfld + 2
-                  jfld_hts = jfld + 3
+                  jfld_ait(jbdy)  = jfld + 1
+                  jfld_htit(jbdy) = jfld + 2
+                  jfld_htst(jbdy) = jfld + 3
                   jfld     = jfld + 3
                   ALLOCATE( dta_bdy(jbdy)%a_i(nblen(1),jpl) )