Changeset 7025 for branches/2016/dev_r6999_CONFIGMAN_1/NEMOGCM/TOOLS/SIREN/src/boundary.f90

Timestamp:

2016-10-13T11:08:38+02:00 (8 years ago)

Author:

jpaul

Message:

see ticket #1781

File:

• branches/2016/dev_r6999_CONFIGMAN_1/NEMOGCM/TOOLS/SIREN/src/boundary.f90 (modified) (2 diffs)

branches/2016/dev_r6999_CONFIGMAN_1/NEMOGCM/TOOLS/SIREN/src/boundary.f90

@@ -1529 +1529 @@
    !> @author J.Paul 
    !> @date November, 2013 - Initial Version 
-   ! 
+   !> @date June, 2016
+   !> - Bug fix: take into account that boundaries are compute on T point,
+   !>   but expressed on U,V point
+   !>
    !> @param[inout] td_bdy boundary structure 
    !> @param[in] td_var    variable structure 
@@ -1552 +1555 @@
       il_max(jp_west )=td_var%t_dim(2)%i_len
  
-      il_maxindex(jp_north)=td_var%t_dim(2)%i_len-ip_ghost
+      ! index expressed on U,V point, move on T point.
+      il_maxindex(jp_north)=td_var%t_dim(2)%i_len-ip_ghost+1
       il_maxindex(jp_south)=td_var%t_dim(2)%i_len-ip_ghost
-      il_maxindex(jp_east )=td_var%t_dim(1)%i_len-ip_ghost
+      il_maxindex(jp_east )=td_var%t_dim(1)%i_len-ip_ghost+1
       il_maxindex(jp_west )=td_var%t_dim(1)%i_len-ip_ghost