NESTING TOOLS AGRIF read updated bathymetry for creating agrif restarts

[jharlass]

Context

Use bathy_updated from base model to correctly compute masks (partial steps) for nest restart creation.

Changes

agrif_create_restart.f90

(a) Revision vs. (b) Working

@@ -113 +113 @@
   !       
   CALL set_child_name(parent_coordinate_file,Childcoordinates)   
   CALL set_child_name(parent_meshmask_file,Childbathy) 
   CALL set_child_name(parent_bathy_meter,Childbathymeter)   
+  ! use bathy_updated for base model ONLY, if bathy_update=true in namelist.input
+  IF (bathy_update) THEN
+    parent_bathy_meter=updated_parent_file
+  END IF
+
   ! create this file
   !
@@ -164 +170 @@
      !
   ENDIF
   !
+  WRITE(*,*) 'Reading Bathymetry: ',TRIM(parent_bathy_meter)
+  CALL Read_Ncdf_var('Bathymetry',parent_bathy_meter,G0%bathy_meter)
+  
+  WRITE(*,*) 'Reading mbathy    : ',TRIM(parent_meshmask_file)
   CALL Init_mask(parent_meshmask_file,G0,x,y)
+  WRITE(*,*)
+  WRITE(*,*) 'Reading Bathymetry: ',TRIM(Childbathymeter)
+  CALL Read_Ncdf_var('Bathymetry',Childbathymeter,G1%bathy_meter)
+  WRITE(*,*) 'Reading mbathy    : ',TRIM(childbathy)
   CALL Init_mask(childbathy,G1,1,1)
 
   G0%tmask = 1.   

[jchanut]

Hi,

Could you please explain why you need the updated parent bathymetry ? Input restart matches the original (i.e. not updated domain) so that I don't understand why you advise this. Maybe I missed something ?
Thanks for your help,
Jerome

[jharlass]

You are right and I can not entirely recall the reason. From my current perspective it is indeed better to use bathy_meter.nc to treat the input fields equally. We came across this issue when adding the partial step support in the creation process of all masks (Init_mask). When a nested configuration starts and the first update takes place immediately, base values will be overwritten by nest values in any way. So either either bathy_updated would not fit not fit to the base restart or the base restart will not fit to 1_bathy_meter.
Ticket can be closed and proposed change dismissed.

[jchanut]

Thanks Jan. I close the ticket.