Ignore:
Timestamp:
2018-06-19T17:09:31+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/trunk/src/OCE/BDY/bdydta.F90

    r9657 r9810  
    5454   INTEGER ::   nice_cat                      ! number of categories in the input file 
    5555   INTEGER ::   jfld_hti, jfld_hts, jfld_ai   ! indices of ice thickness, snow thickness and concentration in bf structure 
     56   INTEGER, DIMENSION(jp_bdy) :: jfld_htit, jfld_htst, jfld_ait 
    5657#endif 
    5758 
     
    351352#if defined key_si3 
    352353               IF( cn_ice(jbdy) /= 'none' .AND. nn_ice_dta(jbdy) == 1 ) THEN 
     354                  jfld_hti = jfld_htit(jbdy) 
     355                  jfld_hts = jfld_htst(jbdy) 
     356                  jfld_ai  = jfld_ait(jbdy) 
    353357                  IF( nice_cat == 1 ) THEN ! case input cat = 1 
    354358                     CALL ice_var_itd ( bf(jfld_hti)%fnow(:,1,1), bf(jfld_hts)%fnow(:,1,1), bf(jfld_ai)%fnow(:,1,1), & 
     
    812816                  dta_bdy(jbdy)%h_s => bf(jfld)%fnow(:,1,:) 
    813817               ELSE                        ! case input cat = 1 OR (/=1 and /=jpl) 
    814                   jfld_ai  = jfld + 1 
    815                   jfld_hti = jfld + 2 
    816                   jfld_hts = jfld + 3 
     818                  jfld_ait(jbdy)  = jfld + 1 
     819                  jfld_htit(jbdy) = jfld + 2 
     820                  jfld_htst(jbdy) = jfld + 3 
    817821                  jfld     = jfld + 3 
    818822                  ALLOCATE( dta_bdy(jbdy)%a_i(nblen(1),jpl) ) 
Note: See TracChangeset for help on using the changeset viewer.