Changeset 11514


Ignore:
Timestamp:
2019-09-09T15:38:47+02:00 (22 months ago)
Author:
clem
Message:

debug ponds with bdy (fraction instead of concentration)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdyice.F90

    r11510 r11514  
    222222               sz_i(ji,jj,:,jl) = sz_i(ib,jb,:,jl) 
    223223               ! 
    224                SELECT CASE( jpbound ) 
    225                   ! 
    226                CASE( 0 )   ! velocity is inward 
    227                   oa_i(ji,jj,jl) = rice_age(jbdy) * a_i(ji,jj,jl)     ! age 
     224               ! ice age and ponds 
     225               IF( jpbound == 0 ) THEN  ! velocity is inward 
     226                  oa_i(ji,jj,jl) = rice_age(jbdy) * a_i(ji,jj,jl) 
    228227                  IF( ln_pnd_CST .OR. ln_pnd_H12 ) THEN 
    229                      a_ip(ji,jj,jl) = rice_apnd(jbdy)                 ! pond concentration 
    230                      v_ip(ji,jj,jl) = rice_apnd(jbdy)*rice_hpnd(jbdy) ! pond volume 
     228                     a_ip_frac(ji,jj,jl) = rice_apnd(jbdy) 
     229                     a_ip     (ji,jj,jl) = rice_apnd(jbdy) * a_i (ji,jj,jl) 
     230                     v_ip     (ji,jj,jl) = rice_hpnd(jbdy) * a_ip(ji,jj,jl) 
    231231                  ENDIF 
    232                CASE( 1 )   ! velocity is outward 
    233                   oa_i(ji,jj,jl) = oa_i(ib,jb,jl)                     ! age 
     232               ELSEIF( jpbound == 1 ) THEN  ! velocity is outward 
     233                  oa_i(ji,jj,jl) = oa_i(ib,jb,jl) 
    234234                  IF( ln_pnd_CST .OR. ln_pnd_H12 ) THEN 
    235                      a_ip(ji,jj,jl) = a_ip(ib,jb,jl)                  ! pond concentration 
    236                      v_ip(ji,jj,jl) = v_ip(ib,jb,jl)                  ! pond volume 
     235                     a_ip_frac(ji,jj,jl) = a_ip_frac(ib,jb,jl) 
     236                     a_ip     (ji,jj,jl) = a_ip     (ib,jb,jl) 
     237                     v_ip     (ji,jj,jl) = v_ip     (ib,jb,jl) 
    237238                  ENDIF 
    238                END SELECT 
     239               ENDIF 
    239240               ! 
    240241               IF( nn_icesal == 1 ) THEN     ! if constant salinity 
     
    272273               t_s (ji,jj,:,jl) = rt0 
    273274               t_i (ji,jj,:,jl) = rt0  
     275 
     276               a_ip_frac(ji,jj,jl) = 0._wp 
     277               h_ip     (ji,jj,jl) = 0._wp 
     278               a_ip     (ji,jj,jl) = 0._wp 
     279               v_ip     (ji,jj,jl) = 0._wp 
    274280                
    275281               IF( nn_icesal == 1 ) THEN     ! if constant salinity 
Note: See TracChangeset for help on using the changeset viewer.