anchor.diff:    If ( varSize(1) /= SIZE(varValues,1) .or. &
          varSize(2) /= SIZE(varValues,2) .or. &
          varSize(3) /= SIZE(varValues,3) ) Then
  -     Print*, "Error: var size mismatches"
  +     Print*, "[err]: var size mismatches"
        Write(*,"(2(A,I8))") "dim(1)=",varSize(1), ", dim_input(1)=", SIZE(varValues,1)
        Write(*,"(2(A,I8))") "dim(2)=",varSize(2), ", dim_input(2)=", SIZE(varValues,2)
        Write(*,"(2(A,I8))") "dim(3)=",varSize(3), ", dim_input(3)=", SIZE(varValues,3)
  
positive.code: ! Open NetCDf file 
    errStatus = NF90_Open( TRIM(fileName), NF90_NOWRITE, ncId)
    If ( errStatus /= NF90_NOERR ) Then
       Print*, "Error: fail to get the file id:",TRIM(fileName)
       Stop
    EndIf
  
  ! Inquire varId
    errStatus = NF90_Inq_Varid( ncId, TRIM(varName), varId )
    If ( errStatus /= NF90_NOERR ) Then
       Print*, "Error: fail to get the var id: ", TRIM(varName)
       Stop
    EndIf 
  
  
negative.code: ! Open NetCDf file 
    errStatus = NF90_Open( TRIM(fileName), NF90_NOWRITE, ncId)
    If ( errStatus /= NF90_NOERR ) Then
       Print*, "[err]: fail to get the file id:",TRIM(fileName)
       Stop
    EndIf
  
  ! Inquire varId
    errStatus = NF90_Inq_Varid( ncId, TRIM(varName), varId )
    If ( errStatus /= NF90_NOERR ) Then
       Print*, "[err]: fail to get the var id: ", TRIM(varName)
       Stop
    EndIf 
  
  