IF ( GGL90TKEbottom .EQ. UNSET_RL ) THEN

      GGL90TKEbottom       = UNSET_RL

      IF ( GM_Visbeck_maxSlope .EQ. UNSET_RL )

      GM_Visbeck_maxSlope = UNSET_RL

       IF ( temp_addMass .EQ. UNSET_RL ) THEN

       IF ( salt_addMass .EQ. UNSET_RL ) THEN

      IF ( ICEFRONTsaltTransCoeff .EQ. UNSET_RL )

      ICEFRONTsaltTransCoeff   = UNSET_RL

      IF (minKPPhbl .EQ. UNSET_RL) THEN

      minKPPhbl = UNSET_RL

        layers_G(k) = UNSET_RL

        IF ( layers_G(k) .EQ. UNSET_RL ) THEN

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

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

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

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

      IF ( RiLimit .EQ. UNSET_RL ) THEN

      RiLimit      = UNSET_RL

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

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

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

        PTRACERS_diffKr(iTracer)=UNSET_RL

        PTRACERS_EvPrRn(iTracer)=UNSET_RL

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

      SaltPlumeCriterion  = UNSET_RL

       IF ( SaltPlumeCriterion .EQ. UNSET_RL ) THEN

      IF ( SEAICE_deltaTevp .NE. UNSET_RL ) THEN

      SEAICE_deltaTevp   = UNSET_RL

      ICE2WATR           = UNSET_RL

      SEAICE_drag_south       = UNSET_RL

      SEAICE_waterDrag_south  = UNSET_RL

      SEAICE_dryIceAlb_south  = UNSET_RL

      SEAICE_wetIceAlb_south  = UNSET_RL

      SEAICE_drySnowAlb_south = UNSET_RL

      SEAICE_wetSnowAlb_south = UNSET_RL

      HO_south                = UNSET_RL

      SEAICE_gamma_t     = UNSET_RL

      SEAICE_gamma_t_frz = UNSET_RL

      SEAICE_availHeatFrac = UNSET_RL

      SEAICE_availHeatFracFrz = UNSET_RL

      SEAICE_sensHeat    = UNSET_RL

      SEAICE_latentWater = UNSET_RL

      SEAICE_latentIce   = UNSET_RL

      IF (ICE2WATR.EQ.UNSET_RL) ICE2WATR = SEAICE_rhoIce*recip_rhoConst

      IF (SEAICE_drag_south       .EQ. UNSET_RL)

      IF (SEAICE_waterDrag_south  .EQ. UNSET_RL)

      IF (SEAICE_dryIceAlb_south  .EQ. UNSET_RL)

      IF (SEAICE_wetIceAlb_south  .EQ. UNSET_RL)

      IF (SEAICE_drySnowAlb_south .EQ. UNSET_RL)

      IF (SEAICE_wetSnowAlb_south .EQ. UNSET_RL)

      IF (HO_south                .EQ. UNSET_RL)

      IF ( SEAICE_deltaTevp .NE. UNSET_RL ) SEAICEuseEVP = .TRUE.

      IF ( SEAICE_availHeatFrac .NE. UNSET_RL .AND.

     &     SEAICE_gamma_t       .NE. UNSET_RL ) THEN

      IF ( SEAICE_availHeatFracFrz .NE. UNSET_RL .AND.

     &     SEAICE_gamma_t_frz      .NE. UNSET_RL ) THEN

      IF ( SEAICE_availHeatFrac .EQ. UNSET_RL ) THEN

      IF ( SEAICE_gamma_t .NE. UNSET_RL .AND.

      IF ( SEAICE_availHeatFracFrz .EQ. UNSET_RL ) THEN

      ELSEIF ( SEAICE_gamma_t_frz .NE. UNSET_RL ) THEN

      IF ( SEAICE_sensHeat    .NE. UNSET_RL ) THEN

      IF ( SEAICE_latentWater .NE. UNSET_RL ) THEN

      IF ( SEAICE_latentIce   .NE. UNSET_RL ) THEN

      IF ( selectAddFluid.NE.0 .AND. temp_addMass.NE.UNSET_RL ) THEN

      IF ( selectAddFluid.NE.0 .AND. salt_addMass.NE.UNSET_RL ) THEN

      IF ( SHELFICEsaltTransCoeff .EQ. UNSET_RL )

      SHELFICEsaltTransCoeff   = UNSET_RL

       IF (temp_EvPrRn.NE.UNSET_RL) THEN

       IF (salt_EvPrRn.NE.UNSET_RL) THEN

        IF (temp_EvPrRn.NE.UNSET_RL) THEN

        IF (salt_EvPrRn.NE.UNSET_RL) THEN

        IF (temp_EvPrRn.NE.UNSET_RL) THEN

        IF (salt_EvPrRn.NE.UNSET_RL) THEN

      hNewIceMax = UNSET_RL

      thSIce_dtTemp  = UNSET_RL

      tauRelax_MxL_salt = UNSET_RL

      IF ( hNewIceMax .EQ. UNSET_RL ) hNewIceMax = hiMax

      IF ( thSIce_dtTemp .EQ. UNSET_RL ) thSIce_dtTemp=thSIce_deltaT

      IF ( tauRelax_MxL_salt .EQ. UNSET_RL ) THEN

      IF ( latBandClimRelax.EQ.UNSET_RL ) THEN

      IF ( latBandClimRelax.EQ.UNSET_RL ) THEN

         IF ( tAlpha .EQ. UNSET_RL ) tAlpha = 2.  _d -4

         IF ( sBeta  .EQ. UNSET_RL ) sBeta  = 7.4 _d -4

      IF ( thetaMin .NE. UNSET_RL ) THEN

      IF ( phiMin .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 ( xgOrigin .EQ. UNSET_RL ) xgOrigin = 0.

      IF ( ygOrigin .EQ. UNSET_RL ) THEN

      IF ( latBandClimRelax .EQ. UNSET_RL) THEN

       IF ( delRc(k).EQ.UNSET_RL ) THEN

       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 ) THEN

      IF (temp_addMass .EQ. UNSET_RL) temp_addMass = temp_EvPrRn

      IF (salt_addMass .EQ. UNSET_RL) salt_addMass = salt_EvPrRn

      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

      thetaMin        = UNSET_RL

      phiMin          = UNSET_RL

      viscAhW  = UNSET_RL

      viscA4W  = UNSET_RL

      viscAhD  = UNSET_RL

      viscAhZ  = UNSET_RL

      viscA4D  = UNSET_RL

      viscA4Z  = UNSET_RL

      viscAz   = UNSET_RL

      viscAp   = UNSET_RL

      viscAr   = UNSET_RL

      diffKzT  = UNSET_RL

      diffKpT  = UNSET_RL

      diffKrT  = UNSET_RL

      diffKzS  = UNSET_RL

      diffKpS  = UNSET_RL

      diffKrS  = UNSET_RL

       viscArNr(k)  = 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

      implicitNHPress     = UNSET_RL

      IF ( gBaro .EQ. UNSET_RL ) gBaro=gravity

      IF ( rhoConst .EQ. UNSET_RL ) rhoConst=rhoNil

      IF ( rhoConstFresh .EQ. UNSET_RL ) rhoConstFresh=rhoConst

      IF ( implicitNHPress.EQ.UNSET_RL )

      IF ( omega .EQ. UNSET_RL ) THEN

      IF (atm_Rd .EQ. UNSET_RL) THEN

      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 ( viscArNr(k).EQ.UNSET_RL ) THEN

      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 ) THEN

      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 ) THEN

      IF (diffKrBLEQsurf .EQ. UNSET_RL) diffKrBLEQsurf = diffKrBL79surf

      IF (diffKrBLEQdeep .EQ. UNSET_RL) diffKrBLEQdeep = diffKrBL79deep

      IF (diffKrBLEQscl  .EQ. UNSET_RL) diffKrBLEQscl  = diffKrBL79scl

      IF (diffKrBLEQHo   .EQ. UNSET_RL) diffKrBLEQHo   = diffKrBL79Ho

      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

      epsAB_CD          = UNSET_RL

      latBandClimRelax  = UNSET_RL

      IF ( tauThetaClimRelax3Dim .NE. UNSET_RL ) THEN

      IF ( tauSaltClimRelax3Dim .NE. UNSET_RL ) THEN

        IF ( epsAB_CD .EQ. UNSET_RL ) epsAB_CD = abEps

       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 ( hybSigmFile.EQ.' ' .AND. rSigmaBnd.EQ.UNSET_RL ) THEN

        IF ( delX(i).EQ.UNSET_RL ) THEN

        IF ( delY(j).EQ.UNSET_RL ) THEN

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

         IF (tRef(k).NE.UNSET_RL) THEN

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

         IF (sRef(k).NE.UNSET_RL) THEN

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

         IF (PTRACERS_EvPrRn(iTrc).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 (temp_EvPrRn.NE.UNSET_RL) THEN

        IF (salt_EvPrRn.NE.UNSET_RL) THEN

        IF (temp_EvPrRn.NE.UNSET_RL) THEN

        IF (salt_EvPrRn.NE.UNSET_RL) THEN

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

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

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

      diffKrBLEQsurf      = UNSET_RL

      diffKrBLEQdeep      = UNSET_RL

      diffKrBLEQscl       = UNSET_RL

      diffKrBLEQHo        = UNSET_RL

      temp_EvPrRn = UNSET_RL

      temp_addMass = UNSET_RL

      salt_addMass = UNSET_RL

      atm_Rd = UNSET_RL

      alph_AB           = UNSET_RL

      beta_AB           = UNSET_RL

      rSigmaBnd           = UNSET_RL

       delRc(k)           = UNSET_RL

      xgOrigin            = UNSET_RL

      ygOrigin            = UNSET_RL

       delX(i)            = UNSET_RL

       delY(j)            = UNSET_RL

      _RL     UNSET_RL

      PARAMETER ( UNSET_RL     = 1.234567D5 )

      undef = UNSET_RL

      undef = UNSET_RL

         IF ( misvalFlt(listId) .NE. UNSET_RL ) THEN

        timePhase(l)  = UNSET_RL

        missing_value(l)     = UNSET_RL

          levels(k,l) = UNSET_RL

        stat_phase(l) = UNSET_RL

        misvalFlt(n) = UNSET_RL

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

         IF ( missing_value(l) .NE. UNSET_RL )

         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

      undef = UNSET_RL

      siobNperiod        = UNSET_RL

      siobSperiod        = UNSET_RL

      siobEperiod        = UNSET_RL

      siobWperiod        = UNSET_RL

      IF(siobNperiod    .EQ.UNSET_RL) siobNperiod     = obcsNperiod

      IF(siobSperiod    .EQ.UNSET_RL) siobSperiod     = obcsSperiod

      IF(siobEperiod    .EQ.UNSET_RL) siobEperiod     = obcsEperiod

      IF(siobWperiod    .EQ.UNSET_RL) siobWperiod     = obcsWperiod