Refactor source codes with 2D advection test to evalaute numerical er…

…rors
ywkawai · Jan 16, 2025 · 3a6e25a · 3a6e25a
1 parent 6130acf
commit 3a6e25a
Show file tree

Hide file tree

Showing 5 changed files with 262 additions and 243 deletions.
diff --git a/sample/advect2d/Makefile b/sample/advect2d/Makefile
@@ -9,7 +9,7 @@ TOPDIR      = $(abspath ../..)
 SAMPLE_DIR  = ../
 
 BINNAME = test_advect2d
-OBJS = 
+OBJS = mod_advect2d_numerror.o
 
 RUNCONF         = test.conf
 TPROC           = 1
@@ -20,3 +20,8 @@ include $(SAMPLE_DIR)/Makefile.common
 
 vis:
 	bash ./visualize/visualize.sh
+
+clean_analysis:
+	rm -f analysis
+
+$(BUILD_DIR)/test_advect2d.o : $(BUILD_DIR)/mod_advect2d_numerror.o
diff --git a/sample/advect2d/test.conf b/sample/advect2d/test.conf
@@ -18,26 +18,33 @@
   ! InitShapeName      = 'top-hat', 
   ! InitShapeParams    = 0.5D0, 0.5D0, 0.15D0, 0.15D0,   
   !-
-  InitCond_GalerkinProjFlag = .false., 
   !** Type of advection velocity **********************
   ! VelTypeName      = 'constant', 
   ! VelTypeParams    = 1.0D0, 1.0D0, 0.0D0, 0.0D0, 
+  ! Do_NumErrorAnalysis = .true., 
   !- 
    VelTypeName      = 'rigid-body-rot', 
    VelTypeParams    = 0.5D0, 0.5D0, 1.0D0, 0.0D0,   
+   Do_NumErrorAnalysis = .true., 
   !-
   ! VelTypeName      = 'swirling', 
   ! VelTypeParams    = 0.0D0, 0.0D0, 5.0D0, 0.0D0,   
+  ! Do_NumErrorAnalysis = .false., 
   !----------------------------------------------------
   TINTEG_SCHEME_TYPE = 'ERK_SSP_3s3o',  
-  nstep_eval_error = 50, 
+  InitGPMatPolyOrder = 7,
 /
 &PARAM_TIME
   TIME_DURATION                 = 5.0D0, 
   TIME_DURATION_UNIT            = 'SEC', 
   TIME_DT                       = 0.0005D0, 
   TIME_DT_UNIT                  = 'SEC', 
 /
+&PARAM_ADVECT2D_NUMERROR
+  PolyOrderErrorCheck = 6,
+  LOG_OUT_BASENAME    = 'LOG_NUMERROR', 
+  LOG_STEP_INTERVAL   = 100, 
+/
 &PARAM_FILE_HISTORY
  FILE_HISTORY_DEFAULT_BASENAME  = "history",
  FILE_HISTORY_DEFAULT_TINTERVAL = 0.1D0,