the constant-in-time base-state densities "rhobf" and "rhobh" are created as time-dependent vars in the output netCDF file

[slayoo]

Here's the relevant part of ncdump output:

netcdf profiles.001 {                                                           
dimensions:                                                                     
        time = UNLIMITED ; // (6 currently)                                     
        zt = 80 ;                                                               
        zm = 80 ;                                                               
variables:                                                                      
        float rhobf(time, zt) ;                                                 
                rhobf:longname = "Full level density" ;                         
                rhobf:units = "kg/m^3" ;                                        
                rhobf:_FillValue = -999.f ;                                     
        float rhobh(time, zm) ;                                                 
                rhobh:longname = "Half level density" ;                         
                rhobh:units = "kg/m^3" ;                                        
                rhobh:_FillValue = -999.f ;