- Timestamp:
- 2020-01-27T15:31:53+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/ISF/isfhdiv.F90
r12077 r12340 24 24 25 25 PUBLIC isf_hdiv 26 !! * Substitutions 27 # include "do_loop_substitute.h90" 26 28 27 29 CONTAINS … … 97 99 ! 98 100 ! update divergence at each level affected by ice shelf top boundary layer 99 DO jj = 1,jpj 100 DO ji = 1,jpi 101 ikt = ktop(ji,jj) 102 ikb = kbot(ji,jj) 103 ! level fully include in the ice shelf boundary layer 104 DO jk = ikt, ikb - 1 105 phdiv(ji,jj,jk) = phdiv(ji,jj,jk) + zhdiv(ji,jj) 106 END DO 107 ! level partially include in ice shelf boundary layer 108 phdiv(ji,jj,ikb) = phdiv(ji,jj,ikb) + zhdiv(ji,jj) * pfrac(ji,jj) 101 DO_2D_11_11 102 ikt = ktop(ji,jj) 103 ikb = kbot(ji,jj) 104 ! level fully include in the ice shelf boundary layer 105 DO jk = ikt, ikb - 1 106 phdiv(ji,jj,jk) = phdiv(ji,jj,jk) + zhdiv(ji,jj) 109 107 END DO 110 END DO 108 ! level partially include in ice shelf boundary layer 109 phdiv(ji,jj,ikb) = phdiv(ji,jj,ikb) + zhdiv(ji,jj) * pfrac(ji,jj) 110 END_2D 111 111 ! 112 112 END SUBROUTINE isf_hdiv_mlt
Note: See TracChangeset
for help on using the changeset viewer.