AKv and AKt in initial condition file

[rtoste]

Hi everyone,

I've just download the new code (svn 1077) and I've been failing to run a 4dvar application if it's not starting from a previous 4dvar output file.
It seems now AKv and AKt are needed in the initial condition file.
This is okay when using dai.nc or a time instant from the previous assimilation cycle but what should I do to run the first cycle?
I'm trying to use the first initial condition based on CMEMS analysis (it works for svn 1005)

A common initial condition file (same for a free run) returns this type of error:

Code: Select all

 WRT_INI_NF90     - NLM: Writing initial state fields,                   2021-07-03 00:00:00.00
                      (Grid 01, t = 8949.0000, File: ini_ciclo.nc, Rec=0002, Index=1)
                   - free-surface
                      (Min = -4.36839044E-01 Max =  1.45649564E+00)
                   - vertically integrated u-momentum component
                      (Min = -3.42652351E-01 Max =  3.63797605E-01)
                   - vertically integrated v-momentum component
                      (Min = -3.30271959E-01 Max =  3.82356495E-01)
                   - u-momentum component
                      (Min = -4.36231226E-01 Max =  6.44137383E-01)
                   - v-momentum component
                      (Min = -7.14611053E-01 Max =  5.83789170E-01)
                   - potential temperature
                      (Min =  0.00000000E+00 Max =  2.66033096E+01)
                   - salinity
                      (Min =  0.00000000E+00 Max =  3.76456642E+01)
 Found Error: **   Line: 1881     Source: ROMS/Modules/mod_netcdf.F, netcdf_inq_varid

NETCDF_INQ_VARID - error while inquiring ID for variable:  AKv
                    in input file:  ini_ciclo.nc
                    call from:  ROMS/Utility/wrt_ini.F, wrt_ini_nf90
                    NetCDF: Variable not found                                                      
 Found Error: **   Line: 424      Source: ROMS/Utility/wrt_ini.F, wrt_ini_nf90

 WRT_INI_NF90 - error while writing variable: AKv
At line 363 of file wrt_ini.f90 (unit = 6, file = 'stdout')
Fortran runtime error: Expected INTEGER for item 2 in formatted transfer, got CHARACTER
(/,' WRT_INI_NF90 - error while writing variable: ',a,             /,16x,'into i

Thanks!

[jivica]

What I would do in your place is:
a) spinup your NLM mode to some physical state (i.e. 1 week and save all you need into history or restart file) and use that for your DA init
b) if you don't want to spinup then start from 0 - no mixing state, it will mix quite fast

Cheers
I

[rtoste]

Thanks for your suggestion, Ivica.
Cheers