Custom Query (964 matches)

I suggest 3 updates to the estaury test so that it is more consistent with the Ocean Modelling 2005 paper.

1) ROMS/External/ocean_estuary_test.in: Change Lm from 50 to 200, and N from 10 to 20.

2) ROMS/External/sediment_estuary_test.in: Change the erosion rate from 1e-5 to 1e-4.

3) ROMS/Functionals/ana_sediment.h: Create it's own settings and change init bed thickness to 1 mm.

Theses files are updated in my branch.

Correction in swanout1.F to allow table output from SWAN during a coupled run.

1) Correct file Waves/SWAN/Src/swanout1.F near lines 302, replace

        IF (RTYPE(1:2) .EQ. 'SP') THEN
          IVTYPE=500
        ELSE
          IVTYPE=CORQ%IVTYP(1)
        END IF

with

        IF (RTYPE(1:3) .EQ. 'BLK') THEN
          IVTYPE=CORQ%IVTYP(1)
        ELSE
          IVTYPE=500
        END IF

2) Modify file Waves/SWAN/External/swan_inlet_test.in to add the line:

TABLE 'point1' HEADER 'point1.table' HSIGN DIR WLEN OUTPUT 20000101.000000 1 HR

This will allow a test for this output option.

This is a complete update for ROMS private storage arrays that are automatically allocated when entering the _tile routines. These array needed to be redisigned to allow the different types of nesting: refining, composite, and mosaics. These arrays where temporarily changed to global dimensions in revision 168 (​src:ticket:142). This is fine in distributed-memory applications but the private storage is much larger than needed in shared-memory and serial with partitions applications. This is not a bug but the model required larger memory in such applications. Perhaps, you didn't noticed...

A new variable is not introduce to tile.h to determine the size of the private storage arrays for each nested grid:

!
!  Set horizontal starting and ending indices for automatic private storage
!  arrays.
!
      IminS=BOUNDS(ng)%Istr(tile)-3
      ImaxS=BOUNDS(ng)%Iend(tile)+3
      JminS=BOUNDS(ng)%Jstr(tile)-3
      JmaxS=BOUNDS(ng)%Jend(tile)+3

So now each tile routines have the following additional arguments: IminS, ImaxS, JminS, and JmaxS, for example:

      CALL xxxxx_tile (ng, tile,                                        &
     &                 LBi, UBi, LBj, UBj,                              &
     &                 IminS, ImaxS, JminS, JmaxS,                      &
     &                 ...)

The values of IminS, ImaxS, JminS, and JmaxS will change in the future to allow for the different size of the contact areas. However, the change will be limited to tile.h. Notice that now the numerous additions and substractions per time-step from the CPP tokens in globaldefs.h are eliminated:

#define PRIVATE_1D_SCRATCH_ARRAY IminS:ImaxS
#define PRIVATE_2D_SCRATCH_ARRAY IminS:ImaxS,JminS:JmaxS

These definitions are keept for compatability but not used in the new version.

See the following forum ​post for more information. Many thanks to John Warner for remind me about this issue.