Changeset 10160 for utils/tools/NESTING/src/agrif_readwrite.f90

Timestamp:

2018-09-28T16:21:53+02:00 (6 years ago)

Author:

clem

Message:

forgotten commits on nesting tool

File:

• utils/tools/NESTING/src/agrif_readwrite.f90 (modified) (5 diffs)

utils/tools/NESTING/src/agrif_readwrite.f90

@@ -502 +502 @@
      REAL*8 ,              DIMENSION(N)     ::   e3t_1d, e3w_1d, gdept_1d
      REAL*8 , ALLOCATABLE, DIMENSION(:,:)   ::   ff_t, ff_f
-     INTEGER, ALLOCATABLE, DIMENSION(:,:)   ::   bottom_level, top_level
+     INTEGER, ALLOCATABLE, DIMENSION(:,:)   ::   bottom_level, top_level, zbathy
      REAL*8 , ALLOCATABLE, DIMENSION(:,:,:) ::   e3t_0, e3u_0, e3v_0, e3f_0, e3w_0, e3uw_0, e3vw_0
      !
@@ -511 +511 @@
      ! allocate needed arrays for domain_cfg
      ALLOCATE( ff_t(nx,ny), ff_f(nx,ny) )
-     ALLOCATE( bottom_level(nx,ny), top_level(nx,ny) )
+     ALLOCATE( bottom_level(nx,ny), top_level(nx,ny), zbathy(nx,ny) )
      ALLOCATE( e3t_0(nx,ny,N), e3u_0 (nx,ny,N), e3v_0 (nx,ny,N), e3f_0(nx,ny,N),  &
         &      e3w_0(nx,ny,N), e3uw_0(nx,ny,N), e3vw_0(nx,ny,N) )
@@ -567 +567 @@
      top_level(  nx,1:ny) = 0
 
+     zbathy(:,:) = Grid%bathy_meter
+     zbathy(1:nx,1   ) = 0.
+     zbathy(1:nx,  ny) = 0.
+     zbathy(1   ,1:ny) = 0.
+     zbathy(  nx,1:ny) = 0.     
+     
      ! -------------------
      ! write domain_cfg.nc
@@ -629 +635 @@
      CALL write_ncdf_var('bottom_level',dimnames(1:2),name,bottom_level,'integer')
      CALL write_ncdf_var('top_level'   ,dimnames(1:2),name,top_level   ,'integer')
-
-     CALL write_ncdf_var('bathy_meter' ,dimnames(1:2),name,Grid%bathy_meter,'float')
+     CALL write_ncdf_var('bathy_meter' ,dimnames(1:2),name,zbathy      ,'float')
 
      ! some attributes      
@@ -660 +665 @@
      !
      DEALLOCATE( ff_t, ff_f )
-     DEALLOCATE( bottom_level, top_level )
+     DEALLOCATE( bottom_level, top_level, zbathy )
      DEALLOCATE( e3t_0, e3u_0, e3v_0, e3f_0, e3w_0, e3uw_0, e3vw_0 )
      !