- Timestamp:
- 2014-11-27T17:13:38+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90
r4901 r4902 692 692 693 693 IF( partfun_swi == 0 ) THEN !--- Linear formulation (Thorndike et al., 1975) 694 DO jl = 0, ice_cat_bounds(1,2) ! only undeformed ice participates694 DO jl = 0, jpl 695 695 DO jj = 1, jpj 696 696 DO ji = 1, jpi … … 715 715 Gsum(:,:,jl) = EXP( -Gsum(:,:,jl) * astari ) * zdummy 716 716 END DO !jl 717 DO jl = 0, ice_cat_bounds(1,2)717 DO jl = 0, jpl 718 718 athorn(:,:,jl) = Gsum(:,:,jl-1) - Gsum(:,:,jl) 719 719 END DO … … 897 897 CALL wrk_alloc( jpi, jpj, afrft, arft1, arft2, virft, vsrft, esrft, smrft, oirft1, oirft2 ) 898 898 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 ) 901 901 902 902 ! Conservation check … … 1193 1193 !------------------------------------------------------------------------------- 1194 1194 ! jl1 looping 1-jpl 1195 DO jl2 = ice_cat_bounds(1,1), ice_cat_bounds(1,2)1195 DO jl2 = 1, jpl 1196 1196 ! over categories to which ridged ice is transferred 1197 1197 !CDIR NODEP … … 1238 1238 END DO ! jl2 (new ridges) 1239 1239 1240 DO jl2 = ice_cat_bounds(1,1), ice_cat_bounds(1,2)1240 DO jl2 = 1, jpl 1241 1241 1242 1242 !CDIR NODEP … … 1302 1302 CALL wrk_dealloc( jpi, jpj, afrft, arft1, arft2, virft, vsrft, esrft, smrft, oirft1, oirft2 ) 1303 1303 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 ) 1306 1306 ! 1307 1307 END SUBROUTINE lim_itd_me_ridgeshift
Note: See TracChangeset
for help on using the changeset viewer.