Custom Query (986 matches)

I updated several scripts in the matlab/4dvar and matlab/utility directories of the repository to enhance their capabilities. This stuff is complicated and the devil is in the details.

We need for these scripts to be flexible. I assume that more fine tunning is required when the users start using them next week at the 4D-Var workshop.

I added the grid_Lm_Mm_N global attribute to the 4D-Var observation metadata. This is needed in quality control and binning into super observations.

I updated the file wc13_obs.nc in the test/WC13/Data directory of the repository. The file now includes the grid_Lm_Mm_N global attribute with integer values:

// global attributes:
...
                :grid_Lm_Mm_N = 54, 53, 30 ;

This attribute is not really needed inside ROMS, but it is added as a good example of 4D-Var observation metadata. This attribute is needed in quality control and binning or the of the observation data during pre-processing.

It is better to add a global attribute than variables or dimensions to the 4D-Var observation file because of backward compatibility.

Lines 2233 and 2353 at set_avg.F

2231 TIDES(ng)%vbar_tide(i,j,nk)=TIDES(ng)%vbar_tide(i,j,    & 
2232 &                                                 nk)+  &    
2233 &                           OCEAN(ng)%zeta(i,j,Kout)*   & 
2234 &                           tide_harmonics(nk) 

2351 TIDES(ng)%vbar_tide(i,j,nk)=TIDES(ng)%vbar_tide(i,j,    & 
2352 &                                               nk)+    & 
2353 &                           OCEAN(ng)%zeta(i,j,Kout)*   & 
235  &                           tide_harmonics(nk) 

Should be replaced by:

&                           OCEAN(ng)%vbar(i,j,Kout)*   & 

This error was introduced on revision 462 and results in wrong values to vbar_detide when AVERAGE_DETIDE is used.

There is also a small typo on line 170 at wrt_tides.F

      IF (Aout(idu3dD,ng)) THEN 

that should be replaced by:

       IF (Aout(idv2dD,ng)) THEN