New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4902 for branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90 – NEMO

Ignore:
Timestamp:
2014-11-27T17:13:38+01:00 (9 years ago)
Author:
cetlod
Message:

2014/dev_CNRS_2014 : Merge in the trunk changes between 4728 and 4879, see ticket #1415

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

    r4901 r4902  
    692692 
    693693      IF( partfun_swi == 0 ) THEN       !--- Linear formulation (Thorndike et al., 1975) 
    694          DO jl = 0, ice_cat_bounds(1,2)       ! only undeformed ice participates 
     694         DO jl = 0, jpl     
    695695            DO jj = 1, jpj  
    696696               DO ji = 1, jpi 
     
    715715            Gsum(:,:,jl) = EXP( -Gsum(:,:,jl) * astari ) * zdummy 
    716716         END DO !jl 
    717          DO jl = 0, ice_cat_bounds(1,2) 
     717         DO jl = 0, jpl 
    718718             athorn(:,:,jl) = Gsum(:,:,jl-1) - Gsum(:,:,jl) 
    719719         END DO 
     
    897897      CALL wrk_alloc( jpi, jpj,             afrft, arft1, arft2, virft, vsrft, esrft, smrft, oirft1, oirft2 ) 
    898898      CALL wrk_alloc( jpi, jpj, jpl,        aicen_init, vicen_init, vsnwn_init, esnwn_init, smv_i_init, oa_i_init ) 
    899       CALL wrk_alloc( jpi, jpj, jkmax,      eirft, erdg1, erdg2, ersw ) 
    900       CALL wrk_alloc( jpi, jpj, jkmax, jpl, eicen_init ) 
     899      CALL wrk_alloc( jpi, jpj, nlay_i+1,      eirft, erdg1, erdg2, ersw ) 
     900      CALL wrk_alloc( jpi, jpj, nlay_i+1, jpl, eicen_init ) 
    901901 
    902902      ! Conservation check 
     
    11931193         !------------------------------------------------------------------------------- 
    11941194         !        jl1 looping 1-jpl 
    1195          DO jl2  = ice_cat_bounds(1,1), ice_cat_bounds(1,2)  
     1195         DO jl2  = 1, jpl  
    11961196            ! over categories to which ridged ice is transferred 
    11971197!CDIR NODEP 
     
    12381238         END DO                 ! jl2 (new ridges)             
    12391239 
    1240          DO jl2 = ice_cat_bounds(1,1), ice_cat_bounds(1,2)  
     1240         DO jl2 = 1, jpl  
    12411241 
    12421242!CDIR NODEP 
     
    13021302      CALL wrk_dealloc( jpi, jpj,             afrft, arft1, arft2, virft, vsrft, esrft, smrft, oirft1, oirft2 ) 
    13031303      CALL wrk_dealloc( jpi, jpj, jpl,        aicen_init, vicen_init, vsnwn_init, esnwn_init, smv_i_init, oa_i_init ) 
    1304       CALL wrk_dealloc( jpi, jpj, jkmax,      eirft, erdg1, erdg2, ersw ) 
    1305       CALL wrk_dealloc( jpi, jpj, jkmax, jpl, eicen_init ) 
     1304      CALL wrk_dealloc( jpi, jpj, nlay_i+1,      eirft, erdg1, erdg2, ersw ) 
     1305      CALL wrk_dealloc( jpi, jpj, nlay_i+1, jpl, eicen_init ) 
    13061306      ! 
    13071307   END SUBROUTINE lim_itd_me_ridgeshift 
Note: See TracChangeset for help on using the changeset viewer.