Changeset 14000
- Timestamp:
- 2020-12-02T15:03:17+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_oce_interp.F90
r13937 r14000 831 831 SUBROUTINE interptsn( ptab, i1, i2, j1, j2, k1, k2, n1, n2, before ) 832 832 !!---------------------------------------------------------------------- 833 !! *** ROUTINE interptsn ***834 !!----------------------------------------------------------------------835 833 REAL(wp), DIMENSION(i1:i2,j1:j2,k1:k2,n1:n2), INTENT(inout) :: ptab 836 834 INTEGER , INTENT(in ) :: i1, i2, j1, j2, k1, k2, n1, n2 … … 902 900 903 901 ! Intermediate grid: 904 DO jk = 1, N_in 905 h_in_i(jk) = e3t0_parent(ji,jj,jk) * & 906 & (1._wp + ptab(ji,jj,k2,n2)/(ht0_parent(ji,jj)*ssmask(ji,jj) + 1._wp - ssmask(ji,jj))) 907 END DO 908 z_in_i(1) = 0.5_wp * h_in_i(1) 909 DO jk=2,N_in 910 z_in_i(jk) = z_in_i(jk-1) + 0.5_wp * ( h_in_i(jk) + h_in_i(jk-1) ) 911 END DO 912 z_in_i(1:N_in) = z_in_i(1:N_in) - ptab(ji,jj,k2,n2) 902 IF ( l_vremap ) THEN 903 DO jk = 1, N_in 904 h_in_i(jk) = e3t0_parent(ji,jj,jk) * & 905 & (1._wp + ptab(ji,jj,k2,n2)/(ht0_parent(ji,jj)*ssmask(ji,jj) + 1._wp - ssmask(ji,jj))) 906 END DO 907 z_in_i(1) = 0.5_wp * h_in_i(1) 908 DO jk=2,N_in 909 z_in_i(jk) = z_in_i(jk-1) + 0.5_wp * ( h_in_i(jk) + h_in_i(jk-1) ) 910 END DO 911 z_in_i(1:N_in) = z_in_i(1:N_in) - ptab(ji,jj,k2,n2) 912 ENDIF 913 913 914 914 ! Output (Child) grid: … … 1080 1080 ! h_in(jk) = ptab(ji,jj,jk,2)/(e2u(ji,jj)*zrhoy) 1081 1081 !ENDIF 1082 h_in(jk) = e3u0_parent(ji,jj,jk) 1083 zhtot = zhtot + h_in(jk) 1082 IF ( l_vremap ) THEN 1083 h_in(jk) = e3u0_parent(ji,jj,jk) 1084 zhtot = zhtot + h_in(jk) 1085 ENDIF 1084 1086 IF( h_in(jk) .GT. 0. ) THEN 1085 1087 tabin(jk) = ptab(ji,jj,jk,1)/(e2u(ji,jj)*zrhoy*h_in(jk)) … … 1203 1205 ! h_in(jk) = ptab(ji,jj,jk,2)/(e1v(ji,jj)*zrhox) 1204 1206 !ENDIF 1205 h_in(jk) = e3v0_parent(ji,jj,jk) 1206 zhtot = zhtot + h_in(jk) 1207 IF (l_vremap) THEN 1208 h_in(jk) = e3v0_parent(ji,jj,jk) 1209 zhtot = zhtot + h_in(jk) 1210 ENDIF 1207 1211 IF( h_in(jk) .GT. 0. ) THEN 1208 1212 tabin(jk) = ptab(ji,jj,jk,1)/(e1v(ji,jj)*zrhox*h_in(jk))
Note: See TracChangeset
for help on using the changeset viewer.