anchor.diff:      
       real (kind = x_precision), intent(in),  dimension(n_cell) :: x, Omega
       real (kind = x_precision), intent(in),  dimension(n_cell) :: T_in, S_in
  -    type(adim_state),          intent(in)                     :: state_0
       type(state),               intent(out)                    :: state_out
   
       call get_parameters(para)
  
positive.code:     enddo
    end subroutine compute_variables
    !--------------------------------------------------------------
    subroutine dim_adim(mode,state_0,state_in)
      use mod_constants
      use mod_read_parameters
      implicit none
      !select 0 or 1 to adimension or dimension your state 
      integer         , intent(in)                 :: mode
      type(adim_state), intent(inout)              :: state_0
      type(state)     , intent(inout)              :: state_in
  	
      select case(mode)
  
negative.code:     end select
    end subroutine dim_adim
   !--------------------------------------------------------------
    subroutine init_variable_0()
      use mod_constants
      use mod_read_parameters
      !Compute the initial adimention parameters
      real(kind = x_precision)         :: omega_max, rs, c2
      type(parameters)                 :: para
  
      c2 = c**2
      
  