Changeset 13782
- Timestamp:
- 2020-11-11T18:41:55+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13333_KERNEL-08_techene_gm_HPG_SPG/src/NST/agrif_oce_update.F90
r13286 r13782 886 886 REAL(wp), DIMENSION(i1:i2,j1:j2), INTENT(inout) :: tabres 887 887 LOGICAL , INTENT(in ) :: before 888 !! 889 REAL(wp), DIMENSION(jpi,jpj) :: zpgu ! 2D workspace 888 890 !! 889 891 INTEGER :: ji, jj, jk … … 905 907 ! 906 908 ! Update "now" 3d velocities: 907 spgu(ji,jj) = 0._wp909 zpgu(ji,jj) = 0._wp 908 910 DO jk=1,jpkm1 909 spgu(ji,jj) = spgu(ji,jj) + e3u(ji,jj,jk,Kmm_a) * uu(ji,jj,jk,Kmm_a)911 zpgu(ji,jj) = zpgu(ji,jj) + e3u(ji,jj,jk,Kmm_a) * uu(ji,jj,jk,Kmm_a) 910 912 END DO 911 913 ! 912 zcorr = (tabres(ji,jj) - spgu(ji,jj)) * r1_hu(ji,jj,Kmm_a)914 zcorr = (tabres(ji,jj) - zpgu(ji,jj)) * r1_hu(ji,jj,Kmm_a) 913 915 DO jk=1,jpkm1 914 916 uu(ji,jj,jk,Kmm_a) = uu(ji,jj,jk,Kmm_a) + zcorr * umask(ji,jj,jk) … … 925 927 ! 926 928 ! Correct "before" velocities to hold correct bt component: 927 spgu(ji,jj) = 0.e0929 zpgu(ji,jj) = 0.e0 928 930 DO jk=1,jpkm1 929 spgu(ji,jj) = spgu(ji,jj) + e3u(ji,jj,jk,Kbb_a) * uu(ji,jj,jk,Kbb_a)931 zpgu(ji,jj) = zpgu(ji,jj) + e3u(ji,jj,jk,Kbb_a) * uu(ji,jj,jk,Kbb_a) 930 932 END DO 931 933 ! 932 zcorr = uu_b(ji,jj,Kbb_a) - spgu(ji,jj) * r1_hu(ji,jj,Kbb_a)934 zcorr = uu_b(ji,jj,Kbb_a) - zpgu(ji,jj) * r1_hu(ji,jj,Kbb_a) 933 935 DO jk=1,jpkm1 934 936 uu(ji,jj,jk,Kbb_a) = uu(ji,jj,jk,Kbb_a) + zcorr * umask(ji,jj,jk) … … 953 955 REAL(wp), DIMENSION(i1:i2,j1:j2), INTENT(inout) :: tabres 954 956 LOGICAL , INTENT(in ) :: before 957 ! 958 REAL(wp), DIMENSION(jpi,jpj) :: zpgv ! 2D workspace 955 959 ! 956 960 INTEGER :: ji, jj, jk … … 971 975 ! 972 976 ! Update "now" 3d velocities: 973 spgv(ji,jj) = 0.e0977 zpgv(ji,jj) = 0.e0 974 978 DO jk=1,jpkm1 975 spgv(ji,jj) = spgv(ji,jj) + e3v(ji,jj,jk,Kmm_a) * vv(ji,jj,jk,Kmm_a)979 zpgv(ji,jj) = zpgv(ji,jj) + e3v(ji,jj,jk,Kmm_a) * vv(ji,jj,jk,Kmm_a) 976 980 END DO 977 981 ! 978 zcorr = (tabres(ji,jj) - spgv(ji,jj)) * r1_hv(ji,jj,Kmm_a)982 zcorr = (tabres(ji,jj) - zpgv(ji,jj)) * r1_hv(ji,jj,Kmm_a) 979 983 DO jk=1,jpkm1 980 984 vv(ji,jj,jk,Kmm_a) = vv(ji,jj,jk,Kmm_a) + zcorr * vmask(ji,jj,jk) … … 991 995 ! 992 996 ! Correct "before" velocities to hold correct bt component: 993 spgv(ji,jj) = 0.e0997 zpgv(ji,jj) = 0.e0 994 998 DO jk=1,jpkm1 995 spgv(ji,jj) = spgv(ji,jj) + e3v(ji,jj,jk,Kbb_a) * vv(ji,jj,jk,Kbb_a)999 zpgv(ji,jj) = zpgv(ji,jj) + e3v(ji,jj,jk,Kbb_a) * vv(ji,jj,jk,Kbb_a) 996 1000 END DO 997 1001 ! 998 zcorr = vv_b(ji,jj,Kbb_a) - spgv(ji,jj) * r1_hv(ji,jj,Kbb_a)1002 zcorr = vv_b(ji,jj,Kbb_a) - zpgv(ji,jj) * r1_hv(ji,jj,Kbb_a) 999 1003 DO jk=1,jpkm1 1000 1004 vv(ji,jj,jk,Kbb_a) = vv(ji,jj,jk,Kbb_a) + zcorr * vmask(ji,jj,jk)
Note: See TracChangeset
for help on using the changeset viewer.