Changeset 14019

Timestamp:

2020-12-02T18:50:19+01:00 (4 years ago)

Author:

jchanut

Message:

#2222, use right thicknesses for initial state interpolation from parent

File:

• NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_oce_interp.F90 (modified) (2 diffs)

NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_oce_interp.F90

@@ -1079 +1079 @@
                      IF ( l_vremap ) THEN
                         h_in(jk) = e3u0_parent(ji,jj,jk) 
-                        zhtot = zhtot + h_in(jk)
+                     ELSE
+                        IF (jk==N_in) THEN
+                           h_in(jk) = hu0_parent(ji,jj) + ptab(ji,jj,k2,2) - zhtot
+                        ELSE
+                           h_in(jk) = ptab(ji,jj,jk,2)/(e2u(ji,jj)*zrhoy) 
+                        ENDIF
                      ENDIF
+                     zhtot = zhtot + h_in(jk)
                      IF( h_in(jk) .GT. 0. ) THEN
                      tabin(jk) = ptab(ji,jj,jk,1)/(e2u(ji,jj)*zrhoy*h_in(jk))
@@ -1204 +1210 @@
                      IF (l_vremap) THEN
                         h_in(jk) = e3v0_parent(ji,jj,jk)
-                        zhtot = zhtot + h_in(jk)
+                     ELSE
+                        IF (jk==N_in) THEN
+                           h_in(jk) = hv0_parent(ji,jj) + ptab(ji,jj,k2,2) - zhtot
+                        ELSE
+                           h_in(jk) = ptab(ji,jj,jk,2)/(e1v(ji,jj)*zrhox) 
+                        ENDIF
                      ENDIF
+                     zhtot = zhtot + h_in(jk)
                      IF( h_in(jk) .GT. 0. ) THEN
                        tabin(jk) = ptab(ji,jj,jk,1)/(e1v(ji,jj)*zrhox*h_in(jk))