From 4e45265aa4cd2982adab890baa933cfefdc44f19 Mon Sep 17 00:00:00 2001 From: Jeff Candy Date: Tue, 19 Nov 2024 14:19:48 -0800 Subject: [PATCH] Fixed forgotten sbeta --- cgyro/src/cgyro_write_initdata.f90 | 1 + f2py/pygacode/cgyro/data.py | 2 ++ 2 files changed, 3 insertions(+) diff --git a/cgyro/src/cgyro_write_initdata.f90 b/cgyro/src/cgyro_write_initdata.f90 index bdb5e14f8..44f32b358 100644 --- a/cgyro/src/cgyro_write_initdata.f90 +++ b/cgyro/src/cgyro_write_initdata.f90 @@ -227,6 +227,7 @@ subroutine cgyro_write_initdata do is=1,n_species write (io,fmtstr) sdlnndr(is) write (io,fmtstr) sdlntdr(is) + write (io,fmtstr) sbeta(is) enddo close(io) diff --git a/f2py/pygacode/cgyro/data.py b/f2py/pygacode/cgyro/data.py index 721a38b79..acbe0b0b0 100644 --- a/f2py/pygacode/cgyro/data.py +++ b/f2py/pygacode/cgyro/data.py @@ -496,9 +496,11 @@ def getgrid(self): # Added 3 May 2024 self.sdlnndr = np.zeros(ns) self.sdlntdr = np.zeros(ns) + self.sbeta = np.zeros(ns) for i in range(ns): self.sdlnndr[i],p = self.eget(data,p) self.sdlntdr[i],p = self.eget(data,p) + self.sbeta[i],p = self.eget(data,p) if p == -1: print('WARNING: (getgrid) Data format outdated. Please run cgyro -t')