Changeset 2977 for branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/OPA_SRC/LBC
- Timestamp:
- 2011-10-22T15:46:41+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/OPA_SRC/LBC/cla.F90
r2715 r2977 387 387 DO ji = mi0(161), mi1(161) 388 388 DO jk = 1, jpkm1 ! surf inflow + reciculation (from Gulf of Aden) 389 t a(ji,jj,jk) = ta(ji,jj,jk) - hdiv_161_88_kt(jk) * tn(ji,jj,jk)390 sa(ji,jj,jk) = sa(ji,jj,jk) - hdiv_161_88_kt(jk) * sn(ji,jj,jk)389 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - hdiv_161_88_kt(jk) * tsn(ji,jj,jk,jp_tem) 390 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - hdiv_161_88_kt(jk) * tsn(ji,jj,jk,jp_sal) 391 391 END DO 392 392 END DO … … 395 395 DO ji = mi0(161), mi1(161) 396 396 jk = 21 ! deep outflow + recirulation (combined flux) 397 t a(ji,jj,jk) = ta(ji,jj,jk) + hdiv_161_88(20) * tn(ji ,jj+1,20) & ! upper recirculation from Gulf of Aden398 & + hdiv_161_88(21) * t n(ji ,jj+1,21) & ! deep recirculation from Gulf of Aden399 & + hdiv_160_89(16) * t n(ji-1,jj+2,16) ! deep inflow from Red sea400 sa(ji,jj,jk) = sa(ji,jj,jk) + hdiv_161_88(20) * sn(ji ,jj+1,20) &401 & + hdiv_161_88(21) * sn(ji ,jj+1,21) &402 & + hdiv_160_89(16) * sn(ji-1,jj+2,16)397 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) + hdiv_161_88(20) * tsn(ji ,jj+1,20,jp_tem) & ! upper recirculation from Gulf of Aden 398 & + hdiv_161_88(21) * tsn(ji ,jj+1,21,jp_tem) & ! deep recirculation from Gulf of Aden 399 & + hdiv_160_89(16) * tsn(ji-1,jj+2,16,jp_tem) ! deep inflow from Red sea 400 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) + hdiv_161_88(20) * tsn(ji ,jj+1,20,jp_sal) & 401 & + hdiv_161_88(21) * tsn(ji ,jj+1,21,jp_sal) & 402 & + hdiv_160_89(16) * tsn(ji-1,jj+2,16,jp_sal) 403 403 END DO 404 404 END DO … … 406 406 DO ji = mi0(160), mi1(160) 407 407 DO jk = 1, 14 ! surface inflow (from Gulf of Aden) 408 t a(ji,jj,jk) = ta(ji,jj,jk) - hdiv_160_89_kt(jk) * tn(ji+1,jj-1,jk)409 sa(ji,jj,jk) = sa(ji,jj,jk) - hdiv_160_89_kt(jk) * sn(ji+1,jj-1,jk)410 END DO 411 ! 412 t a(ji,jj,16) = ta(ji,jj,16) - hdiv_160_89(jk) * tn(ji,jj,jk)413 sa(ji,jj,16) = sa(ji,jj,16) - hdiv_160_89(jk) * sn(ji,jj,jk)408 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - hdiv_160_89_kt(jk) * tsn(ji+1,jj-1,jk,jp_tem) 409 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - hdiv_160_89_kt(jk) * tsn(ji+1,jj-1,jk,jp_sal) 410 END DO 411 ! ! deep outflow (from Red sea) 412 tsa(ji,jj,16,jp_tem) = tsa(ji,jj,16,jp_tem) - hdiv_160_89(16) * tsn(ji,jj,16,jp_tem) 413 tsa(ji,jj,16,jp_sal) = tsa(ji,jj,16,jp_sal) - hdiv_160_89(16) * tsn(ji,jj,16,jp_sal) 414 414 END DO 415 415 END DO … … 577 577 DO ji = mi0(139), mi1(139) 578 578 DO jk = 1, jpkm1 ! surf inflow + mid. & bottom reciculation (from Atlantic) 579 t a(ji,jj,jk) = ta(ji,jj,jk) - hdiv_139_101_kt(jk) * tn(ji,jj,jk)580 sa(ji,jj,jk) = sa(ji,jj,jk) - hdiv_139_101_kt(jk) * sn(ji,jj,jk)579 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - hdiv_139_101_kt(jk) * tsn(ji,jj,jk,jp_tem) 580 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - hdiv_139_101_kt(jk) * tsn(ji,jj,jk,jp_sal) 581 581 END DO 582 582 END DO … … 586 586 DO ji = mi0(139), mi1(139) 587 587 DO jk = 15, 20 ! middle reciculation (Atl 101 -> Atl 102) (div <0) 588 t a(ji,jj,jk) = ta(ji,jj,jk) - hdiv_139_102(jk) * tn(ji,jj-1,jk) ! middle Atlantic recirculation589 sa(ji,jj,jk) = sa(ji,jj,jk) - hdiv_139_102(jk) * sn(ji,jj-1,jk)588 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - hdiv_139_102(jk) * tsn(ji,jj-1,jk,jp_tem) ! middle Atlantic recirculation 589 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - hdiv_139_102(jk) * tsn(ji,jj-1,jk,jp_sal) 590 590 END DO 591 591 ! ! upper & bottom Atl. reciculation (Atl 101 -> Atl 102) - (div <0) 592 592 ! ! deep Med flow (Med 102 -> Atl 102) - (div <0) 593 t a(ji,jj,22) = ta(ji,jj,22) + hdiv_141_102(21) * tn(ji+2,jj ,21) & ! deep Med flow594 & + hdiv_139_101(21) * t n(ji ,jj-1,21) & ! upper Atlantic recirculation595 & + hdiv_139_101(22) * t n(ji ,jj-1,22) ! bottom Atlantic recirculation596 sa(ji,jj,22) = sa(ji,jj,22) + hdiv_141_102(21) * sn(ji+2,jj ,21) &597 & + hdiv_139_101(21) * sn(ji ,jj-1,21) &598 & + hdiv_139_101(22) * sn(ji ,jj-1,22)593 tsa(ji,jj,22,jp_tem) = tsa(ji,jj,22,jp_tem) + hdiv_141_102(21) * tsn(ji+2,jj,21,jp_tem) & ! deep Med flow 594 & + hdiv_139_101(21) * tsn(ji,jj-1,21,jp_tem) & ! upper Atlantic recirculation 595 & + hdiv_139_101(22) * tsn(ji,jj-1,22,jp_tem) ! bottom Atlantic recirculation 596 tsa(ji,jj,22,jp_sal) = tsa(ji,jj,22,jp_sal) + hdiv_141_102(21) * tsn(ji+2,jj,21,jp_sal) & 597 & + hdiv_139_101(21) * tsn(ji,jj-1,21,jp_sal) & 598 & + hdiv_139_101(22) * tsn(ji,jj-1,22,jp_sal) 599 599 END DO 600 600 END DO … … 602 602 DO ji = mi0(141), mi1(141) 603 603 DO jk = 1, 14 ! surface flow from Atlantic to Med sea 604 t a(ji,jj,jk) = ta(ji,jj,jk) - hdiv_141_102_kt(jk) * tn(ji-2,jj-1,jk)605 sa(ji,jj,jk) = sa(ji,jj,jk) - hdiv_141_102_kt(jk) * sn(ji-2,jj-1,jk)604 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - hdiv_141_102_kt(jk) * tsn(ji-2,jj-1,jk,jp_tem) 605 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - hdiv_141_102_kt(jk) * tsn(ji-2,jj-1,jk,jp_sal) 606 606 END DO 607 607 ! ! deeper flow from Med sea to Atlantic 608 t a(ji,jj,21) = ta(ji,jj,21) - hdiv_141_102(21) * tn(ji,jj,21)609 sa(ji,jj,21) = sa(ji,jj,21) - hdiv_141_102(21) * sn(ji,jj,21)608 tsa(ji,jj,21,jp_tem) = tsa(ji,jj,21,jp_tem) - hdiv_141_102(21) * tsn(ji,jj,21,jp_tem) 609 tsa(ji,jj,21,jp_sal) = tsa(ji,jj,21,jp_sal) - hdiv_141_102(21) * tsn(ji,jj,21,jp_sal) 610 610 END DO 611 611 END DO … … 707 707 DO ji = mi0(172), mi1(172) 708 708 DO jk = 1, 8 ! surface inflow (Indian ocean to Persian Gulf) (div<0) 709 t a(ji,jj,jk) = ta(ji,jj,jk) - hdiv_172_94(jk) * tn(ji,jj,jk)710 sa(ji,jj,jk) = sa(ji,jj,jk) - hdiv_172_94(jk) * sn(ji,jj,jk)709 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - hdiv_172_94(jk) * tsn(ji,jj,jk,jp_tem) 710 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - hdiv_172_94(jk) * tsn(ji,jj,jk,jp_sal) 711 711 END DO 712 712 DO jk = 16, 18 ! deep outflow (Persian Gulf to Indian ocean) (div>0) 713 t a(ji,jj,jk) = ta(ji,jj,jk) - hdiv_172_94(jk) * t_171_94_hor(jk)714 sa(ji,jj,jk) = sa(ji,jj,jk) - hdiv_172_94(jk) * s_171_94_hor(jk)713 tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - hdiv_172_94(jk) * t_171_94_hor(jk) 714 tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - hdiv_172_94(jk) * s_171_94_hor(jk) 715 715 END DO 716 716 END DO
Note: See TracChangeset
for help on using the changeset viewer.