Custom Query (964 matches)

Finally, I finished coding the 4DVAR adjustment for the open boundaries (ADJUST_BOUNDARY). This is the final piece missing in the 4DVAR data assimilation (strong and weak constraint) algorithms. We will be able to adjust initial condition, boundary conditions, surface forcing, and interior model error. I have been working on this new option for several weeks. This was technically tricky because of the nature of the boundary arrays which are only processed by tiles adjacent to the open boundaries. However, all the tiles participate in the horizontal convolutions. I rewrote the convolutions 5 times until I found an efficient design.

There are several new files associated with this option:

Nonlinear/bc_bry2d.F
Nonlinear/bc_bry3d.F
Nonlinear/conv_bry2d.F
Nonlinear/conv_bry3d.F

Adjoint/ad_bc_bry2d.F
Adjoint/ad_bc_bry3d.F
Adjoint/ad_conv_bry2d.F
Adjoint/ad_conv_bry3d.F

Tangent/tl_conv_bry2d.F
Tangent/tl_conv_bry3d.F

Utility/nf_fread2d_bry.F
Utility/nf_fread3d_bry.F
Utility/nf_fwrite2d_bry.F
Utility/nf_fwrite3d_bry.F

Data/ROMS/CDL/s4dvar_std_b.cdl
Data/ROMS/CDL/s4dvar_std_f.cdl
Data/ROMS/CDL/s4dvar_std_i.cdl
Data/ROMS/CDL/s4dvar_std_m.cdl

Obviously, we will fine tunning this algorithms in the coming weeks. We will let you know when we are done. Notice that the changes to several routines do not affect previous capabilities due to the modularity of ROMS algorithms.

This new algorithm is similar to the surface forcing adjustment. We can specify in input script ocean.in how often the open boundaries are adjusted. See parameter NOBC:

!  NOBC        Number of time-steps between 4DVAR adjustment of open boundary
!              fields. In strong constraint 4DVAR, it is possible to adjust
!              open boundaries at other time intervals in addition to initial
!              time.  This parameter is used to store the appropriate number
!              of open boundary records in the output history NetCDF files:
!              1+NTIMES/NOBC records. NOBC must be a factor of NTIMES or
!              greater than NTIMES. If NOBC > NTIMES, only one record is
!              stored in the NetCDF files and the adjustment is for constant
!              forcing with constant correction. This parameter is only
!              relevant in 4DVAR when activating ADJUST_BOUNDARY.

Many thanks to Andy Moore for his help in the design of this new option and for coding the variational differentiation of the adjoint-based state variables associated with the open boundaries.

The boundary arrays are now strictly defined and allocated when they are really needed. There are several internal options to do so:

EAST_FSOBC, WEST_FSOBC, SOUTH_FSOBC, NORTH_FSOBC
EAST_M2OBC, WEST_M2OBC, SOUTH_M2OBC, NORTH_M2OBC
EAST_M3OBC, WEST_M3OBC, SOUTH_M3OBC, NORTH_M3OBC
EAST_TOBC,  WEST_TOBC,  SOUTH_TOBC,  NORTH_TOBC

These C-preprocessing options are defined in globaldefs.h and there is not need for the user to manually activate them. However, it is wise to use them in the analytical routines ana_fsobc.h, ana_m2obc.h, ana_m3obc.h, and ana_tobc.h. Sometimes anlytical expressions are coded but not used when playing with various set-ups in the test problems.

There are couple of messages in the forum reporting compilation problems with unused analytical expressions. Using the above options in such expressions solve such compilation errors.

The EcoSim was updated to use the vertical semi-Lagrangian sinking flux for fecal material available in ROMS. This generic sinking is identical to that found in other ROMS ecosystem and sediment models. It also includes a simple biological sedimentation (BIO_SEDIMENT option) which remineralizes particular nitrogen and returns it into the dissolved nitrate pool.

Many thanks to Bronwyn Cahill for implementing and testing these algorithms.