- Timestamp:
- 2016-11-28T17:04:10+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/FLO/floblk.F90
r5836 r7351 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.