CD: Fortran bugfix for malloc error w/ old MAS files

- instead of deallocate, use nullify to clear pointer assignment to the
  input python B data when adding a phi point and updating the B data in mapfl
- add option in testing routine dipole_field to emulate dropping a point

Author: cooperdowns

src/mapfl/mapfl.f

@@ -48,8 +48,8 @@ module ident
 c
 c
       character(*), parameter :: cname='MAPFL'
-      character(*), parameter :: cvers='2.1.0'
-      character(*), parameter :: cdate='04/16/2024'
+      character(*), parameter :: cvers='2.1.0_mapflpy'
+      character(*), parameter :: cdate='02/09/2026'
 c
       end module
 c#######################################################################
@@ -1546,8 +1546,8 @@ subroutine set_btype (b)
         f(:,:,n3+1)=b%r%f(:,:,1)
         z(1:n3)=b%r%scales(3)%f(:)
         z(n3+1)=b%r%scales(3)%f(1)+twopi
-        deallocate (b%r%f)
-        deallocate (b%r%scales(3)%f)
+        nullify (b%r%f)
+        nullify (b%r%scales(3)%f)
         b%r%dims(3)=n3+1
         b%r%f=>f
         b%r%scales(3)%f=>z
@@ -1561,8 +1561,8 @@ subroutine set_btype (b)
         f(:,:,n3+1)=b%t%f(:,:,1)
         z(1:n3)=b%t%scales(3)%f(:)
         z(n3+1)=b%t%scales(3)%f(1)+twopi
-        deallocate (b%t%f)
-        deallocate (b%t%scales(3)%f)
+        nullify (b%t%f)
+        nullify (b%t%scales(3)%f)
         b%t%dims(3)=n3+1
         b%t%f=>f
         b%t%scales(3)%f=>z
@@ -1576,8 +1576,8 @@ subroutine set_btype (b)
         f(:,:,n3+1)=b%p%f(:,:,1)
         z(1:n3)=b%p%scales(3)%f(:)
         z(n3+1)=b%p%scales(3)%f(1)+twopi
-        deallocate (b%p%f)
-        deallocate (b%p%scales(3)%f)
+        nullify (b%p%f)
+        nullify (b%p%scales(3)%f)
         b%p%dims(3)=n3+1
         b%p%f=>f
         b%p%scales(3)%f=>z

tests/utils.py

@@ -30,7 +30,8 @@ def dipole_field(lon: float = 180.,
                  r1: float = 10.0,
                  r_resolution: int = 71,
                  t_resolution: int = 91,
-                 p_resolution: int = 181
+                 p_resolution: int = 181,
+                 remove_phi_point: bool = False,
                  ) -> Tuple[NDArray[float], NDArray[float], NDArray[float], NDArray[float], NDArray[float], NDArray[float]]:
     """Generate a 3D dipole field in spherical coordinates.
 
@@ -54,6 +55,8 @@ def dipole_field(lon: float = 180.,
         number of theta points. Default is 91.
     p_resolution: int
         number of phi points. Default is 181.
+    remove_phi_point: bool
+        Option to delete the last phi point to emulate an "old" MAS file
 
     Returns
     -------
@@ -107,6 +110,12 @@ def dipole_field(lon: float = 180.,
 
     br, bt, bp = cvtosv(bx, by, bz, t3d, p3d)
 
+    if remove_phi_point:
+        br = br[:-1,:,:]
+        bt = bt[:-1,:,:]
+        bp = bp[:-1,:,:]
+        p = p[:-1]
+
     return br, bt, bp, r, t, p