using a gotm.yaml provided variable - fixed_grid - to allow for non-time varying z, zi

Author: bderembl

src/gotm/gotm.F90

@@ -83,6 +83,7 @@ module gotm
    use gotm_fabm, only: fabm_airp, fabm_calendar_date, fabm_julianday
    use gotm_fabm_input,only: configure_gotm_fabm_input, init_gotm_fabm_input
 #endif
+   use register_all_variables,only: fixed_grid
 
    IMPLICIT NONE
    private
@@ -273,6 +274,9 @@ subroutine initialize_gotm()
                    minimum=0._rk, default=0._rk)
    call branch%get(grid_file, 'file', 'path to file with layer thicknesses', &
                    default='')
+   call branch%get(fixed_grid, 'fixed_grid', &
+   'whether the vertical grid is fixed in time (enables static z coordinate in output)', &
+   default=.false.)
 #if 0
    twig => branch%get_child('adaptation')
    call twig%get(c1ad, 'c1ad', 'weighting factor for buoyancy frequency', '-', &

src/gotm/register_all_variables.F90

@@ -21,6 +21,7 @@ module register_all_variables
 !
 ! !PUBLIC DATA MEMBERS:
    type (type_field_manager), public, target :: fm
+   logical, public                           :: fixed_grid
 !
 ! !REVISION HISTORY:
 !  Original author(s): Karsten Bolding & Jorn Bruggeman
@@ -472,10 +473,10 @@ subroutine register_meanflow_variables(nlev)
       call fm%register('Af', 'm^2', 'hypsograph at grid interfaces', standard_name='??', dimensions=(/id_dim_z/), data1d=Af(1:nlev), category='column_structure')
    end if
 #endif
-   call fm%register('z', 'm', 'depth (center)', standard_name='??', dimensions=(/id_dim_z/), data1d=z(1:nlev), coordinate_dimension=id_dim_z,category='column_structure',field=field)
+   call fm%register('z', 'm', 'depth (center)', standard_name='??', dimensions=(/id_dim_z/), no_default_dimensions=fixed_grid, data1d=z(1:nlev), coordinate_dimension=id_dim_z,category='column_structure',field=field)
    call field%attributes%set('positive', 'up')
    call field%attributes%set('axis', 'Z')
-   call fm%register('zi', 'm', 'depth (interface)', standard_name='??', dimensions=(/id_dim_zi/), data1d=zi(0:nlev), coordinate_dimension=id_dim_zi,category='column_structure',field=field)
+   call fm%register('zi', 'm', 'depth (interface)', standard_name='??', dimensions=(/id_dim_zi/), no_default_dimensions=fixed_grid, data1d=zi(0:nlev), coordinate_dimension=id_dim_zi,category='column_structure',field=field)
    call field%attributes%set('positive', 'up')
    call field%attributes%set('axis', 'Z')
    call fm%register('h', 'm', 'layer thickness', standard_name='cell_thickness', dimensions=(/id_dim_z/), data1d=h(1:nlev),category='column_structure',part_of_state=.true.)