Changeset 1007

Timestamp:

01/17/20 13:01:35 (4 years ago)

Author:

dubos

Message:

devel : introduced plugin-based physics (TBC)

File:

• codes/icosagcm/devel/src/physics/physics.f90 (modified) (4 diffs)

codes/icosagcm/devel/src/physics/physics.f90

@@ -9 +9 @@
   INTEGER, PARAMETER :: phys_none=0, phys_column=1, &
        phys_HS94=3, phys_LB2012=4, &
-       phys_DCMIP=11, phys_DCMIP2016=12, &
+       phys_DCMIP=11, phys_DCMIP2016=12, phys_plugin=13, &
        phys_lmdz_generic=21, phys_external=22
   INTEGER :: phys_type
@@ -21 +21 @@
 !$OMP THREADPRIVATE(physics_type)
 
-  PUBLIC :: physics, init_physics, zero_du_phys
+  INTERFACE
+     SUBROUTINE plugin_init_physics
+     END SUBROUTINE plugin_init_physics
+
+     SUBROUTINE plugin_physics
+     END SUBROUTINE plugin_physics
+  END INTERFACE
+
+  PROCEDURE(plugin_init_physics), POINTER :: init_physics_plugin => NULL()
+  PROCEDURE(plugin_physics), POINTER      :: physics_plugin => NULL()
+
+  PUBLIC :: physics, init_physics, init_physics_plugin, physics_plugin, zero_du_phys
 
 CONTAINS
@@ -78 +89 @@
           phys_type = phys_LB2012
           CALL init_phys_venus       
+       CASE ('plugin')
+          IF(.NOT.ASSOCIATED(init_physics_plugin)) THEN
+             PRINT *, 'FATAL : physics = plugin selected by init_physics_plugin not set by driver program. Aborting.'
+             STOP
+          END IF
+          phys_type = phys_plugin
+          CALL init_physics_plugin
        CASE DEFAULT
           IF(is_mpi_root) PRINT*, 'init_physics : Bad selector for variable physics <',&
                TRIM(physics_type), '> options are <none>, <held_suarez>, <Lebonnois2012>,', &
-               '<dcmip>, <dcmip2016>, <phys_lmdz_generic>, <phys_external>'
+               '<dcmip>, <dcmip2016>, <plugin>, <phys_lmdz_generic>, <phys_external>'
           STOP
        END SELECT
@@ -247 +265 @@
     CASE(phys_LB2012)
        IF (is_omp_level_master) CALL full_physics_venus
+    CASE(phys_plugin)
+       IF (is_omp_level_master) CALL physics_plugin
     CASE DEFAULT
        IF(is_master) PRINT *,'Internal error : illegal value of phys_type', phys_type