New stochastic physics API

iangrooms · Aug 20, 2024 · 496d219 · 496d219
1 parent 978b8e6
commit 496d219
Showing 1 changed file with 12 additions and 6 deletions.
diff --git a/config_src/external/stochastic_physics/stochastic_physics.F90 b/config_src/external/stochastic_physics/stochastic_physics.F90
@@ -15,14 +15,20 @@ module stochastic_physics
 contains
 
 !> Initializes the stochastic physics perturbations.
-subroutine init_stochastic_physics_ocn(delt, geoLonT, geoLatT, nx, ny, nz, pert_epbl_in, do_sppt_in, &
-                                       do_skeb_in,mpiroot, mpicomm, iret)
+subroutine init_stochastic_physics_ocn(delt, geoLonT, geoLatT, nxT, nyT, nz, &
+                                             geoLonB, geoLatB, nxB, nyB, &
+                                             pert_epbl_in, do_sppt_in, &
+                                             do_skeb_in, mpiroot, mpicomm, iret)
   real,    intent(in)    :: delt !< timestep in seconds between calls to run_stochastic_physics_ocn [s]
-  integer, intent(in)    :: nx   !< number of gridpoints in the x-direction of the compute grid
-  integer, intent(in)    :: ny   !< number of gridpoints in the y-direction of the compute grid
+  integer, intent(in)    :: nxT  !< number of T gridpoints in the x-direction of the compute grid
+  integer, intent(in)    :: nyT  !< number of T gridpoints in the y-direction of the compute grid
   integer, intent(in)    :: nz   !< number of gridpoints in the z-direction of the compute grid
-  real,    intent(in)    :: geoLonT(nx,ny) !< Longitude in degrees
-  real,    intent(in)    :: geoLatT(nx,ny) !< Latitude in degrees
+  real,    intent(in)    :: geoLonT(nxT,nyT) !< Longitude of T points in degrees
+  real,    intent(in)    :: geoLatT(nxT,nyT) !< Latitude of T points in degrees
+  integer, intent(in)    :: nxB  !< number of B gridpoints in the x-direction of the compute grid
+  integer, intent(in)    :: nyB  !< number of B gridpoints in the y-direction of the compute grid
+  real,    intent(in)    :: geoLonB(nxB,nyB) !< Longitude of B points in degrees
+  real,    intent(in)    :: geoLatB(nxB,nyB) !< Latitude of B points in degrees
   logical, intent(in)    :: pert_epbl_in !< logical flag, if true generate random pattern for ePBL perturbations
   logical, intent(in)    :: do_sppt_in   !< logical flag, if true generate random pattern for SPPT perturbations
   logical, intent(in)    :: do_skeb_in   !< logical flag, if true generate random pattern for SKEB perturbations