Changeset 11625 for NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/NST/agrif_oce_sponge.F90

Timestamp:

2019-10-01T15:18:44+02:00 (5 years ago)

Author:

jchanut

Message:

#2222, add initialization to 0 of tracer open boundary data with vertical interpolation + various neutral optimizations

File:

• NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/NST/agrif_oce_sponge.F90 (modified) (4 diffs)

NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/NST/agrif_oce_sponge.F90

@@ -329 +329 @@
          !
 # if defined key_vertical
-         tabres_child(:,:,:,:) = 0.
          DO jj=j1,j2
             DO ji=i1,i2
+               tabres_child(ji,jj,:,:) = 0._wp 
                N_in = 0
                DO jk=k1,k2 !k2 = jpk of parent grid
@@ -345 +345 @@
                   h_out(jk) = e3t_b(ji,jj,jk) !Child grid scale factors. Could multiply by e1e2t here instead of division above
                ENDDO
-               IF (N_in > 0) THEN
+               IF (N_in*N_out > 0) THEN
                   CALL reconstructandremap(tabin(1:N_in,1:jpts),h_in(1:N_in),tabres_child(ji,jj,1:N_out,1:jpts),h_out(1:N_out),N_in,N_out,jpts)
                ENDIF
@@ -438 +438 @@
       !
       IF( before ) THEN
-         DO jk=1,jpkm1
+         DO jk=k1,k2
             DO jj=j1,j2
                DO ji=i1,i2
@@ -591 +591 @@
       
       IF( before ) THEN 
-         DO jk=1,jpkm1
+         DO jk=k1,k2
             DO jj=j1,j2
                DO ji=i1,i2