Changeset 6445 for branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp
- Timestamp:
- 2016-04-07T17:37:50+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/NEMO/NST_SRC/agrif_opa_interp.F90
r6425 r6445 38 38 PUBLIC Agrif_Init_InterpScales 39 39 ! VERTICAL REFINEMENT END 40 PUBLIC interpun, interpvn , interpun2d, interpvn2d40 PUBLIC interpun, interpvn 41 41 PUBLIC interptsn, interpsshn 42 42 PUBLIC interpunb, interpvnb, interpub2b, interpvb2b … … 1023 1023 1024 1024 1025 SUBROUTINE interpun2d(ptab,i1,i2,j1,j2,before)1026 !!---------------------------------------------1027 !! *** ROUTINE interpun ***1028 !!---------------------------------------------1029 !1030 INTEGER, INTENT(in) :: i1,i2,j1,j21031 REAL(wp), DIMENSION(i1:i2,j1:j2), INTENT(inout) :: ptab1032 LOGICAL, INTENT(in) :: before1033 !1034 INTEGER :: ji,jj1035 REAL(wp) :: ztref1036 REAL(wp) :: zrhoy1037 !!---------------------------------------------1038 !1039 ztref = 1.1040 1041 IF (before) THEN1042 DO jj=j1,j21043 DO ji=i1,MIN(i2,nlci-1)1044 ptab(ji,jj) = e2u(ji,jj) * ((gcx(ji+1,jj) - gcx(ji,jj))/e1u(ji,jj))1045 END DO1046 END DO1047 ELSE1048 zrhoy = Agrif_Rhoy()1049 DO jj=j1,j21050 laplacu(i1:i2,jj) = ztref * (ptab(i1:i2,jj)/(zrhoy*e2u(i1:i2,jj))) !*umask(i1:i2,jj,1)1051 END DO1052 ENDIF1053 !1054 END SUBROUTINE interpun2d1055 1056 1025 1057 1026 SUBROUTINE interpvn(ptab,i1,i2,j1,j2,k1,k2,m1,m2,before,nb,ndir) … … 1059 1028 !! *** ROUTINE interpvn *** 1060 1029 !!---------------------------------------------------------------------- 1061 INTEGER , INTENT(in ) :: i1, i2, j1, j2, k1, k21062 REAL(wp), DIMENSION(i1:i2,j1:j2,k1:k2), INTENT(inout) :: ptab1063 LOGICAL , INTENT(in ) :: before1064 1030 ! 1065 1031 INTEGER, INTENT(in) :: i1,i2,j1,j2,k1,k2,m1,m2
Note: See TracChangeset
for help on using the changeset viewer.