Changeset 9811 for NEMO/releases


Ignore:
Timestamp:
2018-06-19T17:15:12+02:00 (2 years ago)
Author:
clem
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/release-3.6/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

    r8523 r9811  
    5858   LOGICAL :: ll_bdylim3                  ! determine whether ice input is lim2 (F) or lim3 (T) type 
    5959   INTEGER :: jfld_hti, jfld_hts, jfld_ai ! indices of ice thickness, snow thickness and concentration in bf structure 
     60   INTEGER, DIMENSION(jp_bdy) :: jfld_htit, jfld_htst, jfld_ait ! indices of ice thickness, snow thickness and concentration in bf structure 
    6061#endif 
    6162 
     
    380381#if defined key_lim3 
    381382               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) 
     383                  jfld_hti = jfld_htit(ib_bdy) 
     384                  jfld_hts = jfld_htst(ib_bdy) 
     385                  jfld_ai  = jfld_ait(ib_bdy) 
    382386                CALL lim_var_itd ( bf(jfld_hti)%fnow(:,1,1), bf(jfld_hts)%fnow(:,1,1), bf(jfld_ai)%fnow(:,1,1), & 
    383387                                  & dta_bdy(ib_bdy)%ht_i,     dta_bdy(ib_bdy)%ht_s,     dta_bdy(ib_bdy)%a_i     ) 
     
    878882                  dta_bdy(ib_bdy)%ht_s => bf(jfld)%fnow(:,1,:) 
    879883               ELSE ! case input is lim2 type 
    880                   jfld_ai  = jfld + 1 
    881                   jfld_hti = jfld + 2 
    882                   jfld_hts = jfld + 3 
     884                  jfld_ait(ib_bdy)  = jfld + 1 
     885                  jfld_htit(ib_bdy) = jfld + 2 
     886                  jfld_htst(ib_bdy) = jfld + 3 
    883887                  jfld     = jfld + 3 
    884888                  ALLOCATE( dta_bdy(ib_bdy)%a_i (nblen(1),jpl) ) 
Note: See TracChangeset for help on using the changeset viewer.