Changeset 8556


Ignore:
Timestamp:
2017-09-22T11:35:47+02:00 (3 years ago)
Author:
clem
Message:

bug fix

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icedyn_rdgrft.F90

    r8555 r8556  
    501501      !----------------------------------------------------------------- 
    502502 
    503       zaksum(:) = apartf(:,0) 
    504503      zfac = 1._wp / hi_hrft 
     504      zaksum(1:nidx) = apartf(1:nidx,0) 
    505505      ! Transfer function 
    506506      DO jl = 1, jpl !all categories have a specific transfer function 
     
    534534      END WHERE 
    535535       
    536       DO ji = 1, nidx 
    537           
     536      DO ji = 1, nidx   
    538537         ! correction to closing rate and opening if closing rate is excessive 
    539538         !--------------------------------------------------------------------- 
     
    586585      REAL(wp), DIMENSION(jpij) ::   ersw             ! enth of water trapped into ridges 
    587586      REAL(wp), DIMENSION(jpij) ::   zswitch, fvol    ! new ridge volume going to jl2 
    588       ! 
    589       REAL(wp), DIMENSION(jpij,jpl) ::   z1_ai        ! 1 / a 
     587      REAL(wp), DIMENSION(jpij) ::   z1_ai            ! 1 / a 
    590588      ! 
    591589      REAL(wp), DIMENSION(jpij,nlay_i) ::   eirft     ! ice energy of rafting ice 
     
    614612            IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN 
    615613 
    616                z1_ai(ji,jl1) = 1._wp / a_i_2d(ji,jl1) 
     614               z1_ai(ji) = 1._wp / a_i_2d(ji,jl1) 
    617615 
    618616               !-------------------------------------------------------------------- 
     
    628626               ! 2.3) Compute ridging /rafting fractions, make sure afrdg <=1  
    629627               !--------------------------------------------------------------- 
    630                afrdg = airdg1 * z1_ai(ji,jl1) 
    631                afrft = airft1 * z1_ai(ji,jl1) 
     628               afrdg = airdg1 * z1_ai(ji) 
     629               afrft = airft1 * z1_ai(ji) 
    632630 
    633631               ! volume and enthalpy (J/m2, >0) of seawater trapped into ridges 
     
    721719               IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN 
    722720                  ! Compute ridging /rafting fractions 
    723                   afrdg = aridge(ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji,jl1) 
    724                   afrft = araft (ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji,jl1) 
     721                  afrdg = aridge(ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji) 
     722                  afrft = araft (ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji) 
    725723                  ! Compute ridging ice and new ridges for ei 
    726724                  eirdg(ji,jk) = ze_i_2d (ji,jk,jl1) * afrdg + ersw(ji) * r1_nlay_i 
Note: See TracChangeset for help on using the changeset viewer.