IF (PtRACERS_EvPrRn(iTrc).NE.UNSET_RL) THEN

        IF ( PTRACERS_diffKr(iTracer).NE.UNSET_RL ) THEN

        PTRACERS_diffKr(iTracer)=UNSET_RL

        PTRACERS_EvPrRn(iTracer)=UNSET_RL

      IF ( SHELFICEsaltTransCoeff .EQ. UNSET_RL )

      SHELFICEsaltTransCoeff   = UNSET_RL

      _RL    UNSET_RL

      PARAMETER ( UNSET_RL     = 1.234567D5 )

      IF ( RiLimit .EQ. UNSET_RL ) THEN

      RiLimit      = UNSET_RL

        timePhase(l)  = UNSET_RL

          levels(k,l) = UNSET_RL

        stat_phase(l) = UNSET_RL

         IF ( timePhase(l).NE. UNSET_RL ) THEN

         IF ( levels(1,l).NE.UNSET_RL ) THEN

             IF ( levels(k,l).NE.UNSET_RL .AND.

             ELSEIF ( levels(k,l).NE.UNSET_RL ) THEN

         IF ( stat_phase(l).NE. UNSET_RL ) THEN

         if ( tAlpha .eq. UNSET_RL ) tAlpha = 2.  _d -4

         if ( sBeta  .eq. UNSET_RL ) sBeta  = 7.4 _d -4

      IF ( alph_AB.NE.UNSET_RL .OR. beta_AB.NE.UNSET_RL ) THEN

         IF (PtRACERS_EvPrRn(iTrc).NE.UNSET_RL) THEN

         IF (PtRACERS_EvPrRn(iTrc).NE.UNSET_RL) THEN

      IF ( delX(1)   .NE. UNSET_RL ) goptCount = goptCount + 1

      IF ( dxSpacing .NE. UNSET_RL ) goptCount = goptCount + 1

      IF ( dxSpacing .NE. UNSET_RL ) THEN

      IF ( delY(1)   .NE. UNSET_RL ) goptCount = goptCount + 1

      IF ( dySpacing .NE. UNSET_RL ) goptCount = goptCount + 1

      IF ( dySpacing .NE. UNSET_RL ) THEN

      IF ( latBandClimRelax .EQ. UNSET_RL) THEN

      IF (delRc(1).NE.UNSET_RL) setCenterDr=.TRUE.

        IF (delRc(K).EQ.UNSET_RL) setCenterDr = .FALSE.

       IF ( delZ(K) .NE. UNSET_RL ) zCoordInputData = .TRUE.

       IF ( delP(K) .NE. UNSET_RL ) pCoordInputData = .TRUE.

       IF ( delR(K) .NE. UNSET_RL ) rCoordInputData = .TRUE.

       IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delZ(K)

       IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delP(K)

       IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delRDefault(K)

      trac_EvPrRn     = UNSET_RL

      zonal_filt_lat  = UNSET_RL

      gravitySign     = UNSET_RL

      viscAstrain     = UNSET_RL

      viscAtension    = UNSET_RL

      tauThetaClimRelax3Dim = UNSET_RL

      tauSaltClimRelax3Dim  = UNSET_RL

      rkFac           = UNSET_RL

      viscAhW  = UNSET_RL

      viscA4W  = UNSET_RL

      viscAhD  = UNSET_RL

      viscAhZ  = UNSET_RL

      viscA4D  = UNSET_RL

      viscA4Z  = UNSET_RL

      viscAz   = UNSET_RL

      viscAr   = UNSET_RL

      viscAp   = UNSET_RL

      diffKzT  = UNSET_RL

      diffKpT  = UNSET_RL

      diffKrT  = UNSET_RL

      diffKzS  = UNSET_RL

      diffKpS  = UNSET_RL

      diffKrS  = UNSET_RL

       diffKrNrT(k) = UNSET_RL

       diffKrNrS(k) = UNSET_RL

       tRef(k)      = UNSET_RL

       sRef(k)      = UNSET_RL

      gBaro    = UNSET_RL

      rhoConst = UNSET_RL

      omega    = UNSET_RL

      hFacMinDr           = UNSET_RL

      hFacMinDz           = UNSET_RL

      hFacMinDp           = UNSET_RL

      rhoConstFresh  = UNSET_RL

      convertFW2Salt = UNSET_RL

      tAlpha              = UNSET_RL

      sBeta               = UNSET_RL

      IF ( gBaro .EQ. UNSET_RL ) gBaro=gravity

      IF ( rhoConst .EQ. UNSET_RL ) rhoConst=rhoNil

      IF ( rhoConstFresh .EQ. UNSET_RL ) rhoConstFresh=rhoConst

      IF ( omega .EQ. UNSET_RL ) THEN

      IF (atm_Rd .EQ. UNSET_RL) THEN

        IF (tRef(k).EQ.UNSET_RL) tRef(k) = tracerDefault

        IF (sRef(k).EQ.UNSET_RL) sRef(k) = tracerDefault

      IF ( viscAhW .EQ. UNSET_RL ) viscAhW = viscAh

      IF ( viscA4W .EQ. UNSET_RL ) viscA4W = viscA4

      IF ( viscAhD .EQ. UNSET_RL ) viscAhD = viscAh

      IF ( viscAhZ .EQ. UNSET_RL ) viscAhZ = viscAh

      IF ( viscA4D .EQ. UNSET_RL ) viscA4D = viscA4

      IF ( viscA4Z .EQ. UNSET_RL ) viscA4Z = viscA4

      IF ( viscAz .NE. UNSET_RL ) zCoordInputData = .TRUE.

      IF ( viscAp .NE. UNSET_RL ) pCoordInputData = .TRUE.

      IF ( viscAr .NE. UNSET_RL ) rCoordInputData = .TRUE.

      IF ( viscAr .EQ. UNSET_RL )          viscAr = viscAz

      IF ( viscAr .EQ. UNSET_RL )          viscAr = viscAp

      IF ( viscAr .EQ. UNSET_RL )          viscAr = viscArDefault

      IF ( diffKzT .NE. UNSET_RL ) zCoordInputData  = .TRUE.

      IF ( diffKpT .NE. UNSET_RL ) pCoordInputData  = .TRUE.

      IF ( diffKrT .NE. UNSET_RL ) rCoordInputData  = .TRUE.

      IF ( diffKrT .EQ. UNSET_RL )          diffKrT = diffKzT

      IF ( diffKrT .EQ. UNSET_RL )          diffKrT = diffKpT

      IF ( diffKrT .EQ. UNSET_RL )          diffKrT = diffKrTDefault

        IF ( diffKrNrT(k).EQ. UNSET_RL ) diffKrSet = .FALSE.

      IF ( diffKzS .NE. UNSET_RL ) zCoordInputData  = .TRUE.

      IF ( diffKpS .NE. UNSET_RL ) pCoordInputData  = .TRUE.

      IF ( diffKrS .NE. UNSET_RL ) rCoordInputData  = .TRUE.

      IF ( diffKrS .EQ. UNSET_RL )          diffKrS = diffKzS

      IF ( diffKrS .EQ. UNSET_RL )          diffKrS = diffKpS

      IF ( diffKrS .EQ. UNSET_RL )          diffKrS = diffKrSDefault

        IF ( diffKrNrS(k).EQ. UNSET_RL ) diffKrSet = .FALSE.

      IF ( hFacMinDz .NE. UNSET_RL ) zCoordInputData = .TRUE.

      IF ( hFacMinDp .NE. UNSET_RL ) pCoordInputData = .TRUE.

      IF ( hFacMinDr .NE. UNSET_RL ) rCoordInputData = .TRUE.

      IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDz

      IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDp

      IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDrDefault

      IF (convertFW2Salt.EQ.UNSET_RL) THEN

      IF ( zonal_filt_lat .NE. UNSET_RL ) THEN

      IF ( gravitySign .NE. UNSET_RL ) THEN

      IF ( trac_EvPrRn .NE. UNSET_RL ) THEN

      IF ( viscAstrain .NE. UNSET_RL ) THEN

      IF ( viscAtension .NE. UNSET_RL ) THEN

      latBandClimRelax  = UNSET_RL

      IF ( tauThetaClimRelax3Dim .NE. UNSET_RL ) THEN

      IF ( tauSaltClimRelax3Dim .NE. UNSET_RL ) THEN

       delZ(K) = UNSET_RL

       delP(K) = UNSET_RL

       delR(K) = UNSET_RL

      dxSpacing = UNSET_RL

      dySpacing = UNSET_RL

      IF ( rkFac .NE. UNSET_RL ) THEN

         IF (PtRACERS_EvPrRn(iTrc).NE.UNSET_RL) THEN

       IF (temp_EvPrRn.NE.UNSET_RL) THEN

       IF (salt_EvPrRn.NE.UNSET_RL) THEN

      IF ( latBandClimRelax.EQ.UNSET_RL ) THEN

      temp_EvPrRn = UNSET_RL

      atm_Rd = UNSET_RL

      alph_AB           = UNSET_RL

      beta_AB           = UNSET_RL

       delRc(k)           = UNSET_RL

       delX(i)            = UNSET_RL

       delY(j)            = UNSET_RL

      IF ( latBandClimRelax.EQ.UNSET_RL ) THEN

      IF ( GGL90TKEbottom .EQ. UNSET_RL ) THEN

      GGL90TKEbottom       = UNSET_RL