Changes between Version 3 and Version 4 of Developers/RebuildZooms

Timestamp:

2020-08-18T11:54:19+02:00 (4 years ago)

Author:

acc

Comment:

--

Developers/RebuildZooms

@@ -110 +110 @@
  NEMO rebuild completed successfully
 }}}
-This successfully rebuilds the zoom but places it in an otherwise empty global domain.
+This successfully rebuilds the zoom but places it in an otherwise empty global domain:
+
+[[Image(zoom_full.png)]]
 
 Fixing the second issue is trickier. Simply editing the DOMAIN_size_global settings will not suffice because REBUILD_NEMO also uses the DOMAIN_position_first information to place data within the global arrays. Changing the size but not the offset results in Bus errors.
@@ -162 +164 @@
 
 {{{#!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
+--- org.rebuild_nemo.F90        2020-08-12 17:28:17.562148712 +0100
++++ rebuild_nemo.F90    2020-08-18 10:36:12.747885729 +0100
 @@ -70,8 +70,8 @@
     CHARACTER(LEN=50)  :: clibnc ! netcdf library version
@@ -203 +205 @@
 
  !2.1 Set up the output file
-@@ -275,6 +283,29 @@
+@@ -275,6 +283,26 @@
 
     ALLOCATE(global_sizes(rbdims))
@@ -209 +211 @@
 +
 +!2.2.0.1 Override global sizes if zoom attributes are found
++   iszoom = .false.
 +   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 ) )
++       ! Need both zoom_i and zoom_j attributes to determine zoom size
++       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.
++           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 ) )
++       ELSE
++           iszoom = .false.
++       ENDIF
 +   ENDIF
 +
@@ -233 +232 @@
 
 
-@@ -773,6 +804,9 @@
+@@ -773,6 +801,9 @@
                 halo_start(2) = 0
                 di=rebuild_dims(1)
@@ -244 +243 @@
  !3.3.1 Generate local domain interior sizes from local_sizes and halo sizes
 }}}
+
+This modified rebuild_nemo will, if the correct attributes are found, rebuild just the zoom region. It will also ignore any mismatch between the DOMAIN_number_total attribute and the number of files specified by the user; trusting that the user has correctly specified the true number.
+[[Image(zoom_only.png)]]
 
 == Notes for possibly tackling the problem at source ==