diff --git a/force/geo.f90 b/force/geo.f90
index a85961d1beb923c60bee9581912ac6fcd7651a9f..a9c749841f8b5c45c2ded3fccec92ebf56d119da 100644
--- a/force/geo.f90
+++ b/force/geo.f90
@@ -238,6 +238,7 @@ CONTAINS
           input_temp%l_f = input%l_f
           input_temp%gw_neigd = dimension_temp%neigd
           div(:) = MIN(kpts_temp%nmop(:),1)
+          stars_temp%gmax = stars_temp%gmaxInit
           CALL w_inpXML(atoms_new,obsolete_temp,vacuum_temp,input_temp,stars_temp,sliceplot_temp,&
                         banddos_temp,cell_temp,sym_temp,xcpot_temp,noco_temp,jij_temp,oneD_temp,hybrid_temp,&
                         kpts_temp,kpts_temp%nmop,kpts_temp%l_gamma,noel_temp,namex_temp,relcor_temp,a1_temp,a2_temp,a3_temp,&
diff --git a/global/types.F90 b/global/types.F90
index 45c30b49285dfc035a2697956ca99b5f4440944a..911c2dc1dba05c8ad9adce7a9ebd97b6a0b21813 100644
--- a/global/types.F90
+++ b/global/types.F90
@@ -327,6 +327,7 @@
 
         !max-length of star
         REAL :: gmax
+        REAL :: gmaxInit
         !no of 3d-stars
         !INTEGER :: nq3
         INTEGER :: ng3
diff --git a/io/r_inpXML.F90 b/io/r_inpXML.F90
index dfa0f5e5b40a66cee46f146384383bd34e0d64d9..a19df8ed229cee05672bc6a7c18c1c3123fcde49 100644
--- a/io/r_inpXML.F90
+++ b/io/r_inpXML.F90
@@ -326,6 +326,7 @@ SUBROUTINE r_inpXML(&
    IF(numberNodes.EQ.1) THEN
       xcpot%gmaxxc = evaluateFirstOnly(xmlGetAttributeValue(xPathA))
    END IF
+   stars%gmaxInit = stars%gmax
 
    xPathA = '/fleurInput/calculationSetup/cutoffs/@numbands'
    numberNodes = xmlGetNumberOfNodes(xPathA)