- Timestamp:
- 2015-12-16T10:25:22+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC/FLO/floblk.F90
r5836 r6060 23 23 PUBLIC flo_blk ! routine called by floats.F90 24 24 25 !! * Substitutions26 # include "domzgr_substitute.h90"27 25 !!---------------------------------------------------------------------- 28 26 !! NEMO/OPA 3.3 , NEMO Consortium (2010) … … 118 116 119 117 ! compute the transport across the mesh where the float is. 120 !!bug (gm) change e3t into fse3. but never checked121 zsurfx(1) = e2u(iiloc(jfl)-1,ijloc(jfl) ) * fse3u(iiloc(jfl)-1,ijloc(jfl) ,-ikl(jfl))122 zsurfx(2) = e2u(iiloc(jfl) ,ijloc(jfl) ) * fse3u(iiloc(jfl) ,ijloc(jfl) ,-ikl(jfl))123 zsurfy(1) = e1v(iiloc(jfl) ,ijloc(jfl)-1) * fse3v(iiloc(jfl) ,ijloc(jfl)-1,-ikl(jfl))124 zsurfy(2) = e1v(iiloc(jfl) ,ijloc(jfl) ) * fse3v(iiloc(jfl) ,ijloc(jfl) ,-ikl(jfl))118 !!bug (gm) change e3t into e3. but never checked 119 zsurfx(1) = e2u(iiloc(jfl)-1,ijloc(jfl) ) * e3u_n(iiloc(jfl)-1,ijloc(jfl) ,-ikl(jfl)) 120 zsurfx(2) = e2u(iiloc(jfl) ,ijloc(jfl) ) * e3u_n(iiloc(jfl) ,ijloc(jfl) ,-ikl(jfl)) 121 zsurfy(1) = e1v(iiloc(jfl) ,ijloc(jfl)-1) * e3v_n(iiloc(jfl) ,ijloc(jfl)-1,-ikl(jfl)) 122 zsurfy(2) = e1v(iiloc(jfl) ,ijloc(jfl) ) * e3v_n(iiloc(jfl) ,ijloc(jfl) ,-ikl(jfl)) 125 123 126 124 ! for a isobar float zsurfz is put to zero. The vertical velocity will be zero too. 127 125 zsurfz = e1e2t(iiloc(jfl),ijloc(jfl)) 128 zvol = zsurfz * fse3t(iiloc(jfl),ijloc(jfl),-ikl(jfl))126 zvol = zsurfz * e3t_n(iiloc(jfl),ijloc(jfl),-ikl(jfl)) 129 127 130 128 !
Note: See TracChangeset
for help on using the changeset viewer.