Opened 10 years ago

Closed 10 years ago

#650 closed bug (Fixed)

Corrected bugs in inp_par.F

Reported by: arango Owned by: arango
Priority: major Milestone: Release ROMS/TOMS 3.7
Component: Nonlinear Version: 3.7
Keywords: Cc:

Description

Several reported bugs were corrected:

  • inp_par.F: Corrected the Rad+Nud boundary conditions for passive tracer when the compact specification is used and input data is required. Many thanks to Mark Hadfield for reporting this bug and a solution. We just needed to add:
    #ifdef SOLVE3D
    !
    !  If processing tracers and last standard input entry (Icont=0), set
    !  unspecified tracer values to the last tracer entry.
    !
            IF ((iTrcStr.gt.0).and.(iTrcEnd.gt.0)) THEN
              IF ((Icont.eq.0).and.(ifield.lt.isTvar(iTrcEnd))) THEN
                DO i=ifield+1,isTvar(iTrcEnd)
                  DO ibry=1,4
                    ...
                    S(ibry,i,igrid)%acquire   = S(ibry,ifield,igrid)%acquire
                  END DO
                  ic=ic+1
                END DO
              END IF
            END IF
    #endif
    
  • inp_par.F: Corrected the momentum LnudgeM3CLM switch. Many thanks to Xueming zhu and Mark Hadfield for reporting this. We just need to have:
    !
    !  Set internal switch for nudging to climatology fields.
    !
    #ifdef SOLVE3D
            IF (LnudgeM3CLM(ng).or.                                         &
         &      LnudgeM3CLM(ng).or.                                         &
         &      ANY(LnudgeTCLM(:,ng))) THEN
              Lnudging(ng)=.TRUE.
            END IF
    #else
            IF (LnudgeM2CLM(ng)) THEN
              Lnudging(ng)=.TRUE.
            END IF
    #endif
    
  • def_rst.F: Corrected typo when defining wetdry_psi mask. This does not affect solution. Many thanks to Kate Hedstrom for reporting this. We just need to have:
            Aval(5)=REAL(Iinfo(1,idPwet,ng),r8)
    

The other changes are for readability.

Change History (1)

comment:1 by arango, 10 years ago

Resolution: Fixed
Status: newclosed
Note: See TracTickets for help on using tickets.