Custom Query (986 matches)

Dear ROMS-Team,

ROMS 3.0 Version 145

All ASCII input files (./Roms/External/*.in) specify in the command lines the time references 0000-01-01 00:00 for the input time_ref = -1 and time_ref = 0, while in all ./Roms/Utitliy/def_*.F-files the reference is 0001-01-01 00:00 resulting in a difference of one year.

Please correct the inconsistency. Please keep in mind that, for example, the ROMS-Ocean Modeling Discussion thread “Model Time Clocks” refers to the reference the year 0000. Would you also correct this thread if it is necessary.

Thanks !!

Christian

The 4DVAR surface forcing adjustment was revised to include minimization at other times in addition to initialization. A new input parameter (nSFF) was introduced to allow minimization at other times during strong constraint data assimilation. Here, nSFF is the number of time-steps between adjustment of surface fields. This required several changes to the IO NetCDF files involved in 4DVAR data assimilation. All the changes are internal to ROMS.

Several new arrays with an additional dimension (nFrec) are introduced during the minimization. The number of forcing record to process are:

nFrec = 1 + ntimes / nSFF      (integer operation)

nSFF is either a multiple of ntimes or greater than ntimes. Therefore, there are three posibilities:

(1) If  NSFF > NTIMES,       nFrec = 1  (constant adjustment)
(2) If  NSFF = NTIMES,       nFrec = 2  (initial and final)
(3) If  NSFF < NTIMES,       nFrec > 2  (adjustment every NSFF steps)

We will continue looking and testing this capability.

I added the capability to screen for large and unphysical values for ocean currents and density anomaly in diagnostics file diag.F. Two new variables were added to mod_scalars.F:

        real(r8) :: max_speed = 20.0_r8         ! m/s
        real(r8) :: max_rho = 200.0_r8          ! kg/m3

containing the threshold values that will trigger the model to stop (exit_flag=1) because of blowing-up. It will stop execution before of filling output restart file with NaN.

Recall that the routine diag is called once per baroclinic time-step if ninfo=1. However if the model blow-ups in the barotropic time step loop (step2d), you will still get a lot of NaN.

Many thanks to Richard Schmalz for suggesting this capability.