Changes between Version 1 and Version 2 of Developers/RebuildZooms

Timestamp:

2020-08-12T21:08:19+02:00 (4 years ago)

Author:

acc

Comment:

--

Developers/RebuildZooms

@@ -158 +158 @@
 * DOMAIN_size_zoom_i and DOMAIN_size_zoom_j are used in place of DOMAIN_size_global
 * The ibegin and jbegin offsets are subtracted from the DOMAIN_position_first values when deciding where to place values into the output array.
+
+The following changes to rebuild_nemo.F90 achieve the required result:
+
+{{{#!diff
+--- REBUILD_NEMO/src/org.rebuild_nemo.F90       2020-08-12 17:28:17.562148712 +0100
++++ REBUILD_NEMO/src/rebuild_nemo.F90   2020-08-12 20:01:36.044158059 +0100
+@@ -70,8 +70,8 @@
+    CHARACTER(LEN=50)  :: clibnc ! netcdf library version
+
+    INTEGER :: ndomain, ifile, ndomain_file, nslicesize, deflate_level
+-   INTEGER :: ncid, outid, idim, istop
+-   INTEGER :: natts, attid, xtype, varid, rbdims
++   INTEGER :: ncid, outid, idim, istop, istat
++   INTEGER :: natts, attid, xtype, varid, rbdims, rbdims1
+    INTEGER :: jv, ndims, nvars, dimlen, dimids(4)
+    INTEGER :: dimid, unlimitedDimId, di, dj, dr
+    INTEGER :: nmax_unlimited, nt, ntslice
+@@ -92,6 +92,8 @@
+    INTEGER, DIMENSION(2) :: halo_start, halo_end, local_sizes
+    INTEGER, DIMENSION(2) :: idomain, jdomain, rdomain, start_pos
+    INTEGER :: ji, jj, jk, jl, jr
++   INTEGER :: ni_zoom, nj_zoom, ni_off, nj_off
++   LOGICAL :: iszoom
+    INTEGER :: nargs                 ! number of arguments
+    INTEGER, EXTERNAL :: iargc
+
+@@ -252,10 +254,16 @@
+
+    CALL check_nf90( nf90_get_att( ncid, nf90_global, 'DOMAIN_number_total', ndomain_file ) )
+    IF( ndomain /= ndomain_file ) THEN
+-      WRITE(numerr,*) 'ERROR! : number of files to rebuild in file does not agree with namelist'
++      WRITE(numerr,*) 'WARNING! : number of files to rebuild in file does not agree with namelist'
+       WRITE(numerr,*) 'Attribute DOMAIN_number_total is : ', ndomain_file
+       WRITE(numerr,*) 'Number of files specified in namelist is: ', ndomain
+-      STOP 2
++      istat = nf90_inquire_attribute( ncid, nf90_global, 'DOMAIN_size_zoom_i', xtype, rbdims, attid )
++      IF ( istat == nf90_noerr ) THEN
++          WRITE(numerr,*) 'This looks like a zoom region so I will assume you know what you are doing'
++      ELSE
++          WRITE(numerr,*) 'This is a potentially fatal error'
++          STOP 2
++      ENDIF
+    ENDIF
+
+ !2.1 Set up the output file
+@@ -275,6 +283,29 @@
+
+    ALLOCATE(global_sizes(rbdims))
+    CALL check_nf90( nf90_get_att( ncid, nf90_global, 'DOMAIN_size_global', global_sizes ) )
++
++!2.2.0.1 Override global sizes if zoom attributes are found
++   istat = nf90_inquire_attribute( ncid, nf90_global, 'DOMAIN_size_zoom_i', xtype, rbdims1, attid )
++   IF ( istat == nf90_noerr ) THEN
++       CALL check_nf90( nf90_get_att( ncid, nf90_global, 'DOMAIN_size_zoom_i', ni_zoom ) )
++       iszoom = .true.
++   ELSE
++       iszoom = .false.
++   ENDIF
++   istat = nf90_inquire_attribute( ncid, nf90_global, 'DOMAIN_size_zoom_j', xtype, rbdims1, attid )
++   IF ( istat == nf90_noerr ) THEN
++       CALL check_nf90( nf90_get_att( ncid, nf90_global, 'DOMAIN_size_zoom_j', nj_zoom ) )
++       iszoom = .true.
++   ELSE
++       iszoom = .false.
++   ENDIF
++   IF( iszoom ) THEN
++       global_sizes(1) = ni_zoom
++       global_sizes(2) = nj_zoom
++       CALL check_nf90( nf90_get_att( ncid, nf90_global, 'ibegin', ni_off ) )
++       CALL check_nf90( nf90_get_att( ncid, nf90_global, 'jbegin', nj_off ) )
++   ENDIF
++
+    IF (l_verbose) WRITE(numout,*) 'Size of global arrays: ', global_sizes
+
+
+@@ -773,6 +804,9 @@
+                halo_start(2) = 0
+                di=rebuild_dims(1)
+                dj=3-di
++            ELSEIF ( iszoom ) THEN
++               start_pos(di) = start_pos(di) - ni_off
++               start_pos(dj) = start_pos(dj) - nj_off
+             ENDIF
+
+ !3.3.1 Generate local domain interior sizes from local_sizes and halo sizes
+}}}
 
 == Notes for possibly tackling the problem at source ==