Ocean Modeling Discussion

I just ran into the same issue, I occurred in the same place pointed out by powellb, also using a biological model, so that more than 2 tracers are defined and MT>2. ! ! Determine metadata for quadratic tracer averages. ! SELECT CASE (Vinfo(1)) CASE ('tracer2') varid=varid-1 DO i=1,MT varid=varid+1 ...

I created a script called run_roms_da that may be useful for those performing ROMS data assimilation and running multiple consecutive 4DVar cycles. This script has been created a while ago and is now used by several people in our UCSC data assimilation group. It has become more sophisticated over ti...

I am currently working on a project generating tangent linear and adjoint models without the need for tangent linear or adjoint code. When comparing my results with the "reference" results based on the code in the repository, a large discrepancy between the results led me to discover a bug...

I just encountered an error when I tried to compile the most recent ROMS version with IS4DVAR and BGQC defined. There is an unterminated if-clause in obs_write.F: # ifndef WEAK_CONSTRAINT ! ! Set observation scale (ObsScale). The scale factor is used ! for screenning of the observations. This scale ...

There is a bug in the current version (820) of the W4DPSAS code. It appears to be relatively minor but may point to a bigger problem. The error message I get is: NETCDF_CLOSE - error while writing global attribute: bio_file file: wc12_fwd_36524_001.nc call from: w4dpsas_ocean.h, ROMS_run The error o...

For the sake of completeness, if one makes changes to the super_obs routine in the way suggested in the post above, one extra change is required ass well. In summary, the original code % Compute the index in each dimension of the grid cell in which the % observation is located. Xbin = 1.0 + floor((V...

The super_obs matlab function (matlab/4dvar/super_obs.m in the repository) is meant to merge observations of the same type together, if they are located very closely in both time and space. However, it may merge observations in a somewhat arbitrary fashion and does not guarantee the removal of two o...

There is a bug in the NEMURO code that causes 3 of the parameters to be reported incorrectly. These parameters get read in correctly but are written to the log file incorrectly. So this is a relatively minor bug, yet it can be quite confusing to the user. The incorrect lines in "ROMS/Nonlinear/...

The adjoint variables ad_ExpAtt and ad_Itop in the NPZD iron model are not initialized properly: For ad_ExpAtt, the debugger message is: The variable 'ad_biology_mod_mp_ad_biology_tile_$AD_EXPATT' is being used without being defined A "ad_ExpAtt=0.0_r8" and "ad_Itop=0.0_r8" are m...

There is a bug in the 4DVAR code when dealing with observations that have an obs_depth in between 0 and 1. I kept running into problems (negative Ritz values) with a particular model run until I switched on the debug mode which told me about an out of range subscript. The problem occurs in extract_o...

Great! In ROMS 696, I noticed the first bug was still present but having just updated to the latest version 697 I can confirm it is fixed. Thanks

In ROMS: The current version (ROMS 696) of npzd_iron_inp.h contains 2 errors that cause a segmentation fault and a compilation failure, respectively. In line 209 and line 219: There is a missing input argument in function call to load_lbc (this causes a segmentation fault while reading in the biolog...

While attempting to start a 4DVAR data assimilation cycle I encountered the following error message: INQUIRE - unable to find requested variable: zeta in file: boundaryfile.nc where boundaryfile.nc is the file containing the boundary conditions for the simulation. The error message made me wonder wh...

In the data set I downloaded from the ECMWF website the values for e (evaporation) are negative while those for tp (total precipitation) are positive. This makes sense since the description mentions that all downward fluxes are positive. In this case e+tp is the net downward flux. The values for e-t...

I noticed that here is an issue with the way swflux is computed above: ... swflux (cm day-1) (e - tp) * 100 / (3/24); 0.125 day step According to the ECMWF website all downward fluxes are positive (http://www.ecmwf.int/publications/manuals/d/gribapi/param/detail/format=mars/pid=182/) and indeed by l...

If I am not mistaken, it seems like the code needs to check if the file it is supposed to write to exists to distinguish between the case where one simply wants to continue a run and has already got some (potentially unfinished) output files present in the output directory and the case where one wan...

It seems to me that the error occurs whenever multiple snapshots are being written into one file and the starting time of the model is not a multiple of ndefhis. This was not a problem in older versions (e.g. Version 126 2007-11-08 01:44:08Z). Setting nrrec=0 results in a similar, seemingly related ...

This concerns a fairly technical bug that only few ROMS users might encounter, it affects the current version (Version 403 2009-10-01 21:17:45Z) and maybe older versions. To reproduce this error a restart from a previous file is required. Then, for example, a change in the variable ndefHIS (in betwe...

In the new version the error has disappeared. Thanks for the quick fix!

Yes, there is a river_time variable in the river file. Like I mentioned, it runs from -1096 to 933 and the appropriate variables are river_time dependent.

This is a bug concerning the current version of ROMS (Version 396 2009-09-11 18:53:38Z) but it also affects earlier versions at least dating back to September 1st. The error I get when trying to run ROMS is: GET_CYCLE - starting time for variable: river is greater than current model time. TMIN = 1.0...