Custom Query (964 matches)

Corrected couple of bugs:

• get_nudgcoef.F: Corrected index idM3nc when cheching fro 3D momenutum nudging coefficients.

       IF (LnudgeM3CLM(ng)) THEN
         IF (.not.find_string(var_name,n_var,Vname(1,idM3nc),            &
      &                       NUD(ng)%Vid(idM3nc))) THEN
           IF (Master) WRITE (stdout,10) TRIM(Vname(1,idM3nc)),          &
      &                                  TRIM(ncname)
           exit_flag=2
           RETURN
         END IF
       END IF
  

  Many thank to Ryan Holmes for ​reporting this typo.

• def_rst.F: Corrected the output dimension for wetting and drying variables wetdry_mask_u and wetdry_mask_v. We need to have instead:

  !
  !  Define wet/dry mask on U-points.
  !
          Vinfo( 1)=Vname(1,idUwet)
          Vinfo( 2)=Vname(2,idUwet)
          Vinfo( 3)=Vname(3,idUwet)
          Vinfo( 9)='land'
          Vinfo(10)='water'
          Vinfo(14)=Vname(4,idUwet)
          Vinfo(16)=Vname(1,idtime)
          Vinfo(22)='coordinates'
          Aval(5)=REAL(Iinfo(1,idUwet,ng),r8)
          status=def_var(ng, iNLM, RST(ng)%ncid, RST(ng)%Vid(idUwet),     &
       &                 NF_FOUT, nvd3, su2dgrd, Aval, Vinfo, ncname,     &
       &                 SetFillVal = .FALSE.)
          IF (exit_flag.ne.NoError) RETURN
  !
  !  Define wet/dry mask on V-points.
  !
          Vinfo( 1)=Vname(1,idVwet)
          Vinfo( 2)=Vname(2,idVwet)
          Vinfo( 3)=Vname(3,idVwet)
          Vinfo(14)=Vname(4,idVwet)
          Vinfo(16)=Vname(1,idtime)
          Vinfo( 9)='land'
          Vinfo(10)='water'
          Vinfo(22)='coordinates'
          Aval(5)=REAL(Iinfo(1,idVwet,ng),r8)
          status=def_var(ng, iNLM, RST(ng)%ncid, RST(ng)%Vid(idVwet),     &
       &                 NF_FOUT, nvd3, sv2dgrd, Aval, Vinfo, ncname,     &
       &                 SetFillVal = .FALSE.)
          IF (exit_flag.ne.NoError) RETURN
  
  

  Many thanks to Lyon Lanerolle for bringing this to my attention.

WET_DRY and AVERAGES_DETIDE options fail to compile together. The solution is very simple:

1) Edit the file set_avg.F

2) Add the usage of mod_grid module in the set_detide_tile subroutine as in the other subroutines in that file. So, just add:

# ifdef WET_DRY
USE mod_grid
# endif

But it would be nice if this correction is added to ROMS files in the repository.

I commented this in the ​forum.

Cheers, Alexis

editmask.m throws some errors due to conflicts with new or incomplete grid files when variable hraw has not been loaded.

These can be solved by making a few changes to read_mask.m

A grid file may have a variable hraw that has not been filled with data. This will cause an error in editmask. I suggest test for the case hraw has an empty dimension:

case {Vname.hraw}
  got.hraw=true;
  if any(V.Variables(n).Size==0)
     got.hraw=false;
  end

read_mask presently only checks for spherical in the form of a character, and always returns false even when spherical=1 (numeric). Make this robust and compatible with both conventions with:

 if (got.spher)
    spher=nc_read(ncfile,Vname.spher);
    if ischar(spher)
       if strcmpi(spher,'T')
          spherical=true;
       end
    else
       spherical=spher;
    end
  end