- Timestamp:
- 2016-11-21T09:55:07+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/LIM_SRC_2/limdmp_2.F90
r5836 r7277 144 144 145 145 resto_ice(:,:,:) = 0._wp 146 ! Re-calculate the North and South boundary restoring term147 ! because those boundaries may change with the prescribed zoom area.148 146 ! 149 147 irelax = 16 ! width of buffer zone with respect to close boundary … … 156 154 ! REM: if there is no ice in the model and in the data, 157 155 ! no restoring even with non zero resto_ice 158 DO jj = mj0( jpjzoom - 1 + 1), mj1(jpjzoom -1 +irelax)159 zreltim = zdmpmin + zfactor * ( mjg(jj) - jpjzoom + 1)156 DO jj = mj0(1), mj1( irelax) 157 zreltim = zdmpmin + zfactor * mjg(jj) 160 158 resto_ice(:,jj,:) = 1._wp / ( zreltim * 86400._wp ) 161 159 END DO 162 160 163 161 ! North boundary restoring term 164 DO jj = mj0(jpj zoom -1 + jpjglo - irelax), mj1(jpjzoom - 1 +jpjglo)165 zreltim = zdmpmin + zfactor * (jpjglo - ( mjg(jj) - jpjzoom + 1))162 DO jj = mj0(jpjglo - irelax), mj1(jpjglo) 163 zreltim = zdmpmin + zfactor * (jpjglo - mjg(jj)) 166 164 resto_ice(:,jj,:) = 1.e0 / ( zreltim * 86400 ) 167 165 END DO
Note: See TracChangeset
for help on using the changeset viewer.