Changeset 14107 for vendors/AGRIF/dev/AGRIF_FILES/modutil.F90

Timestamp:

2020-12-04T18:02:20+01:00 (4 years ago)

Author:

nicolasmartin

Message:

Reintegration of dev_r12970_AGRIF_CMEMS to AGRIF/dev

File:

• vendors/AGRIF/dev/AGRIF_FILES/modutil.F90 (modified) (6 diffs)

vendors/AGRIF/dev/AGRIF_FILES/modutil.F90

@@ -108 +108 @@
 !---------------------------------------------------------------------------------------------------
 end subroutine Agrif_Step_Child
+!===================================================================================================
+!
+!===================================================================================================
+!  subroutine Agrif_Step_Childs
+!
+!> Apply 'procname' to each child grids of the current grid
+!---------------------------------------------------------------------------------------------------
+!     **************************************************************************
+!!!   Subroutine Agrif_Step_Childs
+!     **************************************************************************
+!
+      Subroutine Agrif_Step_Childs(procname)
+!
+    procedure(step_proc)    :: procname     !< subroutine to call on each grid
+!     Pointer argument
+      Type(Agrif_Grid),pointer   :: g        ! Pointer on the current grid
+!
+
+!
+!     Local pointer
+      Type(Agrif_pgrid),pointer  :: parcours ! Pointer for the recursive
+                                             ! procedure
+!
+      g => Agrif_Curgrid
+      
+      parcours => g % child_list % first
+!
+!     Recursive procedure for the time integration of the grid hierarchy      
+      Do while (associated(parcours))
+!
+!       Instanciation of the variables of the current grid
+        Call Agrif_Instance(parcours % gr)
+
+!     One step on the current grid
+
+         Call procname ()
+        parcours => parcours % next
+      enddo
+   
+      If (associated(g % child_list % first)) Call Agrif_Instance (g)
+      Return
+      End Subroutine Agrif_Step_Childs
 !===================================================================================================
 !
@@ -538 +580 @@
 !===================================================================================================
 !
+!===================================================================================================
+!
 !
 !===================================================================================================
@@ -587 +631 @@
 #ifdef AGRIF_MPI
     else
-#endif    
 ! Continue only if the grid has defined sequences of child integrations.
     if ( .not. associated(save_grid % child_seq) ) return
@@ -610 +653 @@
 !
     enddo
-#ifdef AGRIF_MPI
     endif
 #endif 
@@ -700 +742 @@
 subroutine Agrif_Init_Grids ( procname1, procname2 )
 !---------------------------------------------------------------------------------------------------
-    procedure(typdef_proc), optional   :: procname1 !< (Default: Agrif_probdim_modtype_def)
+    procedure(typedef_proc), optional   :: procname1 !< (Default: Agrif_probdim_modtype_def)
     procedure(alloc_proc),   optional   :: procname2 !< (Default: Agrif_Allocationcalls)
 !
@@ -717 +759 @@
     nunit = Agrif_Get_Unit()
     open(nunit, file='AGRIF_FixedGrids.in', form='formatted', status="old", ERR=98)
-    if (Agrif_Probdim == 3) then
-       read(nunit,*) is_coarse, rhox, rhoy, rhoz, rhot
-    elseif (Agrif_Probdim == 2) then
-       read(nunit,*) is_coarse, rhox, rhoy, rhot
-    elseif (Agrif_Probdim == 2) then
-       read(nunit,*) is_coarse, rhox, rhot
-    endif
+    read(nunit,*) is_coarse
     if (is_coarse == -1) then
        agrif_coarse = .TRUE.
+       rewind(nunit)
        if (Agrif_Probdim == 3) then
+          read(nunit,*) is_coarse, rhox, rhoy, rhoz, rhot
           coarse_spaceref(1:3)=(/rhox,rhoy,rhoz/)
        elseif (Agrif_Probdim == 2) then
+          read(nunit,*) is_coarse, rhox, rhoy, rhot
           coarse_spaceref(1:2)=(/rhox,rhoy/)
-       elseif (Agrif_Probdim == 2) then
+       elseif (Agrif_Probdim == 1) then
+          read(nunit,*) is_coarse, rhox, rhot
           coarse_spaceref(1:1)=(/rhox/)
        endif