Ill defined namelist variable nb_jpk_bdy

[jamesharle]

Context

Interpolation on-the-fly for BDY: nb_jpk_bdy (number of vertical levels in the input BDY files) is defined as an INTEGER - it should be an array to accommodate multiple open boundaries. It is also define under the nambdy and nambdy_dta namelists!

Analysis

If interpolation on the fly is required for >1 number of open boundary conditions with differing numbers of vertical levels the solution will be incorrect.

Recommendation

Redefine nb_jpk_bdy as an integer array. Update associated references to this variable to reflect the change. Only define nb_jpk_bdy in the namebdy_dta namelist / remove all instances from the nambdy namelist

(Apply changes to the trunk too)

[smasson]

another strange point: nambdy_dta is read within a jpby loop...

cn_dir_array(jbdy) = cn_dir
ln_full_vel_array(jbdy) = ln_full_vel

have the same value for each jbdy

[smasson]

I am working on this routine, I think it needs some cleaning...

[jamesharle]

I'm not too sure about your comment re: same value for ln_full_vel

as ln_full_vel is read in in nambdy_tra (and there could be multiple instances of this namelist block if nb_bdy > 1) ln_full_vel could have different values for each nambdy_tra block defined in the namelist.

The only untidy part is that because there is no REWIND of namelist_cfg in the jbdy loop the code attempts to read nambdy_tra jbdy times, which if nb_bdy > 1 is not great as the nambdy_tra block is only defined once in the namelist_cfg file.

One solution would be to put the REWIND statement inside the loop:

 510d509
<       REWIND(numnam_ref)
513a513
>          REWIND(numnam_ref)

to become:

      ! Read namelists
      ! --------------
      REWIND(numnam_cfg)
      jfld = 0
      DO jbdy = 1, nb_bdy
         REWIND(numnam_ref)
         IF( nn_dta(jbdy) == 1 ) THEN
            READ  ( numnam_ref, nambdy_dta, IOSTAT = ios, ERR = 901)
901         IF( ios /= 0 )   CALL ctl_nam ( ios , 'nambdy_dta in reference namelist', lwp )
            READ  ( numnam_cfg, nambdy_dta, IOSTAT = ios, ERR = 902 )
902         IF( ios >  0 )   CALL ctl_nam ( ios , 'nambdy_dta in configuration namelist', lwp )
            IF(lwm) WRITE( numond, nambdy_dta )

            cn_dir_array(jbdy) = cn_dir
            ln_full_vel_array(jbdy) = ln_full_vel

[jamesharle]

Just an update on the original recommendation:

Only define nb_jpk_bdy in the nambdy_dta namelist / remove all instances from the nambdy namelist

nb_jpk_bdy has be left in the nambdy namelist as it's required in bdyini.F90 to allocate array sizes to read in boundary data. It has been removed from the reading in of NAMELIST/nambdy_tra/ in bdydta.F90

[jamesharle]

In 10933:

Error: Failed to load processor CommitTicketReference

No macro or processor named 'CommitTicketReference' found

[jamesharle]

In 10934:

Error: Failed to load processor CommitTicketReference

No macro or processor named 'CommitTicketReference' found

[jamesharle]

In 10951:

Error: Failed to load processor CommitTicketReference

No macro or processor named 'CommitTicketReference' found

[jamesharle]

In 10952:

Error: Failed to load processor CommitTicketReference

No macro or processor named 'CommitTicketReference' found