Custom Query (964 matches)

Modified netcdf_check_var to check few variables for consistency within a tolerance value. The problem is due to round-off. See the following ​forum posting.

For example, now I have:

          CASE ('theta_b')
            CALL netcdf_get_fvar (ng, IDmod, ncname, var_name(i),       &
     &                            FvarS)
            IF (exit_flag.ne.NoError) RETURN

            IF (ABS(theta_b(ng)-FvarS).gt.RoundOff) THEN
              IF (Master) WRITE (stdout,20) TRIM(var_name(i)),          &
     &                                      FvarS, theta_b(ng),         &
     &                                      TRIM(ncname)
              exit_flag=5
              EXIT
            END IF

Many thanks to Mark Hadfield for reporting this problem.

In ROMS/Modules/mod_fourdvar.F, subroutine initialize_fourdvar, the parameter IniVal is declared inside an #ifdef block (lines 503–505):

# ifdef OBSERVATIONS
      real(r8), parameter :: IniVal = 0.0_r8
# endif

but later (lines 758-774) it is used in code that is not conditional on OBSERVATIONS being defined.

PS: I found this while building the CHANNEL_NECK test case.

I am releasing few Matlab pre- and post-processing scripts for ROMS. This by no means are all the scripts that we have but it is a good start. These scripts are well documented. There are a lot of Matlab scripts out there and it is nearly impossible for us to check them for correctness. The ones released here are the official ones. That is, the ones that I maintain. I have been using these scripts for several years.

Currently, the structure of the matlab repository is as follows:

   matlab
             bathymetry
             coastlines
             landmask
             seawater
             utility

Please check the following ​forum message for more details.